De Novo Genes.

Although the majority of annotated new genes in a given genome appear to have arisen from duplication-related mechanisms, recent studies have shown that genes can also originate de novo from ancestrally nongenic sequences. Investigating de novo-originated genes offers rich opportunities to understand the origin and functions of new genes, their regulatory mechanisms, and the associated evolutionary processes. Such studies have uncovered unexpected and intriguing facets of gene origination, offering novel perspectives on the complexity of the genome and gene evolution. In this review, we provide an overview of the research progress in this field, highlight recent advancements, identify key technical and conceptual challenges, and underscore critical questions that remain to be addressed.

Evolution and genetics of accessory gland transcriptome divergence between Drosophila melanogaster and D. simulans.

Studies of allele-specific expression in interspecific hybrids have provided important insights into gene-regulatory divergence and hybrid incompatibilities. Many such investigations in Drosophila have used transcriptome data from complex mixtures of many tissues or from gonads, however, regulatory divergence may vary widely among species, sexes, and tissues. Thus, we lack sufficiently broad sampling to be confident about the general biological principles of regulatory divergence. Here, we seek to fill some of these gaps in the literature by characterizing regulatory evolution and hybrid misexpression in a somatic male sex organ, the accessory gland, in F1 hybrids between Drosophila melanogaster and D. simulans. The accessory gland produces seminal fluid proteins, which play an important role in male and female fertility and may be subject to adaptive divergence due to male-male or male-female interactions. We find that trans differences are relatively more abundant than cis, in contrast to most of the interspecific hybrid literature, though large effect-size trans differences are rare. Seminal fluid protein genes have significantly elevated levels of expression divergence and tend to be regulated through both cis and trans divergence. We find limited misexpression (over- or underexpression relative to both parents) in this organ compared to most other Drosophila studies. As in previous studies, male-biased genes are overrepresented among misexpressed genes and are much more likely to be underexpressed. ATAC-Seq data show that chromatin accessibility is correlated with expression differences among species and hybrid allele-specific expression. This work identifies unique regulatory evolution and hybrid misexpression properties of the accessory gland and suggests the importance of tissue-specific allele-specific expression studies.

Evolution of secondary cell number and position in the Drosophila accessory gland.

In animals with internal fertilization, males transfer gametes and seminal fluid during copulation, both of which are required for successful reproduction. In Drosophila and other insects, seminal fluid is produced in the paired accessory gland (AG), the ejaculatory duct, and the ejaculatory bulb. The D. melanogaster AG has emerged as an important model system for this component of male reproductive biology. Seminal fluid proteins produced in the Drosophila AG are required for proper storage and use of sperm by the females, and are also critical for establishing and maintaining a suite of short- and long-term postcopulatory female physiological responses that promote reproductive success. The Drosophila AG is composed of two main cell types. The majority of AG cells, which are referred to as main cells, are responsible for production of many seminal fluid proteins. A minority of cells, about 4%, are referred to as secondary cells. These cells, which are restricted to the distal tip of the D. melanogaster AG, may play an especially important role in the maintenance of the long-term female post-mating response. Many studies of Drosophila AG evolution have suggested that the proteins produced in the gland evolve quickly, as does the transcriptome. Here, we investigate the evolution of secondary cell number and position in the AG in a collection of eight species spanning the entire history of the Drosophila genus. We document a heretofore underappreciated rapid evolutionary rate for both number and position of these specialized AG cells, raising several questions about the developmental, functional, and evolutionary significance of this variation.

Male-derived transcripts isolated from the mated female reproductive tract in Drosophila melanogaster.

In species with internal fertilization, sperm, and seminal fluid are transferred from male to female during mating. While both sperm and seminal fluid contain various types of molecules, including RNA, the role of most of these molecules in the coordination of fertilization or in other possible functions is poorly understood. In Drosophila, exosomes from the accessory gland, which produces seminal fluid, are transferred to females, but their potential cargoes have not been described. Moreover, while the RNA composition of sperm has been described in several mammalian species, little work on this problem has occurred in Drosophila. Here we use single nucleotide polymorphism differences between males and females from a set of highly inbred lines of D. melanogaster, and transcriptome data from the female reproductive tract, sperm, testis, and accessory gland, to investigate the potential origin, male vs female, RNA molecules isolated from 3 female reproductive tract organs, the seminal receptacle and spermatheca, which store sperm, and the parovaria, which does not. We find that mated females carry male-derived transcripts from many genes, including those that are markers of the accessory gland and known seminal fluid proteins. Our observations also support the idea that intact sperm transcripts can be isolated from the female sperm storage organs.

