A survey of ovary-, testis-, and soma-biased gene expression in Drosophila melanogaster adults

Testis-specific H2B.W1 disrupts nucleosome integrity by reducing DNA-histone interactions

Multiple testis-specific histone variants are involved in the dynamic chromatin transitions during spermatogenesis. H2B.W1 (previously called H2BFWT) is an H2B variant specific to primate testis with hitherto unclear functions, although its single-nucleotide polymorphisms (SNPs) are closely associated with male non-obstructive infertility. Here, we found that H2B.W1 is only expressed in the mid-late spermatogonia stages, and H2B.W1 nucleosomes are defined by a more flexible structure originating from weakened interactions between histones and DNA. Furthermore, one of its SNPs, H2B.W1-H100R, which is associated with infertility, further destabilizes the nucleosomes and increases the nucleosome unwrapping rate by interfering with the R100 and H4 K91/R92 interaction. Our results suggest that destabilizing H2B.W1 containing nucleosomes might change the chromatin structure of spermatogonia, and that H2B.W1-H100R enhances the nucleosome-destabilizing effects, leading to infertility.

Deficiency of ValRS-m Causes Male Infertility in Drosophila melanogaster

Drosophila spermatogenesis involves the renewal of germline stem cells, meiosis of spermatocytes, and morphological transformation of spermatids into mature sperm. We previously demonstrated that Ocnus (ocn) plays an essential role in spermatogenesis. The ValRS-m (Valyl-tRNA synthetase, mitochondrial) gene was down-regulated in ocn RNAi testes. Here, we found that ValRS-m-knockdown induced complete sterility in male flies. The depletion of ValRS-m blocked mitochondrial behavior and ATP synthesis, thus inhibiting the transition from spermatogonia to spermatocytes, and eventually, inducing the accumulation of spermatogonia during spermatogenesis. To understand the intrinsic reason for this, we further conducted transcriptome-sequencing analysis for control and ValRS-m-knockdown testes. The differentially expressed genes (DEGs) between these two groups were selected with a fold change of ≥2 or ≤1/2. Compared with the control group, 4725 genes were down-regulated (dDEGs) and 2985 genes were up-regulated (uDEGs) in the ValRS-m RNAi group. The dDEGs were mainly concentrated in the glycolytic pathway and pyruvate metabolic pathway, and the uDEGs were primarily related to ribosomal biogenesis. A total of 28 DEGs associated with mitochondria and 6 meiosis-related genes were verified to be suppressed when ValRS-m was deficient. Overall, these results suggest that ValRS-m plays a wide and vital role in mitochondrial behavior and spermatogonia differentiation in Drosophila.

Evolution of sex-biased genes in Drosophila species with neo-sex chromosomes: Potential contribution to reducing the sexual conflict

An advantage of sex chromosomes may be the potential to reduce sexual conflict because they provide a basis for selection to operate separately on females and males. However, evaluating the relationship between sex chromosomes and sexual conflict is challenging owing to the difficulty in measuring sexual conflict and substantial divergence between species with and without sex chromosomes. We therefore examined sex-biased gene expression as a proxy for sexual conflict in three sets of Drosophila species with and without young sex chromosomes, the so-called neo-sex chromosomes. In all sets, we detected more sex-biased genes in the species with neo-sex chromosomes than in the species without neo-sex chromosomes in larvae, pupae, and adult somatic tissues but not in gonads. In particular, many unbiased genes became either female- or male-biased after linkage to the neo-sex chromosomes in larvae, despite the low sexual dimorphism. For example, genes involved in metabolism, a key determinant for the rate of development in many animals, were enriched in the genes that acquired sex-biased expression on the neo-sex chromosomes at the larval stage. These genes may be targets of sexually antagonistic selection (i.e., large size and rapid development are selected for in females but selected against in males). These results indicate that acquiring neo-sex chromosomes may have contributed to a reduction in sexual conflict, particularly at the larval stage, in Drosophila..

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.

Origin of functional de novo genes in humans from "hopeful monsters"

For a long time, it was believed that new genes arise only from modifications of preexisting genes, but the discovery of de novo protein-coding genes that originated from noncoding DNA regions demonstrates the existence of a "motherless" origination process for new genes. However, the features, distributions, expression profiles, and origin modes of these genes in humans seem to support the notion that their origin is not a purely "motherless" process; rather, these genes arise preferentially from genomic regions encoding preexisting precursors with gene-like features. In such a case, the gene loci are typically not brand new. In this short review, we will summarize the definition and features of human de novo genes and clarify their process of origination from ancestral non-coding genomic regions. In addition, we define the favored precursors, or "hopeful monsters," for the origin of de novo genes and present a discussion of the functional significance of these young genes in brain development and tumorigenesis in humans. This article is categorized under: RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution.

Sex matters in preclinical research

International Women's Day 2024 has a theme of inclusion. As publishers of preclinical research, we aim to show how inclusion of females in research advances scientific rigor and improves treatment reliability. Sexual reproduction is key to all life across the plant and animal kingdoms. Biological sex takes many forms that are morphologically differentiated during development: stamens versus pistils in plants; color and plumage in birds; fallopian tubes versus vas deferens in mammals; and differences in size, for instance, males are smaller in the fruit fly Drosophila melanogaster. Physical differences may be obvious, but many traits may be more obscure, including hormonal, physiological and metabolic factors. These traits have a big influence on disease and responses to treatment. Thus, we call for improved inclusion, analysis and reporting of sex as a biological variable in preclinical animal modeling research.

Integrating massive RNA-seq data to elucidate transcriptome dynamics in Drosophila melanogaster

The volume of ribonucleic acid (RNA)-seq data has increased exponentially, providing numerous new insights into various biological processes. However, due to significant practical challenges, such as data heterogeneity, it is still difficult to ensure the quality of these data when integrated. Although some quality control methods have been developed, sample consistency is rarely considered and these methods are susceptible to artificial factors. Here, we developed MassiveQC, an unsupervised machine learning-based approach, to automatically download and filter large-scale high-throughput data. In addition to the read quality used in other tools, MassiveQC also uses the alignment and expression quality as model features. Meanwhile, it is user-friendly since the cutoff is generated from self-reporting and is applicable to multimodal data. To explore its value, we applied MassiveQC to Drosophila RNA-seq data and generated a comprehensive transcriptome atlas across 28 tissues from embryogenesis to adulthood. We systematically characterized fly gene expression dynamics and found that genes with high expression dynamics were likely to be evolutionarily young and expressed at late developmental stages, exhibiting high nonsynonymous substitution rates and low phenotypic severity, and they were involved in simple regulatory programs. We also discovered that human and Drosophila had strong positive correlations in gene expression in orthologous organs, revealing the great potential of the Drosophila system for studying human development and disease.

Sex-biased gene expression at single-cell resolution: cause and consequence of sexual dimorphism

Gene expression differences between males and females are thought to be key for the evolution of sexual dimorphism, and sex-biased genes are often used to study the molecular footprint of sex-specific selection. However, gene expression is often measured from complex aggregations of diverse cell types, making it difficult to distinguish between sex differences in expression that are due to regulatory rewiring within similar cell types and those that are simply a consequence of developmental differences in cell-type abundance. To determine the role of regulatory versus developmental differences underlying sex-biased gene expression, we use single-cell transcriptomic data from multiple somatic and reproductive tissues of male and female guppies, a species that exhibits extensive phenotypic sexual dimorphism. Our analysis of gene expression at single-cell resolution demonstrates that nonisometric scaling between the cell populations within each tissue and heterogeneity in cell-type abundance between the sexes can influence inferred patterns of sex-biased gene expression by increasing both the false-positive and false-negative rates. Moreover, we show that, at the bulk level, the subset of sex-biased genes that are the product of sex differences in cell-type abundance can significantly confound patterns of coding-sequence evolution. Taken together, our results offer a unique insight into the effects of allometry and cellular heterogeneity on perceived patterns of sex-biased gene expression and highlight the power of single-cell RNA-sequencing in distinguishing between sex-biased genes that are the result of regulatory change and those that stem from sex differences in cell-type abundance, and hence are a consequence rather than a cause of sexual dimorphism.

