Divergence and introgression among the virilis group of Drosophila

Speciation with gene flow is now widely regarded as common. However, the frequency of introgression between recently diverged species and the evolutionary consequences of gene flow are still poorly understood. The virilis group of Drosophila contains around a dozen species that are geographically widespread and show varying levels of pre-zygotic and post-zygotic isolation. Here, we utilize de novo genome assemblies and whole-genome sequencing data to resolve phylogenetic relationships and describe patterns of introgression and divergence across the group. We suggest that the virilis group consists of three, rather than the traditional two, subgroups. We found evidence of pervasive phylogenetic discordance caused by ancient introgression events between distant lineages within the group, and much more recent gene flow between closely-related species. When assessing patterns of genome-wide divergence in species pairs across the group, we found no consistent genomic evidence of a disproportionate role for the X chromosome. Some genes undergoing rapid sequence divergence across the group were involved in chemical communication and may be related to the evolution of sexual isolation. We suggest that gene flow between closely-related species has potentially had an impact on lineage-specific adaptation and the evolution of reproductive barriers. Our results show how ancient and recent introgression confuse phylogenetic reconstruction, and suggest that shared variation can facilitate adaptation and speciation.

Evolutionary Dynamics of the Pericentromeric Heterochromatin in Drosophila virilis and Related Species

Pericentromeric heterochromatin in Drosophila generally consists of repetitive DNA, forming the environment associated with gene silencing. Despite the expanding knowledge of the impact of transposable elements (TEs) on the host genome, little is known about the evolution of pericentromeric heterochromatin, its structural composition, and age. During the evolution of the Drosophilidae, hundreds of genes have become embedded within pericentromeric regions yet retained activity. We investigated a pericentromeric heterochromatin fragment found in D. virilis and related species, describing the evolution of genes in this region and the age of TE invasion. Regardless of the heterochromatic environment, the amino acid composition of the genes is under purifying selection. However, the selective pressure affects parts of genes in varying degrees, resulting in expansion of gene introns due to TEs invasion. According to the divergence of TEs, the pericentromeric heterochromatin of the species of virilis group began to form more than 20 million years ago by invasions of retroelements, miniature inverted repeat transposable elements (MITEs), and Helitrons. Importantly, invasions into the heterochromatin continue to occur by TEs that fall under the scope of piRNA silencing. Thus, the pericentromeric heterochromatin, in spite of its ability to induce silencing, has the means for being dynamic, incorporating the regions of active transcription.

ebony affects pigmentation divergence and cuticular hydrocarbons in Drosophila americana and D. novamexicana

1AbstractDrosophila pigmentation has been a fruitful model system for understanding the genetic and developmental mechanisms underlying phenotypic evolution. For example, prior work has shown that divergence of the tan gene contributes to pigmentation differences between two members of the virilis group: Drosophila novamexicana, which has a light yellow body color, and D. americana, which has a dark brown body color. Quantitative trait locus (QTL) mapping and expression analysis has suggested that divergence of the ebony gene might also contribute to pigmentation differences between these two species. Here, we directly test this hypothesis by using CRISPR/Cas9 genome editing to generate ebony null mutants in D. americana and D. novamexicana and then using reciprocal hemizygosity testing to compare the effects of each species’ ebony allele on pigmentation. We find that divergence of ebony does indeed contribute to the pigmentation divergence between species, with effects on both the overall body color as well as a difference in pigmentation along the dorsal abdominal midline. Motivated by recent work in D. melanogaster, we also used the ebony null mutants to test for effects of ebony on cuticular hydrocarbon (CHC) profiles. We found that ebony affects CHC abundance in both species, but does not contribute to qualitative differences in the CHC profiles between these two species. Additional transgenic resources for working with D. americana and D. novamexicana, such as white mutants of both species and yellow mutants in D. novamexicana, were generated in the course of this work and are also described. Taken together, this study advances our understanding of loci contributing to phenotypic divergence and illustrates how the latest genome editing tools can be used for functional testing in non-model species.

Genetic variation of the nuclear sequences of mitochondrial origin associated with retrotransposon Tv1 insertions in Drosophila species of the virilis group

Mitochondrial DNA sequences integrated into chromosomes are a promising object for designing genetic markers for studies of phylogenesis and genomic instability. Mitochondrial genomes of D. virilis and other Drosophila species of the virilis group contain (AT)nmicrosatellites in the spacer region between the atp6 and cox3 genes, and this microsatellite sequence is one of the hallmarks of the virilis group. The nuclear genome of D. virilis contains many extended fragments of mitochondrial DNA, which in total are several times longer than the mitochondrial genome. These nuclear sequences of mitochondrial origin contain all types of mitochondrial sequences, including mitochondrial genes and the aforementioned microsatellite sequence. The presence of the (AT)nmicrosatellite allows insertion of retrotransposon Tv1, which can transpose into the (AT)n microsatellite in a site-specific manner. The Tv1 insertion into (AT)n, close to the atp6 or cox3 pseudogenes produces a unique sequence. This sequence is formed by retrotransposon Tv1 and pseudogenes atp6 or cox3. This unique sequence can be detected in the genome by a PCR-based method. We applied this method to the detection and analysis of the nucleotide variability of the pseudogenes atp6 and cox3 associated with Tv1 insertions in a D. virilis cell culture and in the genomes of four Drosophila species of the virilis group: D. virilis, D. montana, D. borealis, and D. lacicola. We discovered new events of mitochondrial sequence transfer to the nucleus in the transplanted cell culture of D. virilis, and new Tv1 insertions, having emerged during the passage of this cell line were detected in the genome of the D. virilis transplanted cell culture. We found atp6 and cox3 pseudogenes associated with insertions of retrotransposon Tv1 in the nuclear genomes of four Drosophila species from the virilis group. These chimeric sequences proved to be species-specific. The age of the Tv1 insertion into the atp6 and cox3 pseudogenes is estimated at 1.50 Ma for D. virilis, 1.31 Ma for D. lacicola, and 1.56 Ma for D. borealis. A specific situation was revealed for D. montana, in which Tv1 insertions with nearly identical 5' and 3' long terminal repeats (LTRs) were present in accessions of flies from Europe and Asia. The age of this insertion was about 300 thousand years, and the insertion was absent from the D. montana fly line from North America.