Identifying candidate de novo genes expressed in the somatic female reproductive tract of Drosophila melanogaster.

Most eukaryotic genes have been vertically transmitted to the present from distant ancestors. However, variable gene number across species indicates that gene gain and loss also occurs. While new genes typically originate as products of duplications and rearrangements of preexisting genes, putative de novo genes-genes born out of ancestrally nongenic sequence-have been identified. Previous studies of de novo genes in Drosophila have provided evidence that expression in male reproductive tissues is common. However, no studies have focused on female reproductive tissues. Here we begin addressing this gap in the literature by analyzing the transcriptomes of 3 female reproductive tract organs (spermatheca, seminal receptacle, and parovaria) in 3 species-our focal species, Drosophila melanogaster-and 2 closely related species, Drosophila simulans and Drosophila yakuba, with the goal of identifying putative D. melanogaster-specific de novo genes expressed in these tissues. We discovered several candidate genes, located in sequence annotated as intergenic. Consistent with the literature, these genes tend to be short, single exon, and lowly expressed. We also find evidence that some of these genes are expressed in other D. melanogaster tissues and both sexes. The relatively small number of intergenic candidate genes discovered here is similar to that observed in the accessory gland, but substantially fewer than that observed in the testis.

Identifying candidate de novo genes expressed in the somatic female reproductive tract of Drosophila melanogaster.

Most eukaryotic genes have been vertically transmitted to the present from distant ancestors. However, variable gene number across species indicates that gene gain and loss also occurs. While new genes typically originate as products of duplications and rearrangements of pre-existing genes, putative de novo genes - genes born out of previously non-genic sequence - have been identified. Previous studies of de novo genes in Drosophila have provided evidence that expression in male reproductive tissues is common. However, no studies have focused on female reproductive tissues. Here we begin addressing this gap in the literature by analyzing the transcriptomes of three female reproductive tract organs (spermatheca, seminal receptacle, and parovaria) in three species - our focal species, D. melanogaster - and two closely related species, D. simulans and D. yakuba , with the goal of identifying putative D. melanogaster -specific de novo genes expressed in these tissues. We discovered several candidate genes, which, consistent with the literature, tend to be short, simple, and lowly expressed. We also find evidence that some of these genes are expressed in other D. melanogaster tissues and both sexes. The relatively small number of candidate genes discovered here is similar to that observed in the accessory gland, but substantially fewer than that observed in the testis.

Selection and geography shape male reproductive tract transcriptomes in Drosophila melanogaster.

Transcriptome analysis of several animal clades suggests that male reproductive tract gene expression evolves quickly. However, the factors influencing the abundance and distribution of within-species variation, the ultimate source of interspecific divergence, are poorly known. Drosophila melanogaster, an ancestrally African species that has recently spread throughout the world and colonized the Americas in the last roughly 100 years, exhibits phenotypic and genetic latitudinal clines on multiple continents, consistent with a role for spatially varying selection in shaping its biology. Nevertheless, geographic expression variation in the Americas is poorly described, as is its relationship to African expression variation. Here, we investigate these issues through the analysis of two male reproductive tissue transcriptomes [testis and accessory gland (AG)] in samples from Maine (USA), Panama, and Zambia. We find dramatic differences between these tissues in differential expression between Maine and Panama, with the accessory glands exhibiting abundant expression differentiation and the testis exhibiting very little. Latitudinal expression differentiation appears to be influenced by the selection of Panama expression phenotypes. While the testis shows little latitudinal expression differentiation, it exhibits much greater differentiation than the accessory gland in Zambia vs American population comparisons. Expression differentiation for both tissues is non-randomly distributed across the genome on a chromosome arm scale. Interspecific expression divergence between D. melanogaster and D. simulans is discordant with rates of differentiation between D. melanogaster populations. Strongly heterogeneous expression differentiation across tissues and timescales suggests a complex evolutionary process involving major temporal changes in the way selection influences expression evolution in these organs.

