The evolutionary causes and consequences of sex-biased gene expression

Influence of Sex and Gender on Musculoskeletal Conditions and How They Are Reported

ABSTRACT

There is increasing evidence that musculoskeletal tissues are differentiallys regulated by sex hormones in males and females. The influence of sex hormones, in addition to other sex-based differences such as in anatomical alignment and immune-system function, impact the prevalence and severity of disease as well as the types of injuries that affect the musculoskeletal system and the outcomes of prevention measures and treatment. Literature specifically addressing sex differences related to the musculoskeletal system is limited, underscoring the imperative for both basic and clinical research on this topic. This review highlights areas of research that have implications for bone and cartilage health, including growth and development, sports injuries, osteoarthritis, osteoporosis, and bone frailty. It is clear that important aspects of the musculoskeletal system have been understudied. Consideration of how sex hormone therapy will affect musculoskeletal tissues in prepuberty, during puberty, and in adults is vital, yet little is known. The purpose of this article is to foster awareness and interest in advancing our understanding of how sex differences influence orthopaedic practice.

Identification of reproductive sex-biased gene expression in Asparagopsis taxiformis (lineage 6) gametophytes

The sub-tropical red seaweed Asparagopsis taxiformis is of significant interest due to its ability to store halogenated compounds, including bromoform, which can mitigate methane production in ruminants. Significant scale-up of aquaculture production of this seaweed is required; however, relatively little is known about the molecular mechanisms that control fundamental physiological processes, including the regulatory factors that determine sexual dimorphism in gametophytes. In this study, we used comparative RNA-sequencing analysis between different morphological parts of mature male and female A. taxiformis (lineage 6) gametophytes that resulted in greater number of sex-biased gene expression in tips (containing the reproductive structures for both sexes), compared with the somatic main axis and rhizomes. Further comparative RNA-seq against immature tips was used to identify 62 reproductive sex-biased genes (59 male-biased, 3 female-biased). Of the reproductive male-biased genes, 46% had an unknown function, while others were predicted to be regulatory factors and enzymes involved in signaling. We found that bromoform content obtained from female samples (8.5 ± 1.0 mg·g-1 dry weight) was ~10% higher on average than that of male samples (6.5 ± 1.0 mg·g-1 dry weight), although no significant difference was observed (p > 0.05). There was also no significant difference in the marine bromoform biosynthesis locus gene expression. In summary, our comparative RNA-sequencing analysis provides a first insight into the potential molecular factors relevant to gametogenesis and sexual differentiation in A. taxiformis, with potential benefits for identification of sex-specific markers.

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.

Positive Selection Drives cis-regulatory Evolution Across the Threespine Stickleback Y Chromosome

Allele-specific gene expression evolves rapidly on heteromorphic sex chromosomes. Over time, the accumulation of mutations on the Y chromosome leads to widespread loss of gametolog expression, relative to the X chromosome. It remains unclear if expression evolution on degrading Y chromosomes is primarily driven by mutations that accumulate through processes of selective interference, or if positive selection can also favor the down-regulation of coding regions on the Y chromosome that contain deleterious mutations. Identifying the relative rates of cis-regulatory sequence evolution across Y chromosomes has been challenging due to the limited number of reference assemblies. The threespine stickleback (Gasterosteus aculeatus) Y chromosome is an excellent model to identify how regulatory mutations accumulate on Y chromosomes due to its intermediate state of divergence from the X chromosome. A large number of Y-linked gametologs still exist across 3 differently aged evolutionary strata to test these hypotheses. We found that putative enhancer regions on the Y chromosome exhibited elevated substitution rates and decreased polymorphism when compared to nonfunctional sites, like intergenic regions and synonymous sites. This suggests that many cis-regulatory regions are under positive selection on the Y chromosome. This divergence was correlated with X-biased gametolog expression, indicating the loss of expression from the Y chromosome may be favored by selection. Our findings provide evidence that Y-linked cis-regulatory regions exhibit signs of positive selection quickly after the suppression of recombination and allow comparisons with recent theoretical models that suggest the rapid divergence of regulatory regions may be favored to mask deleterious mutations on the Y chromosome.

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.

Comparative Genomics of Sex, Chromosomes, and Sex Chromosomes in Caenorhabditis elegans and Other Nematodes

The nematode phylum has evolved a remarkable diversity of reproductive modes, including the repeated emergence of asexuality and hermaphroditism across divergent clades. The species-richness and small genome size of nematodes make them ideal systems for investigating the genome-wide causes and consequences of such major transitions. The availability of functional annotations for most Caenorhabditis elegans genes further allows the linking of patterns of gene content evolution with biological processes. Such gene-centric studies were recently complemented by investigations of chromosome evolution that made use of the first chromosome-scale genome assemblies outside the Caenorhabditis genus. This review highlights recent comparative genomic studies of reproductive mode evolution addressing the hybrid origin of asexuality and the parallel gene loss following the emergence of hermaphroditism. It further summarizes ongoing efforts to characterize ancient linkage blocks called Nigon elements, which form central units of chromosome evolution. Fusions between Nigon elements have been demonstrated to impact recombination and speciation. Finally, multiple recent fusions between autosomal and the sex-linked Nigon element reveal insights into the dynamic evolution of sex chromosomes across various timescales.

Divergent outcomes of genetic conflict on the UV sex chromosomes of Marchantia polymorpha and Ceratodon purpureus

