Beyond Haldane's rule: Sex-biased hybrid dysfunction for all modes of sex determination

Haldane's rule occupies a special place in biology as one of the few 'rules' of speciation, with empirical support from hundreds of species. And yet, its classic purview is restricted taxonomically to the subset of organisms with heteromorphic sex chromosomes. I propose explicit acknowledgement of generalized hypotheses about Haldane's rule that frame sex bias in hybrid dysfunction broadly and irrespective of the sexual system. The consensus view of classic Haldane's rule holds that sex-biased hybrid dysfunction across taxa is a composite phenomenon that requires explanations from multiple causes. Testing of the multiple alternative hypotheses for Haldane's rule is, in many cases, applicable to taxa with homomorphic sex chromosomes, environmental sex determination, haplodiploidy, and hermaphroditism. Integration of a variety of biological phenomena about hybrids across diverse sexual systems, beyond classic Haldane's rule, will help to derive a more general understanding of the contributing forces and mechanisms that lead to predictable sex biases in evolutionary divergence and speciation.

Sexual size dimorphism in mammals is associated with changes in the size of gene families related to brain development

In mammals, sexual size dimorphism often reflects the intensity of sexual selection, yet its connection to genomic evolution remains unexplored. Gene family size evolution can reflect shifts in the relative importance of different molecular functions. Here, we investigate the associate between brain development gene repertoire to sexual size dimorphism using 124 mammalian species. We reveal significant changes in gene family size associations with sexual size dimorphism. High levels of dimorphism correlate with an expansion of gene families enriched in olfactory sensory perception and a contraction of gene families associated with brain development functions, many of which exhibited particularly high expression in the human adult brain. These findings suggest a relationship between intense sexual selection and alterations in gene family size. These insights illustrate the complex interplay between sexual dimorphism, gene family size evolution, and their roles in mammalian brain development and function, offering a valuable understanding of mammalian genome evolution.

Transcriptomic and proteomic elucidation of Z chromosome dosage compensation in Helicoverpa armigera

Transcriptomic data have been used to study sex chromosome dosage compensation (SCDC) in approximately 10 Lepidoptera ZW species, yielding a consensus compensation pattern of Z≈ ZZ < AA . $$ \approx \mathrm{ZZ}<\mathrm{AA}. $$ It remains unclear whether this compensation pattern holds when examining more Lepidoptera ZW species and/or using proteomic data to analyse SCDC. Here we combined transcriptomic and proteomic data as well as transcriptional level of six individual Z genes to reveal the SCDC pattern in Helicoverpa armigera, a polyphagous lepidopteran pest of economic importance. Transcriptomic analysis showed that the Z chromosome expression of H. armigera was balanced between male and female but substantially reduced relative to autosome expression, exhibiting an SCDC pattern of Z≈ ZZ < AA $$ \approx \mathrm{ZZ}<\mathrm{AA} $$ . When using H. amigera midgut proteomic data, the SCDC pattern of this species changed from Z≈ ZZ < AA $$ \approx \mathrm{ZZ}<\mathrm{AA} $$ at transcriptomic level to Z = ZZ = AA at the proteomic level. RT-qPCR analysis of transcript abundance of six Z genes found that compensation for each Z gene could vary from no compensation to overcompensation, depending on the individual genes and tissues tested. These results demonstrate for the first time the existence of a translational compensation mechanism, which is operating in addition to a translational mechanism, such as has been reported in other lepidopteran species. And the transcriptional compensation mechanism functions to accomplish Z chromosome dosage balance between the sexes (M = F on the Z chromosome), whereas the translation compensation mechanism operates to achieve dosage compensation between Z chromosome and autosome (Z = AA).

Blubber gene expression and cortisol concentrations reveal changing physiological stress in a Southern ocean sentinel species

The health of migratory eastern Australian humpback whales (Megaptera novaeangliae) can reflect the condition of their remote polar foraging environments. This study used gene expression (LEP, LEPR, ADIQ, AhR, TNF-α, HSP-70), blubber hormone concentrations (cortisol, testosterone), and photogrammetric body condition to assess this sentinel species during a period of unprecedented changes to anthropogenic activity and natural processes. The results revealed higher cortisol concentrations in 2020 compared to 2021, suggesting a decline in physiological stress between years. Additionally, metabolic transcripts LEPR, and AhR, which is also linked to xenobiotic metabolism, were upregulated during the 2020 southbound migration. These differences suggest that one or more environmental stressors were reduced between 2020 and 2021, with upregulated AhR possibly indicating a Southern Ocean pollutant declined between the years. This research confirms a Southern Ocean-wide decrease in whale stress during the study period and informs efforts to identify key stressors on Antarctic marine ecosystems.

Gonadal transcriptome analysis of sex-biased gene and genome-wide investigation of dmrt gene family in Acanthogobius ommaturus

Acanthogobius ommaturus is one of the largest goby fish, and widely distributed in the Northwestern Pacific as an annual benthic fish. This study aims to report the gonadal transcriptome of A. ommaturus and identify differentially expressed genes (DEGs) between sexes. A total of 5460 (27.94 %) DEGs were detected from genome, with 3301 (16.89 %) biased towards males and 2159 (11.05 %) towards females. Analysis of 76 known vertebrate sex-related genes revealed multiple key genes, including the male-biased genes dmrt1 (Doublesex and Mab-3 related transcription factor 1) and amh (Anti-Mullerian Hormone), and the female-biased genes foxl2 (Forkhead Box L2) and cyp19a1a (Cytochrome P450 Aromatase 19 Subfamily A1). Furthermore, a genome-wide gene family analysis focused on the most significantly differentially expressed male-biased gene, dmrt1, was conducted using the chromosomal-level genome. Six Aodmrt genes were identified and subjected to phylogenetic and protein interaction network analyses. To validate the expression pattern, quantitative real-time PCR (qRT-PCR) was performed and compared with gonadal transcriptome data. The results showed that only dmrt1 exhibited significant male-bias, while the expression levels and sex differences of other dmrt genes in the gonads were inconclusive. Interestingly, the other dmrt genes displayed higher expression levels in other tissues, suggesting currently unknown functions. In conclusion, this study provides valuable genetic information contributing to the understanding of the sex determination mechanism of A. ommaturus and bony fish.

Reference genes to study the sex-biased expression of genes regulating Drosophila metabolism

Sex is an important variable in biology. Notable differences have been observed between male and female Drosophila in regulation of metabolism, in response to nutritional challenges, and in phenotypes relevant for obesity and metabolic disorders. The differences between males and females can be expected to result from differences in gene expression. We observed that expression levels of reference genes commonly used for normalization of qRT-PCR results such as GAPDH, β-actin, and 18SrRNA, show prominent sexual dimorphism. Since this will impact relative expression and conclusions related to that, we performed a systematic analysis of candidate reference genes with the objective of identifying reference genes with stable expression in male and female Drosophila. These reference genes (LamCa, βTub60D and βTub97EF) were then used to assess sex-specific differences in expression of metabolism associated genes. Additionally, we evaluated the utility of these reference genes following a nutritional challenge and showed that LamCa and βtub97EF are stably expressed between sexes and under different nutritional conditions and are thus suitable as reference genes. Our results highlight the importance of evaluating the stability of reference genes when sex-specific differences in gene expression are studied, and identify structural genes as a category worth exploring as reference genes in other species. Finally, we also uncovered hitherto unknown sexually dimorphic expression of a number of metabolism-associated genes, information of interest to others working in the field of metabolic disorders.

Navigating sex and sex roles: deciphering sex-biased gene expression in a species with sex-role reversal (Syngnathus typhle)

