Patterns of Genomic Differentiation in the Drosophila nasuta Species Complex

Metazoa-level USCOs as markers in species delimitation and classification

Metazoa-level universal single-copy orthologs (mzl-USCOs) are universally applicable markers for DNA taxonomy in animals that can replace or supplement single-gene barcodes. Previously, mzl-USCOs from target enrichment data were shown to reliably distinguish species. Here, we tested whether USCOs are an evenly distributed, representative sample of a given metazoan genome and therefore able to cope with past hybridization events and incomplete lineage sorting. This is relevant for coalescent-based species delimitation approaches, which critically depend on the assumption that the investigated loci do not exhibit autocorrelation due to physical linkage. Based on 239 chromosome-level assembled genomes, we confirmed that mzl-USCOs are genetically unlinked for practical purposes and a representative sample of a genome in terms of reciprocal distances between USCOs on a chromosome and of distribution across chromosomes. We tested the suitability of mzl-USCOs extracted from genomes for species delimitation and phylogeny in four case studies: Anopheles mosquitos, Drosophila fruit flies, Heliconius butterflies and Darwin's finches. In almost all instances, USCOs allowed delineating species and yielded phylogenies that corresponded to those generated from whole genome data. Our phylogenetic analyses demonstrate that USCOs may complement single-gene DNA barcodes and provide more accurate taxonomic inferences. Combining USCOs from sources that used different versions of ortholog reference libraries to infer marker orthology may be challenging and, at times, impact taxonomic conclusions. However, we expect this problem to become less severe as the rapidly growing number of reference genomes provides a better representation of the number and diversity of organismal lineages.

Decoupled evolution of the Sex Peptide gene family and Sex Peptide Receptor in Drosophilidae

Across internally fertilising species, males transfer ejaculate proteins that trigger wide-ranging changes in female behaviour and physiology. Much theory has been developed to explore the drivers of ejaculate protein evolution. The accelerating availability of high-quality genomes now allows us to test how these proteins are evolving at fine taxonomic scales. Here, we use genomes from 264 species to chart the evolutionary history of Sex Peptide (SP), a potent regulator of female post-mating responses in Drosophila melanogaster. We infer that SP first evolved in the Drosophilinae subfamily and has followed markedly different evolutionary trajectories in different lineages. Outside of the Sophophora-Lordiphosa, SP exists largely as a single-copy gene with independent losses in several lineages. Within the Sophophora-Lordiphosa, the SP gene family has repeatedly and independently expanded. Up to seven copies, collectively displaying extensive sequence variation, are present in some species. Despite these changes, SP expression remains restricted to the male reproductive tract. Alongside, we document considerable interspecific variation in the presence and morphology of seminal microcarriers that, despite the critical role SP plays in microcarrier assembly in D. melanogaster, appear to be independent of changes in the presence/absence or sequence of SP. We end by providing evidence that SP's evolution is decoupled from that of its receptor, SPR, in which we detect no evidence of correlated diversifying selection. Collectively, our work describes the divergent evolutionary trajectories that a novel gene has taken following its origin and finds a surprisingly weak coevolutionary signal between a supposedly sexually antagonistic protein and its receptor.

Low Endosymbiont Incidence in Drosophila Species Across Peninsula Thailand

Arthropods are known to harbor several endosymbionts, such as Cardinium, Rickettsia, Spiroplasma, and Wolbachia. Wolbachia, for example, are the most widespread known endosymbionts in the world, which are found in about half of all arthropod species. To increase their transmission, these endosymbionts must manipulate their hosts in several ways such as cytoplasmic incompatibility and male killing. In tropical regions, endosymbiont diversity has not been studied exhaustively. Here, we checked four endosymbionts, including Cardinium, Rickettsia, Spiroplasma, and Wolbachia, in eleven Drosophila species found in Thai Peninsula. The Wolbachia strain wRi-like was found in all populations of Drosophila ananassae and Drosophila simulans. Furthermore, we found two new strains, wMalA and wMalB, in two populations of Drosophila malerkotliana. Besides Wolbachia, we did not find any of the above endosymbionts in all fly species. This work reveals the hidden diversity of endosymbionts in Drosophila and is the first exhaustive study on Drosophila in the region.

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 12 species that are geographically widespread and show varying levels of prezygotic and postzygotic isolation. Here, we use 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. Some genes undergoing rapid sequence divergence across the group were involved in chemical communication and desiccation tolerance, and may be related to the evolution of sexual isolation and adaptation. We found evidence of pervasive phylogenetic discordance caused by ancient introgression events between distant lineages within the group, and 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 as has been found in other systems. Our results show how ancient and recent introgressions confuse phylogenetic reconstruction, but may play an important role during early radiation of a group.

