Inferring the demographic history and rate of adaptive substitution in Drosophila

An important goal of population genetics is to determine the forces that have shaped the pattern of genetic variation in natural populations. We developed a maximum likelihood method that allows us to infer demographic changes and detect recent positive selection (selective sweeps) in populations of varying size from DNA polymorphism data. Applying this approach to single nucleotide polymorphism data at more than 250 noncoding loci on the X chromosome of Drosophila melanogaster from an (ancestral) African population and a (derived) European, we found that the African population expanded about 60,000 y ago and that the European population split off from the African lineage about 15,800 y ago, thereby suffering a severe population size bottleneck. We estimated that about 160 beneficial mutations (with selection coefficients s between 0.05% and 0.5%) were fixed in the euchromatic portion of the X in the African population since population size expansion, and about 60 mutations (with s around 0.5%) in the diverging European lineage.     

The genetic basis of adaptive pigmentation variation in Drosophila melanogaster

In a broad survey of Drosophila melanogaster population samples, levels of abdominal pigmentation were found to be highly variable and geographically differentiated. A strong positive correlation was found between dark pigmentation and high altitude, suggesting adaptation to specific environments. DNA sequence polymorphism at the candidate gene ebony revealed a clear association with the pigmentation of homozygous third chromosome lines. The darkest lines sequenced had nearly identical haplotypes spanning 14.5 kb upstream of the protein-coding exons of ebony. Thus, natural selection may have elevated the frequency of an allele that confers dark abdominal pigmentation by influencing the regulation of ebony.     

Selection, recombination and demographic history in Drosophila miranda

Selection, recombination, and the demographic history of a species can all have profound effects on genomewide patterns of variability. To assess the impact of these forces in the genome of Drosophila miranda, we examine polymorphism and divergence patterns at 62 loci scattered across the genome. In accordance with recent findings in D. melanogaster, we find that noncoding DNA generally evolves more slowly than synonymous sites, that the distribution of polymorphism frequencies in noncoding DNA is significantly skewed toward rare variants relative to synonymous sites, and that long introns evolve significantly slower than short introns or synonymous sites. These observations suggest that most noncoding DNA is functionally constrained and evolving under purifying selection. However, in contrast to findings in the D. melanogaster species group, we find little evidence of adaptive evolution acting on either coding or noncoding sequences in D. miranda. Levels of linkage disequilibrium (LD) in D. miranda are comparable to those observed in D. melanogaster, but vary considerably among chromosomes. These patterns suggest a significantly lower rate of recombination on autosomes, possibly due to the presence of polymorphic autosomal inversions and/or differences in chromosome sizes. All chromosomes show significant departures from the standard neutral model, including too much heterogeneity in synonymous site polymorphism relative to divergence among loci and a general excess of rare synonymous polymorphisms. These departures from neutral equilibrium expectations are discussed in the context of nonequilibrium models of demography and selection.     

Genomic analysis of adaptive differentiation in Drosophila melanogaster

Drosophila melanogaster shows clinal variation along latitudinal transects on multiple continents for several phenotypes, allozyme variants, sequence variants, and chromosome inversions. Previous investigation suggests that many such clines are due to spatially varying selection rather than demographic history, but the genomic extent of such selection is unknown. To map differentiation throughout the genome, we hybridized DNA from temperate and subtropical populations to Affymetrix tiling arrays. The dense genomic sampling of variants and low level of linkage disequilibrium in D. melanogaster enabled identification of many small, differentiated regions. Many regions are differentiated in parallel in the United States and Australia, strongly supporting the idea that they are influenced by spatially varying selection. Genomic differentiation is distributed nonrandomly with respect to gene function, even in regions differentiated on only one continent, providing further evidence for the role of selection. These data provide candidate genes for phenotypes known to vary clinally and implicate interesting new processes in genotype-by-environment interactions, including chorion proteins, proteins regulating meiotic recombination and segregation, gustatory and olfactory receptors, and proteins affecting synaptic function and behavior. This portrait of differentiation provides a genomic perspective on adaptation and the maintenance of variation through spatially varying selection.     

Indirect thermal selection in Drosophila melanogaster and adaptive consequences

Summary Short-term indirect selection in Drosophila melanogaster for heat-sensitivity and heat resistance resulted in two strains, one heat sensitive and another heat resistant, and correlated responses were found for the rate of heat shock protein synthesis, behavioral patterns (asymmetrical sexual isolation) and fitness components (fecundity, fertility, viability, developmental time), as well as for several enzyme activities (MDH, G-6-PDH, ADH, ACHE). These responses associated with temperature selection may reflect the effects of differential inbreeding depression caused by homozygosity of temperature sensitive mutations with different pleiotropic effects. Selection even of a very short duration can induce significant adaptive and evolutionary changes.     