Sexual dimorphism, the divergence in morphological traits between males and females of the same species, is often accompanied by sex-biased gene expression. However, the majority of research has focused on species with conventional sex roles, where females have the highest energy burden with both egg production and parental care, neglecting the diversity of reproductive roles found in nature. We investigated sex-biased gene expression in Syngnathus typhle, a sex-role reversed species with male pregnancy, allowing us to separate two female traits: egg production and parental care. Using RNA sequencing, we examined gene expression across organs (brain, head kidney and gonads) at various life stages, encompassing differences in age, sex and reproductive status. While some gene groups were more strongly associated with sex roles, such as stress resistance and immune defence, others were driven by biological sex, such as energy and lipid storage regulation in an organ- and age-specific manner. By investigating how genes regulate and are regulated by changing reproductive roles and resource allocation in a model system with an unconventional life-history strategy, we aim to better understand the importance of sex and sex role in regulating gene expression patterns, broadening the scope of this discussion to encompass a wide range of organisms.

Accurate Bayesian inference of sex chromosome karyotypes and sex-linked scaffolds from low-depth sequencing data

The identification of sex-linked scaffolds and the genetic sex of individuals, i.e. their sex karyotype, is a fundamental step in population genomic studies. If sex-linked scaffolds are known, single individuals may be sexed based on read counts of next-generation sequencing data. If both sex-linked scaffolds as well as sex karyotypes are unknown, as is often the case for non-model organisms, they have to be jointly inferred. For both cases, current methods rely on arbitrary thresholds, which limits their power for low-depth data. In addition, most current methods are limited to euploid sex karyotypes (XX and XY). Here we develop BeXY, a fully Bayesian method to jointly infer the posterior probabilities for each scaffold to be autosomal, X- or Y-linked and for each individual to be any of the sex karyotypes XX, XY, X0, XXX, XXY, XYY and XXYY. If the sex-linked scaffolds are known, it also identifies autosomal trisomies and estimates the sex karyotype posterior probabilities for single individuals. As we show with downsampling experiments, BeXY has higher power than all existing methods. It accurately infers the sex karyotype of ancient human samples with as few as 20,000 reads and accurately infers sex-linked scaffolds from data sets of just a handful of samples or with highly imbalanced sex ratios, also in the case of low-quality reference assemblies. We illustrate the power of BeXY by applying it to both whole-genome shotgun and target enrichment sequencing data of ancient and modern humans, as well as several non-model organisms.

Sex-biased gene regulation varies across human populations as a result of adaptive evolution

OBJECTIVES

Studies of human sexual dimorphism and gender disparities in health focus on ostensibly universal molecular sex differences, such as sex chromosomes and circulating hormone levels, while ignoring the extraordinary diversity in biology, behavior, and culture acquired by different human populations over their unique evolutionary histories.

MATERIALS AND METHODS

Using RNA-Seq data and whole genome sequences from 1000G and HGDP, we investigate variation in sex-biased gene expression across 11 human populations and test whether population-level variation in sex-biased expression may have resulted from adaptive evolution in regions containing sex-specific regulatory variants.

RESULTS

We find that sex-biased gene expression in humans is highly variable, mostly population-specific, and demonstrates between population reversals. Expression quantitative trait locus mapping reveals sex-specific regulatory regions with evidence of recent positive natural selection, suggesting that variation in sex-biased expression may have evolved as an adaptive response to ancestral environments experienced by human populations.

DISCUSSION

These results indicate that sex-biased gene expression is more flexible than previously thought and is not generally shared among human populations. Instead, molecular phenotypes associated with sex depend on complex interactions between population-specific molecular evolution and physiological responses to contemporary socioecologies.

Positive selection and relaxed purifying selection contribute to rapid evolution of male-biased genes in a dioecious flowering plant

Sex-biased genes offer insights into the evolution of sexual dimorphism. Sex-biased genes, especially those with male bias, show elevated evolutionary rates of protein sequences driven by positive selection and relaxed purifying selection in animals. Although rapid sequence evolution of sex-biased genes and evolutionary forces have been investigated in animals and brown algae, less is known about evolutionary forces in dioecious angiosperms. In this study, we separately compared the expression of sex-biased genes between female and male floral buds and between female and male flowers at anthesis in dioecious Trichosanthes pilosa (Cucurbitaceae). In floral buds, sex-biased gene expression was pervasive, and had significantly different roles in sexual dimorphism such as physiology. We observed higher rates of sequence evolution for male-biased genes in floral buds compared to female-biased and unbiased genes. Male-biased genes under positive selection were mainly associated with functions to abiotic stress and immune responses, suggesting that high evolutionary rates are driven by adaptive evolution. Additionally, relaxed purifying selection may contribute to accelerated evolution in male-biased genes generated by gene duplication. Our findings, for the first time in angiosperms, suggest evident rapid evolution of male-biased genes, advance our understanding of the patterns and forces driving the evolution of sexual dimorphism in dioecious plants.

Gender specific SNP markers in Coscinium fenestratum (Gaertn.) Colebr. for resource augmentation

BACKGROUND

Unregulated extraction of highly traded medicinal plant species results in drastic decline of the natural resources and alters viable sex ratio of populations. Conservation and long-term survival of such species, require gender specific restoration programs to ensure reproductive success. However, it is often difficult to differentiate sex of individuals before reaching reproductive maturity. C. fenestratum is one of the medicinally important and overexploited dioecious woody liana, with a reproductive maturity of 15 years. Currently, no information is available to identify sex of C. fenestratum in seedling stage while augmenting the resources. Thus, the current study envisages to utilize transcriptomics approach for gender differentiation which is imperative for undertaking viable resource augmentation programmes.

METHODS AND RESULTS

Gender specific SNPs with probable role in sexual reproduction/sex determination was located using comparative transcriptomics approach (sampling male and female individuals), alongside gene ontology and annotation. Nine sets of primers were synthesized from 7 transcripts (involved in sexual reproduction/other biological process) containing multiple SNP variants. Out of the nine primer pairs, only one SNP locus with no available information of its role in reproduction, showed consistent and accurate results (males-heterozygous and females-homozygous), in the analyzed 40 matured individuals of known sexes. Thus validated the efficiency of this SNP marker in differentiating male and female individuals.

CONCLUSIONS

The study could identify SNPs linked to the loci with apparent role in gender differentiation. This SNP marker can be used for early sexing of seedlings for in-situ conservation and resource augmentation of C. fenestratum in Kerala, India.

Sex-limited experimental evolution drives transcriptomic divergence in a hermaphrodite

The evolution of gonochorism from hermaphroditism is linked with the formation of sex chromosomes, as well as the evolution of sex-biased and sex-specific gene expression to allow both sexes to reach their fitness optimum. There is evidence that sexual selection drives the evolution of male-biased gene expression in particular. However, previous research in this area in animals comes from either theoretical models or comparative studies of already old sex chromosomes. We therefore investigated changes in gene expression under 3 different selection regimes for the simultaneous hermaphrodite Macrostomum lignano subjected to sex-limited experimental evolution (i.e. selection for fitness via eggs, sperm, or a control regime allowing both). After 21 and 22 generations of selection for male-specific or female-specific fitness, we characterized changes in whole-organism gene expression. We found that female-selected lines had changed the most in their gene expression. Although annotation for this species is limited, gene ontology term and Kyoto Encyclopedia of Genes and Genomes pathway analyses suggest that metabolic changes (e.g. biosynthesis of amino acids and carbon metabolism) are an important adaptive component. As predicted, we found that the expression of genes previously identified as testis-biased candidates tended to be downregulated in the female-selected lines. We did not find any significant expression differences for previously identified candidates of other sex-specific organs, but this may simply reflect that few transcripts have been characterized in this way. In conclusion, our experiment suggests that changes in testis-biased gene expression are important in the early evolution of sex chromosomes and gonochorism.

Sex-biased gene content is associated with sex chromosome turnover in Danaini butterflies