Neo-sex chromosome evolution shapes sex-dependent asymmetrical introgression barrier

It is increasingly recognized that sex chromosomes are not only the battlegrounds between sexes but also the Great Walls fencing off introgression between diverging lineages. Here we dissect the multifaceted roles of sex chromosomes using experimental evolution, whole-genome resequencing, and theoretical modeling, taking advantage of hybrid populations between a Drosophila sister species pair in the early stage of speciation that have different sex chromosome systems. Our work sheds light onto the complex roles of neo-sex chromosome evolution in creating a sex-dependent asymmetrical introgression barrier at a species boundary, and we show how diverse population genetic forces act in concert to explain observed patterns of introgression across the genome.

Sex chromosomes play a special role in the evolution of reproductive barriers between species. Here we describe conflicting roles of nascent sex chromosomes on patterns of introgression in an experimental hybrid swarm. Drosophila nasuta and Drosophila albomicans are recently diverged, fully fertile sister species that have different sex chromosome systems. The fusion between an autosome (Muller CD) with the ancestral X and Y gave rise to neo-sex chromosomes in D. albomicans, while Muller CD remains unfused in D. nasuta. We found that a large block containing overlapping inversions on the neo-sex chromosome stood out as the strongest barrier to introgression. Intriguingly, the neo-sex chromosome introgression barrier is asymmetrical and sex-dependent. Female hybrids showed significant D. albomicans–biased introgression on Muller CD (neo-X excess), while males showed heterosis with excessive (neo-X, D. nasuta Muller CD) genotypes. We used a population genetic model to dissect the interplay of sex chromosome drive, heterospecific pairing incompatibility between the neo-sex chromosomes and unfused Muller CD, neo-Y disadvantage, and neo-X advantage in generating the observed sex chromosome genotypes in females and males. We show that moderate neo-Y disadvantage and D. albomicans specific meiotic drive are required to observe female-specific D. albomicans–biased introgression in this system, together with pairing incompatibility and neo-X advantage. In conclusion, this hybrid swarm between a young species pair sheds light onto the multifaceted roles of neo-sex chromosomes in a sex-dependent asymmetrical introgression barrier at a species boundary.

Dynamics and Impacts of Transposable Element Proliferation in the Drosophila nasuta Species Group Radiation

Transposable element (TE) mobilization is a constant threat to genome integrity. Eukaryotic organisms have evolved robust defensive mechanisms to suppress their activity, yet TEs can escape suppression and proliferate, creating strong selective pressure for host defense to adapt. This genomic conflict fuels a never-ending arms race that drives the rapid evolution of TEs and recurrent positive selection of genes involved in host defense; the latter has been shown to contribute to postzygotic hybrid incompatibility. However, how TE proliferation impacts genome and regulatory divergence remains poorly understood. Here, we report the highly complete and contiguous (N50 = 33.8-38.0 Mb) genome assemblies of seven closely related Drosophila species that belong to the nasuta species group-a poorly studied group of flies that radiated in the last 2 My. We constructed a high-quality de novo TE library and gathered germline RNA-seq data, which allowed us to comprehensively annotate and compare TE insertion patterns between the species, and infer the evolutionary forces controlling their spread. We find a strong negative association between TE insertion frequency and expression of genes nearby; this likely reflects survivor bias from reduced fitness impact of TEs inserting near lowly expressed, nonessential genes, with limited TE-induced epigenetic silencing. Phylogenetic analyses of insertions of 147 TE families reveal that 53% of them show recent amplification in at least one species. The most highly amplified TE is a nonautonomous DNA element (Drosophila INterspersed Element; DINE) which has gone through multiple bouts of expansions with thousands of full-length copies littered throughout each genome. Across all TEs, we find that TEs expansions are significantly associated with high expression in the expanded species consistent with suppression escape. Thus, whereas horizontal transfer followed by the invasion of a naïve genome has been highlighted to explain the long-term survival of TEs, our analysis suggests that evasion of host suppression of resident TEs is a major strategy to persist over evolutionary times. Altogether, our results shed light on the heterogenous and context-dependent nature in which TEs affect gene regulation and the dynamics of rampant TE proliferation amidst a recently radiated species group.

