Widespread introgression across a phylogeny of 155 Drosophila genomes

How Important Is Variation in Extrinsic Reproductive Isolation to the Process of Speciation?

The strength of reproductive isolation (RI) between two or more lineages during the process of speciation can vary by the ecological conditions. However, most speciation research has been limited to studying how ecologically dependent RI varies among a handful of broadly categorized environments. Very few studies consider the variability of RI and its effects on speciation at finer scales-that is, within each environment due to spatial or temporal environmental heterogeneity. Such variation in RI across time and/or space may inhibit speciation through leaky reproductive barriers or promote speciation by facilitating reinforcement. To investigate this overlooked aspect of speciation research, we conducted a literature review of existing studies of variation in RI in the field and then conducted individual-based simulations to examine how variation in hybrid fitness across time and space affects the degree of gene flow. Our simulations indicate that the presence of variation in hybrid fitness across space and time often leads to an increase in gene flow compared to scenarios where hybrid fitness remains static. This observation can be attributed to the convex relationship between the degree of gene flow and the strength of selection on hybrids. Our simulations also show that the effect of variation in RI on facilitating gene flow is most pronounced when RI, on average, is relatively low. This finding suggests that it could serve as an important mechanism to explain why the completion of speciation is often challenging. While direct empirical evidence documenting variation in extrinsic RI is limited, we contend that it is a prevalent yet underexplored phenomenon. We support this argument by proposing common scenarios in which RI is likely to exhibit variability and thus influence the process of speciation.

The Meaning and Measure of Concordance Factors in Phylogenomics

As phylogenomic datasets have grown in size, researchers have developed new ways to measure biological variation and to assess statistical support for specific branches. Larger datasets have more sites and loci and therefore less sampling variance. While we can more accurately measure the mean signal in these datasets, lower sampling variance is often reflected in uniformly high measures of branch support-such as the bootstrap and posterior probability-limiting their utility. Larger datasets have also revealed substantial biological variation in the topologies found across individual loci, such that the single species tree inferred by most phylogenetic methods represents a limited summary of the data for many purposes. In contrast to measures of statistical support, the degree of underlying topological variation among loci should be approximately constant regardless of the size of the dataset. "Concordance factors" (CFs) and similar statistics have therefore become increasingly important tools in phylogenetics. In this review, we explain why CFs should be thought of as descriptors of topological variation rather than as measures of statistical support, and argue that they provide important information about the predictive power of the species tree not contained in measures of support. We review a growing suite of statistics for measuring concordance, compare them in a common framework that reveals their interrelationships, and demonstrate how to calculate them using an example from birds. We also discuss how measures of topological variation might change in the future as we move beyond estimating a single "tree of life" toward estimating the myriad evolutionary histories underlying genomic variation.

Deep Ancestral Introgressions between Ovine Species Shape Sheep Genomes via Argali-Mediated Gene Flow

Previous studies revealed extensive genetic introgression between Ovis species, which affects genetic adaptation and morphological traits. However, the exact evolutionary scenarios underlying the hybridization between sheep and allopatric wild relatives remain unknown. To address this problem, we here integrate the reference genomes of several ovine and caprine species: domestic sheep, argali, bighorn sheep, snow sheep, and domestic goats. Additionally, we use 856 whole genomes representing 169 domestic sheep populations and their six wild relatives: Asiatic mouflon, urial, argali, snow sheep, thinhorn sheep, and bighorn sheep. We implement a comprehensive set of analyses to test introgression among these species. We infer that the argali lineage originated ∼3.08 to 3.35 Mya and hybridized with the ancestor of Pachyceriforms (e.g. bighorn sheep and snow sheep) at ∼1.56 Mya. Previous studies showed apparent introgression from North American Pachyceriforms into the Bashibai sheep, a Chinese native sheep breed, despite of their wide geographic separation. We show here that, in fact, the apparent introgression from the Pachyceriforms into Bashibai can be explained by the old introgression from Pachyceriforms into argali and subsequent recent introgression from argali into Bashibai. Our results illustrate the challenges of estimating complex introgression histories and provide an example of how indirect and direct introgression can be distinguished.

Towards Reliable Detection of Introgression in the Presence of Among-Species Rate Variation

The role of interspecific hybridization has recently seen increasing attention, especially in the context of diversification dynamics. Genomic research has now made it abundantly clear that both hybridization and introgression-the exchange of genetic material through hybridization and backcrossing-are far more common than previously thought. Besides cases of ongoing or recent genetic exchange between taxa, an increasing number of studies report "ancient introgression"- referring to results of hybridization that took place in the distant past. However, it is not clear whether commonly used methods for the detection of introgression are applicable to such old systems, given that most of these methods were originally developed for analyses at the level of populations and recently diverged species, affected by recent or ongoing genetic exchange. In particular, the assumption of constant evolutionary rates, which is implicit in many commonly used approaches, is more likely to be violated as evolutionary divergence increases. To test the limitations of introgression detection methods when being applied to old systems, we simulated thousands of genomic datasets under a wide range of settings, with varying degrees of among-species rate variation and introgression. Using these simulated datasets, we showed that some commonly applied statistical methods, including the D-statistic and certain tests based on sets of local phylogenetic trees, can produce false-positive signals of introgression between divergent taxa that have different rates of evolution. These misleading signals are caused by the presence of homoplasies occurring at different rates in different lineages. To distinguish between the patterns caused by rate variation and genuine introgression, we developed a new test that is based on the expected clustering of introgressed sites along the genome and implemented this test in the program Dsuite.

Contrast and Genomic Characterisation of Ancient and Recent Interspecific Introgression Between Deeply Diverged Moustache Toads (Leptobrachium)

Recent genomic analyses have provided new insights into the process of interspecific introgression and its consequences on species evolution. Most recent studies, however, focused on hybridization between recently radiated species, with few examining the genomic outcomes of ancient hybridization across deeply diverged species. Using whole genome data of moustache toads (Leptobrachium), we identified signals of three hybridization events among nine species that diverged at the Eocene. An ancient introgression from L. leishanense to the ancestral branch (C1) of L. liui introduced adaptive variants. The highly introgressed regions include genes with important functions in odorant detection and immune responses. These genes are preserved in all three descendent populations of L. liui_C1, and these regions likely have been positively selected over a long filtering process. A recent introgression occurred from L. huashen to L. tengchongense, with the introgressed regions being mostly neutral. Furthermore, one F1 hybrid individual was detected between sympatric L. ailaonicum and L. promustache. The signals of introgression largely disappeared after removing the hybrid individual, indicating an occasional hybridization but minimal introgression. Further examination of highly divergent but low introgressed genomic regions revealed both pre-mating isolation and genetic incompatibility as potential mechanisms of resisting introgression and maintaining species boundaries. Additionally, no large X-effect was found in these introgression events. Hybridization between deeply diverged amphibian species may be common, but detectable introgressions are likely less so, with recent introgression being mostly neutral and the rare ancient one potentially adaptive. Our findings complement recent genomic work, and together they provide a better understanding of the genomic characteristics of interspecific introgression and its significance in species adaptation and evolution.

Ancestral neural circuits potentiate the origin of a female sexual behavior in Drosophila

Courtship interactions are remarkably diverse in form and complexity among species. How neural circuits evolve to encode new behaviors that are functionally integrated into these dynamic social interactions is unknown. Here we report a recently originated female sexual behavior in the island endemic Drosophila species D. santomea, where females signal receptivity to male courtship songs by spreading their wings, which in turn promotes prolonged songs in courting males. Copulation success depends on this female signal and correlates with males' ability to adjust his singing in such a social feedback loop. Functional comparison of sexual circuitry across species suggests that a pair of descending neurons, which integrates male song stimuli and female internal state to control a conserved female abdominal behavior, drives wing spreading in D. santomea. This co-option occurred through the refinement of a pre-existing, plastic circuit that can be optogenetically activated in an outgroup species. Combined, our results show that the ancestral potential of a socially-tuned key circuit node to engage the wing motor circuit facilitates the expression of a new female behavior in appropriate sensory and motivational contexts. More broadly, our work provides insights into the evolution of social behaviors, particularly female behaviors, and the underlying neural mechanisms.

Phylogenomic Discordance is Driven by Wide-Spread Introgression and Incomplete Lineage Sorting During Rapid Species Diversification Within Rattlesnakes (Viperidae: Crotalus and Sistrurus)