Human demographic history: refining the recent African origin model

Recent studies of large portions of the human genome support a recent origin of modern humans from an African stock after a bottleneck of moderate size followed by a range expansion out of Africa. Under this simple scenario, patterns of molecular diversity suggest that balancing selection could be more prevalent than positive selection in coding regions.     

The recent spread of a vertically transmitted virus through populations of Drosophila melanogaster

The sigma virus is a vertically transmitted pathogen that commonly infects natural populations of Drosophila melanogaster. This virus is the only known host-specific pathogen of D. melanogaster, and so offers a unique opportunity to study the genetics of Drosophila-viral interactions in a natural system. To elucidate the population genetic processes that operate in sigma virus populations, we collected D. melanogaster from 10 populations across three continents. We found that the sigma virus had a prevalence of 0-15% in these populations. Compared to other RNA viruses, we found that levels of viral genetic diversity are very low across Europe and North America. Based on laboratory measurements of the viral substitution rate, we estimate that most European and North American viral isolates shared a common ancestor approximately 200 years ago. We suggest two explanations for this: the first is that D. melanogaster has recently acquired the sigma virus; the second is that a single viral type has recently swept through D. melanogaster populations. Furthermore, in contrast to Drosophila populations, we find that the sigma viral populations are highly structured. This is surprising for a vertically transmitted pathogen that has a similar migration rate to its host. We suggest that the low structure in the viral populations can be explained by the smaller effective population size of the virus.     

High rate of recent transposable element-induced adaptation in Drosophila melanogaster

Although transposable elements (TEs) are known to be potent sources of mutation, their contribution to the generation of recent adaptive changes has never been systematically assessed. In this work, we conduct a genome-wide screen for adaptive TE insertions in Drosophila melanogaster that have taken place during or after the spread of this species out of Africa. We determine population frequencies of 902 of the 1,572 TEs in Release 3 of the D. melanogaster genome and identify a set of 13 putatively adaptive TEs. These 13 TEs increased in population frequency sharply after the spread out of Africa. We argue that many of these TEs are in fact adaptive by demonstrating that the regions flanking five of these TEs display signatures of partial selective sweeps. Furthermore, we show that eight out of the 13 putatively adaptive elements show population frequency heterogeneity consistent with these elements playing a role in adaptation to temperate climates. We conclude that TEs have contributed considerably to recent adaptive evolution (one TE-induced adaptation every 200-1,250 y). The majority of these adaptive insertions are likely to be involved in regulatory changes. Our results also suggest that TE-induced adaptations arise more often from standing variants than from new mutations. Such a high rate of TE-induced adaptation is inconsistent with the number of fixed TEs in the D. melanogaster genome, and we discuss possible explanations for this discrepancy.     

Effects of recent population bottlenecks on reconstructing the demographic history of prairie-chickens