In species with separate sexes, the genome must produce two distinct developmental programs. Sexually dimorphic development may be controlled by either sex-limited loci or biased expression of loci transmitted through both sexes. Variation in the gene content of sex-limited chromosomes demonstrates that eukaryotic species differ markedly in the roles of these two mechanisms in governing sexual dimorphism. The bryophyte model systems Marchantia polymorpha and Ceratodon purpureus provide a particularly striking contrast. Although both species possess a haploid UV sex chromosome system, in which females carry a U chromosome and males carry a V, M. polymorpha relies on biased autosomal expression, while in C. purpureus, sex-linked genes drive dimorphism. Framing these genetic architectures as divergent outcomes of genetic conflict highlights comparative genomic analyses to better understand the evolution of sexual dimorphism.

Sex-biased gene expression across mammalian organ development and evolution

Sexually dimorphic traits are common among mammals and are specified during development through the deployment of sex-specific genetic programs. Because little is known about these programs, we investigated them using a resource of gene expression profiles in males and females throughout the development of five organs in five mammals (human, mouse, rat, rabbit, and opossum) and a bird (chicken). We found that sex-biased gene expression varied considerably across organs and species and was often cell-type specific. Sex differences increased abruptly around sexual maturity instead of increasing gradually during organ development. Finally, sex-biased gene expression evolved rapidly at the gene level, with differences between organs in the evolutionary mechanisms used, but more slowly at the cellular level, with the same cell types being sexually dimorphic across species.

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.

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.

Domestication affects sex-biased gene expression evolution in the duck

Genes with sex-biased expression are thought to underlie sexually dimorphic phenotypes and are therefore subject to different selection pressures in males and females. Many authors have proposed that sexual conflict leads to the evolution of sex-biased expression, which allows males and females to reach separate phenotypic and fitness optima. The selection pressures associated with domestication may cause changes in population architectures and mating systems, which in turn can alter their direction and strength. We compared sex-biased expression and genetic signatures in wild and domestic ducks (Anas platyrhynchos), and observed changes of sexual selection and identified the genomic divergence affected by selection forces. The extent of sex-biased expression in both sexes is positively correlated with the level of both d _N /d _S and nucleotide diversity. This observed changing pattern may mainly be owing to relaxed genetic constraints. We also demonstrate a clear link between domestication and sex-biased evolutionary rate in a comparative framework. Decreased polymorphism and evolutionary rate in domesticated populations generally matched life-history phenotypes known to experience artificial selection. Taken together, our work suggests the important implications of domestication in sex-biased evolution and the roles of artificial selection and sexual selection for shaping the diversity and evolutionary rate of the genome.

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.

Detection of positive selection acting on protein surfaces at the whole-genome scale in the human malaria parasite Plasmodium falciparum

The host-parasite evolutionary arms race is a fundamental process with medical implications. During this process, the host develops parasite resistance, and the parasite develops host immune evasion strategies. Thus, this process accelerates relevant protein evolution. This study test hypothesizes that proteins subject to sequence evolution structural constraints play a crucial role and that these constraints hinder the modification of such proteins in this process. These hypotheses were tested using Plasmodium falciparum model and evaluated protein structures predicted for the entire proteome by the AlphaFold method. Based on dN/dS test results and P. falciparum and P. reichenowi comparisons, the presented approach identified proteins subject to purifying selection acting on the whole sequence and buried residues (dN < dS) and positive selection on nonburied residues. Of the 26 proteins, some known antigens (ring-exported protein 3, RAP protein, erythrocyte binding antigen-140, and protein P47) targeted by the host immune system are promising vaccine candidates. The set also contained 11 enzymes, including FIKK kinase, which modifies host proteins. This set was compared with genes for which the dN/dS test suggested that positive selection acts on the whole gene (i.e., dN > dS). The present study found that such genes encode enzymes and antigenic vaccine candidates less frequently than genes for which evolution is not subject to selection constraints and positive selection acts on only exposed residues. The analysis was repeated comparing P. falciparum with P. alderi, which is more distantly related. The study discusses the potential implications of the presented methodology for rational vaccine design and the parasitology and evolutionary biology fields.

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.

Genomic architecture and sexually dimorphic expression underlying immunity in the red mason bee, Osmia bicornis

Insect pollinators provide crucial ecosystem services yet face increasing environmental pressures. The challenges posed by novel and reemerging pathogens on bee health means we need to improve our understanding of the immune system, an important barrier to infections and disease. Despite the importance of solitary bees, which are ecologically relevant, our understanding of the genomic basis and molecular mechanisms underlying their immune potential, and how intrinsic and extrinsic factors may influence it is limited. To improve our understanding of the genomic architecture underlying immunity of a key solitary bee pollinator, we characterized putative immune genes of the red mason bee, Osmia bicornis. In addition, we used publicly available RNA-seq datasets to determine how sexes differ in immune gene expression and splicing but also how pesticide exposure may affect immune gene expression in females. Through comparative genomics, we reveal an evolutionarily conserved set of more than 500 putative immune-related genes. We found genome-wide patterns of sex-biased gene expression, with greater enrichment of immune-related processes among genes with higher constitutive expression in males than females. Our results also suggest an up-regulation of immune-related genes in response to exposure to two common neonicotinoids, thiacloprid and imidacloprid. Collectively, our study provides important insights into the gene repertoire, regulation and expression differences in the sexes of O. bicornis, as well as providing additional support for how neonicotinoids can affect immune gene expression, which may affect the capacity of solitary bees to respond to pathogenic threats.

Effects of biological sex and pregnancy in experimental autoimmune encephalomyelitis: It's complicated