-Phylogenomics allows us to uncover the historical signal of evolutionary processes through time and estimate phylogenetic networks accounting for these signals. Insight from genome-wide data further allows us to pinpoint the contributions to phylogenetic signal from hybridization, introgression, and ancestral polymorphism across the genome. Here, we focus on how these processes have contributed to phylogenetic discordance among rattlesnakes (genera Crotalus and Sistrurus), a group for which there are numerous conflicting phylogenetic hypotheses based on a diverse array of molecular datasets and analytical methods. We address the instability of the rattlesnake phylogeny using genomic data generated from transcriptomes sampled from nearly all known species. These genomic data, analyzed with coalescent and network-based approaches, reveal numerous instances of rapid speciation where individual gene trees conflict with the species tree. Moreover, the evolutionary history of rattlesnakes is dominated by incomplete speciation and frequent hybridization, both of which have likely influenced past interpretations of phylogeny. We present a new framework in which the evolutionary relationships of this group can only be understood in light of genome-wide data and network-based analytical methods. Our data suggest that network radiations, like those seen within the rattlesnakes, can only be understood in a phylogenomic context, necessitating similar approaches in our attempts to understand evolutionary history in other rapidly radiating species.

Function and evolution of Ir52 receptors in mate detection in Drosophila

Identifying a suitable mating partner is an ancient and critical biological problem. How a fruit fly distinguishes a fly of the same species from flies of innumerable related species remains unclear. We analyze the Ir52 receptors, expressed in taste neurons on the fly legs and encoded by a cluster of genes. We find that the cluster shows dynamic evolution, rapidly expanding and contracting over evolutionary time. We develop a novel in vivo expression system and find that Ir52 receptors respond differently to pheromone extracts of different fly species. The receptors are activated by some compounds and inhibited by others, with different receptors showing distinct response profiles. Circuit mapping shows that Ir52 neurons are pre-synaptic to sexually dimorphic neurons that overlap with neurons acting in courtship behavior. Our results support a model in which Ir52 receptors detect information about the species of a potential mating partner.

Dollo's law of irreversibility in the post-genomic age

Dollo's law of irreversibility argues that evolution cannot revert to earlier states. It has remained controversial ever since its inception in the 19th century. Enabled by advances in phylogenomics and functional genomics, recent studies show that there are very likely some cases of 'breaking Dollo's law'. As post-genomic research grows from showing patterns to revealing processes, new emphasis is needed on the molecular mechanisms by which Dollo's law might be broken. Shifting the argument from 'if it happened' to 'how it happened' will provide richer understanding of organismal and evolutionary biology. Motivated by case studies and novel avenues to test trait loss and regain, we outline a set of alternative hypotheses to be evaluated and what the outcomes tell us about evolution.

Females drive postmating reproductive trait evolution across Drosophila species, but not via remating rate

While traits that contribute to premating sexual interactions are known to be wildly diverse, much less is known about the diversity of postmating (especially female) reproductive traits and the mechanisms shaping this diversity. To assess the rate, pattern, and potential drivers of postmating reproductive trait evolution, we analyzed male and female traits across up to 30 Drosophila species within a phylogenetic comparative framework. In addition to postmating reproductive morphology (e.g., sperm length, reproductive tract length and mass), we also quantified mating behaviors including female remating rate-a common proxy for the strength of postmating sexual selection. We found evidence for strong coevolution between male and female postmating traits (specifically sperm length and sperm storage organ size). However, remating rate was not associated with the rate of evolution or exaggeration of either male or female postmating reproductive morphology, once phylogenetic relatedness was accounted for. We infer that female-mediated and intersexual selection predominantly drive the evolution of our postmating morphological traits, including via divergent male and female interests in controlling paternity. In comparison, remating rate has a complex and likely secondary role in shaping this evolution, in part because this trait can be both a driver and a product of postmating selection.

Prezygotic barriers effectively limit hybridization in a rapid evolutionary radiation

Hybridization is increasingly recognized as an important evolutionary process across the tree of life. In many clades, phylogenomic approaches have permitted unparalleled insight into the extent and frequency of hybridization. However, we continue to lack a deep understanding of the factors that limit and shape patterns of hybridization, especially in evolutionary radiations. In this study, we characterized patterns of introgression across Costus (Costaceae), a young evolutionary radiation of tropical understory plants that maintain widespread interfertility despite exhibiting strong prezygotic reproductive isolation. We analyzed a phylogenomic dataset of 756 genes from 54 Costus species using multiple complementary approaches - D-statistics, gene-tree-based tests, and phylogenetic network analyses - to detect and characterize introgression events throughout the evolutionary history of the radiation. Our results identified a moderate number of introgression events, including a particularly ancient, well-supported event spanning one of the deepest divergences in the clade. Most introgression events occurred between taxa or ancestral lineages that shared the same pollination syndrome (bee-pollinated or hummingbird-pollinated). These findings suggest that prezygotic barriers, including pollinator specialization, have been key to the balance between introgression and reproductive isolation in Costus.

Adhesive and mechanical properties of the glue produced by 25 Drosophila species

Drosophila glue, a bioadhesive produced by fly larvae to attach themselves to a substrate for several days, has recently gained attention for its peculiar adhesive and mechanical properties. Although Drosophila glue production was described more than 50 years ago, a general survey of the adhesive and mechanical properties of this proteinaceous gel across Drosophila species is lacking. To measure adhesion, we present here a protocol that is robust to variations in protocol parameters, pupal age and calculation methods. We find that the glue, which covers the entire pupal surface, increases the animal rigidity and plasticity when bound to a glass slide. Our survey of pupal adhesion in 25 Drosophilidae species reveals a wide range of phenotypes, from species that produce no or little glue and adhere little, to species that produce high amounts of glue and adhere strongly. One species, D. hydei, stands out from the rest and emerges as a promising model for the development of future bioadhesives, as it has the highest detachment force per glue area and produces relatively large amounts of glue relative to its size. We also observe that species that invest more in glue tend to live in more windy and less rainy climates, suggesting that differences in pupal adhesion properties across species are shaped by ecological factors. Our present survey provides a basis for future biomimetic studies based on Drosophila glue.

Whole genome phylogenomics helps to resolve the phylogenetic position of the Zygothrica genus group (Diptera, Drosophilidae) and the causes of previous incongruences

Incomplete Lineage Sorting (ILS) and introgression are among the two main factors causing incongruence between gene and species trees. Advances in phylogenomic studies have allowed us to overcome most of these issues, providing reliable phylogenetic hypotheses while revealing the underlying evolutionary scenario. Across the last century, many incongruent phylogenetic reconstructions were recovered for Drosophilidae, employing a limited sampling of genetic markers or species. In these studies, the monophyly and the phylogenetic positioning of the Zygothrica genus group stood out as one of the most controversial questions. Thus, here, we addressed these issues using a phylogenomic approach, while accessing the influence of ILS and introgressions on the diversification of these species and addressing the spatio-temporal scenario associated with their evolution. For this task, the genomes of nine specimens from six Neotropical species belonging to the Zygothrica genus group were sequenced and evaluated in a phylogenetic framework encompassing other 39 species of Drosophilidae. Nucleotide and amino acid sequences recovered for a set of 2,534 single-copy genes by BUSCO were employed to reconstruct maximum likelihood (ML) concatenated and multi-species coalescent (MSC) trees. Likelihood mapping, quartet sampling, and reticulation tests were employed to infer the level and causes of incongruence. Lastly, a penalized-likelihood molecular clock strategy with fossil calibrations was performed to infer divergence times. Taken together, our results recovered the subdivision of Drosophila into six different lineages, one of which clusters species of the Zygothrica genus group (except for H. duncani). The divergence of this lineage was dated to Oligocene ∼ 31 Mya and seems to have occurred in the same timeframe as other key diversification within Drosophila. According to the concatenated and MSC strategies, this lineage is sister to the clade joining Drosophila (Siphlodora) with the Hawaiian Drosophila and Scaptomyza. Likelihood mapping, quartet sampling, reticulation reconstructions as well as introgression tests revealed that this lineage was the target of several hybridization events involving the ancestors of different Drosophila lineages. Thus, our results generally show introgression as a major source of previous incongruence. Nevertheless, the similar diversification times recovered for several of the Neotropical Drosophila lineages also support the scenario of multiple and simultaneous diversifications taking place at the base of Drosophilidae phylogeny, at least in the Neotropics.