Sex chromosomes play an outsized role in adaptation and speciation, and thus deserve particular attention in evolutionary genomics. In particular, fusions between sex chromosomes and autosomes can produce neo-sex chromosomes, which offer important insights into the evolutionary dynamics of sex chromosomes. Here, we investigate the evolutionary origin of the previously reported Danaus neo-sex chromosome within the tribe Danaini. We assembled and annotated genomes of Tirumala septentrionis (subtribe Danaina), Ideopsis similis (Amaurina), Idea leuconoe (Euploeina) and Lycorea halia (Itunina) and identified their Z-linked scaffolds. We found that the Danaus neo-sex chromosome resulting from the fusion between a Z chromosome and an autosome corresponding to the Melitaea cinxia chromosome (McChr) 21 arose in a common ancestor of Danaina, Amaurina and Euploina. We also identified two additional fusions as the W chromosome further fused with the synteny block McChr31 in I. similis and independent fusion occurred between ancestral Z chromosome and McChr12 in L. halia. We further tested a possible role of sexually antagonistic selection in sex chromosome turnover by analysing the genomic distribution of sex-biased genes in I. leuconoe and L. halia. The autosomes corresponding to McChr21 and McChr31 involved in the fusions are significantly enriched in female- and male-biased genes, respectively, which could have hypothetically facilitated fixation of the neo-sex chromosomes. This suggests a role of sexual antagonism in sex chromosome turnover in Lepidoptera. The neo-Z chromosomes of both I. leuconoe and L. halia appear fully compensated in somatic tissues, but the extent of dosage compensation for the ancestral Z varies across tissues and species.

A novel RT-qPCR health assay reveals differential expression of stress and immunoregulatory genes between the seasonal migrations of humpback whales (Megaptera novaeangliae)

Health information is essential for the conservation management of whale species. However, assessing the health of free-ranging whales is challenging as samples are primarily limited to skin and blubber tissue. Reverse transcription quantitative polymerase chain reaction (RT-qPCR) offers a method to measure health from blubber RNA, providing insights into energetic status, stress and immune activity. To identify changes in health, natural differences in baseline gene expression linked to an individual's sex, reproductive status and life-history stage must first be quantified. This study aimed to establish baseline gene expression indices of health in migrating humpback whales (Megaptera novaeangliae). To do this, we developed an assay to quantify seven health-related gene transcripts (Leptin, Leptin Receptor, Adiponectin, Aryl Hydrocarbon Receptor, Tumour Necrosis Factor-α, Interleukin-6, Heat Shock Protein-70) and used Bayesian mixed effect models to assess differential baseline expression based on sex, lactation status and migration stage (northbound to and southbound from the annual breeding grounds). Results showed no significant contribution of sex to differential baseline expression. However, lactating individuals exhibited downregulated AhR and HSP-70 compared to non-lactating conspecifics. Additionally, southbound individuals demonstrated significantly upregulated HSP-70 and downregulated TNF-alpha, suggesting a relationship between these inflammation-linked transcripts and migratory fasting. Our results suggest that baseline differences due to migratory stage and lactation status should be considered in health applications of this assay. Future monitoring efforts can use our baseline measurements to better understand how gene expression is tied to population-level impacts, such as reduced prey availability or migratory stressors.

Bivalves as Emerging Model Systems to Study the Mechanisms and Evolution of Sex Determination: A Genomic Point of View

Bivalves are a diverse group of molluscs that have recently attained a central role in plenty of biological research fields, thanks to their peculiar life history traits. Here, we propose that bivalves should be considered as emerging model systems also in sex-determination (SD) studies, since they would allow to investigate: 1) the transition between environmental and genetic SD, with respect to different reproductive backgrounds and sexual systems (from species with strict gonochorism to species with various forms of hermaphroditism); 2) the genomic evolution of sex chromosomes (SCs), considering that no heteromorphic SCs are currently known and that homomorphic SCs have been identified only in a few species of scallops; 3) the putative role of mitochondria at some level of the SD signaling pathway, in a mechanism that may resemble the cytoplasmatic male sterility of plants; 4) the evolutionary history of SD-related gene (SRG) families with respect to other animal groups. In particular, we think that this last topic may lay the foundations for expanding our understanding of bivalve SD, as our current knowledge is quite fragmented and limited to a few species. As a matter of fact, tracing the phylogenetic history and diversity of SRG families (such as the Dmrt, Sox, and Fox genes) would allow not only to perform more targeted functional experiments and genomic analyses, but also to foster the possibility of establishing a solid comparative framework.

The genome and sex-dependent responses to temperature in the common yellow butterfly, Eurema hecabe

BACKGROUND

Lepidoptera (butterflies and moths) is one of the most geographically widespread insect orders in the world, and its species play important and diverse ecological and applied roles. Climate change is one of the biggest challenges to biodiversity this century, and lepidopterans are vulnerable to climate change. Temperature-dependent gene expression differences are of relevance under the ongoing climate crisis. However, little is known about how climate affects gene expression in lepidopterans and the ecological consequences of this, particularly with respect to genes with biased expression in one of the sexes. The common yellow butterfly, Eurema hecabe (Family Pieridae), is one of the most geographically widespread lepidopterans that can be found in Asia, Africa, and Australia. Nevertheless, what temperature-dependent effects there may be and whether the effects differ between the sexes remain largely unexplored.

RESULTS

Here, we generated high-quality genomic resources for E. hecabe along with transcriptomes from eight developmental stages. Male and female butterflies were subjected to varying temperatures to assess sex-specific gene expression responses through mRNA and microRNA transcriptomics. We find that there are more temperature-dependent sex-biased genes in females than males, including genes that are involved in a range of biologically important functions, highlighting potential ecological impacts of increased temperatures. Further, by considering available butterfly data on sex-biased gene expression in a comparative genomic framework, we find that the pattern of sex-biased gene expression identified in E. hecabe is highly species-specific, rather than conserved across butterfly species, suggesting that sex-biased gene expression responses to climate change are complex in butterflies.

CONCLUSIONS

Our study lays the foundation for further understanding of differential responses to environmental stress in a widespread lepidopteran model and demonstrates the potential complexity of sex-specific responses of lepidopterans to climate change.

Post-insemination sexual selection in males indirectly masculinizes the female transcriptome

Sex-specific regulation of gene expression is the most plausible way for generating sexually differentiated phenotypes from an essentially shared genome. However, since genetic material is shared, sex-specific selection in one sex can have an indirect response in the other sex. From a gene expression perspective, this tethered response can move one sex away from their wildtype expression state and impact potentially many gene regulatory networks. Here, using experimental evolution in the model nematode Caenorhabditis elegans , we explore the coupling of direct sexual selection on males with the transcriptomic response in females over microevolutionary timescales to uncover the extent to which post-insemination reproductive traits share a genetic basis between the sexes. We find that differential gene expression is driven by female ancestral or evolved generation alone and that male generation has no impact on changes in gene expression. Almost all differentially expressed genes were downregulated in evolved females. Moreover, 80% of these gene were located on the X chromosome and have wildtype female-biased expression profiles. Changes in gene expression profiles were likely driven through trans -acting pathways that are shared between the sexes. We found no evidence that the core dosage compensation machinery was impacted by experimental evolution. Together these data suggest masculinization of the female transcriptome driven by direct selection on male sperm competitive ability. Our results indicate that on short evolutionary timescales sexual selection can generate sexual conflict in expression space.

LAY SUMMARY

