Justin B. Lack

The Drosophila Genome Nexus: A Population Genomic Resource of 623 Drosophila melanogaster Genomes, Including 197 from a Single Ancestral Range Population

Hundreds of wild-derived Drosophila melanogaster genomes have been published, but rigorous comparisons across data sets are precluded by differences in alignment methodology. The most common approach to reference-based genome assembly is a single round of alignment followed by quality filtering and variant detection. We evaluated variations and extensions of this approach and settled on an assembly strategy that utilizes two alignment programs and incorporates both substitutions and short indels to construct an updated reference for a second round of mapping prior to final variant detection. Utilizing this approach, we reassembled published D. melanogaster population genomic data sets and added unpublished genomes from several sub-Saharan populations. Most notably, we present aligned data from phase 3 of the Drosophila Population Genomics Project (DPGP3), which provides 197 genomes from a single ancestral range population of D. melanogaster (from Zambia). The large sample size, high genetic diversity, and potentially simpler demographic history of the DPGP3 sample will make this a highly valuable resource for fundamental population genetic research. The complete set of assemblies described here, termed the Drosophila Genome Nexus, presently comprises 623 consistently aligned genomes and is publicly available in multiple formats with supporting documentation and bioinformatic tools. This resource will greatly facilitate population genomic analysis in this model species by reducing the methodological differences between data sets.

A Thousand Fly Genomes: An Expanded Drosophila Genome Nexus

The Drosophila Genome Nexus is a population genomic resource that provides D. melanogaster genomes from multiple sources. To facilitate comparisons across data sets, genomes are aligned using a common reference alignment pipeline which involves two rounds of mapping. Regions of residual heterozygosity, identity-by-descent, and recent population admixture are annotated to enable data filtering based on the user’s needs. Here, we present a significant expansion of the Drosophila Genome Nexus, which brings the current data object to a total of 1,121 wild-derived genomes. New additions include 305 previously unpublished genomes from inbred lines representing six population samples in Egypt, Ethiopia, France, and South Africa, along with another 193 genomes added from recently-published data sets. We also provide an aligned D. simulans genome to facilitate divergence comparisons. This improved resource will broaden the range of population genomic questions that can addressed from multi-population allele frequencies and haplotypes in this model species. The larger set of genomes will also enhance the discovery of functionally relevant natural variation that exists within and between populations.

Spatially varying selection shapes life history clines among populations of Drosophila melanogaster from sub-Saharan Africa

Clines in life history traits, presumably driven by spatially varying selection, are widespread. Major latitudinal clines have been observed, for example, in Drosophila melanogaster, an ancestrally tropical insect from Africa that has colonized temperate habitats on multiple continents. Yet, how geographic factors other than latitude, such as altitude or longitude, affect life history in this species remains poorly understood. Moreover, most previous work has been performed on derived European, American and Australian populations, but whether life history also varies predictably with geography in the ancestral Afro‐tropical range has not been investigated systematically. Here, we have examined life history variation among populations of D. melanogaster from sub‐Saharan Africa. Viability and reproductive diapause did not vary with geography, but body size increased with altitude, latitude and longitude. Early fecundity covaried positively with altitude and latitude, whereas lifespan showed the opposite trend. Examination of genetic variance–covariance matrices revealed geographic differentiation also in trade‐off structure, and QST‐FST analysis showed that life history differentiation among populations is likely shaped by selection. Together, our results suggest that geographic and/or climatic factors drive adaptive phenotypic differentiation among ancestral African populations and confirm the widely held notion that latitude and altitude represent parallel gradients.

Decanalization of wing development accompanied the evolution of large wings in high-altitude Drosophila

