Polymorphism and Divergence of Novel Gene Expression Patterns in Drosophila melanogaster

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 new 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. After validating these cross-species cell type classifications using RNA fluorescence in situ hybridization (FISH) and imaging, we then investigate the influence of genome organization on the molecular evolution of spermatogenesis vis-a-vis transcriptional bursting. We first show altering transcriptional burst size contributes to pre-meiotic transcription and altering bursting frequency contributes to post-meiotic expression. We then report global differences in autosomal vs. X chromosomal transcription may 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 neighboring genes' transcriptional bursting properties during key stages of spermatogenesis.

Transcriptome analysis of different stages of testis development in Portunus trituberculatus

The swimming crab (Portunus trituberculatus) is an important marine economic species, however its artificial breeding yield is relatively low. Currently, the main challenge faced by the swimming crab seed industry is the reliance on wild populations for seed cultivation, which results in unstable yield and quality, affecting the healthy development of the crab farming industry to some extent. The quality of germplasm resources depends on the quality of gametes, and the quality of sperm depends on the orderly genetic regulation process of spermatogenesis. Therefore, elucidating the genetic regulatory mechanisms of spermatogenesis is of great significance for improving the germplasm resources of P. trituberculatus. To gain a deeper understanding of this process, we conducted a comparative transcriptome study on the testis of the swimming crab at different developmental stages. This study aims to identify key genes that regulate testicular development. We performed paraffin section identification on the testicular tissue of male crabs and conducted transcriptome analysis on the testicular tissue at five different developmental stages and somatic cells. Through differential expression analysis, we screened a total of 31,788 differentially expressed genes (DEGs) from stages I to VI. Through Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, we found that these DEGs were significantly enriched in 15 pathways, including important functional pathways such as the adrenergic signaling pathway, HIF-1 signaling pathway, and TGF-β signaling pathway. GO analysis results showed that calcium ion homeostasis and cell skeleton-related activities were significantly enriched in stage II. Further protein-protein interaction network analysis revealed 68 hub genes, including 13 eukaryotic initiation factors, 6 Ras superfamily members, and 6 genes related to cell division. In addition, genes such as Actin, Myosin, and Nup50 consistently showed high expression at all developmental stages, while genes related to calcium ion homeostasis, such as CaM, significantly increased in expression during stage II. Hsp90 and apoptosis-related genes had higher expression in stage IV, while Smad4 had higher expression in stage V. These results suggest that stage II of the swimming crab sperm development may be a critical period for spermatogenesis, and stage IV may be an important period for regulating sperm quality and quantity. This study not only provides a foundation for further research on the molecular mechanisms of testicular development and spermatogenesis in the swimming crab but also offers theoretical support for improving breeding yield, which has significant practical application value.

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 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.

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.

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.

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.

Transcriptomic and functional genetic evidence for distinct ecophysiological responses across complex life cycle stages

Organisms with complex life cycles demonstrate a remarkable ability to change their phenotypes across development, presumably as an evolutionary adaptation to developmentally variable environments. Developmental variation in environmentally sensitive performance, and thermal sensitivity in particular, has been well documented in holometabolous insects. For example, thermal performance in adults and juvenile stages exhibit little genetic correlation (genetic decoupling) and can evolve independently, resulting in divergent thermal responses. Yet, we understand very little about how this genetic decoupling occurs. We tested the hypothesis that genetic decoupling of thermal physiology is driven by fundamental differences in physiology between life stages, despite a potentially conserved Cellular Stress Response. We used RNAseq to compare transcript expression in response to a cold stressor in Drosophila melanogaster larvae and adults and used RNAi (RNA interference) to test whether knocking down nine target genes differentially affected larval and adult cold tolerance. Transcriptomic responses of whole larvae and adults during and following exposure to -5°C were largely unique both in identity of responding transcripts and in temporal dynamics. Further, we analyzed the tissue-specificity of differentially-expressed transcripts from FlyAtlas 2 data, and concluded that stage-specific differences in transcription were not simply driven by differences in tissue composition. In addition, RNAi of target genes resulted in largely stage-specific and sometimes sex-specific effects on cold tolerance. The combined evidence suggests that thermal physiology is largely stage-specific at the level of gene expression, and thus natural selection may be acting on different loci during the independent thermal adaptation of different life stages.