Sex-biased gene expression in nutrient-sensing pathways

Differences in lifespan between males and females are found across many taxa and may be determined, at least in part, by differential responses to diet. Here we tested the hypothesis that the higher dietary sensitivity of female lifespan is mediated by higher and more dynamic expression in nutrient-sensing pathways in females. We first reanalysed existing RNA-seq data, focusing on 17 nutrient-sensing genes with reported lifespan effects. This revealed, consistent with the hypothesis, a dominant pattern of female-biased gene expression, and among sex-biased genes there tended to be a loss of female-bias after mating. We then tested directly the expression of these 17 nutrient-sensing genes in wild-type third instar larvae, once-mated 5- and 16-day-old adults. This confirmed sex-biased gene expression and showed that it was generally absent in larvae, but frequent and stable in adults. Overall, the findings suggest a proximate explanation for the sensitivity of female lifespan to dietary manipulations. We suggest that the contrasting selective pressures to which males and females are subject create differing nutritional demands and requirements, resulting in sex differences in lifespan. This underscores the potential importance of the health impacts of sex-specific dietary responses.

The evolution of ovary-biased gene expression in Hawaiian Drosophila

With detailed data on gene expression accessible from an increasingly broad array of species, we can test the extent to which our developmental genetic knowledge from model organisms predicts expression patterns and variation across species. But to know when differences in gene expression across species are significant, we first need to know how much evolutionary variation in gene expression we expect to observe. Here we provide an answer by analyzing RNAseq data across twelve species of Hawaiian Drosophilidae flies, focusing on gene expression differences between the ovary and other tissues. We show that over evolutionary time, there exists a cohort of ovary specific genes that is stable and that largely corresponds to described expression patterns from laboratory model Drosophila species. Our results also provide a demonstration of the prediction that, as phylogenetic distance increases, variation between species overwhelms variation between tissue types. Using ancestral state reconstruction of expression, we describe the distribution of evolutionary changes in tissue-biased expression, and use this to identify gains and losses of ovary-biased expression across these twelve species. We then use this distribution to calculate the evolutionary correlation in expression changes between genes, and demonstrate that genes with known interactions in D. melanogaster are significantly more correlated in their evolution than genes with no or unknown interactions. Finally, we use this correlation matrix to infer new networks of genes that share evolutionary trajectories, and we present these results as a dataset of new testable hypotheses about genetic roles and interactions in the function and evolution of the Drosophila ovary.

Stochastic Modeling of Gene Expression Evolution Uncovers Tissue- and Sex-Specific Properties of Expression Evolution in the Drosophila Genus

Gene expression evolution is typically modeled with the stochastic Ornstein-Uhlenbeck (OU) process. It has been suggested that the estimation of within-species variations using replicated data can increase the predictive power of such models, but this hypothesis has not been fully tested. We developed EvoGeneX, a computationally efficient implementation of the OU-based method that models within-species variation. Using extensive simulations, we show that modeling within-species variations and appropriate selection of species improve the performance of the model. Further, to facilitate a comparative analysis of expression evolution, we introduce a formal measure of evolutionary expression divergence for a group of genes using the rate and the asymptotic level of divergence. With these tools in hand, we performed the first-ever analysis of the evolution of gene expression across different body-parts, species, and sexes of the Drosophila genus. We observed that genes with adaptive expression evolution tend to be body-part specific, whereas the genes with constrained evolution tend to be shared across body-parts. Among the neutrally evolving gene expression patterns, gonads in both sexes have higher expression divergence relative to other tissues and the male gonads have even higher divergence than the female gonads. Among the evolutionarily constrained genes, the gonads show different divergence patterns, where the male gonads, and not the female gonads, show less constrained divergence than other body-parts. Finally, we show interesting examples of adaptive expression evolution, including adaptation of odor binding proteins.

Comprehensive Transcriptome Analysis Reveals Sex-Specific Alternative Splicing Events in Zebrafish Gonads

Alternative splicing is an important way of regulating gene functions in eukaryotes. Several key genes involved in sex determination and gonadal differentiation, such as nr5a1 and ddx4, have sex-biased transcripts between males and females, suggesting a potential regulatory role of alternative splicing in gonads. Currently, the sex-specific alternative splicing events and genes have not been comprehensively studied at the genome-wide level in zebrafish. In this study, through global splicing analysis on three independent sets of RNA-seq data from matched zebrafish testes and ovaries, we identified 120 differentially spliced genes shared by the three datasets, most of which haven’t been reported before. Functional enrichment analysis showed that the GO terms of mRNA processing, mRNA metabolism and microtubule-based process were strongly enriched. The testis- and ovary-biased alternative splicing genes were identified, and part of them (tp53bp1, tpx2, mapre1a, kif2c, and ncoa5) were further validated by RT-PCR. Sequence characteristics analysis suggested that the lengths, GC contents, and splice site strengths of the alternative exons or introns may have different influences in different types of alternative splicing events. Interestingly, we identified an unexpected high proportion (over 70%) of non-frameshift exon-skipping events, suggesting that in these cases the two protein isoforms derived from alternative splicing may both have functions. Furthermore, as a representative example, we found that the alternative splicing of ncoa5 causes the loss of a conserved RRM domain in the short transcript predominantly produced in testes. Our study discovers novel sex-specific alternative splicing events and genes with high reliabilities in zebrafish testes and ovaries, which would provide attractive targets for follow-up studies to reveal the biological significances of alternative splicing events and genes in sex determination and gonadal differentiation.

The chromosomal distribution of sex-biased microRNAs in Drosophila is non-adaptive

Genes are often differentially expressed between males and females. In Drosophila melanogaster, the analysis of sex-biased microRNAs (short noncoding regulatory molecules) has revealed striking differences with protein-coding genes. Mainly, the X chromosome is enriched in male-biased microRNA genes, although it is depleted of male-biased protein-coding genes. The paucity of male-biased genes in the X chromosome is generally explained by an evolutionary process called demasculinization. I suggest that the excess of male-biased microRNAs in the X chromosome is due to high rates of de novo emergence of microRNAs (mostly in other neighboring microRNAs), a tendency of novel microRNAs in the X chromosome to be expressed in testis, and to a lack of a demasculinization process. To test this hypothesis, I analyzed the expression profile of microRNAs in males, females, and gonads in D. pseudoobscura, in which an autosome translocated into the X chromosome effectively becoming part of a sex chromosome (neo-X). I found that the pattern of sex-biased expression is generally conserved between D. melanogaster and D. pseudoobscura. Also, orthologous microRNAs in both species conserve their chromosomal location, indicating that there is no evidence of demasculinization or other interchromosomal movement of microRNAs. Drosophila pseudoobscura-specific microRNAs in the neo-X chromosome tend to be male-biased and particularly expressed in testis. In summary, the apparent paradox resulting from male-biased protein-coding genes depleted in the X chromosome and an enrichment in male-biased microRNAs is consistent with different evolutionary dynamics between coding genes and short RNAs.

Comparative transcriptional analysis uncovers molecular processes in early and mature somatic cyst cells of Drosophila testes

