Genomic differentiation between temperate and tropical Australian populations of Drosophila melanogaster

Determining the genetic basis of environmental adaptation is a central problem of evolutionary biology. This issue has been fruitfully addressed by examining genetic differentiation between populations that are recently separated and/or experience high rates of gene flow. A good example of this approach is the decades-long investigation of selection acting along latitudinal clines in Drosophila melanogaster. Here we use next-generation genome sequencing to reexamine the well-studied Australian D. melanogaster cline. We find evidence for extensive differentiation between temperate and tropical populations, with regulatory regions and unannotated regions showing particularly high levels of differentiation. Although the physical genomic scale of geographic differentiation is small--on the order of gene sized--we observed several larger highly differentiated regions. The region spanned by the cosmopolitan inversion polymorphism In(3R)P shows higher levels of differentiation, consistent with the major difference in allele frequencies of Standard and In(3R)P karyotypes in temperate vs. tropical Australian populations. Our analysis reveals evidence for spatially varying selection on a number of key biological processes, suggesting fundamental biological differences between flies from these two geographic regions.     

Genomic imprinting and position-effect variegation in Drosophila melanogaster

Genomic imprinting is a phenomenon in which the expression of a gene or chromosomal region depends on the sex of the individual transmitting it. The term imprinting was first coined to describe parent-specific chromosome behavior in the dipteran insect Sciara and has since been described in many organisms, including other insects, plants, fish, and mammals. In this article we describe a mini-X chromosome in Drosophila melanogaster that shows genomic imprinting of at least three closely linked genes. The imprinting of these genes is observed as mosaic silencing when the genes are transmitted by the male parent, in contrast to essentially wild-type expression when the same genes are maternally transmitted. We show that the imprint is due to the sex of the parent rather than to a conventional maternal effect, differential mitotic instability of the mini-X chromosome, or an allele-specific effect. Finally, we have examined the effects of classical modifiers of position-effect variegation on the maintenance and the establishment of the imprint. Factors that modify position-effect variegation alter the somatic expression but not the establishment of the imprint. This suggests that chromatin structure is important in maintenance of the imprint, but a separate mechanism may be responsible for its initiation.     

Genomic imprinting absent in Drosophila melanogaster adult females

Genomic imprinting occurs when expression of an allele differs based on the sex of the parent that transmitted the allele. In D. melanogaster, imprinting can occur, but its impact on allelic expression genome-wide is unclear. Here, we search for imprinted genes in D. melanogaster using RNA-seq to compare allele-specific expression between pools of 7- to 10-day-old adult female progeny from reciprocal crosses. We identified 119 genes with allelic expression consistent with imprinting, and these genes showed significant clustering within the genome. Surprisingly, additional analysis of several of these genes showed that either genomic heterogeneity or high levels of intrinsic noise caused imprinting-like allelic expression. Consequently, our data provide no convincing evidence of imprinting for D. melanogaster genes in their native genomic context. Elucidating sources of false-positive signals for imprinting in allele-specific RNA-seq data, as done here, is critical given the growing popularity of this method for identifying imprinted genes.     

Genomic heterogeneity of background substitutional patterns in Drosophila melanogaster

Mutation is the underlying force that provides the variation upon which evolutionary forces can act. It is important to understand how mutation rates vary within genomes and how the probabilities of fixation of new mutations vary as well. If substitutional processes across the genome are heterogeneous, then examining patterns of coding sequence evolution without taking these underlying variations into account may be misleading. Here we present the first rigorous test of substitution rate heterogeneity in the Drosophila melanogaster genome using almost 1500 nonfunctional fragments of the transposable element DNAREP1_DM. Not only do our analyses suggest that substitutional patterns in heterochromatic and euchromatic sequences are different, but also they provide support in favor of a recombination-associated substitutional bias toward G and C in this species. The magnitude of this bias is entirely sufficient to explain recombination-associated patterns of codon usage on the autosomes of the D. melanogaster genome. We also document a bias toward lower GC content in the pattern of small insertions and deletions (indels). In addition, the GC content of noncoding DNA in Drosophila is higher than would be predicted on the basis of the pattern of nucleotide substitutions and small indels. However, we argue that the fast turnover of noncoding sequences in Drosophila makes it difficult to assess the importance of the GC biases in nucleotide substitutions and small indels in shaping the base composition of noncoding sequences.     

Genomic organization of gypsy chromatin insulators in Drosophila melanogaster