Sexual selection drives the evolution of some of the most dramatic phenotypic differences between the sexes. Such sexual dimorphism is so common across multicellular organisms that we often overlook how remarkable it is for shared genetic material to create numerous and complex sex differences. At an evolutionary level, sexual dimorphism furthers the opportunity for sex-specific selection to optimize the fitness of a given sex. As a consequence, sex-specific selection, such as sexual selection, can have an indirect evolutionary response in the other sex due to genetic associations created by the sexes sharing the same genome. This correlated evolutionary response can create sexual conflict by shifting a sex away from their fitness optimum. At the functional level, sexual dimorphism is generated is through sex-specific regulation of gene expression. Bridging the evolutionary response to sexual selection with the evolution of sex-specific gene regulation during post-mating interactions has proved challenging. We previously used experimental evolution to increase male fertility by directly selecting for increased sperm competitive ability. In this study, we examined the effect of this direct selection on males on gene expression patterns in females. Differential gene expression was determined by whether a female was ancestral or evolved generation, indicating that gene expression changes were an evolved response due to indirect selection on females. Significantly differentially expressed genes were downregulated in evolved females. These genes tended to be female-biased in wildtype individuals and located on the X chromosome. The downregulation of X-linked genes suggests expression levels in females equal to or lower than that in males. Together these results indicate a less female-like transcriptome after experimental evolution. This supports a sexual conflict scenario by which direct sexual selection on males indirectly masculinizes the female transcriptome over short evolutionary timescales.

Evaluating the role of sexual antagonism in the evolution of sex chromosomes: new data from fish

The recent increase in available molecular and genomic data for diverse taxa helps to shed new light on long-standing theories. Research into sex chromosome evolution has particularly benefited from a growing number of studies of fish, motivated by their highly diverse mechanisms of sex determination. Sexual antagonism is regularly cited as an influential force in sex chromosome emergence; however, this so far proves difficult to demonstrate. In this review, we highlight recent developments in the investigation of sexual antagonism in sex chromosome research in fish. We find strong emphasis placed on study-organism specific genomic features and patterns of recombination, rather than evidence for a comprehensive role of sexual antagonism. In this light, we discuss the alternative models of sex chromosome evolution. We conclude that fish represents a key resource for further research, provided attention is given to species-specific effects while simultaneously integrating comparative studies across taxa for a vital and comprehensive understanding of sex chromosome evolution and investigation of proposed models.

The roles of sexual selection and sexual conflict in shaping patterns of genome and transcriptome variation

Sexual dimorphism is one of the most prevalent, and often the most extreme, examples of phenotypic variation within species, and arises primarily from genomic variation that is shared between females and males. Many sexual dimorphisms arise through sex differences in gene expression, and sex-biased expression is one way that a single, shared genome can generate multiple, distinct phenotypes. Although many sexual dimorphisms are expected to result from sexual selection, and many studies have invoked the possible role of sexual selection to explain sex-specific traits, the role of sexual selection in the evolution of sexually dimorphic gene expression remains difficult to differentiate from other forms of sex-specific selection. In this Review, we propose a holistic framework for the study of sex-specific selection and transcriptome evolution. We advocate for a comparative approach, across tissues, developmental stages and species, which incorporates an understanding of the molecular mechanisms, including genomic variation and structure, governing gene expression. Such an approach is expected to yield substantial insights into the evolution of genetic variation and have important applications in a variety of fields, including ecology, evolution and behaviour.

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.

Evolutionary dynamics of sex-biased gene expression in a young XY system: insights from the brown alga genus Fucus

Sex-biased gene expression is considered to be an underlying cause of sexually dimorphic traits. Although the nature and degree of sex-biased expression have been well documented in several animal and plant systems, far less is known about the evolution of sex-biased genes in more distant eukaryotic groups. Here, we investigate sex-biased gene expression in two brown algal dioecious species, Fucus serratus and Fucus vesiculosus, where male heterogamety (XX/XY) has recently emerged. We find that in contrast to evolutionary distant plant and animal lineages, male-biased genes do not experience high turnover rates, but instead reveal remarkable conservation of bias and expression levels between the two species, suggesting their importance in sexual differentiation. Genes with consistent male bias were enriched in functions related to gamete production, along with sperm competition and include three flagellar proteins under positive selection. We present one of the first reports, outside of the animal kingdom, showing that male-biased genes display accelerated rates of coding sequence evolution compared with female-biased or unbiased genes. Our results imply that evolutionary forces affect male and female sex-biased genes differently on structural and regulatory levels, resulting in unique properties of differentially expressed transcripts during reproductive development in Fucus algae.

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.

Sexual dimorphism in chronic respiratory diseases

Sex differences in susceptibility, severity, and progression are prevalent for various diseases in multiple organ systems. This phenomenon is particularly apparent in respiratory diseases. Asthma demonstrates an age-dependent pattern of sexual dimorphism. However, marked differences between males and females exist in other pervasive conditions such as chronic obstructive pulmonary disease (COPD) and lung cancer. The sex hormones estrogen and testosterone are commonly considered the primary factors causing sexual dimorphism in disease. However, how they contribute to differences in disease onset between males and females remains undefined. The sex chromosomes are an under-investigated fundamental form of sexual dimorphism. Recent studies highlight key X and Y-chromosome-linked genes that regulate vital cell processes and can contribute to disease-relevant mechanisms. This review summarises patterns of sex differences in asthma, COPD and lung cancer, highlighting physiological mechanisms causing the observed dimorphism. We also describe the role of the sex hormones and present candidate genes on the sex chromosomes as potential factors contributing to sexual dimorphism in disease.

Multi-tissue RNA-Seq Analysis and Long-read-based Genome Assembly Reveal Complex Sex-specific Gene Regulation and Molecular Evolution in the Manila Clam

The molecular factors and gene regulation involved in sex determination and gonad differentiation in bivalve molluscs are unknown. It has been suggested that doubly uniparental inheritance (DUI) of mitochondria may be involved in these processes in species such as the ubiquitous and commercially relevant Manila clam, Ruditapes philippinarum. We present the first long-read-based de novo genome assembly of a Manila clam, and a RNA-Seq multi-tissue analysis of 15 females and 15 males. The highly contiguous genome assembly was used as reference to investigate gene expression, alternative splicing, sequence evolution, tissue-specific co-expression networks, and sexual contrasting SNPs. Differential expression (DE) and differential splicing (DS) analyses revealed sex-specific transcriptional regulation in gonads, but not in somatic tissues. Co-expression networks revealed complex gene regulation in gonads, and genes in gonad-associated modules showed high tissue specificity. However, male gonad-associated modules showed contrasting patterns of sequence evolution and tissue specificity. One gene set was related to the structural organization of male gametes and presented slow sequence evolution but high pleiotropy, whereas another gene set was enriched in reproduction-related processes and characterized by fast sequence evolution and tissue specificity. Sexual contrasting SNPs were found in genes overrepresented in mitochondrial-related functions, providing new candidates for investigating the relationship between mitochondria and sex in DUI species. Together, these results increase our understanding of the role of DE, DS, and sequence evolution of sex-specific genes in an understudied taxon. We also provide resourceful genomic data for studies regarding sex diagnosis and breeding in bivalves.

Masculinizer-induced dosage compensation is achieved by transcriptional downregulation of both copies of Z-linked genes in the silkworm, Bombyx mori

The evolution of dosage compensation produces similar expression of sex-linked and autosomal genes in the heterogametic sex. The silkworm (Bombyx mori), a lepidopteran insect, has a female heterogametic WZ sex determination system. A Z-linked gene, Masculinizer (Masc), is the primary determinant of maleness and dosage compensation in B. mori. However, it remains unknown whether one of the two Z chromosomes is inactivated or both Z chromosomes are suppressed in B. mori males. Hence, we performed transcriptome analysis using hybrids between two B. mori strains and analysed allele-specific expression to distinguish these alternatives. Our analysis revealed that genes on both the maternal and paternal Z chromosomes are transcriptionally upregulated in Masc knocked down males. We therefore conclude that both Z chromosomes are transcriptionally downregulated in B. mori males, similar to the system in Caenorhabditis elegans.

Transcriptome analysis of different life-history stages and screening of male-biased genes in Daphnia sinensis

Background

In the life history of Daphnia, the reproductive mode of parthenogenesis and sexual reproduction alternate in aquatic ecosystem, which are often affected by environmental and genetic factors. Recently, the sex-biased genes are of great significance for clarifying the origin and evolution of reproductive transformation and the molecular regulation mechanism of sex determination in Daphnia. Although some genes on reproductive transition of Daphnia had been researched, molecular mechanism on the maintenance of sexually dimorphic phenotypes of Daphnia are still not well known, including differentially expressed genes in different life-history stages.