The tight interaction between somatic and germline cells is conserved in animal spermatogenesis. The testes of Drosophila melanogaster are the model of choice to identify processes responsible for mature gamete production. However, processes of differentiation and soma-germline interactions occurring in somatic cyst cells are currently understudied. Here we focused on the comparison of transcriptome expression patterns of early and mature somatic cyst cells to find out the developmental changes taking place in them. We employed a FACS-based approach for the isolation of early and mature somatic cyst cells from fly testes, subsequent preparation of RNA-Seq libraries, and analysis of gene differential expression in the sorted cells. We found increased expression of genes involved in cell cycle-related processes in early cyst cells, which is necessary for the proliferation and self-renewal of a crucial population of early cyst cells, cyst stem cells. Genes proposedly required for lamellipodium-like projection organization for proper cyst formation were also detected among the upregulated ones in early cyst cells. Gene Ontology and interactome analyses of upregulated genes in mature cyst cells revealed a striking over-representation of gene categories responsible for metabolic and catabolic cellular processes, as well as genes supporting the energetic state of the cells provided by oxidative phosphorylation that is carried out in mitochondria. Our comparative analyses of differentially expressed genes revealed major peculiarities in early and mature cyst cells and provide novel insight into their regulation, which is important for male fertility.

Transcriptome Profiling Reveals Genes Related to Sex Determination and Differentiation in Sugarcane Borer (Chilo sacchariphagus Bojer)

Simple Summary

Chilo sacchariphagus Bojer is an important sugarcane pest globally. The identification of key genes associated with sex determination and differentiation will provide important basic information for the sterile insect technique control strategy. In this study, the comparative transcriptomic analysis of female and male adults revealed sex-biased gene expression, indicating putative genetic elements of sex determination and differentiation in this species.

Abstract

Chilo sacchariphagus Bojer is an important sugarcane pest globally. Along with genetic modification strategies, the sterile insect technique (SIT) has gained more attention as an environment-friendly method for pest control. The identification of key genes associated with sex determination and differentiation will provide important basic information for this control strategy. As such, the transcriptome sequencing of female and male adults was conducted in order to understand the sex-biased gene expression and molecular basis of sex determination and differentiation in this species. A total of 60,429 unigenes were obtained; among them, 34,847 genes were annotated. Furthermore, 11,121 deferentially expressed genes (DEGs) were identified, of which 8986 were male-biased and 2135 were female-biased genes. The male-biased genes were enriched for carbon metabolism, peptidase activity and transmembrane transport, while the female-biased genes were enriched for the cell cycle, DNA replication, and the MAPK signaling pathway. In addition, 102 genes related to sex-determination and differentiation were identified, including the protein toll, ejaculatory bulb-specific protein, fruitless, transformer-2, sex-lethal, beta-Catenin, sox, gata4, beta-tubulin, cytosol aminopeptidase, seminal fluid, and wnt4. Furthermore, transcription factors such as myb, bhlh and homeobox were also found to be potentially related to sex determination and differentiation in this species. Our data provide new insights into the genetic elements associated with sex determination and differentiation in Chilo sacchariphagus, and identified potential candidate genes to develop pest-control strategies.

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.

Msl3 promotes germline stem cell differentiation in female Drosophila

Gamete formation from germline stem cells (GSCs) is essential for sexual reproduction. However, the regulation of GSC differentiation is incompletely understood. Set2, which deposits H3K36me3 modifications, is required for GSC differentiation during Drosophila oogenesis. We discovered that the H3K36me3 reader Male-specific lethal 3 (Msl3) and histone acetyltransferase complex Ada2a-containing (ATAC) cooperate with Set2 to regulate GSC differentiation in female Drosophila. Msl3, acting independently of the rest of the male-specific lethal complex, promotes transcription of genes, including a germline-enriched ribosomal protein S19 paralog RpS19b. RpS19b upregulation is required for translation of RNA-binding Fox protein 1 (Rbfox1), a known meiotic cell cycle entry factor. Thus, Msl3 regulates GSC differentiation by modulating translation of a key factor that promotes transition to an oocyte fate.

Male-biased protein expression in primordial germ cells, identified through a comparative study of UAS vectors in Drosophila

In Drosophila, three types of UAS vectors (UASt, UASp, and UASz) are currently available for use with the Gal4-UAS system. They have been used successfully in somatic cells and germline cells from ovaries. However, it remains unclear whether they are functional in the germline cells of embryos, larvae, and adult testes. In this study, we found that all three types of UAS vectors were functional in the germline cells of embryos and larvae and that the UASt and UASz vectors were active in the germline of the distal tip region in adult testes. Moreover, we observed that protein expression from the UAS vectors was male-biased in germline cells of late embryos, whereas their respective mRNA expression levels were not. Furthermore, O-propargyl-puromycin (OPP) staining revealed that protein synthesis was male-biased in these germline cells. In addition, GO terms related to translation and ribosomal maturation were significantly enriched in the male germline. These observations show that translational activity is higher in male than in female germline cells. Therefore, we propose that male-biased protein synthesis may be responsible for the sex differences observed in the early germline.

A putative de novo evolved gene required for spermatid chromatin condensation in Drosophila melanogaster

Comparative genomics has enabled the identification of genes that potentially evolved de novo from non-coding sequences. Many such genes are expressed in male reproductive tissues, but their functions remain poorly understood. To address this, we conducted a functional genetic screen of over 40 putative de novo genes with testis-enriched expression in Drosophila melanogaster and identified one gene, atlas, required for male fertility. Detailed genetic and cytological analyses showed that atlas is required for proper chromatin condensation during the final stages of spermatogenesis. Atlas protein is expressed in spermatid nuclei and facilitates the transition from histone- to protamine-based chromatin packaging. Complementary evolutionary analyses revealed the complex evolutionary history of atlas. The protein-coding portion of the gene likely arose at the base of the Drosophila genus on the X chromosome but was unlikely to be essential, as it was then lost in several independent lineages. Within the last ~15 million years, however, the gene moved to an autosome, where it fused with a conserved non-coding RNA and evolved a non-redundant role in male fertility. Altogether, this study provides insight into the integration of novel genes into biological processes, the links between genomic innovation and functional evolution, and the genetic control of a fundamental developmental process, gametogenesis.

A cross-species approach for the identification of Drosophila male sterility genes

Male reproduction encompasses many essential cellular processes and interactions. As a focal point for these events, sperm offer opportunities for advancing our understanding of sexual reproduction at multiple levels during development. Using male sterility genes identified in human, mouse, and fruit fly databases as a starting point, 103 Drosophila melanogaster genes were screened for their association with male sterility by tissue-specific RNAi knockdown and CRISPR/Cas9-mediated mutagenesis. This list included 56 genes associated with male infertility in the human databases, but not found in the Drosophila database, resulting in the discovery of 63 new genes associated with male fertility in Drosophila. The phenotypes identified were categorized into six distinct classes affecting sperm development. Interestingly, the second largest class (Class VI) caused sterility despite apparently normal testis and sperm morphology suggesting that these proteins may have functions in the mature sperm following spermatogenesis. We focused on one such gene, Rack 1, and found that it plays an important role in two developmental periods, in early germline cells or germline stem cells and in spermatogenic cells or sperm. Taken together, many genes are yet to be identified and their role in male reproduction, especially after ejaculation, remains to be elucidated in Drosophila, where a wealth of data from human and other model organisms would be useful.

Sexual Dimorphism in Growth Rate and Gene Expression Throughout Immature Development in Wild Type Chrysomya rufifacies (Diptera: Calliphoridae) Macquart