Limited population structure but signals of recent selection in introduced African Fig Fly (Zaprionus indianus) in North America

Invasive species have devastating consequences for human health, food security, and the environment. Many invasive species adapt to new ecological niches following invasion, but little is known about the early steps of adaptation. Here we examine population genomics of a recently introduced drosophilid in North America, the African Fig Fly, Zaprionus indianus. This species is likely intolerant of subfreezing temperatures and recolonizes temperate environments yearly. We generated a new chromosome-level genome assembly for Z. indianus. Using resequencing of over 200 North American individuals collected over four years in temperate Virginia, plus a single collection from subtropical Florida, we tested for signatures of recolonization, population structure, and adaptation within invasive populations. We show founding populations are sometimes small and contain close genetic relatives, yet temporal population structure and differentiation of populations is mostly absent across recurrent recolonization events. Although we find limited signals of genome-wide spatial or temporal population structure, we identify haplotypes on the X chromosome that are repeatedly differentiated between Virginia and Florida populations. These haplotypes show signatures of natural selection and are not found in African populations. We also find evidence for several large structural polymorphisms segregating within North America populations and show X chromosome evolution in invasive populations is strikingly different from the autosomes. These results show that despite limited population structure, populations may rapidly evolve genetic differences early in an invasion. Further uncovering how these genomic regions influence invasive potential and success in new environments will advance our understanding of how organisms evolve in changing environments.

Molecular phylogeny and species diversity of the genus Dichaetophora Duda and related taxa (Diptera: Drosophilidae)

Our intensive surveys of wild drosophilids in East and Southeast Asia discovered a great species diversity (more than 100 putatively new species) of the genus Dichaetophora, which is currently comprised of 67 formally described species assigned into five species groups, i.e., agbo, tenuicauda, acutissima, sinensis and trilobita. In the present study, we delimited species from a huge amount of samples of Dichaetophora and allied taxa (the genus Mulgravea and the subgenus Dudaica of Drosophila) collected from a wide range of the Oriental and east Palearctic regions. We first sorted all specimens into morpho-species, and representative specimen(s) selected from each morpho-species were subjected to barcoding of COI (the cytochrome c oxidase subunit I gene) sequences. The applied ASAP (Assemble Species by Automatic Partitioning) analysis estimated a total of 166 to 168 MOTUs (molecular operational taxonomic units). Integrating this result with morphological evidence from re-examined, detailed structures of male terminalia, we recognized a total of 144 (109 new and 35 known) species in our sample. Out of them, 83 species representing the supraspecific taxa of Dichaetophora, Mulgravea and Dudaica were selected, along with 33 species from major genera and subgenera of Drosophila in the tribe Drosophilini, as in-group and four species from the tribe Colocasiomyini as out-group for phylogenetic reconstruction based on 12 nuclear gene markers. In the trees constructed by the maximum likelihood and Bayesian inference methods, the three focal taxa (i.e., Dichaetophora, Mulgravea and Dudaica) formed a clade provisionally called the "pan-Dichaetophora". Within this large clade, the agbo, tenuicauda, sinensis and trilobita groups of Dichaetophora, Mulgravea and Dudaica were recovered as monophyletic groups, but Dichaetophora and its acutissima group were regarded as paraphyletic. In addition, two clusters were recognized among ungrouped species of Dichaetophora. Thus, the present study has uncovered some issues concerning the taxonomy of the pan-Dichaetophora. Such issues will be addressed elsewhere in the phylogenetic reclassification of the pan-Dichaetophora, along with descriptions/redescriptions of a large number of new/known species delimited in the present study.

Reticulate evolution: Detection and utility in the phylogenomics era

Phylogenomics has enriched our understanding that the Tree of Life can have network-like or reticulate structures among some taxa and genes. Two non-vertical modes of evolution - hybridization/introgression and horizontal gene transfer - deviate from a strictly bifurcating tree model, causing non-treelike patterns. However, these reticulate processes can produce similar patterns to incomplete lineage sorting or recombination, potentially leading to ambiguity. Here, we present a brief overview of a phylogenomic workflow for inferring organismal histories and compare methods for distinguishing modes of reticulate evolution. We discuss how the timing of coalescent events can help disentangle introgression from incomplete lineage sorting and how horizontal gene transfer events can help determine the relative timing of speciation events. In doing so, we identify pitfalls of certain methods and discuss how to extend their utility across the Tree of Life. Workflows, methods, and future directions discussed herein underscore the need to embrace reticulate evolutionary patterns for understanding the timing and rates of evolutionary events, providing a clearer view of life's history.

Prevalent Fast Evolution of Genes Involved in Heterochromatin Functions

Heterochromatin is a gene-poor and repeat-rich genomic compartment universally found in eukaryotes. Despite its low transcriptional activity, heterochromatin plays important roles in maintaining genome stability, organizing chromosomes, and suppressing transposable elements. Given the importance of these functions, it is expected that genes involved in heterochromatin regulation would be highly conserved. Yet, a handful of these genes were found to evolve rapidly. To investigate whether these previous findings are anecdotal or general to genes modulating heterochromatin, we compile an exhaustive list of 106 candidate genes involved in heterochromatin functions and investigate their evolution over short and long evolutionary time scales in Drosophila. Our analyses find that these genes exhibit significantly more frequent evolutionary changes, both in the forms of amino acid substitutions and gene copy number change, when compared to genes involved in Polycomb-based repressive chromatin. While positive selection drives amino acid changes within both structured domains with diverse functions and intrinsically disordered regions, purifying selection may have maintained the proportions of intrinsically disordered regions of these proteins. Together with the observed negative associations between the evolutionary rate of these genes and the genomic abundance of transposable elements, we propose an evolutionary model where the fast evolution of genes involved in heterochromatin functions is an inevitable outcome of the unique functional roles of heterochromatin, while the rapid evolution of transposable elements may be an effect rather than cause. Our study provides an important global view of the evolution of genes involved in this critical cellular domain and provides insights into the factors driving the distinctive evolution of heterochromatin.

Genomic Diversity Illuminates the Environmental Adaptation of Drosophila suzukii

Biological invasions carry substantial practical and scientific importance and represent natural evolutionary experiments on contemporary timescales. Here, we investigated genomic diversity and environmental adaptation of the crop pest Drosophila suzukii using whole-genome sequencing data and environmental metadata for 29 population samples from its native and invasive range. Through a multifaceted analysis of this population genomic data, we increase our understanding of the D. suzukii genome, its diversity and its evolution, and we identify an appropriate genotype-environment association pipeline for our dataset. Using this approach, we detect genetic signals of local adaptation associated with nine distinct environmental factors related to altitude, wind speed, precipitation, temperature, and human land use. We uncover unique functional signatures for each environmental variable, such as the prevalence of cuticular genes associated with annual precipitation. We also infer biological commonalities in the adaptation to diverse selective pressures, particularly in terms of the apparent contribution of nervous system evolution to enriched processes (ranging from neuron development to circadian behavior) and to top genes associated with all nine environmental variables. Our findings therefore depict a finer-scale adaptive landscape underlying the rapid invasion success of this agronomically important species.

Incomplete lineage sorting and hybridization underlie tree discordance in Petunia and related genera (Petunieae, Solanaceae)

Despite the overarching history of species divergence, phylogenetic studies often reveal distinct topologies across regions of the genome. The sources of these gene tree discordances are variable, but incomplete lineage sorting (ILS) and hybridization are among those with the most biological importance. Petunia serves as a classic system for studying hybridization in the wild. While field studies suggest that hybridization is frequent, the extent of reticulation within Petunia and its closely related genera has never been examined from a phylogenetic perspective. In this study, we used transcriptomic data from 11 Petunia, 16 Calibrachoa, and 10 Fabiana species to illuminate the relationships between these species and investigate whether hybridization played a significant role in the diversification of the clade. We inferred that gene tree discordance within genera is linked to hybridization events along with high levels of ILS due to their rapid diversification. Moreover, network analyses estimated deeper hybridization events between Petunia and Calibrachoa, genera that have different chromosome numbers. Although these genera cannot hybridize at the present time, ancestral hybridization could have played a role in their parallel radiations, as they share the same habitat and life history.