Single-nucleus transcriptomes reveal evolutionary and functional properties of cell types in the Drosophila accessory gland.

Many traits responsible for male reproduction evolve quickly, including gene expression phenotypes in germline and somatic male reproductive tissues. Rapid male evolution in polyandrous species is thought to be driven by competition among males for fertilizations and conflicts between male and female fitness interests that manifest in postcopulatory phenotypes. In Drosophila, seminal fluid proteins secreted by three major cell types of the male accessory gland and ejaculatory duct are required for female sperm storage and use, and influence female postcopulatory traits. Recent work has shown that these cell types have overlapping but distinct effects on female postcopulatory biology, yet relatively little is known about their evolutionary properties. Here, we use single-nucleus RNA-Seq of the accessory gland and ejaculatory duct from Drosophila melanogaster and two closely related species to comprehensively describe the cell diversity of these tissues and their transcriptome evolution for the first time. We find that seminal fluid transcripts are strongly partitioned across the major cell types, and expression of many other genes additionally defines each cell type. We also report previously undocumented diversity in main cells. Transcriptome divergence was found to be heterogeneous across cell types and lineages, revealing a complex evolutionary process. Furthermore, protein adaptation varied across cell types, with potential consequences for our understanding of selection on male postcopulatory traits.

Population biology of accessory gland-expressed de novo genes in Drosophila melanogaster.

Early work on de novo gene discovery in Drosophila was consistent with the idea that many such genes have male-biased patterns of expression, including a large number expressed in the testis. However, there has been little formal analysis of variation in the abundance and properties of de novo genes expressed in different tissues. Here, we investigate the population biology of recently evolved de novo genes expressed in the Drosophila melanogaster accessory gland, a somatic male tissue that plays an important role in male and female fertility and the post mating response of females, using the same collection of inbred lines used previously to identify testis-expressed de novo genes, thus allowing for direct cross tissue comparisons of these genes in two tissues of male reproduction. Using RNA-seq data, we identify candidate de novo genes located in annotated intergenic and intronic sequence and determine the properties of these genes including chromosomal location, expression, abundance, and coding capacity. Generally, we find major differences between the tissues in terms of gene abundance and expression, though other properties such as transcript length and chromosomal distribution are more similar. We also explore differences between regulatory mechanisms of de novo genes in the two tissues and how such differences may interact with selection to produce differences in D. melanogaster de novo genes expressed in the two tissues.

Sex and tissue-specific evolution of developmental plasticity in Drosophila melanogaster.

Developmental plasticity influences the size of adult tissues in insects. Tissues can have unique responses to environmental perturbation during development; however, the prevalence of within species evolution of tissue-specific developmental plasticity remains unclear. To address this, we studied the effects of temperature and nutrition on wing and femur size in D. melanogaster populations from a temperate and tropical region. Wings were more sensitive to temperature, while wings and femurs were equally responsive to nutrition in both populations and sexes. The temperate population was larger under all conditions, except for femurs of starved females. In line with this, we observed greater femur size plasticity in response to starvation in temperate females, leading to differences in sexual dimorphism between populations such that the slope of the reaction norm of sexual dimorphism in the tropical population was double that of the temperate population. Lastly, we observed a significant trend for steeper slopes of reaction norms in temperate than in tropical females, but not in males. These findings highlight that plasticity divergence between populations can evolve heterogeneously across sexes and tissues and that nutritional plasticity can alter sexual dimorphism in D. melanogaster.

Phenotypic coupling of sleep and starvation resistance evolves in D. melanogaster.