Reliability of forensic entomology analyses to produce relevant information to a given case requires an understanding of the underlying arthropod population(s) of interest and the factors contributing to variability. Common traits for analyses are affected by a variety of genetic and environmental factors. One trait of interest in forensic investigations has been species-specific temperature-dependent growth rates. Recent work indicates sexual dimorphism may be important in the analysis of such traits and related genetic markers of age. However, studying sexual dimorphic patterns of gene expression throughout immature development in wild-type insects can be difficult due to a lack of genetic tools, and the limits of most sex-determination mechanisms. Chrysomya rufifacies, however, is a particularly tractable system to address these issues as it has a monogenic sex determination system, meaning females have only a single-sex of offspring throughout their life. Using modified breeding procedures (to ensure single-female egg clutches) and transcriptomics, we investigated sexual dimorphism in development rate and gene expression. Females develop slower than males (9 h difference from egg to eclosion respectively) even at 30°C, with an average egg-to-eclosion time of 225 h for males and 234 h for females. Given that many key genes rely on sex-specific splicing for the development and maintenance of sexually dimorphic traits, we used a transcriptomic approach to identify different expression of gene splice variants. We find that 98.4% of assembled nodes exhibited sex-specific, stage-specific, to sex-by-stage specific patterns of expression. However, the greatest signal in the expression data is differentiation by developmental stage, indicating that sexual dimorphism in gene expression during development may not be investigatively important and that markers of age may be relatively independent of sex. Subtle differences in these gene expression patterns can be detected as early as 4 h post-oviposition, and 12 of these nodes demonstrate homology with key Drosophila sex determination genes, providing clues regarding the distinct sex determination mechanism of C. rufifacies. Finally, we validated the transcriptome analyses through qPCR and have identified five genes that are developmentally informative within and between sexes.

Evolutionary dynamics of sex-biased genes expressed in cricket brains and gonads

Sex-biased gene expression, particularly sex-biased expression in the gonad, has been linked to rates of protein sequence evolution (nonsynonymous to synonymous substitutions, dN/dS) in animals. However, in insects, sex-biased expression studies remain centred on a few holometabolous species. Moreover, other major tissue types such as the brain remain underexplored. Here, we studied sex-biased gene expression and protein evolution in a hemimetabolous insect, the cricket Gryllus bimaculatus. We generated novel male and female RNA-seq data for two sexual tissue types, the gonad and somatic reproductive system, and for two core components of the nervous system, the brain and ventral nerve cord. From a genome-wide analysis, we report several core findings. Firstly, testis-biased genes had accelerated evolution, as compared to ovary-biased and unbiased genes, which was associated with positive selection events. Secondly, although sex-biased brain genes were much less common than for the gonad, they exhibited a striking tendency for rapid protein sequence evolution, an effect that was stronger for the female than male brain. Further, some sex-biased brain genes were linked to sexual functions and mating behaviours, which we suggest may have accelerated their evolution via sexual selection. Thirdly, a tendency for narrow cross-tissue expression breadth, suggesting low pleiotropy, was observed for sex-biased brain genes, suggesting relaxed purifying selection, which we speculate may allow enhanced freedom to evolve adaptive protein functional changes. The findings of rapid evolution of testis-biased genes and male and female-biased brain genes are discussed with respect to pleiotropy, positive selection and the mating biology of this cricket.

A comparative genomic approach using mouse and fruit fly data to discover genes involved in testis function in hymenopterans with a focus on Nasonia vitripennis

Background

Spermatogenesis appears to be a relatively well-conserved process even among distantly related animal taxa such as invertebrates and vertebrates. Although Hymenopterans share many characteristics with other organisms, their complex haplodiploid reproduction system is still relatively unknown. However, they serve as a complementary insect model to Drosophila for studying functional male fertility. In this study, we used a comparative method combining taxonomic, phenotypic data and gene expression to identify candidate genes that could play a significant role in spermatogenesis in hymenopterans.

Results

Of the 546 mouse genes predominantly or exclusively expressed in the mouse testes, 36% had at least one ortholog in the fruit fly. Of these genes, 68% had at least one ortholog in one of the six hymenopteran species we examined. Based on their gene expression profiles in fruit fly testes, 71 of these genes were hypothesized to play a marked role in testis function. Forty-three of these 71 genes had an ortholog in at least one of the six hymenopteran species examined, and their enriched GO terms were related to the G2/M transition or to cilium organization, assembly, or movement. Second, of the 379 genes putatively involved in male fertility in Drosophila, 224 had at least one ortholog in each of the six Hymenoptera species. Finally, we showed that 199 of these genes were expressed in early pupal testis in Nasonia vitripennis; 86 exhibited a high level of expression, and 54 displayed modulated expression during meiosis.

Conclusions

In this study combining phylogenetic and experimental approaches, we highlighted genes that may have a major role in gametogenesis in hymenopterans; an essential prerequisite for further research on functional importance of these genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01825-6.

Impact of male trait exaggeration on sex-biased gene expression and genome architecture in a water strider

BACKGROUND

Exaggerated secondary sexual traits are widespread in nature and often evolve under strong directional sexual selection. Although heavily studied from both theoretical and empirical viewpoints, we have little understanding of how sexual selection influences sex-biased gene regulation during the development of exaggerated secondary sexual phenotypes, and how these changes are reflected in genomic architecture. This is primarily due to the limited availability of representative genomes and associated tissue and sex transcriptomes to study the development of these traits. Here we present the genome and developmental transcriptomes, focused on the legs, of the water strider Microvelia longipes, a species where males exhibit strikingly long third legs compared to females, which they use as weapons.

RESULTS

We generated a high-quality genome assembly with 90% of the sequence captured in 13 scaffolds. The most exaggerated legs in males were particularly enriched in both sex-biased and leg-biased genes, indicating a specific signature of gene expression in association with trait exaggeration. We also found that male-biased genes showed patterns of fast evolution compared to non-biased and female-biased genes, indicative of directional or relaxed purifying selection. By contrast to male-biased genes, female-biased genes that are expressed in the third legs, but not the other legs, are over-represented in the X chromosome compared to the autosomes. An enrichment analysis for sex-biased genes along the chromosomes revealed also that they arrange in large genomic regions or in small clusters of two to four consecutive genes. The number and expression of these enriched regions were often associated with the exaggerated legs of males, suggesting a pattern of common regulation through genomic proximity in association with trait exaggeration.

CONCLUSION

Our findings indicate how directional sexual selection may drive sex-biased gene expression and genome architecture along the path to trait exaggeration and sexual dimorphism.

A Superfamily-wide Activity Atlas of Serine Hydrolases in Drosophila melanogaster

The serine hydrolase (SH) superfamily is, perhaps, one of the largest functional enzyme classes in all forms of life and consists of proteases, peptidases, lipases, and carboxylesterases as representative members. Consistent with the name of this superfamily, all members, without any exception to date, use a nucleophilic serine residue in the enzyme active site to perform hydrolytic-type reactions via a two-step ping-pong mechanism involving a covalent enzyme intermediate. Given the highly conserved catalytic mechanism, this superfamily has served as a classical prototype in the development of several platforms of chemical proteomics techniques, activity-based protein profiling (ABPP), to globally interrogate the functions of its different members in various native, yet complex, biological settings. While ABPP-based proteome-wide activity atlases for SH activities are available in numerous organisms, including humans, to the best of our knowledge, such an analysis for this superfamily is lacking in any insect model. To address this, we initially report a bioinformatics analysis toward the identification and categorization of nonredundant SHs in Drosophila melanogaster. Following up on this in silico analysis, leveraging discovery chemoproteomics, we identify and globally map the full complement of SH activities during various developmental stages and in different adult tissues of Drosophila. Finally, as a proof of concept of the utility of this activity atlas, we highlight sexual dimorphism in SH activities across different tissues in adult D. melanogaster, and we propose new research directions, resources, and tools that this study can provide to the fly community.