A newly evolved gene is essential for efficient sperm entry into eggs in Drosophila melanogaster

New genes arise through a variety of evolutionary processes and provide raw material for adaptation in the face of both natural and sexual selection. De novo evolved genes emerge from previously non-protein-coding DNA sequences, and many such genes are expressed in male reproductive structures. In Drosophila melanogaster, several putative de novo genes have evolved essential roles in spermatogenesis, but whether such genes can also impact sperm function beyond the male has not been investigated. We identified a putative de novo gene, katherine johnson (kj), that is required for high levels of male fertility. Males that do not express kj produce and transfer sperm that are stored normally in females, but sperm from these males enter eggs with severely reduced efficiency. Using a tagged transgenic rescue construct, we observed that KJ protein localizes to the nuclear periphery in various stages of spermatogenesis, but is not detectable in mature sperm. These data suggest that kj exerts an effect on sperm development, the loss of which results in reduced fertilization ability. While previous bioinformatic analyses suggested the kj gene was restricted to the melanogaster group of Drosophila, we identified putative orthologs with conserved synteny, male-biased expression, and predicted protein features across the genus, as well as instances of gene loss in some lineages. Thus, kj potentially arose in the Drosophila common ancestor and subsequently evolved an essential role in D. melanogaster. Our results demonstrate a new aspect of male reproduction that has been shaped by new gene evolution and provide a molecular foothold for further investigating the mechanism of sperm entry into eggs in Drosophila.

Swordtail fish hybrids reveal that genome evolution is surprisingly predictable after initial hybridization

Over the past 2 decades, biologists have come to appreciate that hybridization, or genetic exchange between distinct lineages, is remarkably common-not just in particular lineages but in taxonomic groups across the tree of life. As a result, the genomes of many modern species harbor regions inherited from related species. This observation has raised fundamental questions about the degree to which the genomic outcomes of hybridization are repeatable and the degree to which natural selection drives such repeatability. However, a lack of appropriate systems to answer these questions has limited empirical progress in this area. Here, we leverage independently formed hybrid populations between the swordtail fish Xiphophorus birchmanni and X. cortezi to address this fundamental question. We find that local ancestry in one hybrid population is remarkably predictive of local ancestry in another, demographically independent hybrid population. Applying newly developed methods, we can attribute much of this repeatability to strong selection in the earliest generations after initial hybridization. We complement these analyses with time-series data that demonstrates that ancestry at regions under selection has remained stable over the past approximately 40 generations of evolution. Finally, we compare our results to the well-studied X. birchmanni × X. malinche hybrid populations and conclude that deeper evolutionary divergence has resulted in stronger selection and higher repeatability in patterns of local ancestry in hybrids between X. birchmanni and X. cortezi.

Evidence for gene flow and trait reversal during radiation of Mexican Goodeid fish

Understanding the phylogeographic history of a group and identifying the factors contributing to speciation is an important challenge in evolutionary biology. The Goodeinae are a group of live-bearing fishes endemic to Mexico. Here, we develop genomic resources for species within the Goodeinae and use phylogenomic approaches to characterise their evolutionary history. We sequenced, assembled and annotated the genomes of four Goodeinae species, including Ataeniobius toweri, the only matrotrophic live-bearing fish without a trophotaenia in the group. We estimated timings of species divergence and examined the extent and timing of introgression between the species to assess if this may have occurred during an early radiation, or in more recent episodes of secondary contact. We used branch-site models to detect genome-wide positive selection across Goodeinae, and we specifically asked whether this differs in A. toweri, where loss of placental viviparity has recently occurred. We found evidence of gene flow between geographically isolated species, suggesting vicariant speciation was supplemented by limited post-speciation gene flow, and gene flow may explain previous uncertainties about Goodeid phylogeny. Genes under positive selection in the group are likely to be associated with the switch to live-bearing. Overall, our studies suggest that both volcanism-driven vicariance and changes in reproductive mode influenced radiation in the Goodeinae.

Prevalent fast evolution of genes involved in heterochromatin functions

Heterochromatin is a gene-poor and repeat-rich genomic compartment universally found in eukaryotes. Despite its low transcriptional activity, heterochromatin plays important roles in maintaining genome stability, organizing chromosomes, and suppressing transposable elements (TEs). Given the importance of these functions, it is expected that the genes involved in heterochromatin regulation would be highly conserved. Yet, a handful of these genes were found to evolve rapidly. To investigate whether these previous findings are anecdotal or general to genes modulating heterochromatin, we compile an exhaustive list of 106 candidate genes involved in heterochromatin functions and investigate their evolution over short and long evolutionary time scales in Drosophila. Our analyses find that these genes exhibit significantly more frequent evolutionary changes, both in the forms of amino acid substitutions and gene copy number change, when compared to genes involved in Polycomb-based repressive chromatin. While positive selection drives amino acid changes within both structured domains with diverse functions and intrinsically disordered regions (IDRs), purifying selection may have maintained the proportions of IDRs of these proteins. Together with the observed negative associations between evolutionary rates of these genes and genomic TE abundance, we propose an evolutionary model where the fast evolution of genes involved in heterochromatin functions is an inevitable outcome of the unique functional roles of heterochromatin, while the rapid evolution of TEs may be an effect rather than cause. Our study provides an important global view of the evolution of genes involved in this critical cellular domain and provides insights into the factors driving the distinctive evolution of heterochromatin.

Gene Protein Sequence Evolution Can Predict the Rapid Divergence of Ovariole Numbers in the Drosophila melanogaster Subgroup

Ovaries play key roles in fitness and evolution: they are essential female reproductive structures that develop and house the eggs in sexually reproducing animals. In Drosophila, the mature ovary contains multiple tubular egg-producing structures known as ovarioles. Ovarioles arise from somatic cellular structures in the larval ovary called terminal filaments (TFs), formed by TF cells and subsequently enclosed by sheath (SH) cells. As in many other insects, ovariole number per female varies extensively in Drosophila. At present, however, there is a striking gap of information on genetic mechanisms and evolutionary forces that shape the well-documented rapid interspecies divergence of ovariole numbers. To address this gap, here we studied genes associated with Drosophila melanogaster ovariole number or functions based on recent experimental and transcriptional datasets from larval ovaries, including TFs and SH cells, and assessed their rates and patterns of molecular evolution in five closely related species of the melanogaster subgroup that exhibit species-specific differences in ovariole numbers. From comprehensive analyses of protein sequence evolution (dN/dS), branch-site positive selection, expression specificity (tau), and phylogenetic regressions (phylogenetic generalized least squares), we report evidence of 42 genes that showed signs of playing roles in the genetic basis of interspecies evolutionary change of Drosophila ovariole number. These included the signaling genes upd2 and Ilp5 and extracellular matrix genes vkg and Col4a1, whose dN/dS predicted ovariole numbers among species. Together, we propose a model whereby a set of ovariole-involved gene proteins have an enhanced evolvability, including adaptive evolution, facilitating rapid shifts in ovariole number among Drosophila species.

Ancient and Recent Hybridization in the Oreochromis Cichlid Fishes

Cichlid fishes of the genus Oreochromis (tilapia) are among the most important fish for inland capture fisheries and global aquaculture. Deliberate introductions of non-native species for fisheries improvement and accidental escapees from farms have resulted in admixture with indigenous species. Such hybridization may be detrimental to native biodiversity, potentially leading to genomic homogenization of populations and the loss of important genetic material associated with local adaptation. By contrast, introgression may fuel diversification when combined with ecological opportunity, by supplying novel genetic combinations. To date, the role of introgression in the evolutionary history of tilapia has not been explored. Here we studied both ancient and recent hybridization in tilapia, using whole genome resequencing of 575 individuals from 23 species. We focused on Tanzania, a natural hotspot of tilapia diversity, and a country where hybridization between exotic and native species in the natural environment has been previously reported. We reconstruct the first genome-scale phylogeny of the genus and reveal prevalent ancient gene flow across the Oreochromis phylogeny. This has likely resulted in the hybrid speciation of one species, O. chungruruensis. We identify multiple cases of recent hybridization between native and introduced species in the wild, linked to the use of non-native species in both capture fisheries improvement and aquaculture. This has potential implications for both conservation of wild populations and the development of the global tilapia aquaculture industry.