Research gaps and new insights in the evolution of Drosophila seminal fluid proteins

While the striking effects of seminal fluid proteins (SFPs) on females are fairly conserved among Diptera, most SFPs lack detectable homologues among the SFP repertoires of phylogenetically distant species. How such a rapidly changing proteome conserves functions across taxa is a fascinating question. However, this and other pivotal aspects of SFPs' evolution remain elusive because discoveries on these proteins have been mainly restricted to the model Drosophila melanogaster. Here, we provide an overview of the current knowledge on the inter-specific divergence of the SFP repertoire in Drosophila and compile the increasing amount of relevant genomic information from multiple species. Capitalizing on the accumulated knowledge in D. melanogaster, we present novel sets of high-confidence SFP candidates and transcription factors presumptively involved in regulating the expression of SFPs. We also address open questions by performing comparative genomic analyses that failed to support the existence of many conserved SFPs shared by most dipterans and indicated that gene co-option is the most frequent mechanism accounting for the origin of Drosophila SFP-coding genes. We hope our update establishes a starting point to integrate further data and thus widen the understanding of the intricate evolution of these proteins.

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.

Molecular Evolution across Mouse Spermatogenesis

Genes involved in spermatogenesis tend to evolve rapidly, but we lack a clear understanding of how protein sequences and patterns of gene expression evolve across this complex developmental process. We used fluorescence-activated cell sorting (FACS) to generate expression data for early (meiotic) and late (postmeiotic) cell types across 13 inbred strains of mice (Mus) spanning ∼7 My of evolution. We used these comparative developmental data to investigate the evolution of lineage-specific expression, protein-coding sequences, and expression levels. We found increased lineage specificity and more rapid protein-coding and expression divergence during late spermatogenesis, suggesting that signatures of rapid testis molecular evolution are punctuated across sperm development. Despite strong overall developmental parallels in these components of molecular evolution, protein and expression divergences were only weakly correlated across genes. We detected more rapid protein evolution on the X chromosome relative to the autosomes, whereas X-linked gene expression tended to be relatively more conserved likely reflecting chromosome-specific regulatory constraints. Using allele-specific FACS expression data from crosses between four strains, we found that the relative contributions of different regulatory mechanisms also differed between cell types. Genes showing cis-regulatory changes were more common late in spermatogenesis, and tended to be associated with larger differences in expression levels and greater expression divergence between species. In contrast, genes with trans-acting changes were more common early and tended to be more conserved across species. Our findings advance understanding of gene evolution across spermatogenesis and underscore the fundamental importance of developmental context in molecular evolutionary studies.

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.

Adaptive mechanisms of plant specialized metabolism connecting chemistry to function

As sessile organisms, plants evolved elaborate metabolic systems that produce a plethora of specialized metabolites as a means to survive challenging terrestrial environments. Decades of research have revealed the genetic and biochemical basis for a multitude of plant specialized metabolic pathways. Nevertheless, knowledge is still limited concerning the selective advantages provided by individual and collective specialized metabolites to the reproductive success of diverse host plants. Here we review the biological functions conferred by various classes of plant specialized metabolites in the context of the interaction of plants with their surrounding environment. To achieve optimal multifunctionality of diverse specialized metabolic processes, plants use various adaptive mechanisms at subcellular, cellular, tissue, organ and interspecies levels. Understanding these mechanisms and the evolutionary trajectories underlying their occurrence in nature will ultimately enable efficient bioengineering of desirable metabolic traits in chassis organisms.