BACKGROUND: One hypothesis for the function of sleep is that it serves as a mechanism to conserve energy. Recent studies have suggested that increased sleep can be an adaptive mechanism to improve survival under food deprivation in Drosophila melanogaster. To test the generality of this hypothesis, we compared sleep and its plastic response to starvation in a temperate and tropical population of Drosophila melanogaster. RESULTS: We found that flies from the temperate population were more starvation resistant, and hypothesized that they would engage in behaviors that are considered to conserve energy, including increased sleep and reduced movement. Surprisingly, temperate flies slept less and moved more when they were awake compared to tropical flies, both under fed and starved conditions, therefore sleep did not correlate with population-level differences in starvation resistance. In contrast, total sleep and percent change in sleep when starved were strongly positively correlated with starvation resistance within the tropical population, but not within the temperate population. Thus, we observe unexpectedly complex relationships between starvation and sleep that vary both within and across populations. These observations falsify the simple hypothesis of a straightforward relationship between sleep and energy conservation. We also tested the hypothesis that starvation is correlated with metabolic phenotypes by investigating stored lipid and carbohydrate levels, and found that stored metabolites partially contributed towards variation starvation resistance. CONCLUSIONS: Our findings demonstrate that the function of sleep under starvation can rapidly evolve on short timescales and raise new questions about the physiological correlates of sleep and the extent to which variation in sleep is shaped by natural selection.

Polymorphism and Divergence of Novel Gene Expression Patterns in Drosophila melanogaster.

Transcriptomes may evolve by multiple mechanisms, including the evolution of novel genes, the evolution of transcript abundance, and the evolution of cell, tissue, or organ expression patterns. Here, we focus on the last of these mechanisms in an investigation of tissue and organ shifts in gene expression in Drosophila melanogaster. In contrast to most investigations of expression evolution, we seek to provide a framework for understanding the mechanisms of novel expression patterns on a short population genetic timescale. To do so, we generated population samples of D. melanogaster transcriptomes from five tissues: accessory gland, testis, larval salivary gland, female head, and first-instar larva. We combined these data with comparable data from two outgroups to characterize gains and losses of expression, both polymorphic and fixed, in D. melanogaster We observed a large number of gain- or loss-of-expression phenotypes, most of which were polymorphic within D. melanogaster Several polymorphic, novel expression phenotypes were strongly influenced by segregating cis-acting variants. In support of previous literature on the evolution of novelties functioning in male reproduction, we observed many more novel expression phenotypes in the testis and accessory gland than in other tissues. Additionally, genes showing novel expression phenotypes tend to exhibit greater tissue-specific expression. Finally, in addition to qualitatively novel expression phenotypes, we identified genes exhibiting major quantitative expression divergence in the D. melanogaster lineage.

Patterns of transposable element variation and clinality in Drosophila.

Natural populations often exist in spatially diverse environments and may experience variation in the strength and targets of natural selection across their ranges. Drosophila provides an excellent opportunity to study the effects of spatially varying selection in natural populations, as both Drosophila melanogaster and Drosophila simulans live across a wide range of environments in North America. Here, we characterize patterns of variation in transposable elements (TEs) from six populations of D. melanogaster and nine populations of D. simulans sampled from multiple latitudes across North America. We find a nearly twofold excess of TEs in D. melanogaster relative to D. simulans, with this difference largely driven by TEs segregating at the lowest and highest allele frequencies. We find no effect of latitude on either total TE abundance or average TE allele frequencies in either species. Moreover, we show that, as a class of mutations, the most common patterns of TE variation do not coincide with the sampled latitudinal gradient, nor are they consistent with local adaptation acting on environmental differences found in the most extreme latitudes. We also do not find a cline in ancestry for North American D. melanogaster-for either TEs or single nucleotide polymorphisms-suggesting a limited role for demography in shaping patterns of TE variation. Though we find little evidence for widespread clinality among TEs in Drosophila, this does not necessarily imply a limited role for TEs in adaptation. We discuss the need for improved models of adaptation to large-scale environmental heterogeneity, and how these might be applied to TEs.

Functional Analysis of a Putative Target of Spatially Varying Selection in the Menin1 Gene of Drosophila melanogaster.