Chromatin insulators have been implicated in the regulation of higher-order chromatin structure and may function to compartmentalize the eukaryotic genome into independent domains of gene expression. To test this possibility, we used biochemical and computational approaches to identify gypsy-like genomic-binding sites for the Suppressor of Hairy-wing [Su(Hw)] protein, a component of the gypsy insulator. EMSA and FISH analyses suggest that these are genuine Su(Hw)-binding sites. In addition, functional tests indicate that genomic Su(Hw)-binding sites can inhibit enhancer-promoter interactions and thus function as bona fide insulators. The insulator strength is dependent on the genomic location of the transgene and the number of Su(Hw)-binding sites, with clusters of two to three sites showing a stronger effect than individual sites. These clusters of Su(Hw)-binding sites are located mostly in intergenic regions or in introns of large genes, an arrangement that fits well with their proposed role in the formation of chromatin domains. Taken together, these data suggest that genomic gypsy-like insulators may provide a means for the compartmentalization of the genome within the nucleus.     

Genomic signatures of experimental adaptive radiation in Drosophila

Abiotic environmental factors play a fundamental role in determining the distribution, abundance and adaptive diversification of species. Empowered by new technologies enabling rapid and increasingly accurate examination of genomic variation in populations, researchers may gain new insights into the genomic background of adaptive radiation and stress resistance. We investigated genomic variation across generations of large‐scale experimental selection regimes originating from a single founder population of Drosophila melanogaster, diverging in response to ecologically relevant environmental stressors: heat shock, heat knock down, cold shock, desiccation and starvation. When compared to the founder population, and to parallel unselected controls, there were more than 100,000 single nucleotide polymorphisms (SNPs) displaying consistent allelic changes in response to selective pressures across generations. These SNPs were found in both coding and noncoding sequences, with the highest density in promoter regions, and involved a broad range of functionalities, including molecular chaperoning by heat‐shock proteins. The SNP patterns were highly stressor‐specific despite considerable variation among line replicates within each selection regime, as reflected by a principal component analysis, and co‐occurred with selective sweep regions. Only ~15% of SNPs with putatively adaptive changes were shared by at least two selective regimes, while less than 1% of SNPs diverged in opposite directions. Divergent stressors driving evolution in the experimental system of adaptive radiation left distinct genomic signatures, most pronounced in starvation and heat‐shock selection regimes.     

An analysis of genetic differentiation among assortatively mating Drosophila melanogaster in Zimbabwe

African Drosophila melanogaster populations, and those from Zimbabwe in particular, have attracted much interest recently. African flies differ genetically from 'cosmopolitan' populations and were found to exhibit discriminative mating behaviour against individuals from 'cosmopolitan' populations. It has therefore been proposed that Zimbabwean and some other African populations are in an 'incipient stage of speciation'. However, whether the mating behaviour is an effective barrier against gene flow from other populations, and whether intra-population genetic differentiation has already evolved in sympatry is not known. Here, we took a population-based approach to test whether the well-characterized mating type differences have resulted in a genome-wide differentiation at the population level. Using 122 polymorphic microsatellite loci mapping to the third chromosome, we demonstrate a significant genetic differentiation between Zimbabwean flies differing in their mating behaviour. We also provide evidence to suggest that this difference is unlikely to be attributable to population structure within Zimbabwe. However, the analysis of individual microsatellite loci did not indicate more loci differentiating these two groups than expected by chance. Our data suggest that the 'Z'-'M' mating behaviour is strong enough to result in a small but significant genetic differentiation. Thus, future studies based on a larger population sample of flies characterized for their mating behaviour and using more markers are expected to provide more information on the genetic basis of the mating traits in the Zimbabwe flies.     

Genomic deletions of the Drosophila melanogaster Hsp70 genes

Homologous recombination can produce directed mutations in the genomes of a number of model organisms, including Drosophila melanogaster. One of the most useful applications has been to delete target genes to generate null alleles. In Drosophila, specific gene deletions have not yet been produced by this method. To test whether such deletions could be produced by homologous recombination in D. melanogaster we set out to delete the Hsp70 genes. Six nearly identical copies of this gene, encoding the major heat-shock protein in Drosophila, are found at two separate but closely linked loci. This arrangement has thwarted standard genetic approaches to generate an Hsp70-null fly, making this an ideal test of gene targeting. In this study, ends-out targeting was used to generate specific deletions of all Hsp70 genes, including one deletion that spanned approximately 47 kb. The Hsp70-null flies are viable and fertile. The results show that genomic deletions of varied sizes can be readily generated by homologous recombination in Drosophila.     

Genomic distribution of copia-like elements in laboratory stocks of Drosophila melanogaster