Current methods of DNA sequence analysis attempt to reconstruct historical patterns of population structure and growth from contemporary samples. However, these techniques may be influenced by recent population bottlenecks, which have the potential to eliminate lineages that reveal past changes in demography. One way to examine the performance of these demographic methods is to compare samples from populations before and after recent bottlenecks. We compared estimates of demographic history from populations of greater prairie-chickens (Tympanuchus cupido) before and after recent bottlenecks using four common methods (nested clade analysis [NCA], Tajima's D, mismatch distribution, and MDIV). We found that NCA did not perform well in the presence of bottleneck events, although it did recover some genetic signals associated with increased isolation and the extinction of intermediate populations. The majority of estimates for Tajima's D, including those from bottlenecked populations, were not significantly different from zero, suggesting our data conformed to neutral expectations. In contrast, mismatch distributions including the raggedness index were more likely to identify recently bottlenecked populations with this data set. Estimates of population mutation rate (theta), population divergence time (t), and time to the most recent common ancestor (TMRCA) from MDIV were similar before and after bottlenecks; however, estimates of gene flow (M) were significantly lower in a few cases following a bottleneck. These results suggest that caution should be used when assessing demographic history from contemporary data sets, as recently fragmented and bottlenecked populations may have lost lineages that affect inferences of their demographic history.     

Genetic differentiation and demographic history in Drosophila pachea from the Sonoran Desert

Genetic variation at six microsatellite DNA loci and a segment of the mitochondrial cytochrome oxidase subunit I (COI) locus was used to estimate gene flow, population structure, and demographic history in the cactophilic Drosophila pachea from the Sonoran Desert of North America, a species that shows a strict association with its senita host cactus (genus Lophocereus). For microsatellite analyses, thirteen populations of D. pachea were sampled, five in mainland Mexico and the southwestern USA, and eight on the Baja California (Baja) peninsula, covering essentially the entire range of the species. Analysis of molecular variance (AMOVA) of microsatellite data revealed that populations from both the mainland and the Baja peninsula generally showed little structure, although there were a few exceptions, suggesting some local differentiation and restriction of gene flow within both regions. Pairwise comparisons of F(ST) among each of the mainland and Baja populations showed evidence of both panmixia and population subdivision. AMOVA performed on grouped populations from both the mainland and Baja, however, revealed significant partitioning of genetic variation among the two regions, but no partitioning among localities within each region. Bayesian skyline analyses of the COI data set, consisting of four mainland and seven peninsular populations, revealed population expansions dating to the Pleistocene or late Pliocene in D. pachea from both regions, although regional differences were seen in the estimated timing of the expansions and in changes in effective population size over time.     

Effects of assay conditions in life history experiments with Drosophila melanogaster

Selection experiments with Drosophila have revealed constraints on the simultaneous evolution of life history traits. However, the responses to selection reported by different research groups have not been consistent. Two possible reasons for these inconsistencies are (i) that different groups used different environments for their experiments and (ii) that the selection environments were not identical to the assay environments in which the life history traits were measured. We tested for the effect of the assay environment in life history experiments by measuring a set of Drosophila selection lines in laboratories working on life history evolution with Drosophila in Basel, Groningen, Irvine and London. The lines measured came from selection experiments from each of these laboratories. In each assay environment, we measured fecundity, longevity, development time and body size. The results show that fecundity measurements were particularly sensitive to the assay environment. Differences between assay and selection environment in the same laboratory or differences between assay environments between laboratories could have contributed to the differences in the published results. The other traits measured were less sensitive to the assay environment. However, for all traits there were cases where the measurements in one laboratory suggested that selection had an effect on the trait, whereas in other laboratories no such conclusion would have been drawn. Moreover, we provide good evidence for local adaptation in early fecundity for lines from two laboratories.     

The Drosophila melanogaster Cajal body

Cajal bodies (CBs) are nuclear organelles that are usually identified by the marker protein p80-coilin. Because no orthologue of coilin is known in Drosophila melanogaster, we identified D. melanogaster CBs using probes for other components that are relatively diagnostic for CBs in vertebrate cells. U85 small CB-specific RNA, U2 small nuclear RNA, the survival of motor neurons protein, and fibrillarin occur together in a nuclear body that is closely associated with the nucleolus. Based on its similarity to CBs in other organisms, we refer to this structure as the D. melanogaster CB. Surprisingly, the D. melanogaster U7 small nuclear RNP resides in a separate nuclear body, which we call the histone locus body (HLB). The HLB is invariably colocalized with the histone gene locus. Thus, canonical CB components are distributed into at least two nuclear bodies in D. melanogaster. The identification of these nuclear bodies now permits a broad range of questions to be asked about CB structure and function in a genetically tractable organism.     

The translation factors of Drosophila melanogaster

Synthesis of polypeptides from mRNA (translation) is a fundamental cellular process that is coordinated and catalyzed by a set of canonical ‘translation factors’. Surprisingly, the translation factors of Drosophila melanogaster have not yet been systematically identified, leading to inconsistencies in their nomenclature and shortcomings in functional (Gene Ontology, GO) annotations. Here, we describe the complete set of translation factors in D. melanogaster, applying nomenclature already in widespread use in other species, and revising their functional annotation. The collection comprises 43 initiation factors, 12 elongation factors, 3 release factors and 6 recycling factors, totaling 64 of which 55 are cytoplasmic and 9 are mitochondrial. We also provide an overview of notable findings and particular insights derived from Drosophila about these factors. This catalog, together with the incorporation of the improved nomenclature and GO annotation into FlyBase, will greatly facilitate access to information about the functional roles of these important proteins.     

The generation of cloned Drosophila melanogaster

We report here the first successful use of embryonic nuclear transfer to create viable adult Drosophila melanogaster clones. Given the generation time, cost effectiveness, and relative ease of embryonic nuclear transplant in Drosophila, this method can provide an opportunity to further study the constraints on development imposed by transplanting determined or differentiated nuclei.