While significant effort has been devoted to investigating the potential influence of spatially varying selection on genomic variation, relatively little effort has been devoted to experimental analysis of putative variants or genes experiencing such selection. Previous population genetic work identified an amino acid polymorphism in the Mnn1 gene as one of the most strongly latitudinally differentiated SNPs in the genome of Drosophila melanogaster in the United States and Australia. Here we report the results of our transgenic analysis of this amino acid polymorphism. Genotypes carrying alternative Mnn1 alleles differed in multiple phenotypes in a direction generally consistent with phenotypic differences previously observed along latitudinal clines. These results support inferences from earlier population genomic work that this variant influences fitness, and support the idea that the alleles exhibiting clines may be likely to have pleiotropic effects that are correlated along the axes favored by natural selection.

Genomics of parallel adaptation at two timescales in Drosophila.

Two interesting unanswered questions are the extent to which both the broad patterns and genetic details of adaptive divergence are repeatable across species, and the timescales over which parallel adaptation may be observed. Drosophila melanogaster is a key model system for population and evolutionary genomics. Findings from genetics and genomics suggest that recent adaptation to latitudinal environmental variation (on the timescale of hundreds or thousands of years) associated with Out-of-Africa colonization plays an important role in maintaining biological variation in the species. Additionally, studies of interspecific differences between D. melanogaster and its sister species D. simulans have revealed that a substantial proportion of proteins and amino acid residues exhibit adaptive divergence on a roughly few million years long timescale. Here we use population genomic approaches to attack the problem of parallelism between D. melanogaster and a highly diverged conger, D. hydei, on two timescales. D. hydei, a member of the repleta group of Drosophila, is similar to D. melanogaster, in that it too appears to be a recently cosmopolitan species and recent colonizer of high latitude environments. We observed parallelism both for genes exhibiting latitudinal allele frequency differentiation within species and for genes exhibiting recurrent adaptive protein divergence between species. Greater parallelism was observed for long-term adaptive protein evolution and this parallelism includes not only the specific genes/proteins that exhibit adaptive evolution, but extends even to the magnitudes of the selective effects on interspecific protein differences. Thus, despite the roughly 50 million years of time separating D. melanogaster and D. hydei, and despite their considerably divergent biology, they exhibit substantial parallelism, suggesting the existence of a fundamental predictability of adaptive evolution in the genus.

The Adaptive Significance of Natural Genetic Variation in the DNA Damage Response of Drosophila melanogaster.

Despite decades of work, our understanding of the distribution of fitness effects of segregating genetic variants in natural populations remains largely incomplete. One form of selection that can maintain genetic variation is spatially varying selection, such as that leading to latitudinal clines. While the introduction of population genomic approaches to understanding spatially varying selection has generated much excitement, little successful effort has been devoted to moving beyond genome scans for selection to experimental analysis of the relevant biology and the development of experimentally motivated hypotheses regarding the agents of selection; it remains an interesting question as to whether the vast majority of population genomic work will lead to satisfying biological insights. Here, motivated by population genomic results, we investigate how spatially varying selection in the genetic model system, Drosophila melanogaster, has led to genetic differences between populations in several components of the DNA damage response. UVB incidence, which is negatively correlated with latitude, is an important agent of DNA damage. We show that sensitivity of early embryos to UVB exposure is strongly correlated with latitude such that low latitude populations show much lower sensitivity to UVB. We then show that lines with lower embryo UVB sensitivity also exhibit increased capacity for repair of damaged sperm DNA by the oocyte. A comparison of the early embryo transcriptome in high and low latitude embryos provides evidence that one mechanism of adaptive DNA repair differences between populations is the greater abundance of DNA repair transcripts in the eggs of low latitude females. Finally, we use population genomic comparisons of high and low latitude samples to reveal evidence that multiple components of the DNA damage response and both coding and non-coding variation likely contribute to adaptive differences in DNA repair between populations.

Parallel Evolution of Copy-Number Variation across Continents in Drosophila melanogaster.