Combined transcriptomic, proteomic and genomic analysis identifies reproductive-related proteins and potential modulators of female behaviors in Spodoptera litura

The common cutworm, Spodoptera litura, is a polyandrous moth with high reproductive ability. Sexual reproduction is a unique strategy for survival and reproduction of population in this species. However, to date available information about its reproductive genes is rare. Here, we combined transcriptomics, genomics and proteomics approaches to characterize reproductive-related proteins in S. litura. Illumina sequencing in parallel with the reference genome led to the yields of 12,161 reproductive genes, representing 47.83% of genes annotated in the genome. Further, 524 genes of 19 specific gene families annotated in the genome were detected in reproductive tissues of both sexes, some of which exhibited sex-biased and/or tissue-enriched expression. Of these, manual efforts together with the transcriptome analyses re-annotated 54 odorant binding proteins (OBPs) and 23 chemosensory proteins (CSPs) with an increase of 18 OBPs and one CSP compared to those previously annotated in the genome. Interestingly, at least 35 OBPs and 22 CSPs were transcribed in at least one reproductive tissue, suggestive of their involvement in reproduction. Further proteomic analysis revealed 2381 common proteins between virgin and mated female reproductive systems, 79 of which were differentially expressed. More importantly, 74 proteins exclusive to mated females were identified as transferred relatives, coupled with their specific or high expression in male reproductive systems. Of the transferred proteins, several conserved protein classes across insects were observed including OBPs, serpins, trypsins and juvenile hormone-binding proteins. Our current study has extensively surveyed reproductive genes in S. litura with an emphasis on the roles of OBPs and CSPs in reproduction, and identifies potentially transferred proteins serving as modulators of female post-mating behaviors.

Origin of Sex-Biased Mental Disorders: An Evolutionary Perspective

Sexual dimorphism or sex bias in diseases and mental disorders have two biological causes: sexual selection and sex hormones. We review the role of sexual selection theory and bring together decades of molecular studies on the variation and evolution of sex-biased genes and provide a theoretical basis for the causes of sex bias in disease and health. We present a Sexual Selection-Sex Hormone theory and show that male-driven evolution, including sexual selection, leads to: (1) increased male vulnerability due to negative pleiotropic effects associated with male-driven sexual selection and evolution; (2) increased rates of male-driven mutations and epimutations in response to early fitness gains and at the cost of late fitness; and (3) enhanced female immunity due to antagonistic responses to mutations that are beneficial to males but harmful to females, reducing female vulnerability to diseases and increasing the thresholds for disorders such as autism. Female-driven evolution, such as reproduction-related fluctuation in female sex hormones in association with stress and social condition, has been shown to be associated with increased risk of certain mental disorders such as major depression disorder in women. Bodies have history, cells have memories. An evolutionary framework, such as the Sexual Selection–Sex Hormone theory, provides a historical perspective for understanding how the differences in the sex-biased diseases and mental disorders have evolved over time. It has the potential to direct the development of novel preventive and treatment strategies.

Multi-level analysis of reproduction in an Antarctic midge identifies female and male accessory gland products that are altered by larval stress and impact progeny viability

The Antarctic midge, Belgica antarctica, is a wingless, non-biting midge endemic to Antarctica. Larval development requires at least 2 years, but adults live only 2 weeks. The nonfeeding adults mate in swarms and females die shortly after oviposition. Eggs are suspended in a gel of unknown composition that is expressed from the female accessory gland. This project characterizes molecular mechanisms underlying reproduction in this midge by examining differential gene expression in whole males, females, and larvae, as well as in male and female accessory glands. Functional studies were used to assess the role of the gel encasing the eggs, as well as the impact of stress on reproductive biology. RNA-seq analyses revealed sex- and development-specific gene sets along with those associated with the accessory glands. Proteomic analyses were used to define the composition of the egg-containing gel, which is generated during multiple developmental stages and derived from both the accessory gland and other female organs. Functional studies indicate the gel provides a larval food source as well as a buffer for thermal and dehydration stress. All of these function are critical to juvenile survival. Larval dehydration stress directly reduces production of storage proteins and key accessory gland components, a feature that impacts adult reproductive success. Modeling reveals that bouts of dehydration may have a significant impact on population growth. This work lays a foundation for further examination of reproduction in midges and provides new information related to general reproduction in dipterans. A key aspect of this work is that reproduction and stress dynamics, currently understudied in polar organisms, are likely to prove critical in determining how climate change will alter their survivability.

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.

Sexual dimorphism and sex-biased gene expression in an egg parasitoid species, Anastatus disparis

Background

Differences in the expression of genes present in both sexes are assumed to contribute to sex differences including behavioural, physiological and morphological dimorphisms. For enriching our knowledge of gender differences in an important egg parasitoid wasp, Anastatus disparis (Hymenoptera: Eupelmidae), sex-biased differences in gene expression were investigated using Illumina-based transcriptomic analysis.

Results

A total of 15,812 resulting unigenes were annotated, and a large set of genes accounting for 50.09% of the total showed sex-biased expression and included 630 sex-specific genes. Gene Ontology (GO) enrichment analyses showed that the functional categories associated with sex-biased genes were mainly related to reproduction. In addition, the transcriptome data provided evidence that sex pheromones in A. disparis are produced by the female, and activity of Δ12-desaturases appear to have been replaced by Δ9-desaturases playing roles in sex pheromone production. The large set of sex-biased genes identified in this study provide a molecular background for sexually dimorphic traits such as flyability, longevity, and aggression in this species and suggests candidate venom proteins expressed only in females that could be used for biological control.

Conclusions

This study provides comprehensive insight into sexually dimorphic traits of a parasitoid wasp and can inform future research into the molecular mechanisms underlying such traits and the application of parasitoids to the biological control of pest species.

A role for triglyceride lipase brummer in the regulation of sex differences in Drosophila fat storage and breakdown

Triglycerides are the major form of stored fat in all animals. One important determinant of whole-body fat storage is whether an animal is male or female. Here, we use Drosophila, an established model for studies on triglyceride metabolism, to gain insight into the genes and physiological mechanisms that contribute to sex differences in fat storage. Our analysis of triglyceride storage and breakdown in both sexes identified a role for triglyceride lipase brummer (bmm) in the regulation of sex differences in triglyceride homeostasis. Normally, male flies have higher levels of bmm mRNA both under normal culture conditions and in response to starvation, a lipolytic stimulus. We find that loss of bmm largely eliminates the sex difference in triglyceride storage and abolishes the sex difference in triglyceride breakdown via strongly male-biased effects. Although we show that bmm function in the fat body affects whole-body triglyceride levels in both sexes, in males, we identify an additional role for bmm function in the somatic cells of the gonad and in neurons in the regulation of whole-body triglyceride homeostasis. Furthermore, we demonstrate that lipid droplets are normally present in both the somatic cells of the male gonad and in neurons, revealing a previously unrecognized role for bmm function, and possibly lipid droplets, in these cell types in the regulation of whole-body triglyceride homeostasis. Taken together, our data reveal a role for bmm function in the somatic cells of the gonad and in neurons in the regulation of male–female differences in fat storage and breakdown and identify bmm as a link between the regulation of triglyceride homeostasis and biological sex.

An investigation of the genetic and physiological mechanisms underlying sex differences in fat storage and breakdown in the fruit fly Drosophila identifies previously unrecognized sex- and cell type-specific roles for the conserved triglyceride lipase brummer.

A genomic analysis and transcriptomic atlas of gene expression in Psoroptes ovis reveals feeding- and stage-specific patterns of allergen expression

Background