Evolutionary dynamics of male reproductive genes in the Drosophila virilis subgroup

Postcopulatory sexual selection (PCSS) is a potent evolutionary force that can drive rapid changes of reproductive genes within species, and thus has the potential to generate reproductive incompatibilities between species. Male seminal fluid proteins (SFPs) are major players in postmating interactions, and likely the main targets of PCSS in males. The virilis subgroup of Drosophila exhibits strong interspecific gametic incompatibilities, and can serve as a model to study the genetic basis of PCSS and gametic isolation. However, reproductive genes in this group have not been characterized. Here we use short-read RNA sequencing of male reproductive organs to examine the evolutionary dynamics of reproductive genes in members of the virilis subgroup: D. americana, D. lummei, D. novamexicana, and D. virilis. For each of the three male reproductive organs (accessory glands, ejaculatory bulb, and testes), we identify genes that show strong expression bias in a given tissue relative to the remaining tissues. We find that the majority of male reproductive transcripts are testes-biased, accounting for ~15% of all annotated genes. Ejaculatory bulb-biased transcripts largely code for lipid metabolic enzymes, and contain orthologs of the D. melanogaster ejaculatory bulb protein, Peb-me, which is involved in mating-plug formation. In addition, we identify 71 candidate SFPs, and show that this set of genes has the highest rate of nonsynonymous codon substitution relative to testes- and ejaculatory bulb-biased genes. Furthermore, these SFPs are underrepresented on the X chromosome and are enriched for proteolytic enzymes, which is consistent with SFPs in other insect species. Surprisingly, we find 35 D. melanogaster SFPs with conserved accessory gland expression in the virilis group, suggesting these genes may have conserved reproductive roles in Drosophila. Finally, we show that several of the SFPs that have the highest rate of nonsynonymous codon substitutions reside on the centromeric half of chromosome 2, which contributes to paternal gametic incompatibility between species. Our results suggest that SFPs are under strong selection in the virilis group, and likely play a major role in PCSS and/or gametic isolation.

Multiple Genes Cause Postmating Prezygotic Reproductive Isolation in the Drosophila virilis Group

Understanding the genetic basis of speciation is a central problem in evolutionary biology. Studies of reproductive isolation have provided several insights into the genetic causes of speciation, especially in taxa that lend themselves to detailed genetic scrutiny. Reproductive barriers have usually been divided into those that occur before zygote formation (prezygotic) and after (postzygotic), with the latter receiving a great deal of attention over several decades. Reproductive barriers that occur after mating but before zygote formation [postmating prezygotic (PMPZ)] are especially understudied at the genetic level. Here, I present a phenotypic and genetic analysis of a PMPZ reproductive barrier between two species of the Drosophila virilis group: D. americana and D. virilis. This species pair shows strong PMPZ isolation, especially when D. americana males mate with D. virilis females: ∼99% of eggs laid after these heterospecific copulations are not fertilized. Previous work has shown that the paternal loci contributing to this incompatibility reside on two chromosomes, one of which (chromosome 5) likely carries multiple factors. The other (chromosome 2) is fixed for a paracentric inversion that encompasses nearly half the chromosome. Here, I present two results. First, I show that PMPZ in this species cross is largely due to defective sperm storage in heterospecific copulations. Second, using advanced intercross and backcross mapping approaches, I identify genomic regions that carry genes capable of rescuing heterospecific fertilization. I conclude that paternal incompatibility between D. americana males and D. virilis females is underlain by four or more genes on chromosomes 2 and 5.

A Single Gene Causes an Interspecific Difference in Pigmentation in Drosophila

The genetic basis of species differences remains understudied. Studies in insects have contributed significantly to our understanding of morphological evolution. Pigmentation traits in particular have received a great deal of attention and several genes in the insect pigmentation pathway have been implicated in inter- and intraspecific differences. Nonetheless, much remains unknown about many of the genes in this pathway and their potential role in understudied taxa. Here we genetically analyze the puparium color difference between members of the Virilis group of Drosophila. The puparium of Drosophila virilis is black, while those of D. americana, D. novamexicana, and D. lummei are brown. We used a series of backcross hybrid populations between D. americana and D. virilis to map the genomic interval responsible for the difference between this species pair. First, we show that the pupal case color difference is caused by a single Mendelizing factor, which we ultimately map to an ~11kb region on chromosome 5. The mapped interval includes only the first exon and regulatory region(s) of the dopamine N-acetyltransferase gene (Dat). This gene encodes an enzyme that is known to play a part in the insect pigmentation pathway. Second, we show that this gene is highly expressed at the onset of pupation in light-brown taxa (D. americana and D. novamexicana) relative to D. virilis, but not in the dark-brown D. lummei. Finally, we examine the role of Dat in adult pigmentation between D. americana (heavily melanized) and D. novamexicana (lightly melanized) and find no discernible effect of this gene in adults. Our results demonstrate that a single gene is entirely or almost entirely responsible for a morphological difference between species.