Evolution of assortative mating following selective introgression of pigmentation genes between two Drosophila species

Adaptive introgression is ubiquitous in animals, but experimental support for its role in driving speciation remains scarce. In the absence of conscious selection, admixed laboratory strains of Drosophila asymmetrically and progressively lose alleles from one parental species and reproductive isolation against the predominant parent ceases after 10 generations. Here, we selectively introgressed during 1 year light pigmentation genes of D. santomea into the genome of its dark sibling D. yakuba, and vice versa. We found that the pace of phenotypic change differed between the species and the sexes and identified through genome sequencing common as well as distinct introgressed loci in each species. Mating assays showed that assortative mating between introgressed flies and both parental species persisted even after 4 years (~60 generations) from the end of the selection. Those results indicate that selective introgression of as low as 0.5% of the genome can beget morphologically distinct and reproductively isolated strains, two prerequisites for the delimitation of new species. Our findings hence represent a significant step toward understanding the genome‐wide dynamics of speciation‐through‐introgression.

Widespread introgression across a phylogeny of 155 Drosophila genomes

Genome-scale sequence data have invigorated the study of hybridization and introgression, particularly in animals. However, outside of a few notable cases, we lack systematic tests for introgression at a larger phylogenetic scale across entire clades. Here, we leverage 155 genome assemblies from 149 species to generate a fossil-calibrated phylogeny and conduct multilocus tests for introgression across 9 monophyletic radiations within the genus Drosophila. Using complementary phylogenomic approaches, we identify widespread introgression across the evolutionary history of Drosophila. Mapping gene-tree discordance onto the phylogeny revealed that both ancient and recent introgression has occurred across most of the 9 clades that we examined. Our results provide the first evidence of introgression occurring across the evolutionary history of Drosophila and highlight the need to continue to study the evolutionary consequences of hybridization and introgression in this genus and across the tree of life.

Current and future potential global distribution of the invading species Drosophila nasuta (Diptera: Drosophilidae)

Species distribution modelling has been widely employed to indicate probable areas of invasion and to guide management strategies. Drosophila nasuta is native to Asia and has invaded Africa, islands of the Indian and Pacific Oceans and, more recently, the Americas. This species has been dispersing rapidly in the past decade, dominating the assemblage of drosophilids in numerous invaded territories, especially in protected areas. Here, we model the potential geographic distribution of D. nasuta for the present and two future scenarios. We also determine the environmental variables that most influence its distribution and investigate the risk of invasion in protected areas. Drosophila nasuta has the potential to expand its occurrence, especially on continents that have already been invaded. Variables related to greater rainfall were those that most influenced its distribution. The projections for the two future scenarios revealed a small increase in the distribution of the species compared to the projection for the present. The largest overlaps between the projected areas to be invaded by D. nasuta and territories in protected areas were found for Central and South America. The predictive maps delineated here can assist in the establishment of management plans directed at the conservation of biodiversity.

Rapid Divergence of Key Spermatogenesis Genes in nasuta-Subgroup of Drosophila

