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

Cell cycle variants during Drosophila male accessory gland development

The Drosophila melanogaster male accessory gland (AG) is a functional analog of the mammalian prostate and seminal vesicles containing two secretory epithelial cell types, termed main and secondary cells. This tissue is responsible for making and secreting seminal fluid proteins and other molecules that contribute to successful reproduction. The cells of this tissue are binucleate and polyploid, due to variant cell cycles that include endomitosis and endocycling during metamorphosis. Here, we provide evidence of additional cell cycle variants in this tissue. We show that main cells of the gland are connected by ring canals that form after the penultimate mitosis, and we describe an additional post-eclosion endocycle required for gland maturation that is dependent on juvenile hormone signaling. We present evidence that the main cells of the D. melanogaster AG undergo a unique cell cycle reprogramming throughout organ development that results in step-wise cell cycle truncations culminating in cells containing two octoploid nuclei with under-replicated heterochromatin in the mature gland. We propose this tissue as a model to study developmental and hormonal temporal control of cell cycle variants in terminally differentiating tissues.

Mating system variation and gene expression in the male reproductive tract of Peromyscus mice

Genes involved in reproduction often evolve rapidly at the sequence level due to postcopulatory sexual selection (PCSS) driven by male-male competition and male-female sexual conflict, but the impact of PCSS on gene expression has been under-explored. Further, though multiple tissues contribute to male reproductive success, most studies have focused on the testes. To explore the influence of mating system variation on reproductive tract gene expression in natural populations, we captured adult males from monogamous Peromyscus californicus and polygynandrous P. boylii and P. maniculatus. We generated RNAseq libraries, quantified gene expression in the testis, seminal vesicle, epididymis, and liver, and identified 3627 mating system-associated differentially expressed genes (MS-DEGs), where expression shifted in the same direction in P. maniculatus and P. boylii relative to P. californicus. Gene expression variation was most strongly associated with mating behaviour in the seminal vesicles, where 89% of differentially expressed genes were MS-DEGs, including the key seminal fluid proteins Svs2 and Pate4. We also used published rodent genomes to test for positive and relaxed selection on Peromyscus-expressed genes. Though we did not observe more overlap than expected by chance between MS-DEGs and positively selected genes, 203 MS-DEGs showed evidence of positive selection. Fourteen reproductive genes were under tree-wide positive selection but convergent relaxed selection in P. californicus and Microtus ochrogaster, a distantly related monogamous species. Changes in transcript abundance and gene sequence evolution in association with mating behaviour suggest that male mice may respond to sexual selection intensity by altering aspects of sperm motility, sperm-egg binding and copulatory plug formation.

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.

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

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

The seminal odorant binding protein Obp56g is required for mating plug formation and male fertility in Drosophila melanogaster

In Drosophila melanogaster and other insects, the seminal fluid proteins (SFPs) and male sex pheromones that enter the female with sperm during mating are essential for fertility and induce profound post-mating effects on female physiology. The SFPs in D. melanogaster and other taxa include several members of the large gene family known as odorant binding proteins (Obps). Work in Drosophila has shown that some Obp genes are highly expressed in the antennae and can mediate behavioral responses to odorants, potentially by binding and carrying these molecules to odorant receptors. These observations have led to the hypothesis that the seminal Obps might act as molecular carriers for pheromones or other compounds important for male fertility, though functional evidence in any species is lacking. Here, we used functional genetics to test the role of the seven seminal Obps in D. melanogaster fertility and the post-mating response (PMR). We found that Obp56g is required for male fertility and the induction of the PMR, whereas the other six genes are dispensable. We found males lacking Obp56g fail to form a mating plug in the mated female's reproductive tract, leading to ejaculate loss and reduced sperm storage, likely due to its expression in the male ejaculatory bulb. We also examined the evolutionary history of these seminal Obp genes, as several studies have documented rapid evolution and turnover of SFP genes across taxa. We found extensive lability in gene copy number and evidence of positive selection acting on two genes, Obp22a and Obp51a. Comparative RNAseq data from the male reproductive tract of multiple Drosophila species revealed that Obp56g shows high male reproductive tract expression in a subset of taxa, though conserved head expression across the phylogeny. Together, these functional and expression data suggest that Obp56g may have been co-opted for a reproductive function over evolutionary time.

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

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

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.

The seminal odorant binding protein Obp56g is required for mating plug formation and male fertility in Drosophila melanogaster

In Drosophila melanogaster and other insects, the seminal fluid proteins (SFPs) and male sex pheromones that enter the female with sperm during mating are essential for fertility and induce profound post-mating effects on female physiology and behavior. The SFPs in D. melanogaster and other taxa include several members of the large gene family known as odorant binding proteins (Obps). Previous work in Drosophila has shown that some Obp genes are highly expressed in the antennae and can mediate behavioral responses to odorants, potentially by binding and carrying these molecules to odorant receptors. These observations have led to the hypothesis that the seminal Obps might act as molecular carriers for pheromones or other compounds important for male fertility in the ejaculate, though functional evidence in any species is lacking. Here, we used RNAi and CRISPR/Cas9 generated mutants to test the role of the seven seminal Obps in D. melanogaster fertility and the post-mating response (PMR). We found that Obp56g is required for male fertility and the induction of the PMR, whereas the other six genes had no effect on fertility when mutated individually. Obp56g is expressed in the male's ejaculatory bulb, an important tissue in the reproductive tract that synthesizes components of the mating plug. We found males lacking Obp56g fail to form a mating plug in the mated female's reproductive tract, leading to ejaculate loss and reduced sperm storage. We also examined the evolutionary history of these seminal Obp genes, as several studies have documented rapid evolution and turnover of SFP genes across taxa. We found extensive lability in gene copy number and evidence of positive selection acting on two genes, Obp22a and Obp51a. Comparative RNAseq data from the male reproductive tract of multiple Drosophila species revealed that Obp56g shows high male reproductive tract expression only in species of the melanogaster and obscura groups, though conserved head expression in all species tested. Together, these functional and expression data suggest that Obp56g may have been co-opted for a reproductive function over evolutionary time.

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.