Genetic differentiation across populations that is maintained in the presence of gene flow is a hallmark of spatially varying selection. In Drosophila melanogaster, the latitudinal clines across the eastern coasts of Australia and North America appear to be examples of this type of selection, with recent studies showing that a substantial portion of the D. melanogaster genome exhibits allele frequency differentiation with respect to latitude on both continents. As of yet there has been no genome-wide examination of differentiated copy-number variants (CNVs) in these geographic regions, despite their potential importance for phenotypic variation in Drosophila and other taxa. Here, we present an analysis of geographic variation in CNVs in D. melanogaster. We also present the first genomic analysis of geographic variation for copy-number variation in the sister species, D. simulans, in order to investigate patterns of parallel evolution in these close relatives. In D. melanogaster we find hundreds of CNVs, many of which show parallel patterns of geographic variation on both continents, lending support to the idea that they are influenced by spatially varying selection. These findings support the idea that polymorphic CNVs contribute to local adaptation in D. melanogaster In contrast, we find very few CNVs in D. simulans that are geographically differentiated in parallel on both continents, consistent with earlier work suggesting that clinal patterns are weaker in this species.

Genomic Patterns of Geographic Differentiation in Drosophila simulans.

Geographic patterns of genetic differentiation have long been used to understand population history and to learn about the biological mechanisms of adaptation. Here we present an examination of genomic patterns of differentiation between northern and southern populations of Australian and North American Drosophila simulans, with an emphasis on characterizing signals of parallel differentiation. We report on the genomic scale of differentiation and functional enrichment of outlier SNPs. While, overall, signals of shared differentiation are modest, we find the strongest support for parallel differentiation in genomic regions that are associated with regulation. Comparisons to Drosophila melanogaster yield potential candidate genes involved in local adaptation in both species, providing insight into common selective pressures and responses. In contrast to D. melanogaster, in D. simulans we observe patterns of variation that are inconsistent with a model of temperate adaptation out of a tropical ancestral range, highlighting potential differences in demographic and colonization histories of this cosmopolitan species pair.

Parallel Gene Expression Differences between Low and High Latitude Populations of Drosophila melanogaster and D. simulans.

Gene expression variation within species is relatively common, however, the role of natural selection in the maintenance of this variation is poorly understood. Here we investigate low and high latitude populations of Drosophila melanogaster and its sister species, D. simulans, to determine whether the two species show similar patterns of population differentiation, consistent with a role for spatially varying selection in maintaining gene expression variation. We compared at two temperatures the whole male transcriptome of D. melanogaster and D. simulans sampled from Panama City (Panama) and Maine (USA). We observed a significant excess of genes exhibiting differential expression in both species, consistent with parallel adaptation to heterogeneous environments. Moreover, the majority of genes showing parallel expression differentiation showed the same direction of differential expression in the two species and the magnitudes of expression differences between high and low latitude populations were correlated across species, further bolstering the conclusion that parallelism for expression phenotypes results from spatially varying selection. However, the species also exhibited important differences in expression phenotypes. For example, the genomic extent of genotype × environment interaction was much more common in D. melanogaster. Highly differentiated SNPs between low and high latitudes were enriched in the 3' UTRs and CDS of the geographically differently expressed genes in both species, consistent with an important role for cis-acting variants in driving local adaptation for expression-related phenotypes.

Evidence that natural selection maintains genetic variation for sleep in Drosophila melanogaster.

BACKGROUND: Drosophila melanogaster often shows correlations between latitude and phenotypic or genetic variation on different continents, which suggests local adaptation with respect to a heterogeneous environment. Previous phenotypic analyses of latitudinal clines have investigated mainly physiological, morphological, or life-history traits. Here, we studied latitudinal variation in sleep in D. melanogaster populations from North and Central America. In parallel, we used RNA-seq to identify interpopulation gene expression differences. RESULTS: We found that in D. melanogaster the average nighttime sleep bout duration exhibits a latitudinal cline such that sleep bouts of equatorial populations are roughly twice as long as those of temperate populations. Interestingly, this pattern of latitudinal variation is not observed for any daytime measure of activity or sleep. We also found evidence for geographic variation for sunrise anticipation. Our RNA-seq experiment carried out on heads from a low and high latitude population identified a large number of gene expression differences, most of which were time dependent. Differentially expressed genes were enriched in circadian regulated genes and enriched in genes potentially under spatially varying selection. CONCLUSION: Our results are consistent with a mechanistic and selective decoupling of nighttime and daytime activity. Furthermore, the present study suggests that natural selection plays a major role in generating transcriptomic variation associated with circadian behaviors. Finally, we identified genomic variants plausibly causally associated with the observed behavioral and transcriptomic variation.