Single-fly genome assemblies fill major phylogenomic gaps across the Drosophilidae Tree of Life

Long-read sequencing is driving rapid progress in genome assembly across all major groups of life, including species of the family Drosophilidae, a longtime model system for genetics, genomics, and evolution. We previously developed a cost-effective hybrid Oxford Nanopore (ONT) long-read and Illumina short-read sequencing approach and used it to assemble 101 drosophilid genomes from laboratory cultures, greatly increasing the number of genome assemblies for this taxonomic group. The next major challenge is to address the laboratory culture bias in taxon sampling by sequencing genomes of species that cannot easily be reared in the lab. Here, we build upon our previous methods to perform amplification-free ONT sequencing of single wild flies obtained either directly from the field or from ethanol-preserved specimens in museum collections, greatly improving the representation of lesser studied drosophilid taxa in whole-genome data. Using Illumina Novaseq X Plus and ONT P2 sequencers with R10.4.1 chemistry, we set a new benchmark for inexpensive hybrid genome assembly at US $150 per genome while assembling genomes from as little as 35 ng of genomic DNA from a single fly. We present 183 new genome assemblies for 179 species as a resource for drosophilid systematics, phylogenetics, and comparative genomics. Of these genomes, 62 are from pooled lab strains and 121 from single adult flies. Despite the sample limitations of working with small insects, most single-fly diploid assemblies are comparable in contiguity (>1 Mb contig N50), completeness (>98% complete dipteran BUSCOs), and accuracy (>QV40 genome-wide with ONT R10.4.1) to assemblies from inbred lines. We present a well-resolved multi-locus phylogeny for 360 drosophilid and 4 outgroup species encompassing all publicly available (as of August 2023) genomes for this group. Finally, we present a Progressive Cactus whole-genome, reference-free alignment built from a subset of 298 suitably high-quality drosophilid genomes. The new assemblies and alignment, along with updated laboratory protocols and computational pipelines, are released as an open resource and as a tool for studying evolution at the scale of an entire insect family.

Rapid host switching of Wolbachia and even more rapid turnover of their phages and incompatibility-causing loci

About half of all insect species carry maternally inherited Wolbachia alphaproteobacteria, making Wolbachia the most common endosymbionts known in nature. Often Wolbachia spread to high frequencies within populations due to cytoplasmic incompatibility (CI), a Wolbachia-induced sperm modification caused by prophage-associated genes (cifs) that kill embryos without Wolbachia. Several Wolbachia variants also block viruses, including wMel from Drosophila melanogaster when transinfected into the mosquito Aedes aegypti. CI enables the establishment and stable maintenance of pathogen-blocking wMel in natural Ae. aegypti populations. These transinfections are reducing dengue disease incidence on multiple continents. While it has long been known that closely related Wolbachia occupy distantly related hosts, the timing of Wolbachia host switching and molecular evolution has not been widely quantified. We provide a new, conservative calibration for Wolbachia chronograms based on examples of co-divergence of Wolbachia and their insect hosts. Synthesizing publicly available and new genomic data, we use our calibration to demonstrate that wMel-like variants separated by only about 370,000 years have naturally colonized holometabolous dipteran and hymenopteran insects that diverged approximately 350 million years ago. Data from Wolbachia variants closely related to those currently dominant in D. melanogaster and D. simulans illustrate that cifs are rapidly acquired and lost among Wolbachia genomes, on a time scale of 104-105 years. This turnover occurs with and without the Wovirus prophages that contain them, with closely related cifs found in distantly related phages and distantly related cifs found in closely related phages. We present evidence for purifying selection on CI rescue function and on particular Cif protein domains. Our results quantify the tempo and mode of rapid host switching and horizontal gene transfer that underlie the spread and diversity of Wolbachia sampled from diverse host species. The wMel variants we highlight from hosts in different climates may offer new options for broadening Wolbachia-based biocontrol of diseases and pests.

Comparative analysis of Wolbachia maternal transmission and localization in host ovaries

Many insects and other animals carry microbial endosymbionts that influence their reproduction and fitness. These relationships only persist if endosymbionts are reliably transmitted from one host generation to the next. Wolbachia are maternally transmitted endosymbionts found in most insect species, but transmission rates can vary across environments. Maternal transmission of wMel Wolbachia depends on temperature in natural Drosophila melanogaster hosts and in transinfected Aedes aegypti, where wMel is used to block pathogens that cause human disease. In D. melanogaster, wMel transmission declines in the cold as Wolbachia become less abundant in host ovaries and at the posterior pole plasm (the site of germline formation) in mature oocytes. Here, we assess how temperature affects maternal transmission and underlying patterns of Wolbachia localization across 10 Wolbachia strains diverged up to 50 million years-including strains closely related to wMel-and their natural Drosophila hosts. Many Wolbachia maintain high transmission rates across temperatures, despite highly variable (and sometimes low) levels of Wolbachia in the ovaries and at the developing germline in late-stage oocytes. Identifying strains like closely related wMel-like Wolbachia with stable transmission across variable environmental conditions may improve the efficacy of Wolbachia-based biocontrol efforts as they expand into globally diverse environments.

Comparative analysis of Wolbachia maternal transmission and localization in host ovaries

Many insects and other animals carry microbial endosymbionts that influence their reproduction and fitness. These relationships only persist if endosymbionts are reliably transmitted from one host generation to the next. Wolbachia are maternally transmitted endosymbionts found in most insect species, but transmission rates can vary across environments. Maternal transmission of wMel Wolbachia depends on temperature in natural Drosophila melanogaster hosts and in transinfected Aedes aegypti, where wMel is used to block pathogens that cause human disease. In D. melanogaster, wMel transmission declines in the cold as Wolbachia become less abundant in host ovaries and at the posterior pole plasm (the site of germline formation) in mature oocytes. Here, we assess how temperature affects maternal transmission and underlying patterns of Wolbachia localization across 10 Wolbachia strains diverged up to 50 million years-including strains closely related to wMel-and their natural Drosophila hosts. Many Wolbachia maintain high transmission rates across temperatures, despite highly variable (and sometimes low) levels of Wolbachia in the ovaries and at the developing germline in late-stage oocytes. Identifying strains like closely related wMel-like Wolbachia with stable transmission across variable environmental conditions may improve the efficacy of Wolbachia-based biocontrol efforts as they expand into globally diverse environments.

Rapid response of fly populations to gene dosage across development and generations

Although the effects of genetic and environmental perturbations on multicellular organisms are rarely restricted to single phenotypic layers, our current understanding of how developmental programs react to these challenges remains limited. Here, we have examined the phenotypic consequences of disturbing the bicoid regulatory network in early Drosophila embryos. We generated flies with two extra copies of bicoid, which causes a posterior shift of the network's regulatory outputs and a decrease in fitness. We subjected these flies to EMS mutagenesis, followed by experimental evolution. After only 8-15 generations, experimental populations have normalized patterns of gene expression and increased survival. Using a phenomics approach, we find that populations were normalized through rapid increases in embryo size driven by maternal changes in metabolism and ovariole development. We extend our results to additional populations of flies, demonstrating predictability. Together, our results necessitate a broader view of regulatory network evolution at the systems level.

Detection of Ghost Introgression Requires Exploiting Topological and Branch Length Information

In recent years, the study of hybridization and introgression has made significant progress, with ghost introgression-the transfer of genetic material from extinct or unsampled lineages to extant species-emerging as a key area for research. Accurately identifying ghost introgression, however, presents a challenge. To address this issue, we focused on simple cases involving 3 species with a known phylogenetic tree. Using mathematical analyses and simulations, we evaluated the performance of popular phylogenetic methods, including HyDe and PhyloNet/MPL, and the full-likelihood method, Bayesian Phylogenetics and Phylogeography (BPP), in detecting ghost introgression. Our findings suggest that heuristic approaches relying on site-pattern counts or gene-tree topologies struggle to differentiate ghost introgression from introgression between sampled non-sister species, frequently leading to incorrect identification of donor and recipient species. The full-likelihood method BPP uses multilocus sequence alignments directly-hence taking into account both gene-tree topologies and branch lengths, by contrast, is capable of detecting ghost introgression in phylogenomic datasets. We analyzed a real-world phylogenomic dataset of 14 species of Jaltomata (Solanaceae) to showcase the potential of full-likelihood methods for accurate inference of introgression.