Psoroptic mange, caused by infestation with the ectoparasitic mite, Psoroptes ovis, is highly contagious, resulting in intense pruritus and represents a major welfare and economic concern for the livestock industry Worldwide. Control relies on injectable endectocides and organophosphate dips, but concerns over residues, environmental contamination, and the development of resistance threaten the sustainability of this approach, highlighting interest in alternative control methods. However, development of vaccines and identification of chemotherapeutic targets is hampered by the lack of P. ovis transcriptomic and genomic resources.

Results

Building on the recent publication of the P. ovis draft genome, here we present a genomic analysis and transcriptomic atlas of gene expression in P. ovis revealing feeding- and stage-specific patterns of gene expression, including novel multigene families and allergens. Network-based clustering revealed 14 gene clusters demonstrating either single- or multi-stage specific gene expression patterns, with 3075 female-specific, 890 male-specific and 112, 217 and 526 transcripts showing larval, protonymph and tritonymph specific-expression, respectively. Detailed analysis of P. ovis allergens revealed stage-specific patterns of allergen gene expression, many of which were also enriched in “fed” mites and tritonymphs, highlighting an important feeding-related allergenicity in this developmental stage. Pair-wise analysis of differential expression between life-cycle stages identified patterns of sex-biased gene expression and also identified novel P. ovis multigene families including known allergens and novel genes with high levels of stage-specific expression.

Conclusions

The genomic and transcriptomic atlas described here represents a unique resource for the acarid-research community, whilst the OrcAE platform makes this freely available, facilitating further community-led curation of the draft P. ovis genome.

Sex identification from distinctive gene expression patterns in Antarctic krill (Euphausia superba)

Antarctic krill (Euphausia superba) is a highly abundant keystone species of the Southern Ocean ecosystem, directly connecting primary producers to high-trophic level predators. Sex ratios of krill vary remarkably between swarms and this phenomenon is poorly understood, as identification of krill sex relies on external morphological differences that appear late during development. Sex determination mechanisms in krill are unknown, but could include genetic, environmental or parasitic mechanisms. Similarly, virtually nothing is known about molecular sex differentiation. The krill genome has to date not been sequenced, and due to its enormous size and large amount of repetitive elements, it is currently not feasible to develop sex-specific DNA markers. To produce a reliable molecular marker for sex in krill and to investigate molecular sex differentiation we therefore focused on identifying sex-specific transcriptomic differences. Through transcriptomic analysis, we found large gene expression differences between testes and ovaries and identified three genes exclusively expressed in female whole krill from early juvenile stages onwards. The sex-specific expression of these three genes persisted through sexual regression, although our regressed samples originated from a krill aquarium and may differ from wild-regressed krill. Two slightly male-biased genes did not display sufficient expression differences to clearly differentiate sexes. Based on the expression of the three female-specific genes we developed a molecular test that for the first time allows the unambiguous sex determination of krill samples lacking external sex-specific features from juvenile stages onwards, including the sexually regressed krill we examined.

Patterns and Constraints in the Evolution of Sperm Individualization Genes in Insects, with an Emphasis on Beetles

Gene expression profiles can change dramatically between sexes and sex bias may contribute specific macroevolutionary dynamics for sex-biased genes. However, these dynamics are poorly understood at large evolutionary scales due to the paucity of studies that have assessed orthology and functional homology for sex-biased genes and the pleiotropic effects possibly constraining their evolutionary potential. Here, we explore the correlation of sex-biased expression with macroevolutionary processes that are associated with sex-biased genes, including duplications and accelerated evolutionary rates. Specifically, we examined these traits in a group of 44 genes that orchestrate sperm individualization during spermatogenesis, with both unbiased and sex-biased expression. We studied these genes in the broad evolutionary framework of the Insecta, with a particular focus on beetles (order Coleoptera). We studied data mined from 119 insect genomes, including 6 beetle models, and from 19 additional beetle transcriptomes. For the subset of physically and/or genetically interacting proteins, we also analyzed how their network structure may condition the mode of gene evolution. The collection of genes was highly heterogeneous in duplication status, evolutionary rates, and rate stability, but there was statistical evidence for sex bias correlated with faster evolutionary rates, consistent with theoretical predictions. Faster rates were also correlated with clocklike (insect amino acids) and non-clocklike (beetle nucleotides) substitution patterns in these genes. Statistical associations (higher rates for central nodes) or lack thereof (centrality of duplicated genes) were in contrast to some current evolutionary hypotheses, highlighting the need for more research on these topics.

Testis single-cell RNA-seq reveals the dynamics of de novo gene transcription and germline mutational bias in Drosophila

The testis is a peculiar tissue in many respects. It shows patterns of rapid gene evolution and provides a hotspot for the origination of genetic novelties such as de novo genes, duplications and mutations. To investigate the expression patterns of genetic novelties across cell types, we performed single-cell RNA-sequencing of adult Drosophila testis. We found that new genes were expressed in various cell types, the patterns of which may be influenced by their mode of origination. In particular, lineage-specific de novo genes are commonly expressed in early spermatocytes, while young duplicated genes are often bimodally expressed. Analysis of germline substitutions suggests that spermatogenesis is a highly reparative process, with the mutational load of germ cells decreasing as spermatogenesis progresses. By elucidating the distribution of genetic novelties across spermatogenesis, this study provides a deeper understanding of how the testis maintains its core reproductive function while being a hotbed of evolutionary innovation.

Asymmetric histone inheritance via strand-specific incorporation and biased replication fork movement

Many stem cells undergo asymmetric division to produce a self-renewing stem cell and a differentiating daughter cell. Here we show that, similarly to H3, histone H4 is inherited asymmetrically in Drosophila melanogaster male germline stem cells undergoing asymmetric division. In contrast, both H2A and H2B are inherited symmetrically. By combining super-resolution microscopy and chromatin fiber analyses with proximity ligation assays on intact nuclei, we find that old H3 is preferentially incorporated by the leading strand, whereas newly synthesized H3 is enriched on the lagging strand. Using a sequential nucleoside analog incorporation assay, we detect a high incidence of unidirectional replication fork movement in testes-derived chromatin and DNA fibers. Biased fork movement coupled with a strand preference in histone incorporation would explain how asymmetric old and new H3 and H4 are established during replication. These results suggest a role for DNA replication in patterning epigenetic information in asymmetrically dividing cells in multicellular organisms.

Differentially and Co-expressed Genes in Embryo, Germ-Line and Somatic Tissues of Tribolium castaneum

Transcriptomic studies of Tribolium castaneum have led to significant advances in our understanding of co-regulation and differential expression of genes in development. However, previously used microarray approaches have covered only a subset of known genes. The aim of this study was to investigate gene expression patterns of beetle embryo, germ-line and somatic tissues. We identified 12,302 expressed genes and determined differentially expressed up and down-regulated genes among all samples. For example, 1624 and 3639 genes were differentially increased in expression greater than or equal to twofold change (FDR < 0.01) in testis vs. ovary (virgin female) and ovary vs. embryo (0-5 hr), respectively. Of these, many developmental, somatic and germ-line differentially expressed genes were identified. Furthermore, many maternally deposited transcripts were identified, whose expression either decreased rapidly or persisted during embryogenesis. Genes with the largest change in expression were predominantly decreased during early embryogenesis as compared to ovary or were increased in testis compared to embryo. We also identify zygotic genes induced after fertilization. The genome wide variation in transcript regulation in maternal and zygotic genes could provide additional information on how the anterior posterior axis formation is established in Tribolium embryos as compared to Drosophila Together, our data will facilitate studies of comparative developmental biology as well as help identify candidate genes for identifying cis-elements to drive transgenic constructs.

Sexual Transcription Differences in Brachymeria lasus (Hymenoptera: Chalcididae), a Pupal Parasitoid Species of Lymantria dispar (Lepidoptera: Lymantriidae)