Significance Developmental buffering mechanisms that stabilize phenotypes against perturbations (such as harmful mutations or environmental stress) have previously been inferred. However, it is unclear whether this “canalization” can be maintained when adaptive evolution causes phenotypes and developmental processes to change. Here, we report a loss of canalization accompanying the evolution of larger wings in high-altitude fruit flies (Drosophila melanogaster). We redeploy a classic genetics technique (mutagenesis) to show that wing development of Ethiopian flies is less robust to new mutations and that large wing size is inherited together with decanalized wing development. These results represent the first example, to our knowledge, of adaptive evolution apparently undermining the developmental buffering of a recently evolved trait. Decanalized development, thus, represents a potential “cost of adaptation.” In higher organisms, the phenotypic impacts of potentially harmful or beneficial mutations are often modulated by complex developmental networks. Stabilizing selection may favor the evolution of developmental canalization—that is, robustness despite perturbation—to insulate development against environmental and genetic variability. In contrast, directional selection acts to alter the developmental process, possibly undermining the molecular mechanisms that buffer a trait’s development, but this scenario has not been shown in nature. Here, we examined the developmental consequences of size increase in highland Ethiopian Drosophila melanogaster. Ethiopian inbred strains exhibited much higher frequencies of wing abnormalities than lowland populations, consistent with an elevated susceptibility to the genetic perturbation of inbreeding. We then used mutagenesis to test whether Ethiopian wing development is, indeed, decanalized. Ethiopian strains were far more susceptible to this genetic disruption of development, yielding 26 times more novel wing abnormalities than lowland strains in F2 males. Wing size and developmental perturbability cosegregated in the offspring of between-population crosses, suggesting that genes conferring size differences had undermined developmental buffering mechanisms. Our findings represent the first observation, to our knowledge, of morphological evolution associated with decanalization in the same tissue, underscoring the sensitivity of development to adaptive change.

A Variable Genetic Architecture of Melanic Evolution in Drosophila melanogaster

Unraveling the genetic architecture of adaptive phenotypic divergence is a fundamental quest in evolutionary biology. In Drosophila melanogaster, high-altitude melanism has evolved in separate mountain ranges in sub-Saharan Africa, potentially as an adaptation to UV intensity. We investigated the genetic basis of this melanism in three populations using a new bulk segregant analysis mapping method. We identified 19 distinct QTL regions from nine mapping crosses, with several QTL peaks overlapping between two or all populations, and yet different crosses involving the same melanic population commonly yielded distinct QTL. The strongest QTL often overlapped well-known pigmentation genes, but we typically did not find wide signals of genetic differentiation (FST) between lightly and darkly pigmented populations at these genes. Instead, we found small numbers of highly differentiated SNPs at the probable causative genes. A simulation analysis showed that these patterns of polymorphism were consistent with selection on standing genetic variation. Overall, our results suggest that, even for potentially simpler traits like pigmentation, the complexity of adaptive trait evolution poses important challenges for QTL mapping and population genetic analysis.

Life history evolution and cellular mechanisms associated with increased size in high-altitude Drosophila.

Abstract Understanding the physiological and genetic basis of growth and body size variation has wide‐ranging implications, from cancer and metabolic disease to the genetics of complex traits. We examined the evolution of body and wing size in high‐altitude Drosophila melanogaster from Ethiopia, flies with larger size than any previously known population. Specifically, we sought to identify life history characteristics and cellular mechanisms that may have facilitated size evolution. We found that the large‐bodied Ethiopian flies laid significantly fewer but larger eggs relative to lowland, smaller‐bodied Zambian flies. The highland flies were found to achieve larger size in a similar developmental period, potentially aided by a reproductive strategy favoring greater provisioning of fewer offspring. At the cellular level, cell proliferation was a strong contributor to wing size evolution, but both thorax and wing size increases involved important changes in cell size. Nuclear size measurements were consistent with elevated somatic ploidy as an important mechanism of body size evolution. We discuss the significance of these results for the genetic basis of evolutionary changes in body and wing size in Ethiopian D. melanogaster.

Bat Molecular Phylogenetics: Past, Present, and Future Directions