Results

In this study, four life-history stages of Daphnia sinensis, juvenile female (JF), parthenogenetic female (PF), sexual female (SF) and male (M), were performed for transcriptome, and male-biased genes were screened. A total of 110437 transcripts were obtained and assembled into 22996 unigenes. In the four life-history stages (JF, PF, SF and M), the number of unique unigenes is respectively 2863, 445, 437 and 586, and the number of common unigenes is 9708. The differentially expressed genes (DEGs) between male and other three female stages (M vs JF, M vs PF and M vs SF) were 4570, 4358 and 2855, respectively. GO gene enrichment analysis showed that the up-regulated genes in male were mainly enriched in hydrolase activity and peptidase activity. Thirty-six genes in male were significantly higher expression than in the three female stages, including one Doublesex (Dsx) gene, one laminin gene, five trypsin genes and one serine protease genes, and one chitin synthase gene and two chitinase genes.

Conclusions

Our results showed that thirty-six candidate genes may be as the male-biased genes involving in the maintenance of sexually dimorphic phenotypes. This work will provide a reference for further exploring the functional genes related to sex differentiation in Daphnia species. Moreover, according to previous investigations, we thought that the expression level of functional genes may be related to the life-history stages of organisms, and may be also affected by different Daphnia species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08824-x.

A mammalian messenger RNA sex determination method from humpback whale (Megaptera novaeangliae) blubber biopsies

The large size of free-ranging mysticetes, such as humpback whales (Megaptera novaeangliae), make capture and release health assessments unfeasible for conservation research. However, individual energetic condition or reproductive health may be assessed from the gene expression of remotely biopsied tissue. To do this, researchers must reliably extract RNA and interpret gene expression measurements within the context of an individual's sex. Here, we outline an RNA extraction protocol from blubber tissue and describe a novel mammalian RNA sex determination method. Our method consists of a duplex reverse transcription-quantitative (real-time) polymerase chain reaction (RT-qPCR) with primer sets for a control gene (ACTB) and the X-chromosome inactivation gene (XIST). Products of each RT-qPCR had distinct melting temperature profiles based on the presence (female) or absence (male) of the XIST transcript. Using high-resolution melt analysis, reactions were sorted into one of two clusters (male/female) based on their melting profiles. We validated the XIST method by comparing results with a standard DNA-based method. With adequate quantities of RNA (minimum of approx. 9 ng µl^-1), the XIST sex determination method shows 100% agreement with traditional DNA sex determination. Using the XIST method, future cetacean health studies can interpret gene expression within the context of an individual's sex, all from a single extraction.

Transcriptome analysis of sex-biased gene expression in the spotted-wing Drosophila, Drosophila suzukii (Matsumura)

Sexual dimorphism occurs widely throughout insects and has profound influences on evolutionary path. Sex-biased genes are considered to account for most of phenotypic differences between sexes. In order to explore the sex-biased genes potentially associated with sexual dimorphism and sexual development in Drosophila suzukii, a major devastating and invasive crop pest, we conducted whole-organism transcriptome profiling and sex-biased gene expression analysis on adults of both sexes. We identified transcripts of genes involved in several sex-specific physiological and functional processes, including transcripts involved in sex determination, reproduction, olfaction, and innate immune signals. A total of 11,360 differentially expressed genes were identified in the comparison, and 1,957 differentially expressed genes were female-biased and 4,231 differentially expressed genes were male-biased. The pathway predominantly enriched for differentially expressed genes was related to spliceosome, which might reflect the differences in the alternative splicing mechanism between males and females. Twenty-two sex determination and 16 sex-related reproduction genes were identified, and expression pattern analysis revealed that the majority of genes were differentially expressed between sexes. Additionally, the differences in sex-specific olfactory and immune processes were analyzed and the sex-biased expression of these genes may play important roles in pheromone and odor detection, and immune response. As a valuable dataset, our sex-specific transcriptomic data can significantly contribute to the fundamental elucidation of the molecular mechanisms of sexual dimorphism in fruit flies, and may provide candidate genes potentially useful for the development of genetic sexing strains, an important tool for sterile insect technique applications against this economically important species.

Epigenetic Regulation of Phenotypic Sexual Plasticity Inducing Skewed Sex Ratio in Zebrafish

The plasticity of sexual phenotype in response to environmental conditions results in biased sex ratios, and their variation has an effect on population dynamics. Epigenetic modifications can modulate sex ratio variation in species, where sex is determined by genetic and environmental factors. However, the role of epigenetic mechanisms underlying skewed sex ratios is far from being clear and is still an object of debate in evolutionary developmental biology. In this study, we used zebrafish as a model animal to investigate the effect of DNA methylation on sex ratio variation in sex-biased families in response to environmental temperature. Two sex-biased families with a significant difference in sex ratio were selected for genome-wide DNA methylation analysis using reduced representation bisulfite sequencing (RRBS). The results showed significant genome-wide methylation differences between male-biased and female-biased families, with a greater number of methylated CpG sites in testes than ovaries. Likewise, pronounced differences between testes and ovaries were identified within both families, where the male-biased family exhibited a higher number of methylated sites than the female-biased family. The effect of temperature showed more methylated positions in the high incubation temperature than the control temperature. We found differential methylation of many reproduction-related genes (e.g., sox9a, nr5a2, lhx8a, gata4) and genes involved in epigenetic mechanisms (e.g., dnmt3bb.1, dimt1l, hdac11, h1m) in both families. We conclude that epigenetic modifications can influence the sex ratio variation in zebrafish families and may generate skewed sex ratios, which could have a negative consequence for population fitness in species with genotype-environment interaction sex-determining system under rapid environmental changes.

Dioecy and chromosomal sex determination are maintained through allopolyploid speciation in the plant genus Mercurialis

Polyploidization may precipitate dramatic changes to the genome, including chromosome rearrangements, gene loss, and changes in gene expression. In dioecious plants, the sex-determining mechanism may also be disrupted by polyploidization, with the potential evolution of hermaphroditism. However, while dioecy appears to have persisted through a ploidy transition in some species, it is unknown whether the newly formed polyploid maintained its sex-determining system uninterrupted, or whether dioecy re-evolved after a period of hermaphroditism. Here, we develop a bioinformatic pipeline using RNA-sequencing data from natural populations to demonstrate that the allopolyploid plant Mercurialis canariensis directly inherited its sex-determining region from one of its diploid progenitor species, M. annua, and likely remained dioecious through the transition. The sex-determining region of M. canariensis is smaller than that of its diploid progenitor, suggesting that the non-recombining region of M. annua expanded subsequent to the polyploid origin of M. canariensis. Homeologous pairs show partial sexual subfunctionalization. We discuss the possibility that gene duplicates created by polyploidization might contribute to resolving sexual antagonism.

Sex-specific variation in R-loop formation in Drosophila melanogaster