Experimental autoimmune encephalomyelitis (EAE) can be induced in many animal strains by inoculation with central nervous system antigens and adjuvant or by the passive transfer of lymphocytes reactive with these antigens and is widely used as an animal model for multiple sclerosis (MS). There are reports that female sex and pregnancy affect EAE. Here we review the effects of biological sex and the effects of pregnancy on the clinical features (including disease susceptibility) and pathophysiology of EAE. We also review reports of the possible mechanisms underlying these differences. These include sex-related differences in the immune system and in the central nervous system, the effects of hormones and the sex chromosomes and molecules unique to pregnancy. We also review sex differences in the response to factors that can modify the course of EAE. Our conclusion is that the effects of biological sex in EAE vary amongst animal models and should not be widely extrapolated. In EAE, it is therefore essential that studies looking at the effects of biological sex or pregnancy give full information about the model that is used (i.e. animal strain, sex, the inducing antigen, timing of EAE induction in relation to pregnancy, etc.). In addition, it would be preferable if more than one EAE model were used, to show if any observed effects are generalizable. This is clearly a field that requires further work. However, understanding of the mechanisms of sex differences could lead to greater understanding of EAE, and suggest possible therapies for MS.

Sex differences in interindividual gene expression variability across human tissues

Understanding phenotypic sex differences has long been a goal of biology from both a medical and evolutionary perspective. Although much attention has been paid to mean differences in phenotype between the sexes, little is known about sex differences in phenotypic variability. To gain insight into sex differences in interindividual variability at the molecular level, we analyzed RNA-seq data from 43 tissues from the Genotype-Tissue Expression project (GTEx). Within each tissue, we identified genes that show sex differences in gene expression variability. We found that these sex-differentially variable (SDV) genes are associated with various important biological functions, including sex hormone response, immune response, and other signaling pathways. By analyzing single-cell RNA sequencing data collected from breast epithelial cells, we found that genes with sex differences in gene expression variability in breast tissue tend to be expressed in a cell-type-specific manner. We looked for an association between SDV expression and Graves' disease, a well-known heavily female-biased disease, and found a significant enrichment of Graves' associated genes among genes with higher variability in females in thyroid tissue. This suggests a possible role for SDV expression in sex-biased disease. We then examined the evolutionary constraints acting on genes with sex differences in variability and found that they exhibit evidence of increased selective constraint. Through analysis of sex-biased eQTL data, we found evidence that SDV expression may have a genetic basis. Finally, we propose a simple evolutionary model for the emergence of SDV expression from sex-specific constraints.

Correlation scan: identifying genomic regions that affect genetic correlations applied to fertility traits

Although the genetic correlations between complex traits have been estimated for more than a century, only recently we have started to map and understand the precise localization of the genomic region(s) that underpin these correlations. Reproductive traits are often genetically correlated. Yet, we don't fully understand the complexities, synergism, or trade-offs between male and female fertility. In this study, we used reproductive traits in two cattle populations (Brahman; BB, Tropical Composite; TC) to develop a novel framework termed correlation scan (CS). This framework was used to identify local regions associated with the genetic correlations between male and female fertility traits. Animals were genotyped with bovine high-density single nucleotide polymorphisms (SNPs) chip assay. The data used consisted of ~1000 individual records measured through frequent ovarian scanning for age at first corpus luteum (AGECL) and a laboratory assay for serum levels of insulin growth hormone (IGF1 measured in bulls, IGF1b, or cows, IGF1c). The methodology developed herein used correlations of 500-SNP effects in a 100-SNPs sliding window in each chromosome to identify local genomic regions that either drive or antagonize the genetic correlations between traits. We used Fisher's Z-statistics through a permutation method to confirm which regions of the genome harboured significant correlations. About 30% of the total genomic regions were identified as driving and antagonizing genetic correlations between male and female fertility traits in the two populations. These regions confirmed the polygenic nature of the traits being studied and pointed to genes of interest. For BB, the most important chromosome in terms of local regions is often located on bovine chromosome (BTA) 14. However, the important regions are spread across few different BTA's in TC. Quantitative trait loci (QTLs) and functional enrichment analysis revealed many significant windows co-localized with known QTLs related to milk production and fertility traits, especially puberty. In general, the enriched reproductive QTLs driving the genetic correlations between male and female fertility are the same for both cattle populations, while the antagonizing regions were population specific. Moreover, most of the antagonizing regions were mapped to chromosome X. These results suggest regions of chromosome X for further investigation into the trade-offs between male and female fertility. We compared the CS with two other recently proposed methods that map local genomic correlations. Some genomic regions were significant across methods. Yet, many significant regions identified with the CS were overlooked by other methods.

On how to identify a seminal fluid protein: A commentary on Hurtado et al

Seminal fluid proteins (Sfps) have striking effects on the behaviour and physiology of females in many insects. Some Drosophila melanogaster Sfps are not highly or exclusively expressed in the accessory glands, but derive from, or are additionally expressed in other male reproductive tissues. The full suite of Sfps includes transferred proteins from all male reproductive tissues, regardless of expression level or presence of a signal peptide.

Potential Involvement of ewsr1-w Gene in Ovarian Development of Chinese Tongue Sole, Cynoglossus semilaevis

Simple Summary

Sexual dimorphism is a phenomenon commonly existing in animals. Chinese tongue sole Cynoglossus semilaevis is an economical marine fish with obvious female-biased size dimorphism. So, it is important to explore the molecular mechanism beyond gonadal development for sex control in aquaculture industry. RNA-binding protein Ewing Sarcoma protein-like (ewsr1) gene is important for mouse gonadal development and reproduction, however there are limited studies on this gene in teleost. In this study, two ewsr1 genes were cloned and characterized from C. semilaevis. The ewsr1-w gene, located in W chromosomes, showed female-biased expression during C. semilaevis gonadal development. In addition, knock-down effect and transcriptional regulation of Cs-ewsr1-w further suggested its essential role in ovarian development. This study broadened our understanding on ewsr1 function in teleost, and provided genetic resources for the further development of sex control breeding techniques in C. semilaevis aquaculture.

Abstract