The crosses between closely related Drosophila species usually produce sterile hybrid males with spermatogenesis disrupted at post-meiotic phase, especially in sperm individualization stage than the pre-meiotic stage. This is possibly due to the rapid interspecies divergence of male sex and reproduction-related genes. Here we annotated 11 key spermatogenesis genes in 35 strains of species belonging to nasuta-subgroup of Drosophila, where many interspecies crosses produce sterile males. We characterized the divergence and polymorphism in the protein coding regions by employing gene-wide, codon-wide, and lineage-specific selection analysis to test the mode and strength of selection acting on these genes. Our analysis showed signature of positive selection at bag of marbles (bam) and benign gonial cell neoplasma (bgcn) despite the selection constrains and the absence of endosymbiont infection which could potentially drive rapid divergence due to an arms race while roughex (rux) showed lineage-specific rapid divergence in frontal sheen complex of nasuta-subgroup. cookie monster (comr) showed rapid divergence consistent with the possibility of meiotic arrest observed in sterile hybrids of Drosophila species. Rapid divergence observed at don juan (dj) and Mst98Ca-like was consistent with fused sperm-tail abnormality observed in the hybrids of Drosophila nasuta and Drosophila albomicans. These findings highlight the potential role of rapid nucleotide divergence in bringing about hybrid incompatibility in the form of male sterility; however, additional genetic manipulation studies can widen our understanding of hybrid incompatibilities. Furthermore, our study emphasizes the importance of young species belonging to nasuta-subgroup of Drosophila in studying post-zygotic reproductive isolation mechanisms.

Shared evolutionary trajectories of three independent neo-sex chromosomes in Drosophila

Dosage compensation (DC) on the X Chromosome counteracts the deleterious effects of gene loss on the Y Chromosome. However, DC is not efficient if the X Chromosome also degenerates. This indeed occurs in Drosophila miranda, in which both the neo-Y and the neo-X are under accelerated pseudogenization. To examine the generality of this pattern, we investigated the evolution of two additional neo-sex chromosomes that emerged independently in D. albomicans and D. americana and reanalyzed neo-sex chromosome evolution in D. miranda. Comparative genomic and transcriptomic analyses revealed that the pseudogenization rate on the neo-X is also accelerated in D. albomicans and D. americana although to a lesser extent than in D. miranda. In males, neo-X-linked genes whose neo-Y-linked homologs are pseudogenized tended to be up-regulated more than those whose neo-Y-linked homologs remain functional. Moreover, genes under strong functional constraint and genes highly expressed in the testis tended to remain functional on the neo-X and neo-Y, respectively. Focusing on the D. miranda and D. albomicans neo-sex chromosomes that emerged independently from the same autosome, we further found that the same genes tend to become pseudogenized in parallel on the neo-Y. These genes include Idgf6 and JhI-26, which may be unnecessary or even harmful in males. Our results indicate that neo-sex chromosomes in Drosophila share a common evolutionary trajectory after their emergence, which may prevent sex chromosomes from being an evolutionary dead end.

Highly contiguous assemblies of 101 drosophilid genomes

Over 100 years of studies in Drosophila melanogaster and related species in the genus Drosophila have facilitated key discoveries in genetics, genomics, and evolution. While high-quality genome assemblies exist for several species in this group, they only encompass a small fraction of the genus. Recent advances in long-read sequencing allow high-quality genome assemblies for tens or even hundreds of species to be efficiently generated. Here, we utilize Oxford Nanopore sequencing to build an open community resource of genome assemblies for 101 lines of 93 drosophilid species encompassing 14 species groups and 35 sub-groups. The genomes are highly contiguous and complete, with an average contig N50 of 10.5 Mb and greater than 97% BUSCO completeness in 97/101 assemblies. We show that Nanopore-based assemblies are highly accurate in coding regions, particularly with respect to coding insertions and deletions. These assemblies, along with a detailed laboratory protocol and assembly pipelines, are released as a public resource and will serve as a starting point for addressing broad questions of genetics, ecology, and evolution at the scale of hundreds of species.

De novo assembly, annotation and gene expression profiles of gonads of Cytorace-3, a hybrid lineage of Drosophila nasuta nasuta and D. n. albomicans