R-loops are three-stranded nucleotide structures consisting of a DNA:RNA hybrid and a displaced ssDNA non-template strand. Previous work suggests that R-loop formation is primarily determined by the thermodynamics of DNA:RNA binding, which are governed by base composition (e.g., GC skew) and transcription-induced DNA superhelicity. However, R-loops have been described at genomic locations that lack these properties, suggesting that they may serve other context-specific roles. To better understand the genetic determinants of R-loop formation, we have characterized the Drosophila melanogaster R-loop landscape across strains and between sexes using DNA:RNA immunoprecipitation followed by high-throughput sequencing (DRIP-seq). We find that R-loops are associated with sequence motifs that are G-rich or exhibit G/C skew, as well as highly expressed genes, tRNAs, and small nuclear RNAs, consistent with a role for DNA sequence and torsion in R-loop specification. However, we also find motifs associated with R-loops that are A/T-rich and lack G/C skew as well as a subset of R-loops that are enriched in polycomb-repressed chromatin. Differential enrichment analysis reveals a small number of sex-biased R-loops: while non-differentially enriched and male-enriched R-loops form at similar genetic features and chromatin states and contain similar sequence motifs, female-enriched R-loops form at unique genetic features, chromatin states, and sequence motifs and are associated with genes that show ovary-biased expression. Male-enriched R-loops are most abundant on the dosage-compensated X chromosome, where R-loops appear stronger compared to autosomal R-loops. R-loop-containing genes on the X chromosome are dosage-compensated yet show lower MOF binding and reduced H4K16ac compared to R-loop-absent genes, suggesting that H4K16ac or MOF may attenuate R-loop formation. Collectively, these results suggest that R-loop formation in vivo is not fully explained by DNA sequence and topology and raise the possibility that a distinct subset of these hybrid structures plays an important role in the establishment and maintenance of epigenetic differences between sexes.

Ontogenetic Change in Male Expression of Testosterone-Responsive Genes Contributes to the Emergence of Sex-Biased Gene Expression in Anolis sagrei

Sex differences in gene expression tend to increase with age across a variety of species, often coincident with the development of sexual dimorphism and maturational changes in hormone levels. However, because most transcriptome-wide characterizations of sexual divergence are framed as comparisons of sex-biased gene expression across ages, it can be difficult to determine the extent to which age-biased gene expression within each sex contributes to the emergence of sex-biased gene expression. Using RNAseq in the liver of the sexually dimorphic brown anole lizard (Anolis sagrei), we found that a pronounced increase in sex-biased gene expression with age was associated with a much greater degree of age-biased gene expression in males than in females. This pattern suggests that developmental changes in males, such as maturational increases in circulating testosterone, contribute disproportionately to the ontogenetic emergence of sex-biased gene expression. To test this hypothesis, we used four different experimental contrasts to independently characterize sets of genes whose expression differed as a function of castration and/or treatment with exogenous testosterone. We found that genes that were significantly male-biased in expression or upregulated as males matured tended to be upregulated by testosterone, whereas genes that were female-biased or downregulated as males matured tended to be downregulated by testosterone. Moreover, the first two principal components describing multivariate gene expression indicated that exogenous testosterone reversed many of the feminizing effects of castration on the liver transcriptome of maturing males. Collectively, our results suggest that developmental changes that occur in males contribute disproportionately to the emergence of sex-biased gene expression in the Anolis liver, and that many of these changes are orchestrated by androgens such as testosterone.

Sex determination mechanisms and sex control approaches in aquaculture animals

Aquaculture is one of the most efficient modes of animal protein production and plays an important role in global food security. Aquaculture animals exhibit extraordinarily diverse sexual phenotypes and underlying mechanisms, providing an ideal system to perform sex determination research, one of the important areas in life science. Moreover, sex is also one of the most valuable traits because sexual dimorphism in growth, size, and other economic characteristics commonly exist in aquaculture animals. Here, we synthesize current knowledge of sex determination mechanisms, sex chromosome evolution, reproduction strategies, and sexual dimorphism, and also review several approaches for sex control in aquaculture animals, including artificial gynogenesis, application of sex-specific or sex chromosome-linked markers, artificial sex reversal, as well as gene editing. We anticipate that better understanding of sex determination mechanisms and innovation of sex control approaches will facilitate sustainable development of aquaculture.

Transcriptome analysis in the silkworm Bombyx mori overexpressing piRNA-resistant Masculinizer gene

Dosage compensation is a process that produces a similar expression of sex-linked and autosomal genes. In the silkworm Bombyx mori with a WZ sex-determination system, the expression from the single Z in WZ females matches that of ZZ males due to the suppression of Z-linked genes in males. A primary maleness determinant gene, Masculinizer (Masc), is also required for dosage compensation. In females, silkworm Piwi is complexed with the W chromosome-derived female-specific Feminizer (Fem) PIWI-interacting RNA (piRNA) and cleaves Masc mRNA. When Fem piRNA-resistant Masc cDNA (Masc-R) is overexpressed in both sexes, only female larvae are dead during the larval stage. In this study, transcriptome analysis was performed in neonate larvae to examine the effects of Masc-R overexpression on a global gene expression profile. Z-linked genes were globally repressed in Masc-R-overexpressing females due to force-driven dosage compensation. In contrast, Masc-R overexpression had little effect on the expression of Z-linked genes and the male-specific isoform of B. mori insulin-like growth factor II mRNA-binding protein in males, indicating that excessive Masc expression strengthens neither dosage compensation nor maleness in males. Fourteen genes were differentially expressed between Masc-R-overexpressing and control neonate larvae in both sexes, suggesting Masc functions other than dosage compensation and masculinization.

Chromatin landscape associated with sexual differentiation in a UV sex determination system

In many eukaryotes, such as dioicous mosses and many algae, sex is determined by UV sex chromosomes and is expressed during the haploid phase of the life cycle. In these species, the male and female developmental programs are initiated by the presence of the U- or V-specific regions of the sex chromosomes but, as in XY and ZW systems, sexual differentiation is largely driven by autosomal sex-biased gene expression. The mechanisms underlying the regulation of sex-biased expression of genes during sexual differentiation remain elusive. Here, we investigated the extent and nature of epigenomic changes associated with UV sexual differentiation in the brown alga Ectocarpus, a model UV system. Six histone modifications were quantified in near-isogenic lines, leading to the identification of 16 chromatin signatures across the genome. Chromatin signatures correlated with levels of gene expression and histone PTMs changes in males versus females occurred preferentially at genes involved in sex-specific pathways. Despite the absence of chromosome scale dosage compensation and the fact that UV sex chromosomes recombine across most of their length, the chromatin landscape of these chromosomes was remarkably different to that of autosomes. Hotspots of evolutionary young genes in the pseudoautosomal regions appear to drive the exceptional chromatin features of UV sex chromosomes.

Metabolic Sex Dimorphism of the Brain at the Gene, Cell, and Tissue Level

The palpable observation in the sex bias of disease prevalence in the CNS has fascinated scientists for several generations. Brain sex dimorphism has been visualized by imaging and analytical tools at the tissue, cellular, and molecular levels. Recent work highlighted the specificity of such sex bias in the brain and its subregions, offering a unique lens through which disease pathogenesis can be investigated. The brain is the largest consumer of energy in the body and provides a unique metabolic environment for diverse lineages of cells. Immune cells are increasingly recognized as an integral part of brain physiology, and their function depends on metabolic homeostasis. This review focuses on metabolic sex dimorphism in brain tissue, resident, and infiltrating immune cells. In this context, we highlight the relevance of recent advances in metabolomics and RNA sequencing technologies at the single cell resolution and the development of novel computational approaches.

Faster Rates of Molecular Sequence Evolution in Reproduction-Related Genes and in Species with Hypodermic Sperm Morphologies

Sexual selection drives the evolution of many striking behaviors and morphologies and should leave signatures of selection at loci underlying these phenotypes. However, although loci thought to be under sexual selection often evolve rapidly, few studies have contrasted rates of molecular sequence evolution at such loci across lineages with different sexual selection contexts. Furthermore, work has focused on separate sexed animals, neglecting alternative sexual systems. We investigate rates of molecular sequence evolution in hermaphroditic flatworms of the genus Macrostomum. Specifically, we compare species that exhibit contrasting sperm morphologies, strongly associated with multiple convergent shifts in the mating strategy, reflecting different sexual selection contexts. Species donating and receiving sperm in every mating have sperm with bristles, likely to prevent sperm removal. Meanwhile, species that hypodermically inject sperm lack bristles, potentially as an adaptation to the environment experienced by hypodermic sperm. Combining functional annotations from the model, Macrostomum lignano, with transcriptomes from 93 congeners, we find genus-wide faster sequence evolution in reproduction-related versus ubiquitously expressed genes, consistent with stronger sexual selection on the former. Additionally, species with hypodermic sperm morphologies had elevated molecular sequence evolution, regardless of a gene's functional annotation. These genome-wide patterns suggest reduced selection efficiency following shifts to hypodermic mating, possibly due to higher selfing rates in these species. Moreover, we find little evidence for convergent amino acid changes across species. Our work not only shows that reproduction-related genes evolve rapidly also in hermaphroditic animals, but also that well-replicated contrasts of different sexual selection contexts can reveal underappreciated genome-wide effects.