The evolutionary novelty of insect defensins: from bacterial killing to toxin neutralization

Insect host defense comprises two complementary dimensions, microbial killing-mediated resistance and microbial toxin neutralization-mediated resilience, both jointly providing protection against pathogen infections. Insect defensins are a class of effectors of innate immunity primarily responsible for resistance to Gram-positive bacteria. Here, we report a newly originated gene from an ancestral defensin via genetic deletion following gene duplication in Drosophila virilis, which confers an enhanced resilience to Gram-positive bacterial infection. This gene encodes an 18-mer arginine-rich peptide (termed DvirARP) with differences from its parent gene in its pattern of expression, structure and function. DvirARP specifically expresses in D. virilis female adults with a constitutive manner. It adopts a novel fold with a 310 helix and a two CXC motif-containing loop stabilized by two disulfide bridges. DvirARP exhibits no activity on the majority of microorganisms tested and only a weak activity against two Gram-positive bacteria. DvirARP knockout flies are viable and have no obvious defect in reproductivity but they are more susceptible to the DvirARP-resistant Staphylococcus aureus infection than the wild type files, which can be attributable to its ability in neutralization of the S. aureus secreted toxins. Phylogenetic distribution analysis reveals that DvirARP is restrictedly present in the Drosophila subgenus, but independent deletion variations also occur in defensins from the Sophophora subgenus, in support of the evolvability of this class of immune effectors. Our work illustrates for the first time how a duplicate resistance-mediated gene evolves an ability to increase the resilience of a subset of Drosophila species against bacterial infection.

Naturally occurring viruses of Drosophila reduce offspring number and lifespan

Drosophila remains a pre-eminent insect model system for host-virus interaction, but the host range and fitness consequences of the drosophilid virome are poorly understood. Metagenomic studies have reported approximately 200 viruses associated with Drosophilidae, but few isolates are available to characterize the Drosophila immune response, and most characterization has relied on injection and systemic infection. Here, we use a more natural infection route to characterize the fitness effects of infection and to study a wider range of viruses. We exposed laboratory Drosophila melanogaster to 23 naturally occurring viruses from wild-collected drosophilids. We recorded transmission rates along with two components of female fitness: survival and the lifetime number of adult offspring produced. Nine different viruses transmitted during contact with laboratory D. melanogaster, although for the majority, rates of transmission were less than 20%. Five virus infections led to a significant decrease in lifespan (D. melanogaster Nora virus, D. immigrans Nora virus, Muthill virus, galbut virus and Prestney Burn virus), and three led to a reduction in the total number of offspring. Our findings demonstrate the utility of the Drosophila model for community-level studies of host-virus interactions, and suggest that viral infection could be a substantial fitness burden on wild flies.

Comparative Single Cell Analysis of Transcriptional Bursting Reveals the Role of Genome Organization on de novo Transcript Origination

Spermatogenesis is a key developmental process underlying the origination of newly evolved genes. However, rapid cell type-specific transcriptomic divergence of the Drosophila germline has posed a significant technical barrier for comparative single-cell RNA-sequencing (scRNA-Seq) studies. By quantifying a surprisingly strong correlation between species-and cell type-specific divergence in three closely related Drosophila species, we apply a simple statistical procedure to identify a core set of 198 genes that are highly predictive of cell type identity while remaining robust to species-specific differences that span over 25-30 million years of evolution. We then utilize cell type classifications based on the 198-gene set to show how transcriptional divergence in cell type increases throughout spermatogenic developmental time, contrasting with traditional hourglass models of whole-organism development. With these cross-species cell type classifications, we then investigate the influence of genome organization on the molecular evolution of spermatogenesis vis-a-vis transcriptional bursting. We first demonstrate how mechanistic control of pre-meiotic transcription is achieved by altering transcriptional burst size while post-meiotic control is exerted via altered bursting frequency. We then report how global differences in autosomal vs. X chromosomal transcription likely arise in a developmental stage preceding full testis organogenesis by showing evolutionarily conserved decreases in X-linked transcription bursting kinetics in all examined somatic and germline cell types. Finally, we provide evidence supporting the cultivator model of de novo gene origination by demonstrating how the appearance of newly evolved testis-specific transcripts potentially provides short-range regulation of the transcriptional bursting properties of neighboring genes during key stages of spermatogenesis.

Extreme genome scrambling in marine planktonic Oikopleura dioica cryptic species

Genome structural variations within species are rare. How selective constraints preserve gene order and chromosome structure is a central question in evolutionary biology that remains unsolved. Our sequencing of several genomes of the appendicularian tunicate Oikopleura dioica around the globe reveals extreme genome scrambling caused by thousands of chromosomal rearrangements, although showing no obvious morphological differences between these animals. The breakpoint accumulation rate is an order of magnitude higher than in ascidian tunicates, nematodes, Drosophila, or mammals. Chromosome arms and sex-specific regions appear to be the primary unit of macrosynteny conservation. At the microsyntenic level, scrambling did not preserve operon structures, suggesting an absence of selective pressure to maintain them. The uncoupling of the genome scrambling with morphological conservation in O. dioica suggests the presence of previously unnoticed cryptic species and provides a new biological system that challenges our previous vision of speciation in which similar animals always share similar genome structures.

Pervasive gene flow despite strong and varied reproductive barriers in swordtails

One of the mechanisms that can lead to the formation of new species occurs through the evolution of reproductive barriers. However, recent research has demonstrated that hybridization has been pervasive across the tree of life even in the presence of strong barriers. Swordtail fishes (genus Xiphophorus) are an emerging model system for studying the interface between these barriers and hybridization. We document overlapping mechanisms that act as barriers between closely related species, X. birchmanni and X. cortezi, by combining genomic sequencing from natural hybrid populations, artificial crosses, behavioral assays, sperm performance, and developmental studies. We show that strong assortative mating plays a key role in maintaining subpopulations with distinct ancestry in natural hybrid populations. Lab experiments demonstrate that artificial F1 crosses experience dysfunction: crosses with X. birchmanni females were largely inviable and crosses with X. cortezi females had a heavily skewed sex ratio. Using F2 hybrids we identify several genomic regions that strongly impact hybrid viability. Strikingly, two of these regions underlie genetic incompatibilities in hybrids between X. birchmanni and its sister species X. malinche. Our results demonstrate that ancient hybridization has played a role in the origin of this shared genetic incompatibility. Moreover, ancestry mismatch at these incompatible regions has remarkably similar consequences for phenotypes and hybrid survival in X. cortezi × X. birchmanni hybrids as in X. malinche × X. birchmanni hybrids. Our findings identify varied reproductive barriers that shape genetic exchange between naturally hybridizing species and highlight the complex evolutionary outcomes of hybridization.

Hybrid speciation driven by multilocus introgression of ecological traits

Hybridization allows adaptations to be shared among lineages and may trigger the evolution of new species1,2. However, convincing examples of homoploid hybrid speciation remain rare because it is challenging to demonstrate that hybridization was crucial in generating reproductive isolation3. Here we combine population genomic analysis with quantitative trait locus mapping of species-specific traits to examine a case of hybrid speciation in Heliconius butterflies. We show that Heliconius elevatus is a hybrid species that is sympatric with both parents and has persisted as an independently evolving lineage for at least 180,000 years. This is despite pervasive and ongoing gene flow with one parent, Heliconius pardalinus, which homogenizes 99% of their genomes. The remaining 1% introgressed from the other parent, Heliconius melpomene, and is scattered widely across the H. elevatus genome in islands of divergence from H. pardalinus. These islands contain multiple traits that are under disruptive selection, including colour pattern, wing shape, host plant preference, sex pheromones and mate choice. Collectively, these traits place H. elevatus on its own adaptive peak and permit coexistence with both parents. Our results show that speciation was driven by introgression of ecological traits, and that speciation with gene flow is possible with a multilocus genetic architecture.

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.

Horizontal transfer and the widespread presence of Galileo transposons in Drosophilidae (Insecta: Diptera)