Cytorace-3 is a laboratory evolved hybrid lineage of Drosophila nasuta nasuta males and Drosophila nasuta albomicans females currently passing ~850 generations. To assess interracial hybridization effects on gene expression in Cytorace-3 we profiled the transcriptomes of mature ovaries and testes by employing Illumina sequencing technology and de novo transcriptome assembling strategies. We found 26% of the ovarian, and 14% of testis genes to be differentially expressed in Cytorace-3 relative to the expressed genes in the parental gonadal transcriptomes. About 5% of genes exhibited additive gene expression pattern in the ovary and 3% in the testis, while the remaining genes were misexpressed in Cytorace-3. Nearly 772 of these misexpressed genes in the ovary and 413 in the testis were either over- or under-dominant. Genes following D. n. nasuta dominance was twice (270 genes) than D. n. albomicans dominance (133 genes) in the ovary. In contrast, only 105 genes showed D. n. nasuta dominance and 207 showed D. n. albomicans dominance in testis transcriptome. Of the six expression inheritance patterns, conserved inheritance pattern was predominant for both ovary (73%) and testis (85%) in Cytorace-3. This study is the first to provide an overview of the expression divergence and inheritance patterns of the transcriptomes in an independently evolving distinct hybrid lineage of Drosophila. This recorded expression divergence in Cytorace-3 surpasses that between parental lineages illustrating the strong impact of hybridization driving rapid gene expression changes.

Chromosome-Level Assembly of Drosophila bifasciata Reveals Important Karyotypic Transition of the X Chromosome

The Drosophila obscura species group is one of the most studied clades of Drosophila and harbors multiple distinct karyotypes. Here we present a de novo genome assembly and annotation of D. bifasciata, a species which represents an important subgroup for which no high-quality chromosome-level genome assembly currently exists. We combined long-read sequencing (Nanopore) and Hi-C scaffolding to achieve a highly contiguous genome assembly approximately 193 Mb in size, with repetitive elements constituting 30.1% of the total length. Drosophila bifasciata harbors four large metacentric chromosomes and the small dot, and our assembly contains each chromosome in a single scaffold, including the highly repetitive pericentromeres, which were largely composed of Jockey and Gypsy transposable elements. We annotated a total of 12,821 protein-coding genes and comparisons of synteny with D. athabasca orthologs show that the large metacentric pericentromeric regions of multiple chromosomes are conserved between these species. Importantly, Muller A (X chromosome) was found to be metacentric in D. bifasciata and the pericentromeric region appears homologous to the pericentromeric region of the fused Muller A-AD (XL and XR) of pseudoobscura/affinis subgroup species. Our finding suggests a metacentric ancestral X fused to a telocentric Muller D and created the large neo-X (Muller A-AD) chromosome ∼15 MYA. We also confirm the fusion of Muller C and D in D. bifasciata and show that it likely involved a centromere-centromere fusion.

Ancestral male recombination in Drosophila albomicans produced geographically restricted neo-Y chromosome haplotypes varying in age and onset of decay

Male Drosophila typically have achiasmatic meiosis, and fusions between autosomes and the Y chromosome have repeatedly created non-recombining neo-Y chromosomes that degenerate. Intriguingly, Drosophila nasuta males recombine, but their close relative D. albomicans reverted back to achiasmy after evolving neo-sex chromosomes. Here we use genome-wide polymorphism data to reconstruct the complex evolutionary history of neo-sex chromosomes in D. albomicans and examine the effect of recombination and its cessation on the initiation of neo-Y decay. Population and phylogenomic analyses reveal three distinct neo-Y types that are geographically restricted. Due to ancestral recombination with the neo-X, overall nucleotide diversity on the neo-Y is similar to the neo-X but severely reduced within neo-Y types. Consistently, the neo-Y chromosomes fail to form a monophyletic clade in sliding window trees outside of the region proximal to the fusion. Based on tree topology changes, we inferred the recombination breakpoints that produced haplotypes specific to each neo-Y type. We show that recombination became suppressed at different time points for the different neo-Y haplotypes. Haplotype age correlates with onset of neo-Y decay, and older neo-Y haplotypes show more fixed gene disruption via frameshift indels and down-regulation of neo-Y alleles. Genes are downregulated independently on the different neo-Ys, but are depleted of testes-expressed genes across all haplotypes. This indicates that genes important for male function are initially shielded from degeneration. Our results offer a time course of the early progression of Y chromosome evolution, showing how the suppression of recombination, through the reversal to achiasmy in D. albomicans males, initiates the process of degeneration.