Clustering of Sex-Biased Genes and Transposable Elements in the Genome of the Medaka Fish Oryzias latipes

Although genes with similar expression patterns are sometimes found in the same genomic regions, almost nothing is known about the relative organization in genomes of genes and transposable elements (TEs), which might influence each other at the regulatory level. In this study, we used transcriptomic data from male and female gonads of the Japanese medaka Oryzias latipes to define sexually biased genes and TEs and analyze their relative genomic localization. We identified 20,588 genes expressed in the adult gonads of O. latipes. Around 39% of these genes are differentially expressed between male and female gonads. We further analyzed the expression of TEs using the program SQuIRE and showed that more TE copies are overexpressed in testis than in ovaries (36% vs. 10%, respectively). We then developed a method to detect genomic regions enriched in testis- or ovary-biased genes. This revealed that sex-biased genes and TEs are not randomly distributed in the genome and a part of them form clusters with the same expression bias. We also found a correlation of expression between TE copies and their closest genes, which increases with decreasing intervening distance. Such a genomic organization suggests either that TEs hijack the regulatory sequences of neighboring sexual genes, allowing their expression in germ line cells and consequently new insertions to be transmitted to the next generation, or that TEs are involved in the regulation of sexual genes, and might therefore through their mobility participate in the rewiring of sex regulatory networks.

High rates of evolution preceded shifts to sex-biased gene expression in Leucadendron, the most sexually dimorphic angiosperms

Differences between males and females are usually more subtle in dioecious plants than animals, but strong sexual dimorphism has evolved convergently in the South African Cape plant genus Leucadendron. Such sexual dimorphism in leaf size is expected largely to be due to differential gene expression between the sexes. We compared patterns of gene expression in leaves among 10 Leucadendron species across the genus. Surprisingly, we found no positive association between sexual dimorphism in morphology and the number or the percentage of sex-biased genes (SBGs). Sex bias in most SBGs evolved recently and was species specific. We compared rates of evolutionary change in expression for genes that were sex biased in one species but unbiased in others and found that SBGs evolved faster in expression than unbiased genes. This greater rate of expression evolution of SBGs, also documented in animals, might suggest the possible role of sexual selection in the evolution of gene expression. However, our comparative analysis clearly indicates that the more rapid rate of expression evolution of SBGs predated the origin of bias, and shifts towards bias were depleted in signatures of adaptation. Our results are thus more consistent with the view that sex bias is simply freer to evolve in genes less subject to constraints in expression level.

Harnessing natural variation to identify cis regulators of sex-biased gene expression in a multi-strain mouse liver model

Sex differences in gene expression are widespread in the liver, where many autosomal factors act in tandem with growth hormone signaling to regulate individual variability of sex differences in liver metabolism and disease. Here, we compare hepatic transcriptomic and epigenetic profiles of mouse strains C57BL/6J and CAST/EiJ, representing two subspecies separated by 0.5-1 million years of evolution, to elucidate the actions of genetic factors regulating liver sex differences. We identify 144 protein coding genes and 78 lncRNAs showing strain-conserved sex bias; many have gene ontologies relevant to liver function, are more highly liver-specific and show greater sex bias, and are more proximally regulated than genes whose sex bias is strain-dependent. The strain-conserved genes include key growth hormone-dependent transcriptional regulators of liver sex bias; however, three other transcription factors, Trim24, Tox, and Zfp809, lose their sex-biased expression in CAST/EiJ mouse liver. To elucidate the observed strain specificities in expression, we characterized the strain-dependence of sex-biased chromatin opening and enhancer marks at cis regulatory elements (CREs) within expression quantitative trait loci (eQTL) regulating liver sex-biased genes. Strikingly, 208 of 286 eQTLs with strain-specific, sex-differential effects on expression were associated with a complete gain, loss, or reversal of the sex differences in expression between strains. Moreover, 166 of the 286 eQTLs were linked to the strain-dependent gain or loss of localized sex-biased CREs. Remarkably, a subset of these CREs apparently lacked strain-specific genetic variants yet showed coordinated, strain-dependent sex-biased epigenetic regulation. Thus, we directly link hundreds of strain-specific genetic variants to the high variability in CRE activity and expression of sex-biased genes and uncover underlying genetically-determined epigenetic states controlling liver sex bias in genetically diverse mouse populations.

A partially sex-reversed giant kelp sheds light into the mechanisms of sexual differentiation in a UV sexual system

In UV sexual systems, sex is determined during the haploid phase of the life cycle and males have a V chromosome whereas females have a U chromosome. Previous work in the brown alga Ectocarpus revealed that the V chromosome has a dominant role in male sex determination and suggested that the female developmental programme may occur by 'default'. Here, we describe the identification of a genetically male giant kelp strain presenting phenotypic features typical of a female, despite lacking the U-specific region. The conversion to the female developmental programme is however incomplete, because gametes of this feminized male are unable to produce the sperm-attracting pheromone lamoxirene. We identify the transcriptomic patterns underlying the male and female specific developmental programmes, and show that the phenotypic feminization is associated with both feminization and de-masculinization of gene expression patterns. Importantly, the feminization phenotype was associated with dramatic downregulation of two V-specific genes including a candidate male-determining gene. Our results reveal the transcriptional changes associated with sexual differentiation in a UV system, and contribute to disentangling the role of sex-linked and autosomal gene expression in the initiation of sex-specific developmental programmes. Overall, the data presented here imply that the U-specific region is not required to initiate the female developmental programme, but is critical to produce fully functional eggs, arguing against the idea that female is the 'default' sex in this species.

Sex-biased gene expression and recent sex chromosome turnover

Cichlids are well known for their propensity to radiate generating arrays of morphologically and ecologically diverse species in short evolutionary time. Following this rapid evolutionary pace, cichlids show high rates of sex chromosome turnover. We here studied the evolution of sex-biased gene (SBG) expression in 14 recently diverged taxa of the Lake Tanganyika Tropheini cichlids, which show different XY sex chromosomes. Across species, sex chromosome sequence divergence predates divergence in expression between the sexes. Only one sex chromosome, the oldest, showed signs of demasculinization in gene expression and potentially contribution to the resolution of sexual conflict. SBGs in general showed high rates of turnovers and evolved mostly under drift. Sexual selection did not shape the rapid evolutionary changes of SBGs. Male-biased genes evolved faster than female-biased genes, which seem to be under more phylogenetic constraint. We found a relationship between the degree of sex bias and sequence evolution driven by sequence differences among the sexes. Consistent with other species, strong sex bias towards sex-limited expression contributes to resolving sexual conflict in cichlids. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part II)'.

Intraspecific Transcriptome Variation and Sex-Biased Expression in Anopheles arabiensis

The magnitude and functional patterns of intraspecific transcriptional variation in the anophelines, including those of sex-biased genes underlying sex-specific traits relevant for malaria transmission, remain understudied. As a result, how changes in expression levels drive adaptation in these species is poorly understood. We sequenced the female, male, and larval transcriptomes of three populations of Anopheles arabiensis from Burkina Faso. One-third of the genes were differentially expressed between populations, often involving insecticide resistance-related genes in a sample type-specific manner, and with the females showing the largest number of differentially expressed genes. At the genomic level, the X chromosome appears depleted of differentially expressed genes compared with the autosomes, chromosomes harboring inversions do not exhibit evidence for enrichment of such genes, and genes that are top contributors to functional enrichment patterns of population differentiation tend to be clustered in the genome. Further, the magnitude of variation for the sex expression ratio across populations did not substantially differ between male- and female-biased genes, except for some populations in which male-limited expressed genes showed more variation than their female counterparts. In fact, female-biased genes exhibited a larger level of interpopulation variation than male-biased genes, both when assayed in males and females. Beyond uncovering the extensive adaptive potential of transcriptional variation in An. Arabiensis, our findings suggest that the evolutionary rate of changes in expression levels on the X chromosome exceeds that on the autosomes, while pointing to female-biased genes as the most variable component of the An. Arabiensis transcriptome.