Ewsr1 encodes a protein that acts as a multifunctional molecule in a variety of cellular processes. The full-length of Cs-ewsr1-w and Cs-ewsr1-z were cloned in Chinese tongue sole (Cynoglossus semilaevis). The open reading frame (ORF) of Cs-ewsr1-w was 1,767 bp that encoded 589 amino acids, while Cs-ewsr1-z was 1,794 bp that encoded 598 amino acids. Real-time PCR assays showed that Cs-ewsr1-w exhibited significant female-biased expression and could be hardly detected in male. It has the most abundant expression in ovaries among eight healthy tissues. Its expression in ovary increased gradually from 90 d to 3 y with C. semilaevis ovarian development and reached the peak at 3 y. After Cs-ewsr1-w knockdown with siRNA interference, several genes related to gonadal development including foxl2, sox9b and pou5f1 were down-regulated in ovarian cell line, suggesting the possible participation of Cs-ewsr1-w in C. semilaevis ovarian development. The dual-luciferase reporter assay revealed that the -733/-154 bp Cs-ewsr1-w promoter fragment exhibited strong transcription activity human embryonic kidney (HEK) 293T cell line. The mutation of a MAF BZIP Transcription Factor K (Mafk) binding site located in this fragment suggested that transcription factor Mafk might play an important role in Cs-ewsr1-w basal transcription. Our results will provide clues on the gene expression level, transcriptional regulation and knock-down effect of ewsr1 gene during ovarian development in teleost.

X chromosomes show relaxed selection and complete somatic dosage compensation across Timema stick insect species

Sex chromosomes have evolved repeatedly across the tree of life. As they are present in different copy numbers in males and females, they are expected to experience different selection pressures than the autosomes, with consequences including a faster rate of evolution, increased accumulation of sexually antagonistic alleles and the evolution of dosage compensation. Whether these consequences are general or linked to idiosyncrasies of specific taxa is not clear as relatively few taxa have been studied thus far. Here, we use whole-genome sequencing to identify and characterize the evolution of the X chromosome in five species of Timema stick insects with XX:X0 sex determination. The X chromosome had a similar size (approximately 12% of the genome) and gene content across all five species, suggesting that the X chromosome originated prior to the diversification of the genus. Genes on the X showed evidence of relaxed selection (elevated dN/dS) and a slower evolutionary rate (dN + dS) than genes on the autosomes, likely due to sex-biased mutation rates. Genes on the X also showed almost complete dosage compensation in somatic tissues (heads and legs), but dosage compensation was absent in the reproductive tracts. Contrary to prediction, sex-biased genes showed little enrichment on the X, suggesting that the advantage X-linkage provides to the accumulation of sexually antagonistic alleles is weak. Overall, we found the consequences of X-linkage on gene sequences and expression to be similar across Timema species, showing the characteristics of the X chromosome are surprisingly consistent over 30 million years of evolution.

Sex chromosomes as supergenes of speciation: why amphibians defy the rules?