With the development of techniques for the isolation, amplification, and sequencing of DNA, studies addressing phylogenetic relationships among bats moved in the 1980s and early 1990s from restriction fragment length polymorphisms and DNA hybridization to examination of changes at the level of individual nucleotides via DNA sequence analysis. Coinciding with these molecular advances were increases in computational capacity and the development of sophisticated analytical models of the nucleotide and amino acid substitution processes. Thus, molecular phylogenetics moved primarily from distance- and parsimony-based algorithms to complex optimality criterion, such as maximum likelihood and Bayesian phylogenetics. These advances helped to clarify our understanding of the evolutionary relationships and biogeographic history of bats as well as the evolution of echolocation and flight. We have now entered the age of phylogenomics, where phylogenetic datasets represent genome-scale variation. With the recent explosion in available data and the ever-expanding ability to inexpensively produce massive datasets for any organism, we now have the ability to extend phylogenomic approaches to non-model organisms such as bats, which have been historically neglected in studies of molecular evolution. Bats represent an extraordinary group characterized by evolutionary novelty (i.e., echolocation and powered flight) and adaptation (e.g., seven distinct feeding strategies within the family Phyllostomidae) to an extent arguably paralleled by no other group of mammals. As we move into the era of phylogenomics, it is time for bats to move to the forefront of the study of evolutionary novelty and adaptation, which for mammals has been dominated by rodents due to the mouse and rat models being so ubiquitous in all aspects of biology.

Parallel Evolution of Cold Tolerance Within Drosophila melanogaster

Drosophila melanogaster originated in tropical Africa before expanding into strikingly different temperate climates in Eurasia and beyond. Here, we find elevated cold tolerance in three distinct geographic regions: beyond the well-studied non-African case, we show that populations from the highlands of Ethiopia and South Africa have significantly increased cold tolerance as well. We observe greater cold tolerance in outbred versus inbred flies, but only in populations with higher inversion frequencies. Each cold-adapted population shows lower inversion frequencies than a closely-related warm-adapted population, suggesting that inversion frequencies may decrease with altitude in addition to latitude. Using the FST-based “Population Branch Excess” statistic (PBE), we found only limited evidence for parallel genetic differentiation at the scale of ∼4 kb windows, specifically between Ethiopian and South African cold-adapted populations. And yet, when we looked for single nucleotide polymorphisms (SNPs) with codirectional frequency change in two or three cold-adapted populations, strong genomic enrichments were observed from all comparisons. These findings could reflect an important role for selection on standing genetic variation leading to “soft sweeps”. One SNP showed sufficient codirectional frequency change in all cold-adapted populations to achieve experiment-wide significance: an intronic variant in the synaptic gene Prosap. Another codirectional outlier SNP, at senseless-2, had a strong association with our cold trait measurements, but in the opposite direction as predicted. More generally, proteins involved in neurotransmission were enriched as potential targets of parallel adaptation. The ability to study cold tolerance evolution in a parallel framework will enhance this classic study system for climate adaptation.

Oligogenic Adaptation, Soft Sweeps, and Parallel Melanic Evolution in Drosophila melanogaster

Unraveling the genetic architecture of adaptive phenotypic divergence is a fundamental quest in evolutionary biology. In Drosophila melanogaster, high-altitude melanism has evolved in separate mountain ranges in sub-Saharan Africa, potentially as an adaptation to UV intensity. We investigated the genetic basis of this melanism in three populations using a new bulk segregant analysis mapping method. Although hundreds of genes are known to affect cuticular pigmentation in D. melanogaster, we identified only 19 distinct QTLs from 9 mapping crosses, with several QTL peaks being shared among two or all populations. Surprisingly, we did not find wide signals of genetic differentiation (Fst) between lightly and darkly pigmented populations at these QTLs, in spite of the pronounced phenotypic difference in pigmentation. Instead, we found small numbers of highly differentiated SNPs at the probable causative genes. A simulation analysis showed that these patterns of polymorphism are consistent with selection on standing genetic variation (leading to “soft sweeps“). Our results thus support a role for oligogenic selection on standing genetic variation in driving parallel ecological adaptation.