Sex differences in gene expression have been extensively documented, but little is known about these differences in parasitoid species that are widely applied to control pests. Brachymeria lasus is a solitary parasitoid species and has been evaluated as a potential candidate for release to control Lymantria dispar. In this study, gender differences in B. lasus were investigated using Illumina-based transcriptomic analysis. The resulting 37,453 unigene annotations provided a large amount of useful data for molecular studies of B. lasus. A total of 1416 differentially expressed genes were identified between females and males, and the majority of the sex-biased genes were female biased. Gene Ontology (GO) and Pathway enrichment analyses showed that (1) the functional categories DNA replication, fatty acid biosynthesis, and metabolism were enhanced in females and that (2) the only pathway enriched in males was phototransduction, while the GO subcategories enriched in males were those involved in membrane and ion transport. In addition, thirteen genes involving transient receptor potential (TRP) channels were annotated in B. lasus. We further explored and discussed the functions of TRPs in sensory signaling of light and temperature. In general, this study provides new molecular insights into the biological and sexually dimorphic traits of parasitoids, which may improve the application of these insects to the biological control of pests.

Selection shapes turnover and magnitude of sex-biased expression in Drosophila gonads

BACKGROUND

Sex-biased gene expression is thought to drive the phenotypic differences in males and females in metazoans. Drosophila has served as a primary model for studying male-female differences in gene expression, and its effects on protein sequence divergence. However, the forces shaping evolution of sex-biased expression remain largely unresolved, including the roles of selection and pleiotropy. Research on sex organs in Drosophila, employing original approaches and multiple-species contrasts, provides a means to gain insights into factors shaping the turnover and magnitude (fold-bias) of sex-biased expression.

RESULTS

Here, using recent RNA-seq data, we studied sex-biased gonadal expression in 10,740 protein coding sequences in four species of Drosophila, D. melanogaster, D. simulans, D. yakuba and D. ananassae (5 to 44 My divergence). Using an approach wherein we identified genes with lineage-specific transitions (LSTs) in sex-biased status (amongst testis-biased, ovary-biased and unbiased; thus, six transition types) standardized to the number of genes with the ancestral state (S-LSTs), and those with clade-wide expression bias status, we reveal several key findings. First, the six categorical types of S-LSTs in sex-bias showed disparate rates of turnover, consistent with differential selection pressures. Second, the turnover in sex-biased status was largely unrelated to cross-tissue expression breadth, suggesting pleiotropy does not restrict evolution of sex-biased expression. Third, the fold-sex-biased expression, for both testis-biased and ovary-biased genes, evolved directionally over time toward higher values, a crucial finding that could be interpreted as a selective advantage of greater sex-bias, and sexual antagonism. Fourth, in terms of protein divergence, genes with LSTs to testis-biased expression exhibited weak signals of elevated rates of evolution (than ovary-biased) in as little as 5 My, which strengthened over time. Moreover, genes with clade-wide testis-specific expression (44 My), a status not observed for any ovary-biased genes, exhibited striking acceleration of protein divergence, which was linked to low pleiotropy.

CONCLUSIONS

By studying LSTs and clade-wide sex-biased gonadal expression in a multi-species clade of Drosophila, we describe evidence that interspecies turnover and magnitude of sex-biased expression have been influenced by selection. Further, whilst pleiotropy was not connected to turnover in sex-biased gonadal expression, it likely explains protein sequence divergence.

Sexually dimorphic gene expression and transcriptome evolution provide mixed evidence for a fast-Z effect in Heliconius

Sex chromosomes have different evolutionary properties compared to autosomes due to their hemizygous nature. In particular, recessive mutations are more readily exposed to selection, which can lead to faster rates of molecular evolution. Here, we report patterns of gene expression and molecular evolution for a group of butterflies. First, we improve the completeness of the Heliconius melpomene reference annotation, a neotropical butterfly with a ZW sex determination system. Then, we analyse RNA from male and female whole abdomens and sequence female ovary and gut tissue to identify sex‐ and tissue‐specific gene expression profiles in H. melpomene. Using these expression profiles, we compare (a) sequence divergence and polymorphism; (b) the strength of positive and negative selection; and (c) rates of adaptive evolution, for Z and autosomal genes between two species of Heliconius butterflies, H. melpomene and H. erato. We show that the rate of adaptive substitutions is higher for Z than autosomal genes, but contrary to expectation, it is also higher for male‐biased than female‐biased genes. Additionally, we find no significant increase in the rate of adaptive evolution or purifying selection on genes expressed in ovary tissue, a heterogametic‐specific tissue. Our results contribute to a growing body of literature from other ZW systems that also provide mixed evidence for a fast‐Z effect where hemizygosity influences the rate of adaptive substitutions.

Re-annotation of eight Drosophila genomes

The sequenced genomes of the Drosophila phylogeny are a central resource for comparative work supporting the understanding of the Drosophila melanogaster non-mammalian model system. These have also facilitated evolutionary studies on the selected and random differences that distinguish the thousands of extant species of Drosophila. However, full utility has been hampered by uneven genome annotation. We have generated a large expression profile dataset for nine species of Drosophila and trained a transcriptome assembly approach on D. melanogaster that best matched the extensively curated annotation. We then applied this to the other species to add more than 10000 transcript models per species. We also developed new orthologs to facilitate cross-species comparisons. We validated the new annotation of the distantly related Drosophila grimshawi with an extensive collection of newly sequenced cDNAs. This re-annotation will facilitate understanding both the core commonalities and the species differences in this important group of model organisms, and suggests a strategy for annotating the many forthcoming genomes covering the tree of life.

Variation in Mitochondria-Derived Transcript Levels Associated With DDT Resistance in the 91-R Strain of Drosophila melanogaster (Diptera: Drosophilidae)

The organochloride insecticide dichlorodiphenyltrichloroethane (DDT) and its metabolites can increase cellular levels of reactive oxygen species (ROS), cause mitochondrial dysfunction, and induce apoptosis. The highly DDT-resistant Drosophila melanogaster Meigen 1830 (Drosophila) strain, 91-R, and its susceptible control, 91-C, were used to investigate functional and structural changes among mitochondrial-derived pathways. Resequencing of mitochondrial genomes (mitogenomes) detected no structural differences between 91-R and 91-C, whereas RNA-seq suggested the differential expression of 221 mitochondrial-associated genes. Reverse transcriptase-quantitative PCR validation of 33 candidates confirmed that transcripts for six genes (Cyp12d1-p, Cyp12a4, cyt-c-d, COX5BL, COX7AL, CG17140) were significantly upregulated and two genes (Dif, Rel) were significantly downregulated in 91-R. Among the upregulated genes, four genes are duplicated within the reference genome (cyt-c-d, CG17140, COX5BL, and COX7AL). The predicted functions of the differentially expressed genes, or known functions of closely related genes, suggest that 91-R utilizes existing ROS regulation pathways of the mitochondria to combat increased ROS levels from exposure to DDT. This study represents, to our knowledge, the initial investigation of mitochondrial genome sequence variants and functional adaptations in responses to intense DDT selection and provides insights into potential adaptations of ROS management associated with DDT selection in Drosophila.

Masculinization of gene expression is associated with male quality in Drosophila melanogaster