As reflected by the two rules of speciation (Haldane's rule and the large X-/Z-effect), sex chromosomes are expected to behave like supergenes of speciation: they recombine only in one sex (XX females or ZZ males), supposedly recruit sexually antagonistic genes and evolve faster than autosomes, which can all contribute to pre-zygotic and post-zygotic isolation. While this has been mainly studied in organisms with conserved sex-determining systems and highly differentiated (heteromorphic) sex chromosomes like mammals, birds and some insects, these expectations are less clear in organismal groups where sex chromosomes repeatedly change and remain mostly homomorphic, like amphibians. In this article, we review the proposed roles of sex-linked genes in isolating nascent lineages throughout the speciation continuum and discuss their support in amphibians given current knowledge of sex chromosome evolution and speciation modes. Given their frequent recombination and lack of differentiation, we argue that amphibian sex chromosomes are not expected to become supergenes of speciation, which is reflected by the rarity of empirical studies consistent with a 'large sex chromosome effect' in frogs and toads. The diversity of sex chromosome systems in amphibians has a high potential to disentangle the evolutionary mechanisms responsible for the emergence of sex-linked speciation genes in other organisms. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

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.

The rattlesnake W chromosome: A GC-rich retroelement refugium with retained gene function across ancient evolutionary strata

Sex chromosomes diverge after the establishment of recombination suppression, resulting in differential sex-linkage of genes involved in genetic sex determination and dimorphic traits. This process produces systems of male or female heterogamety wherein the Y and W chromosomes are only present in one sex and are often highly degenerated. Sex-limited Y and W chromosomes contain valuable information about the evolutionary transition from autosomes to sex chromosomes, yet detailed characterizations of the structure, composition, and gene content of sex-limited chromosomes are lacking for many species. In this study, we characterize the female-specific W chromosome of the prairie rattlesnake (Crotalus viridis) and evaluate how recombination suppression and other processes have shaped sex chromosome evolution in ZW snakes. Our analyses indicate that the rattlesnake W chromosome is over 80% repetitive and that an abundance of GC-rich mdg4 elements has driven an overall high degree of GC-richness despite a lack of recombination. The W chromosome is also highly enriched for repeat sequences derived from endogenous retroviruses and likely acts as a “refugium” for these and other retroelements. We annotated 219 putatively functional W-linked genes across at least two evolutionary strata identified based on estimates of sequence divergence between Z and W gametologs. The youngest of these strata is relatively gene-rich, however gene expression across strata suggests retained gene function amidst a greater degree of degeneration following ancient recombination suppression. Functional annotation of W-linked genes indicates a specialization of the W chromosome for reproductive and developmental function since recombination suppression from the Z chromosome.

Sex-limited chromosomes and non-reproductive traits

Sex chromosomes are typically viewed as having originated from a pair of autosomes, and differentiated as the sex-limited chromosome (e.g. Y) has degenerated by losing most genes through cessation of recombination. While often thought that degenerated sex-limited chromosomes primarily affect traits involved in sex determination and sex cell production, accumulating evidence suggests they also influence traits not sex-limited or directly involved in reproduction. Here, we provide an overview of the effects of sex-limited chromosomes on non-reproductive traits in XY, ZW or UV sex determination systems, and discuss evolutionary processes maintaining variation at sex-limited chromosomes and molecular mechanisms affecting non-reproductive traits.

Detecting signatures of selection on gene expression

A substantial amount of phenotypic diversity results from changes in gene expression levels and patterns. Understanding how the transcriptome evolves is therefore a key priority in identifying mechanisms of adaptive change. However, in contrast to powerful models of sequence evolution, we lack a consensus model of gene expression evolution. Furthermore, recent work has shown that many of the comparative approaches used to study gene expression are subject to biases that can lead to false signatures of selection. Here we first outline the main approaches for describing expression evolution and their inherent biases. Next, we bridge the gap between the fields of phylogenetic comparative methods and transcriptomics to reinforce the main pitfalls of inferring selection on expression patterns and use simulation studies to show that shifts in tissue composition can heavily bias inferences of selection. We close by highlighting the multi-dimensional nature of transcriptional variation and identifying major unanswered questions in disentangling how selection acts on the transcriptome.

Evolution of sexual systems, sex chromosomes and sex-linked gene transcription in flatworms and roundworms

Many species with separate male and female individuals (termed ‘gonochorism’ in animals) have sex-linked genome regions. Here, we investigate evolutionary changes when genome regions become completely sex-linked, by analyses of multiple species of flatworms (Platyhelminthes; among which schistosomes recently evolved gonochorism from ancestral hermaphroditism), and roundworms (Nematoda) which have undergone independent translocations of different autosomes. Although neither the evolution of gonochorism nor translocations fusing ancestrally autosomal regions to sex chromosomes causes inevitable loss of recombination, we document that formerly recombining regions show genomic signatures of recombination suppression in both taxa, and become strongly genetically degenerated, with a loss of most genes. Comparisons with hermaphroditic flatworm transcriptomes show masculinisation and some defeminisation in schistosome gonad gene expression. We also find evidence that evolution of sex-linkage in nematodes is accompanied by transcriptional changes and dosage compensation. Our analyses also identify sex-linked genes that could assist future research aimed at controlling some of these important parasites.

Transitions between hermaphroditic and separate sexes are relatively understudied in animals compared to pants. Here, Wang et al. reconstruct the evolution of separate sexes in the flatworms and complex changes of sex chromosomes in the roundworms.

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.

Genome-wide association study reveals the genetic basis of growth trait in yellow catfish with sexual size dimorphism

Sexual size dimorphism has been widely observed in a large number of animals including fish species. Genome-wide association study (GWAS) is a powerful tool to dissect the genetic basis of complex traits, whereas the sex-differences in the genomics of animal complex traits have been ignored in the GWAS analysis. Yellow catfish (Pelteobagrus fulvidraco) is an important aquaculture fish in China with significant sexual size dimorphism. In this study, GWAS was conducted to identify candidate SNPs and genes related to body length (BL) and body weight (BW) in 125 female yellow catfish from a breeding population. In total, one BL-related SNP and three BW-related SNPs were identified to be significantly associated with the traits. Besides, one of these SNPs (Chr15:19195072) was shared in both the BW and BL traits in female yellow catfish, which was further validated in 185 male individuals and located on the exon of stat5b gene. Transgenic yellow catfish and zebrafish that expressed yellow catfish stat5b showed increased growth rate and reduction of sexual size dimorphism. These results not only reveal the genetic basis of growth trait and sexual size dimorphism in fish species, but also provide useful information for the marker-assisted breeding in yellow catfish.

Selection drives convergent gene expression changes during transitions to co-sexuality in haploid sexual systems

Co-sexuality has evolved repeatedly from unisexual (dioicous) ancestors across a wide range of taxa. However, the molecular changes underpinning this important transition remain unknown, particularly in organisms with haploid sexual systems such as bryophytes, red algae and brown algae. Here we explore four independent events of emergence of co-sexuality from unisexual ancestors in brown algal clades to examine the nature, evolution and degree of convergence of gene expression changes that accompany the breakdown of dioicy. The amounts of male versus female phenotypic differences in dioicous species were not correlated with the extent of sex-biased gene expression, in stark contrast to what is observed in animals. Although sex-biased genes exhibited a high turnover rate during brown alga diversification, some of their predicted functions were conserved across species. Transitions to co-sexuality consistently involved adaptive gene expression shifts and rapid sequence evolution, particularly for male-biased genes. Gene expression in co-sexual species was more similar to that in females rather than males of related dioicous species, suggesting that co-sexuality may have arisen from ancestral females. Finally, extensive convergent gene expression changes, driven by selection, were associated with the transition to co-sexuality. Together, our observations provide insights on how co-sexual systems arise from ancestral, haploid UV sexual systems.

Molecular Evolution across Mouse Spermatogenesis

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

Detection of sexually antagonistic transmission distortions in trio datasets

Sexual dimorphisms are widespread in animals and plants, for morphological as well as physiological traits. Understanding the genetic basis of sexual dimorphism and its evolution is crucial for understanding biological differences between the sexes. Genetic variants with sex‐antagonistic effects on fitness are expected to segregate in populations at the early phases of sexual dimorphism emergence. Detecting such variants is notoriously difficult, and the few genome‐scan methods employed so far have limited power and little specificity. Here, we propose a new framework to detect a signature of sexually antagonistic (SA) selection. We rely on trio datasets where sex‐biased transmission distortions can be directly tracked from parents to offspring, and identify signals of SA transmission distortions in genomic regions. We report the genomic location of six candidate regions detected in human populations as potentially under sexually antagonist selection. We find an enrichment of genes associated with embryonic development within these regions. Last, we highlight two candidate regions for SA selection in humans.

Sex differences in immune gene expression in the brain of a small shorebird

Males and females often exhibit differences in behaviour, life histories, and ecology, many of which are typically reflected in their brains. Neuronal protection and maintenance include complex processes led by the microglia, which also interacts with metabolites such as hormones or immune components. Despite increasing interest in sex-specific brain function in laboratory animals, the significance of sex-specific immune activation in the brain of wild animals along with the variables that could affect it is widely lacking. Here, we use the Kentish plover (Charadrius alexandrinus) to study sex differences in expression of immune genes in the brain of adult males and females, in two wild populations breeding in contrasting habitats: a coastal sea-level population and a high-altitude inland population in China. Our analysis yielded 379 genes associated with immune function. We show a significant male-biased immune gene upregulation. Immune gene expression in the brain did not differ in upregulation between the coastal and inland populations. We discuss the role of dosage compensation in our findings and their evolutionary significance mediated by sex-specific survival and neuronal deterioration. Similar expression profiles in the coastal and inland populations suggest comparable genetic control by the microglia and possible similarities in pathogen pressures between habitats. We call for further studies on gene expression of males and females in wild population to understand the implications of immune function for life-histories and demography in natural systems.

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.

Tissue-Specific Expression Pattern in Ancherythroculter nigrocauda, a Sexually Size Dimorphic Fish

Certain members of the Actinopterygii class are known to exhibit sexual dimorphism (SD) that results in major phenotypic differences between male and female fishes of a species. One of the most common differences between the two sexes is in body weight, a factor with a high economic value in aquaculture. In this study, we used RNA sequencing (RNA-seq) to study the liver and brain transcriptomes of Ancherythroculter nigrocauda, a fish exhibiting SD. Females attain about fourfold body weight of males at sexual maturity. Sample clustering showed that both sexes were grouped well with their sex phenotypes. In addition, 2,395 and 457 differentially expressed genes (DEGs) were identified in the liver and brain tissues, respectively. The gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses predicted the association of PPAR signaling, cytochrome P450, and steroid hormone biosynthesis to the differences in sexual size. In addition, weighted gene co-expression network analyses (WGCNA) were conducted, and the green module was identified to be significantly correlated with sexual size dimorphism (SSD). Altogether, these results improve our understanding of the molecular mechanism underlying SSD in A. nigrocauda.

A Rapid Evolving microRNA Cluster Rewires Its Target Regulatory Networks in Drosophila

New miRNAs are evolutionarily important but their functional evolution remains unclear. Here we report that the evolution of a microRNA cluster, mir-972C rewires its downstream regulatory networks in Drosophila. Genomic analysis reveals that mir-972C originated in the common ancestor of Drosophila where it comprises six old miRNAs. It has subsequently recruited six new members in the melanogaster subgroup after evolving for at least 50 million years. Both the young and the old mir-972C members evolved rapidly in seed and non-seed regions. Combining target prediction and cell transfection experiments, we found that the seed and non-seed changes in individual mir-972C members cause extensive target divergence among D. melanogaster, D. simulans, and D. virilis, consistent with the functional evolution of mir-972C reported recently. Intriguingly, the target pool of the cluster as a whole remains relatively conserved. Our results suggest that clustering of young and old miRNAs broadens the target repertoires by acquiring new targets without losing many old ones. This may facilitate the establishment of new miRNAs in existing 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 sex-specific evolutionary interaction between ADCY9 and CETP

Pharmacogenomic studies have revealed associations between rs1967309 in the adenylyl cyclase type 9 (ADCY9) gene and clinical responses to the cholesteryl ester transfer protein (CETP) modulator dalcetrapib, however, the mechanism behind this interaction is still unknown. Here, we characterized selective signals at the locus associated with the pharmacogenomic response in human populations and we show that rs1967309 region exhibits signatures of positive selection in several human populations. Furthermore, we identified a variant in CETP, rs158477, which is in long-range linkage disequilibrium with rs1967309 in the Peruvian population. The signal is mainly seen in males, a sex-specific result that is replicated in the LIMAA cohort of over 3400 Peruvians. Analyses of RNA-seq data further suggest an epistatic interaction on CETP expression levels between the two SNPs in multiple tissues, which also differs between males and females. We also detected interaction effects of the two SNPs with sex on cardiovascular phenotypes in the UK Biobank, in line with the sex-specific genotype associations found in Peruvians at these loci. We propose that ADCY9 and CETP coevolved during recent human evolution due to sex-specific selection, which points toward a biological link between dalcetrapib’s pharmacogene ADCY9 and its therapeutic target CETP.

The obstetrical dilemma hypothesis: there's life in the old dog yet

The term 'obstetrical dilemma' was coined by Washburn in 1960 to describe the trade-off between selection for a larger birth canal, permitting successful passage of a big-brained human neonate, and the smaller pelvic dimensions required for bipedal locomotion. His suggested solution to these antagonistic pressures was to give birth prematurely, explaining the unusual degree of neurological and physical immaturity, or secondary altriciality, observed in human infants. This proposed trade-off has traditionally been offered as the predominant evolutionary explanation for why human childbirth is so challenging, and inherently risky, compared to that of other primates. This perceived difficulty is likely due to the tight fit of fetal to maternal pelvic dimensions along with the convoluted shape of the birth canal and a comparatively low degree of ligamentous flexibility. Although the ideas combined under the obstetrical dilemma hypothesis originated almost a century ago, they have received renewed attention and empirical scrutiny in the last decade, with some researchers advocating complete rejection of the hypothesis and its assumptions. However, the hypothesis is complex because it presently captures several, mutually non-exclusive ideas: (i) there is an evolutionary trade-off resulting from opposing selection pressures on the pelvis; (ii) selection favouring a narrow pelvis specifically derives from bipedalism; (iii) human neonates are secondarily altricial because they are born relatively immature to ensure that they fit through the maternal bony pelvis; (iv) as a corollary to the asymmetric selection pressure for a spacious birth canal in females, humans evolved pronounced sexual dimorphism of pelvic shape. Recently, the hypothesis has been challenged on both empirical and theoretical grounds. Here, we appraise the original ideas captured under the 'obstetrical dilemma' and their subsequent evolution. We also evaluate complementary and alternative explanations for a tight fetopelvic fit and obstructed labour, including ecological factors related to nutrition and thermoregulation, constraints imposed by the stability of the pelvic floor or by maternal and fetal metabolism, the energetics of bipedalism, and variability in pelvic shape. This reveals that human childbirth is affected by a complex combination of evolutionary, ecological, and biocultural factors, which variably constrain maternal pelvic form and fetal growth. Our review demonstrates that it is unwarranted to reject the obstetrical dilemma hypothesis entirely because several of its fundamental assumptions have not been successfully discounted despite claims to the contrary. As such, the obstetrical dilemma remains a tenable hypothesis that can be used productively to guide evolutionary research.

Transitions to asexuality and evolution of gene expression in Artemia brine shrimp

While sexual reproduction is widespread among many taxa, asexual lineages have repeatedly evolved from sexual ancestors. Despite extensive research on the evolution of sex, it is still unclear whether this switch represents a major transition requiring major molecular reorganization, and how convergent the changes involved are. In this study, we investigated the phylogenetic relationship and patterns of gene expression of sexual and asexual lineages of Eurasian Artemia brine shrimp, to assess how gene expression patterns are affected by the transition to asexuality. We find only a few genes that are consistently associated with the evolution of asexuality, suggesting that this shift may not require an extensive overhauling of the meiotic machinery. While genes with sex-biased expression have high rates of expression divergence within Eurasian Artemia, neither female- nor male-biased genes appear to show unusual evolutionary patterns after sexuality is lost, contrary to theoretical expectations.

Epistatic interactions between sex chromosomes and autosomes can affect the stability of sex determination systems

Sex determination (SD) is an essential and ancient developmental process, but the genetic systems that regulate this process are surprisingly variable. Why SD mechanisms vary so much is a longstanding question in evolutionary biology. SD genes are generally located on sex chromosomes which also carry genes that interact epistatically with autosomes to affect fitness. How this affects the evolutionary stability of SD mechanisms is still unknown. Here, we explore how epistatic interactions between a sexually antagonistic (SA) non‐SD gene, located on either an ancestral or novel sex chromosome, and an autosomal gene affect the conditions under which an evolutionary transition to a new SD system occurs. We find that when the SD gene is linked to an ancestral sex‐chromosomal gene which engages in epistatic interactions, epistasis enhances the stability of the sex chromosomes so that they are retained under conditions where transitions would otherwise occur. This occurs both when weaker fitness effects are associated with the ancestral sex chromosome pair or stronger fitness effects associated with a newly evolved SD gene. However, the probability that novel SD genes spread is unaffected if they arise near genes involved in epistasis. This discrepancy occurs because, on autosomes, SA allele frequencies are typically lower than on sex chromosomes. In our model, increased frequencies of these alleles contribute to a higher frequency of epistasis which may therefore more readily occur on sex chromosomes. Because sex chromosome–autosome interactions are abundant and can take several forms, they may play a large role in maintaining sex chromosomes.

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

Sex-Biased Gene Expression and Isoform Profile of Brine Shrimp Artemia franciscana by Transcriptome Analysis

Simple Summary

The brine shrimp Artemia is a promising model organism for ZW sex determination system, but the genes related to sex determination and differentiation of Artemia have not yet been examined in detail. In this study, the first isoform-level transcriptome sequencing was performed on female and male Artemia franciscana. By using PacBio Iso-Seq and Illumina RNA-Seq technologies, we found 39 candidate sex determination genes that showed sex-biased gene expression. The male-biased expressed genes included DMRT1 and Sad genes, which had three and seven isoforms, respectively. Among these, the Sad gene is an ecdysteroid biosynthetic pathway gene associated with arthropod molting and metamorphosis. We propose the importance and the necessity of further research on genes involved in ecdysteroid biosynthesis. These results will contribute to understand sex determination and differentiation of Artemia and other crustaceans having ZW systems.

Abstract

The brine shrimp Artemia has a ZW sex determination system with ZW chromosomes in females and ZZ chromosomes in males. Artemia has been considered a promising model organism for ZW sex-determining systems, but the genes involved in sex determination and differentiation of Artemia have not yet been identified. Here, we conducted transcriptome sequencing of female and male A. franciscana using PacBio Iso-Seq and Illumina RNA-Seq techniques to identify candidate sex determination genes. Among the 42,566 transcripts obtained from Iso-Seq, 23,514 were analyzed. Of these, 2065 (8.8%) were female specific, 2513 (10.7%) were male specific, and 18,936 (80.5%) were co-expressed in females and males. Based on GO enrichment analysis and expression values, we found 10 female-biased and 29 male-biased expressed genes, including DMRT1 and Sad genes showing male-biased expression. Our results showed that DMRT1 has three isoforms with five exons, while Sad has seven isoforms with 2–11 exons. The Sad gene is involved in ecdysteroid signaling related to molting and metamorphosis in arthropods. Further studies on ecdysteroid biosynthetic genes are needed to improve our understanding of Artemia sex determination. This study will provide a valuable resource for sex determination and differentiation studies on Artemia and other crustaceans with ZW systems.

Support for faster and more adaptive Z chromosome evolution in two divergent lepidopteran lineages

The rateof divergence for Z or X chromosomes is usually observed to be greater than autosomes, but the proposed evolutionary causes for this pattern vary, as do empirical results from diverse taxa. Even among moths and butterflies (Lepidoptera), which generally share a single-origin Z chromosome, the handful of available studies give mixed support for faster or more adaptive evolution of the Z chromosome, depending on the species assayed. Here, we examine the molecular evolution of Z chromosomes in two additional lepidopteran species: the Carolina sphinx moth and the monarch butterfly, the latter of which possesses a recent chromosomal fusion yielding a segment of newly Z-linked DNA. We find evidence for both faster and more adaptive Z chromosome evolution in both species, although this effect is strongest in the neo-Z portion of the monarch sex chromosome. The neo-Z is less male-biased than expected of a Z chromosome, and unbiased and female-biased genes drive the signal for adaptive evolution here. Together these results suggest that male-biased gene accumulation and haploid selection have opposing effects on long-term rates of adaptation and may help explain the discrepancies in previous findings as well as the repeated evolution of neo-sex chromosomes in Lepidoptera.

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.

piRNA and Transposon Dynamics in Drosophila: A Female Story

The germlines of metazoans contain transposable elements (TEs) causing genetic instability and affecting fitness. To protect the germline from TE activity, gonads of metazoans produce TE-derived PIWI-interacting RNAs (piRNAs) that silence TE expression. In Drosophila, our understanding of piRNA biogenesis is mainly based on studies of the Drosophila melanogaster female germline. However, it is not known whether piRNA functions are also important in the male germline or whether and how piRNAs are affected by the global genomic context. To address these questions, we compared genome sequences, transcriptomes, and small RNA libraries extracted from entire testes and ovaries of two sister species: D. melanogaster and Drosophila simulans. We found that most TE-derived piRNAs were produced in ovaries and that piRNA pathway genes were strongly overexpressed in ovaries compared with testes, indicating that the silencing of TEs by the piRNA pathway mainly took place in the female germline. To study the relationship between host piRNAs and TE landscape, we analyzed TE genomic features and how they correlate with piRNA production in the two species. In D. melanogaster, we found that TE-derived piRNAs target recently active TEs. In contrast, although Drosophila simulans TEs do not display any features of recent activity, the host still intensively produced silencing piRNAs targeting old TE relics. Together, our results show that the piRNA silencing response mainly takes place in Drosophila ovaries and indicate that the host piRNA response is implemented following a burst of TE activity and could persist long after the extinction of active TE families.

Invasion and rapid adaptation of guppies (Poecilia reticulata) across the Hawaiian Archipelago

How much does natural selection, as opposed to genetic drift, admixture, and gene flow, contribute to the evolution of invasive species following introduction to a new environment? Here we assess how evolution can shape biological invasions by examining population genomic variation in non-native guppies (Poecilia reticulata) introduced to the Hawaiian Islands approximately a century ago. By examining 18 invasive populations from four Hawaiian islands and four populations from the native range in northern South America, we reconstructed the history of introductions and evaluated population structure as well as the extent of ongoing gene flow across watersheds and among islands. Patterns of differentiation indicate that guppies have developed significant population structure, with little natural or human-mediated gene flow having occurred among populations following introduction. Demographic modeling and admixture graph analyses together suggest that guppies were initially introduced to O'ahu and Maui and then translocated to Hawai'i and Kaua'i. We detected evidence for only one introduction event from the native range, implying that any adaptive evolution in introduced populations likely utilized the genetic variation present in the founding population. Environmental association tests accounting for population structure identified loci exhibiting signatures of adaptive variation related to predators and landscape characteristics but not nutrient regimes. When paired with high estimates of effective population sizes and detectable population structure, the presence of environment-associated loci supports the role of natural selection in shaping contemporary evolution of Hawaiian guppy populations. Our findings indicate that local adaptation may engender invasion success, particularly in species with life histories that facilitate rapid evolution. Finally, evidence of low gene flow between populations suggests that removal could be an effective approach to control invasive guppies across the Hawaiian archipelago.

Gene-Level, but Not Chromosome-Wide, Divergence between a Very Young House Fly Proto-Y Chromosome and Its Homologous Proto-X Chromosome

X and Y chromosomes are usually derived from a pair of homologous autosomes, which then diverge from each other over time. Although Y-specific features have been characterized in sex chromosomes of various ages, the earliest stages of Y chromosome evolution remain elusive. In particular, we do not know whether early stages of Y chromosome evolution consist of changes to individual genes or happen via chromosome-scale divergence from the X. To address this question, we quantified divergence between young proto-X and proto-Y chromosomes in the house fly, Musca domestica. We compared proto-sex chromosome sequence and gene expression between genotypic (XY) and sex-reversed (XX) males. We find evidence for sequence divergence between genes on the proto-X and proto-Y, including five genes with mitochondrial functions. There is also an excess of genes with divergent expression between the proto-X and proto-Y, but the number of genes is small. This suggests that individual proto-Y genes, but not the entire proto-Y chromosome, have diverged from the proto-X. We identified one gene, encoding an axonemal dynein assembly factor (which functions in sperm motility), that has higher expression in XY males than XX males because of a disproportionate contribution of the proto-Y allele to gene expression. The upregulation of the proto-Y allele may be favored in males because of this gene's function in spermatogenesis. The evolutionary divergence between proto-X and proto-Y copies of this gene, as well as the mitochondrial genes, is consistent with selection in males affecting the evolution of individual genes during early Y chromosome evolution.