Sex-Biased Gene Expression of Mesobuthus martensii Collected from Gansu Province, China, Reveals Their Different Therapeutic Potentials

The scorpions, named Mesobuthus martensii, commonly called Quanxie (全蝎) in Chinese, have been widely used as one of the animal medicines for more than 1,000 years because of the strong toxicity of their venoms. Meanwhile, scorpions are sexually dimorphic in appearance, and many exhibit traits associated with sex-biased gene expression, including maternal care, mating competition, female mating choices, ecology, and even venom composition and lethality. This study aims to explore the differences in composition of the venom of scorpions of different sex using the method of transcriptomics. Whole de novo transcriptomes were performed on the samples of M. martensii captured from Gansu Province to identify their sex-biased gene expression. The conserved CO-1 sequences of the captured samples matched that of M. martensii. A total of 8,444 (35.15%), 7,636 (31.78%), 8,510 (35.42%), 7,840 (32.63%), 9,980 (41.54%), and 11,829 (49.23%) unigenes were annotated with GO, KEGG, Pfam, Swissprot, eggNOG, and NR databases. Moreover, a total of 43 metalloproteases, 40 potassium channel toxins, 24 phospholipases, 12 defensins, 10 peroxiredoxins, 9 cysteine proteinase inhibitors, 7 serine protease inhibitors, 6 sodium channel toxins, 2 NDBPs, 1 calcium channel toxin, 1 waprin-like peptide, 1 antibacterial peptide, 1 antimicrobial peptide, and 1 anticoagulant peptide were screened out. With the fold change of 2 and 0.5, p value < 0.01, and q value < 0.05 as thresholds, a total of 41 out of 157 (26.11%) toxin-related unigenes had significant differential expression, and this ratio was much higher than the ratio of differentially expressed unigenes out of all annotated ones (8.84%). Of these differentially expressed toxins, 28 were upregulated and occupied the majority, up to 68.30%. The female scorpions showed more upregulated unigenes that annotated with toxins and had the potential to be used as more effective therapeutic drugs. In addition, this method of omics can be further used as a useful way to identify the difference between female and male toxic animals.

Sex-Specific Development in Haplodiploid Honeybee Is Controlled by the Female-Embryo-Specific Activation of Thousands of Intronic LncRNAs

Embryonic development depends on a highly coordinated shift in transcription programs known as the maternal-to-zygotic transition (MZT). It remains unclear how haploid and diploid embryo coordinate their genomic activation and embryonic development during MZT in haplodiploid animals. Here, we applied a single-embryo RNA-seq approach to characterize the embryonic transcriptome dynamics in haploid males vs. diploid females of the haplodiploid insect honeybee (Apis mellifera). We observed typical zygotic genome activation (ZGA) occurred in three major waves specifically in female honeybee embryos; haploid genome activation was much weaker and occurred later. Strikingly, we also observed three waves of transcriptional activation for thousands of long non-coding transcripts (lncRNA), 73% of which are transcribed from intronic regions and 65% were specific to female honeybee embryos. These findings support a model in which introns encode thousands of lncRNAs that are expressed in a diploid-embryo-specific and ZGA-triggered manner that may have potential functions to regulate gene expression during early embryonic development in the haplodiploid insect honeybee.

Joint identification of sex and sex-linked scaffolds in non-model organisms using low depth sequencing data

Being able to assign sex to individuals and identify autosomal and sex-linked scaffolds are essential in most population genomic analyses. Non-model organisms often have genome assemblies at scaffold-level and lack characterization of sex-linked scaffolds. Previous methods to identify sex and sex-linked scaffolds have relied on synteny between the non-model organism and a closely related species or prior knowledge about the sex of the samples to identify sex-linked scaffolds. In the latter case, the difference in depth of coverage between the autosomes and the sex chromosomes are used. Here, we present "sex assignment through coverage" (SATC), a method to assign sex to samples and identify sex-linked scaffolds from next generation sequencing (NGS) data. The method works for species with a homogametic/heterogametic sex determination system and only requires a scaffold-level reference assembly and sampling of both sexes with whole genome sequencing (WGS) data. We use the sequencing depth distribution across scaffolds to jointly identify: (i) male and female individuals, and (ii) sex-linked scaffolds. This is achieved through projecting the scaffold depths into a low-dimensional space using principal component analysis (PCA) and subsequent Gaussian mixture clustering. We demonstrate the applicability of our method using data from five mammal species and a bird species complex. The method is freely available at https://github.com/popgenDK/SATC as R code and a graphical user interface (GUI).

The transcriptomic blueprint of molt in rooster using various tissues from Ginkkoridak (Korean long-tailed chicken)

BACKGROUND

Annual molt is a critical stage in the life cycle of birds. Although the most extensively documented aspects of molt are the renewing of plumage and the remodeling of the reproductive tract in laying hens, in chicken, molt deeply affects various tissues and physiological functions. However, with exception of the reproductive tract, the effect of molt on gene expression across the tissues known to be affected by molt has to date never been investigated. The present study aimed to decipher the transcriptomic effects of molt in Ginkkoridak, a Korean long-tailed chicken. Messenger RNA data available across 24 types of tissue samples (9 males) and a combination of mRNA and miRNA data on 10 males and 10 females blood were used.

RESULTS

The impact of molt on gene expression and gene transcript usage appeared to vary substantially across tissues types in terms of histological entities or physiological functions particularly related to nervous system. Blood was the tissue most affected by molt in terms of differentially expressed genes in both sexes, closely followed by meninges, bone marrow and heart. The effect of molt in blood appeared to differ between males and females, with a more than fivefold difference in the number of down-regulated genes between both sexes. The blueprint of molt in roosters appeared to be specific to tissues or group of tissues, with relatively few genes replicating extensively across tissues, excepted for the spliceosome genes (U1, U4) and the ribosomal proteins (RPL21, RPL23). By integrating miRNA and mRNA data, when chickens molt, potential roles of miRNA were discovered such as regulation of neurogenesis, regulation of immunity and development of various organs. Furthermore, reliable candidate biomarkers of molt were found, which are related to cell dynamics, nervous system or immunity, processes or functions that have been shown to be extensively modulated in response to molt.

CONCLUSIONS

Our results provide a comprehensive description at the scale of the whole organism deciphering the effects of molt on the transcriptome in chicken. Also, the conclusion of this study can be used as a valuable resource in transcriptome analyses of chicken in the future and provide new insights related to molt.

Widespread cryptic variation in genetic architecture between the sexes

The majority of the genome is shared between the sexes, and it is expected that the genetic architecture of most traits is shared as well. This common architecture has been viewed as a major source of constraint on the evolution of sexual dimorphism (SD). SD is nonetheless common in nature, leading to assumptions that it results from differential regulation of shared genetic architecture. Here, we study the effect of thousands of gene knockout mutations on 202 mouse phenotypes to explore how regulatory variation affects SD. We show that many traits are dimorphic to some extent, and that a surprising proportion of knockouts have sex-specific phenotypic effects. Many traits, regardless whether they are monomorphic or dimorphic, harbor cryptic differences in genetic architecture between the sexes, resulting in sexually discordant phenotypic effects from sexually concordant regulatory changes. This provides an alternative route to dimorphism through sex-specific genetic architecture, rather than differential regulation of shared architecture.

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.