The genomic location of five copia-like transposable elements has been compared by the Southern technique in laboratory lines of Oregon R and Canton S strains of Drosophila melanogaster. The results show that extensive rearrangements have taken place in the few decades of separation between the stocks and suggest that transposition occurs at a sizable rate.     

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.     

Rates and Genomic Consequences of Spontaneous Mutational Events in Drosophila melanogaster

Because spontaneous mutation is the source of all genetic diversity, measuring mutation rates can reveal how natural selection drives patterns of variation within and between species. We sequenced eight genomes produced by a mutation-accumulation experiment in Drosophila melanogaster. Our analysis reveals that point mutation and small indel rates vary significantly between the two different genetic backgrounds examined. We also find evidence that ∼2% of mutational events affect multiple closely spaced nucleotides. Unlike previous similar experiments, we were able to estimate genome-wide rates of large deletions and tandem duplications. These results suggest that, at least in inbred lines like those examined here, mutational pressures may result in net growth rather than contraction of the Drosophila genome. By comparing our mutation rate estimates to polymorphism data, we are able to estimate the fraction of new mutations that are eliminated by purifying selection. These results suggest that ∼99% of duplications and deletions are deleterious—making them 10 times more likely to be removed by selection than nonsynonymous mutations. Our results illuminate not only the rates of new small- and large-scale mutations, but also the selective forces that they encounter once they arise.     

Genomic and structural organization of Drosophila melanogaster G elements.

The properties and the genomic organization of G elements, a moderately repeated DNA family of D. melanogaster, are reported. G elements lack terminal repeats, generate target site duplications at the point of insertion and exhibit at one end a stretch of A residues of variable length. In a large number of recombinant clones analyzed G elements occur in tandem arrays, interspersed with specific ribosomal DNA (rDNA) segments. This arrangement results from the insertion of members of the G family within the nontranscribed spacer (NTS) of rDNA units. Similarity of the site of integration of G elements to that of ribosomal DNA insertions suggests that distinct DNA sequences might have been inserted into rDNA through a partly common pathway.     

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.     

The adaptive significance of chromosomal inversion polymorphisms in Drosophila melanogaster

Chromosomal inversions, structural mutations that reverse a segment of a chromosome, cause suppression of recombination in the heterozygous state. Several studies have shown that inversion polymorphisms can form clines or fluctuate predictably in frequency over seasonal time spans. These observations prompted the hypothesis that chromosomal rearrangements might be subject to spatially and/or temporally varying selection. Here, we review what has been learned about the adaptive significance of inversion polymorphisms in the vinegar fly Drosophila melanogaster, the species in which they were first discovered by Sturtevant in 1917. A large body of work provides compelling evidence that several inversions in this system are adaptive; however, the precise selective mechanisms that maintain them polymorphic in natural populations remain poorly understood. Recent advances in population genomics, modelling and functional genetics promise to greatly improve our understanding of this long‐standing and fundamental problem in the near future.     

The recent demographic and adaptive history of Drosophila melanogaster

Population genetic analyses of the past two decades confirmed an earlier hypothesis by L Tsacas and D Lachaise that the cosmopolitan species Drosophila melanogaster has an Afrotropical origin, and that it colonized the rest of the world only very recently. Maximum likelihood analyses based on multilocus data suggest that the putative ancestral African population expanded its size about 60,000 years ago (ya). These demographic changes were accompanied by the fixation of numerous beneficial mutations, as revealed by signatures of positive directional selection in the genome (selective sweeps). The estimated rate of adaptive substitution on the X chromosome is in the order of 10(-11) per nucleotide site per generation. Comparable (but not significantly higher) substitution rates are found in derived populations that colonized new habitats outside Africa, such as in a European population that branched off from the African lineage about 16,000 ya.     

Genomic analysis of Drosophila melanogaster telomeres: full-length copies of HeT-A and TART elements at telomeres

The repetitive nature of heterochromatin hampers its analysis in general genome-sequencing projects. Specific studies are needed to extend the sequence into telomeric and centromeric heterochromatin. Drosophila telomeres lack the telomerase-generated repeats that are characteristic of other eukaryotic chromosomes. Instead, they consist of tandem arrays of HeT-A and TART elements. Herein, we present the genomic organization of the telomeres in the isogenic strain (y; cn bw sp) that was used for the Drosophila melanogaster sequencing project. The data indicate that the canonical features of telomere organization are widely conserved in evolution. In addition, we have identified full-length elements, likely competent elements, for HeT-A and TART.