Galileo is a transposon notoriously involved with inversions in Drosophila buzzatii by ectopic recombination. Although widespread in Drosophila, little is known about this transposon in other lineages of Drosophilidae. Here, the abundance of the canonical Galileo and its evolutionary history in Drosophilidae genomes was estimated and reconstructed across genera within its two subfamilies. Sequences of this transposon were masked in these genomes and their transposase sequences were recovered using BLASTn. Phylogenetic analyses were employed to reconstruct their evolutionary history and compare it to that of host genomes. Galileo was found in nearly all 163 species, however, only 37 harbored nearly complete transposase sequences. In the remaining, Galileo was found highly fragmented. Copies from related species were clustered, however horizontal transfer events were detected between the melanogaster and montium groups of Drosophila, and between the latter and the Lordiphosa genus. The similarity of sequences found in the virilis and willistoni groups of Drosophila was found to be a consequence of lineage sorting. Therefore, the evolution of Galileo is primarily marked by vertical transmission and long-term inactivation, mainly through the deletion of open reading frames. The latter has the potential to lead copies of this transposon to become miniature inverted-repeat transposable elements.

Chromosome-level Subgenome-aware de novo Assembly of Saccharomyces bayanus Provides Insight into Genome Divergence after Hybridization

Interspecies hybridization is prevalent in various eukaryotic lineages and plays important roles in phenotypic diversification, adaption, and speciation. To better understand the changes that occurred in the different subgenomes of a hybrid species and how they facilitated adaptation, we completed chromosome-level de novo assemblies of all 16 pairs chromosomes for a recently formed hybrid yeast, Saccharomyces bayanus strain CBS380 (IFO11022), using Nanopore MinION long-read sequencing. Characterization of S. bayanus subgenomes and comparative analysis with the genomes of its parent species, S. uvarum and S. eubayanus, provide several new insights into understanding genome evolution after a relatively recent hybridization. For instance, multiple recombination events between the two subgenomes have been observed in each chromosome, followed by loss of heterozygosity (LOH) in most chromosomes in nine chromosome pairs. In addition to maintaining nearly all gene content and synteny from its parental genomes, S. bayanus has acquired many genes from other yeast species, primarily through the introgression of S. cerevisiae, such as those involved in the maltose metabolism. In addition, the patterns of recombination and LOH suggest an allotetraploid origin of S. bayanus. The gene acquisition and rapid LOH in the hybrid genome probably facilitated its adaption to maltose brewing environments and mitigated the maladaptive effect of hybridization.

The genome of the blind bee louse fly reveals deep convergences with its social host and illuminates Drosophila origins

Social insects' nests harbor intruders known as inquilines,1 which are usually related to their hosts.2,3 However, distant non-social inquilines may also show convergences with their hosts,4,5 although the underlying genomic changes remain unclear. We analyzed the genome of the wingless and blind bee louse fly Braula coeca, an inquiline kleptoparasite of the western honey bee, Apis mellifera.6,7 Using large phylogenomic data, we confirmed recent accounts that the bee louse fly is a drosophilid8,9 and showed that it had likely evolved from a sap-breeder ancestor associated with honeydew and scale insects' wax. Unlike many parasites, the bee louse fly genome did not show significant erosion or strict reliance on an endosymbiont, likely due to a relatively recent age of inquilinism. However, we observed a horizontal transfer of a transposon and a striking parallel evolution in a set of gene families between the honey bee and the bee louse fly. Convergences included genes potentially involved in metabolism and immunity and the loss of nearly all bitter-tasting gustatory receptors, in agreement with life in a protective nest and a diet of honey, pollen, and beeswax. Vision and odorant receptor genes also exhibited rapid losses. Only genes whose orthologs in the closely related Drosophila melanogaster respond to honey bee pheromone components or floral aroma were retained, whereas the losses included orthologous receptors responsive to the anti-ovarian honey bee queen pheromones. Hence, deep genomic convergences can underlie major phenotypic transitions during the evolution of inquilinism between non-social parasites and their social hosts.

A new hybrid incompatibility locus between Drosophila melanogaster and Drosophila sechellia

Despite the fundamental importance of hybrid incompatibilities to the process of speciation, there are few cases where the evolution and genetic architecture of hybrid incompatibilities are understood. One of the longest studied hybrid incompatibilities causes F1 hybrid male inviability in crosses between Drosophila melanogaster females and males from the Drosophila simulans clade of species-Drosophila simulans, Drosophila mauritiana, and Drosophila sechellia. Here, we discover dramatic differences in the manifestation of this lethal hybrid incompatibility among the D. simulans clade of species. In particular, F1 hybrid males between D. melanogaster and D. sechellia are resistant to hybrid rescue through RNAi knockdown of an essential hybrid incompatibility gene. To understand the genetic basis of this inter-species difference in hybrid rescue, we developed a triple-hybrid mapping method. Our results show that 2 discrete large effect loci and many dispersed small effect changes across the genome underlie D. sechellia aversion to hybrid rescue. The large effect loci encompass a known incompatibility gene Lethal hybrid rescue (Lhr) and previously unknown factor, Sechellia aversion to hybrid rescue (Satyr). These results show that the genetic architecture of F1 hybrid male inviability is overlapping but not identical in the 3 inter-species crosses. Our results raise questions about whether new hybrid incompatibility genes can integrate into an existing hybrid incompatibility thus increasing in complexity over time, or if the continued evolution of genes can gradually strengthen an existing hybrid incompatibility.

A comprehensive evolutionary scenario for the origin and neofunctionalization of the Drosophila speciation gene Odysseus (OdsH)

Odysseus (OdsH) was the first speciation gene described in Drosophila related to hybrid sterility in offspring of mating between Drosophila mauritiana and Drosophila simulans. Its origin is attributed to the duplication of the gene unc-4 in the subgenus Sophophora. By using a much larger sample of Drosophilidae species, we showed that contrary to what has been previously proposed, OdsH origin occurred 62 MYA. Evolutionary rates, expression, and transcription factor-binding sites of OdsH evidence that it may have rapidly experienced neofunctionalization in male sexual functions. Furthermore, the analysis of the OdsH peptide allowed the identification of mutations of D. mauritiana that could result in incompatibility in hybrids. In order to find if OdsH could be related to hybrid sterility, beyond Sophophora, we explored the expression of OdsH in Drosophila arizonae and Drosophila mojavensis, a pair of sister species with incomplete reproductive isolation. Our data indicated that OdsH expression is not atypical in their male-sterile hybrids. In conclusion, we have proposed that the origin of OdsH occurred earlier than previously proposed, followed by neofunctionalization. Our results also suggested that its role as a speciation gene might be restricted to D. mauritiana and D. simulans.

Effect of acetic acid bacteria colonization on oviposition and feeding site choice in Drosophila suzukii and its related species

Unlike many species of Drosophila flies that colonize decaying fruits, Drosophila suzukii lay eggs in ripening fruits. The oviposition and feeding site preferences for bacterial growth were quantified in multiple strains of D. suzukii and its closely related species, D. subpulchrella and D. biarmipes . A continuous degree of preference for oviposition sites with Acetobacter growth both within and across species suggested that the separation in resource usage is notable but not complete among these species. The lack of interspecific differences in feeding site preference for Acetobacter -containing media implied that the oviposition site preferences evolved independently from the feeding site preference.

Evolution of chemosensory tissues and cells across ecologically diverse Drosophilids

Chemosensory tissues exhibit significant between-species variability, yet the evolution of gene expression and cell types underlying this diversity remain poorly understood. To address these questions, we conducted transcriptomic analyses of five chemosensory tissues from six Drosophila species and integrated the findings with single-cell datasets. While stabilizing selection predominantly shapes chemosensory transcriptomes, thousands of genes in each tissue have evolved expression differences. Genes that have changed expression in one tissue have often changed in multiple other tissues but at different past epochs and are more likely to be cell type-specific than unchanged genes. Notably, chemosensory-related genes have undergone widespread expression changes, with numerous species-specific gains/losses including novel chemoreceptors expression patterns. Sex differences are also pervasive, including a D. melanogaster-specific excess of male-biased expression in sensory and muscle cells in its forelegs. Together, our analyses provide new insights for understanding evolutionary changes in chemosensory tissues at both global and individual gene levels.

Higher evolutionary dynamics of gene copy number for Drosophila glue genes located near short repeat sequences

BACKGROUND