The signature of sexual selection has been revealed through the study of differences in patterns of genome-wide gene expression, both between the sexes and between alternative reproductive morphs within a single sex. What remains unclear, however, is whether differences in gene expression patterns between individuals of a given sex consistently map to variation in individual quality. Such a pattern, particularly if found in males, would provide unambiguous evidence that the phenotypic response to sexual selection is shaped through sex-specific alterations to the transcriptome. To redress this knowledge gap, we explored whether patterns of sex-biased gene expression are associated with variation in male reproductive quality in Drosophila melanogaster. We measured two male reproductive phenotypes, and their association with sex-biased gene expression, across a selection of inbred lines from the Drosophila Genetic Reference Panel. Genotypes with higher expression of male-biased genes produced males exhibiting shorter latencies to copulation, and higher capacity to inseminate females. Conversely, female-biased genes tended to show negative associations with these male reproductive traits across genotypes. We uncovered similar patterns, by reanalyzing a published dataset from a second D. melanogaster population. Our results reveal the footprint of sexual selection in masculinising the male transcriptome.

The H3K9 methyltransferase SETDB1 maintains female identity in Drosophila germ cells

The preservation of germ cell sexual identity is essential for gametogenesis. Here we show that H3K9me3-mediated gene silencing is integral to female fate maintenance in Drosophila germ cells. Germ cell specific loss of the H3K9me3 pathway members, the H3K9 methyltransferase SETDB1, WDE, and HP1a, leads to ectopic expression of genes, many of which are normally expressed in testis. SETDB1 controls the accumulation of H3K9me3 over a subset of these genes without spreading into neighboring loci. At phf7, a regulator of male germ cell sexual fate, the H3K9me3 peak falls over the silenced testis-specific transcription start site. Furthermore, H3K9me3 recruitment to phf7 and repression of testis-specific transcription is dependent on the female sex determination gene Sxl. Thus, female identity is secured by an H3K9me3 epigenetic pathway in which Sxl is the upstream female-specific regulator, SETDB1 is the required chromatin writer, and phf7 is one of the critical SETDB1 target genes.

Epigenetic regulation is critical for the maintenance of germ cell identity. Here the authors show that H3K9me3-mediated gene silencing is critical for repression of testis-specific transcription in Drosophila female germ cells, indicating H3K9me3 maintains female germ cell sexual identity.

Reprogramming of regulatory network using expression uncovers sex-specific gene regulation in Drosophila

Gene regulatory networks (GRNs) describe regulatory relationships between transcription factors (TFs) and their target genes. Computational methods to infer GRNs typically combine evidence across different conditions to infer context-agnostic networks. We develop a method, Network Reprogramming using EXpression (NetREX), that constructs a context-specific GRN given context-specific expression data and a context-agnostic prior network. NetREX remodels the prior network to obtain the topology that provides the best explanation for expression data. Because NetREX utilizes prior network topology, we also develop PriorBoost, a method that evaluates a prior network in terms of its consistency with the expression data. We validate NetREX and PriorBoost using the "gold standard" E. coli GRN from the DREAM5 network inference challenge and apply them to construct sex-specific Drosophila GRNs. NetREX constructed sex-specific Drosophila GRNs that, on all applied measures, outperform networks obtained from other methods indicating that NetREX is an important milestone toward building more accurate GRNs.

Ocnus is essential for male germ cell development in Drosophila melanogaster

The ocnus (ocn) gene encodes a protein abundant in the testes, implying its role in testis development. When Drosophila melanogaster is infected with the endosymbiont wMel Wolbachia, which affects the spermatogenesis of its hosts, ocn is downregulated in the third-instar larval testes, suggesting a role of ocn in spermatogenesis. In this study, we knocked down ocn in the testes and found that the hatch rates of embryos derived from ocn-knockdown males were significantly decreased, and 84.38% of the testes were much smaller in comparison to controls. Analysis of the smaller testes showed no germ cells but they had an extended hub. Using RNA-sequencing (RNA-Seq), we identified 69 genes with at least a twofold change (q-value < 5%) in their expression after ocn knockdown; of these, eight testes-specific and three reproduction-related genes were verified to be significantly downregulated using quantitative reverse transcription-PCR. Three genes (orientation disruptor, p24-2 and CG13541) were also significantly downregulated in the presence of Wolbachia. Furthermore, 98 genes were not expressed when ocn was knocked down in testes. These results suggest that ocn plays a crucial role in male germ cell development in Drosophila, possibly by regulating the expression of multiple spermatogenesis-related genes. Our data provide important information to help understand the molecular regulatory mechanisms underlying spermatogenesis.

Tissue Specificity and Dynamics of Sex-Biased Gene Expression in a Common Frog Population with Differentiated, Yet Homomorphic, Sex Chromosomes

Sex-biased genes are central to the study of sexual selection, sexual antagonism, and sex chromosome evolution. We describe a comprehensive de novo assembled transcriptome in the common frog Rana temporaria based on five developmental stages and three adult tissues from both sexes, obtained from a population with karyotypically homomorphic but genetically differentiated sex chromosomes. This allows the study of sex-biased gene expression throughout development, and its effect on the rate of gene evolution while accounting for pleiotropic expression, which is known to negatively correlate with the evolutionary rate. Overall, sex-biased genes had little overlap among developmental stages and adult tissues. Late developmental stages and gonad tissues had the highest numbers of stage- or tissue-specific genes. We find that pleiotropic gene expression is a better predictor than sex bias for the evolutionary rate of genes, though it often interacts with sex bias. Although genetically differentiated, the sex chromosomes were not enriched in sex-biased genes, possibly due to a very recent arrest of XY recombination. These results extend our understanding of the developmental dynamics, tissue specificity, and genomic localization of sex-biased genes.

Identification of General Patterns of Sex-Biased Expression in Daphnia, a Genus with Environmental Sex Determination

Daphnia reproduce by cyclic-parthenogenesis, where phases of asexual reproduction are intermitted by sexual production of diapause stages. This life cycle, together with environmental sex determination, allow the comparison of gene expression between genetically identical males and females. We investigated gene expression differences between males and females in four genotypes of Daphnia magna and compared the results with published data on sex-biased gene expression in two other Daphnia species, each representing one of the major phylogenetic clades within the genus. We found that 42% of all annotated genes showed sex-biased expression in D. magna This proportion is similar both to estimates from other Daphnia species as well as from species with genetic sex determination, suggesting that sex-biased expression is not reduced under environmental sex determination. Among 7453 single copy, one-to-one orthologs in the three Daphnia species, 707 consistently showed sex-biased expression and 675 were biased in the same direction in all three species. Hence these genes represent a core-set of genes with consistent sex-differential expression in the genus. A functional analysis identified that several of them are involved in known sex determination pathways. Moreover, 75% were overexpressed in females rather than males, a pattern that appears to be a general feature of sex-biased gene expression in Daphnia.

The impact of genome variation and diet on the metabolic phenotype and microbiome composition of Drosophila melanogaster

The metabolic phenotype of an organism depends on a complex regulatory network, which integrates the plethora of intrinsic and external information and prioritizes the flow of nutrients accordingly. Given the rise of metabolic disorders including obesity, a detailed understanding of this regulatory network is in urgent need. Yet, our level of understanding is far from completeness and complicated by the discovery of additional layers in metabolic regulation, such as the impact of the microbial community present in the gut on the hosts’ energy storage levels. Here, we investigate the interplay between genome variation, diet and the gut microbiome in the shaping of a metabolic phenotype. For this purpose, we reared a set of fully sequenced wild type Drosophila melanogaster flies under basal and nutritionally challenged conditions and performed metabolic and microbiome profiling experiments. Our results introduce the fly as a model system to investigate the impact of genome variation on the metabolic response to diet alterations and reveal candidate single nucleotide polymorphisms associated with different metabolic traits, as well as metabolite-metabolite and metabolite-microbe correlations. Intriguingly, the dietary changes affected the microbiome composition less than anticipated. These results challenge the current view of a rapidly changing microbiome in response to environmental fluctuations.