During evolution, genes can experience duplications, losses, inversions and gene conversions. Why certain genes are more dynamic than others is poorly understood. Here we examine how several Sgs genes encoding glue proteins, which make up a bioadhesive that sticks the animal during metamorphosis, have evolved in Drosophila species.

RESULTS

We examined high-quality genome assemblies of 24 Drosophila species to study the evolutionary dynamics of four glue genes that are present in D. melanogaster and are part of the same gene family - Sgs1, Sgs3, Sgs7 and Sgs8 - across approximately 30 millions of years. We annotated a total of 102 Sgs genes and grouped them into 4 subfamilies. We present here a new nomenclature for these Sgs genes based on protein sequence conservation, genomic location and presence/absence of internal repeats. Two types of glue genes were uncovered. The first category (Sgs1, Sgs3x, Sgs3e) showed a few gene losses but no duplication, no local inversion and no gene conversion. The second group (Sgs3b, Sgs7, Sgs8) exhibited multiple events of gene losses, gene duplications, local inversions and gene conversions. Our data suggest that the presence of short "new glue" genes near the genes of the latter group may have accelerated their dynamics.

CONCLUSIONS

Our comparative analysis suggests that the evolutionary dynamics of glue genes is influenced by genomic context. Our molecular, phylogenetic and comparative analysis of the four glue genes Sgs1, Sgs3, Sgs7 and Sgs8 provides the foundation for investigating the role of the various glue genes during Drosophila life.

Aneuploidy during development in facultative parthenogenetic Drosophila

From concatenated chromosomes to polyploidization, large-scale genome changes are known to occur in parthenogenetic animals. Here, we report mosaic aneuploidy in larval brains of facultatively parthenogenetic Drosophila. We identified a background of aneuploidy in D. mercatorum strains and found increased levels of aneuploidy in the larval brain tissue of animals arising parthenogenetically versus those arising from sexual reproduction. There is also intra-individual variation in germline-derived aneuploidy within the same strain. To determine if this is a general feature of facultative parthenogenesis in drosophilids, we compared sexually reproduced and parthenogenetic offspring from an engineered facultative parthenogenetic strain of D. melanogaster. In addition to germline-derived aneuploidy, this revealed somatic aneuploidy that increased by up to fourfold in parthenogens compared to sexually reproduced offspring. Therefore, the genetic combination identified in D. mercatorum that causes facultative parthenogenesis in D. melanogaster results in aneuploidy, which indicates that the loss of mitotic control resulting in parthenogenesis causes subsequent genome variation within the parthenogenetic offspring. Our findings challenge the assumption that parthenogenetic offspring are near genetic replicas of their mothers.

Skeleton phylogeny reconstructed with transcriptomes for the tribe Drosophilini (Diptera: Drosophilidae)

The family Drosophilidae is one of the most important model systems in evolutionary biology. Thanks to advances in high-throughput sequencing technology, a number of molecular phylogenetic analyses have been undertaken by using large data sets of many genes and many species sampled across this family. Especially, recent analyses using genome sequences have depicted the family-wide skeleton phylogeny with high confidence. However, the taxon sampling is still insufficient for minor lineages and non-Drosophila genera. In this study, we carried out phylogenetic analyses using a large number of transcriptome-based nucleotide sequences, focusing on the largest, core tribe Drosophilini in the Drosophilidae. In our analyses, some noise factors against phylogenetic reconstruction were taken into account by removing putative paralogy from the datasets and examining the effects of missing data, i.e. gene occupancy and site coverage, and incomplete lineage sorting. The inferred phylogeny has newly resolved the following phylogenetic positions/relationships at the genomic scale: (i) the monophyly of the subgenus Siphlodora including Zaprionus flavofasciatus to be transferred therein; (ii) the paraphyly of the robusta and melanica species groups within a clade comprised of the robusta, melanica and quadrisetata groups and Z. flavofasciatus; (iii) Drosophila curviceps (representing the curviceps group), D. annulipes (the quadrilineata subgroup of the immigrans group) and D. maculinotata clustered into a clade sister to the Idiomyia + Scaptomyza clade, forming together the expanded Hawaiian drosophilid lineage; (iv) Dichaetophora tenuicauda (representing the lineage comprised of the Zygothrica genus group and Dichaetophora) placed as the sister to the clade of the expanded Hawaiian drosophilid lineage and Siphlodora; and (v) relationships of the subgenus Drosophila and the genus Zaprionus as follows: (Zaprionus, (the quadrilineata subgroup, ((D. sternopleuralis, the immigrans group proper), (the quinaria radiation, the tripunctata radiation)))). These results are to be incorporated into the so-far published phylogenomic tree as a backbone (constraint) tree for grafting much more species based on sequences of a limited number of genes. Such a comprehensive, highly confident phylogenetic tree with extensive and dense taxon sampling will provide an essential framework for comparative studies of the Drosophilidae.

IntroUNET: Identifying introgressed alleles via semantic segmentation

A growing body of evidence suggests that gene flow between closely related species is a widespread phenomenon. Alleles that introgress from one species into a close relative are typically neutral or deleterious, but sometimes confer a significant fitness advantage. Given the potential relevance to speciation and adaptation, numerous methods have therefore been devised to identify regions of the genome that have experienced introgression. Recently, supervised machine learning approaches have been shown to be highly effective for detecting introgression. One especially promising approach is to treat population genetic inference as an image classification problem, and feed an image representation of a population genetic alignment as input to a deep neural network that distinguishes among evolutionary models (i.e. introgression or no introgression). However, if we wish to investigate the full extent and fitness effects of introgression, merely identifying genomic regions in a population genetic alignment that harbor introgressed loci is insufficient-ideally we would be able to infer precisely which individuals have introgressed material and at which positions in the genome. Here we adapt a deep learning algorithm for semantic segmentation, the task of correctly identifying the type of object to which each individual pixel in an image belongs, to the task of identifying introgressed alleles. Our trained neural network is thus able to infer, for each individual in a two-population alignment, which of those individual's alleles were introgressed from the other population. We use simulated data to show that this approach is highly accurate, and that it can be readily extended to identify alleles that are introgressed from an unsampled "ghost" population, performing comparably to a supervised learning method tailored specifically to that task. Finally, we apply this method to data from Drosophila, showing that it is able to accurately recover introgressed haplotypes from real data. This analysis reveals that introgressed alleles are typically confined to lower frequencies within genic regions, suggestive of purifying selection, but are found at much higher frequencies in a region previously shown to be affected by adaptive introgression. Our method's success in recovering introgressed haplotypes in challenging real-world scenarios underscores the utility of deep learning approaches for making richer evolutionary inferences from genomic data.

IntroUNET: identifying introgressed alleles via semantic segmentation

A growing body of evidence suggests that gene flow between closely related species is a widespread phenomenon. Alleles that introgress from one species into a close relative are typically neutral or deleterious, but sometimes confer a significant fitness advantage. Given the potential relevance to speciation and adaptation, numerous methods have therefore been devised to identify regions of the genome that have experienced introgression. Recently, supervised machine learning approaches have been shown to be highly effective for detecting introgression. One especially promising approach is to treat population genetic inference as an image classification problem, and feed an image representation of a population genetic alignment as input to a deep neural network that distinguishes among evolutionary models (i.e. introgression or no introgression). However, if we wish to investigate the full extent and fitness effects of introgression, merely identifying genomic regions in a population genetic alignment that harbor introgressed loci is insufficient-ideally we would be able to infer precisely which individuals have introgressed material and at which positions in the genome. Here we adapt a deep learning algorithm for semantic segmentation, the task of correctly identifying the type of object to which each individual pixel in an image belongs, to the task of identifying introgressed alleles. Our trained neural network is thus able to infer, for each individual in a two-population alignment, which of those individual's alleles were introgressed from the other population. We use simulated data to show that this approach is highly accurate, and that it can be readily extended to identify alleles that are introgressed from an unsampled "ghost" population, performing comparably to a supervised learning method tailored specifically to that task. Finally, we apply this method to data from Drosophila, showing that it is able to accurately recover introgressed haplotypes from real data. This analysis reveals that introgressed alleles are typically confined to lower frequencies within genic regions, suggestive of purifying selection, but are found at much higher frequencies in a region previously shown to be affected by adaptive introgression. Our method's success in recovering introgressed haplotypes in challenging real-world scenarios underscores the utility of deep learning approaches for making richer evolutionary inferences from genomic data.