The ontogeny and evolution of sex-biased gene expression in Drosophila melanogaster

Sexual conflict drive in the rapid evolution of new gametogenesis genes

The evolutionary forces underlying the rapid evolution in sequences and functions of new genes remain a mystery. Adaptation by natural selection explains the evolution of some new genes. However, many new genes perform sex-biased functions that have rapidly evolved over short evolutionary time scales, suggesting that new gene evolution may often be driven by conflicting selective pressures on males and females. It is well established that such sexual conflict (SC) plays a central role in maintaining phenotypic and genetic variation within populations, but the role of SC in driving new gene evolution remains essentially unknown. This review explores the connections between SC and new gene evolution through discussions of the concept of SC, the phenotypic and genetic signatures of SC in evolving populations, and the molecular mechanisms by which SC could drive the evolution of new genes. We synthesize recent work in this area with a discussion of the case of Apollo and Artemis, two extremely young genes (<200,000 years) in Drosophila melanogaster, which offered the first empirical insights into the evolutionary process by which SC could drive the evolution of new genes. These new duplicate genes exhibit the hallmarks of sexually antagonistic selection: rapid DNA and protein sequence evolution, essential sex-specific functions in gametogenesis, and complementary sex-biased expression patterns. Importantly, Apollo is essential for male fitness but detrimental to female fitness, while Artemis is essential for female fitness but detrimental to male fitness. These sexually antagonistic fitness effects and complementary changes to expression, sequence, and function suggest that these duplicates were selected for mitigating SC, but that SC has not been fully resolved. Finally, we propose Sexual Conflict Drive as a self-driven model to interpret the rapid evolution of new genes, explain the potential for SC and sexually antagonistic selection to contribute to long-term evolution, and suggest its utility for understanding the rapid evolution of new genes in gametogenesis.

Vascularization underlies differences in sexually selected skin coloration in a wild primate

Male reproductive competition can select for condition-dependent, conspicuous traits that signal some aspect of fighting ability and facilitate assessment of potential rivals. However, the underlying mechanisms that link the signal to a male's current condition are difficult to investigate in wild populations, often requiring invasive experimental manipulation. Here, we use digital photographs and chest skin samples to investigate the mechanisms of a visual signal used in male competition in a wild primate, the red chest patch in geladas (Theropithecus gelada). We analysed photographs collected during natural (n = 144) and anaesthetized conditions (n = 38) to understand variability in male and female chest redness, and we used chest skin biopsies (n = 38) to explore sex differences in gene expression. Male and female geladas showed similar average redness, but males exhibited a wider within-individual range in redness under natural conditions. These sex differences were also reflected at the molecular level, with 10.5% of genes exhibiting significant sex differences in expression. Subadult males exhibited intermediate gene expression patterns between adult males and females, pointing to mechanisms underlying the development of the red chest patch. We found that genes more highly expressed in males were associated with blood vessel development and maintenance but not with androgen or oestrogen activity. Together, our results suggest male gelada redness variability is driven by increased blood vessel branching in the chest skin, providing a potential link between male chest redness and current condition as increased blood circulation to exposed skin could lead to heat loss in the cold, high-altitude environment of geladas.

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.

Recent progress in understanding the genomic architecture of sexual conflict

Genomic conflict between the sexes over shared traits is widely assumed to be resolved through the evolution of sex-biased expression and the subsequent emergence of sexually dimorphic phenotypes. However, while there is support for a broad relationship between genome-wide patterns of expression level and sexual conflict, recent studies suggest that sex differences in the nature and strength of interactions between loci are instead key to conflict resolution. Furthermore, the advent of new technologies for measuring and perturbing expression means we now have much more power to detect genomic signatures of sexual conflict. Here, we review our current understanding of the genomic architecture of sexual conflict in the light of these new studies and highlight the potential for novel approaches to address outstanding knowledge gaps.

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.

Sex-biased expression is associated with chromatin state in D. melanogaster and D. simulans

We propose a new model for the association of chromatin state and sex-bias in expression. We hypothesize enrichment of open chromatin in the sex where we see expression bias (OS) and closed chromatin in the opposite sex (CO). In this study of D. melanogaster and D. simulans head tissue, sex-bias in expression is associated with H3K4me3 (open mark) in males for male-biased genes and in females for female-biased genes in both species. Sex-bias in expression is also largely conserved in direction and magnitude between the two species on the X and autosomes. In male-biased orthologs, the sex-bias ratio is more divergent between species if both species have H3K27me2me3 marks in females compared to when either or neither species has H3K27me2me3 in females. H3K27me2me3 marks in females are associated with male-bias in expression on the autosomes in both species, but on the X only in D. melanogaster . In female-biased orthologs the relationship between the species for the sex-bias ratio is similar regardless of the H3K27me2me3 marks in males. Female-biased orthologs are more similar in the ratio of sex-bias than male-biased orthologs and there is an excess of male-bias in expression in orthologs that gain/lose sex-bias. There is an excess of male-bias in sex-limited expression in both species suggesting excess male-bias is due to rapid evolution between the species. The X chromosome has an enrichment in male-limited H3K4me3 in both species and an enrichment of sex-bias in expression compared to the autosomes.

Ancient homomorphy of molluscan sex chromosomes sustained by reversible sex-biased genes and sex determiner translocation

Contrary to classic theory prediction, sex-chromosome homomorphy is prevalent in the animal kingdom but it is unclear how ancient homomorphic sex chromosomes avoid chromosome-scale degeneration. Molluscs constitute the second largest, Precambrian-originated animal phylum and have ancient, uncharacterized homomorphic sex chromosomes. Here, we profile eight genomes of the bivalve mollusc family of Pectinidae in a phylogenetic context and show 350 million years sex-chromosome homomorphy, which is the oldest known sex-chromosome homomorphy in the animal kingdom, far exceeding the ages of well-known heteromorphic sex chromosomes such as 130-200 million years in mammals, birds and flies. The long-term undifferentiation of molluscan sex chromosomes is potentially sustained by the unexpected intertwined regulation of reversible sex-biased genes, together with the lack of sexual dimorphism and occasional sex chromosome turnover. The pleiotropic constraint of regulation of reversible sex-biased genes is widely present in ancient homomorphic sex chromosomes and might be resolved in heteromorphic sex chromosomes through gene duplication followed by subfunctionalization. The evolutionary dynamics of sex chromosomes suggest a mechanism for 'inheritance' turnover of sex-determining genes that is mediated by translocation of a sex-determining enhancer. On the basis of these findings, we propose an evolutionary model for the long-term preservation of homomorphic sex chromosomes.

Dynamics of sex-biased gene expression during development in the stick insect Timema californicum

Sexually dimorphic phenotypes are thought to arise primarily from sex-biased gene expression during development. Major changes in developmental strategies, such as the shift from hemimetabolous to holometabolous development, are therefore expected to have profound consequences for the dynamics of sex-biased gene expression. However, no studies have previously examined sex-biased gene expression during development in hemimetabolous insects, precluding comparisons between developmental strategies. Here we characterized sex-biased gene expression at three developmental stages in a hemimetabolous stick insect (Timema californicum): hatchlings, juveniles, and adults. As expected, the proportion of sex-biased genes gradually increased during development, mirroring the gradual increase of phenotypic sexual dimorphism. Sex-biased genes identified at early developmental stages were generally consistently male- or female-biased at later stages, suggesting their importance in sexual differentiation. Additionally, we compared the dynamics of sex-biased gene expression during development in T. californicum to those of the holometabolous fly Drosophila melanogaster by reanalyzing publicly available RNA-seq data from third instar larval, pupal and adult stages. In D. melanogaster, 84% of genes were sex-biased at the adult stage (compared to only 20% in T. californicum), and sex-biased gene expression increased abruptly at the adult stage when morphological sexual dimorphism is manifested. Our findings are consistent with the prediction that the dynamics of sex-biased gene expression during development differ extensively between holometabolous and hemimetabolous insect species.

RNA-Seq Analysis Reveals Expression Regulatory Divergence of W-Linked Genes between Two Contrasting Chicken Breeds

Simple Summary

Understanding the mode of gene expression and regulation is essential for understanding the evolutionary process. Many previous studies tried to explain regulatory changes at the autosomal level, but little research has extended these explorations to the field of sex chromosomes due to their complex sex-limit features. Here, we first adopted an innovative method of identifying regulatory divergence of W-linked genes. Compared with cis-regulatory divergence, trans acting genes were more extensive in the W chromosome. We also found that divergent sex specific selection cannot strongly affect the expression evolution of the W chromosome. This insensitivity to selection may be one of the reasons why regulatory divergence is so small between autosomal and sex chromosomes.

Abstract

The regulation of gene expression is a complex process involving organism function and phenotypic diversity, and is caused by cis- and trans- regulation. While prior studies identified the regulatory pattern of the autosome rewiring in hybrids, the role of gene regulation in W sex chromosomes is not clear due to their degradation and sex-limit expression. Here, we developed reciprocal crosses of two chicken breeds, White Leghorn and Cornish Game, which exhibited broad differences in gender-related traits, and assessed the expression of the genes on the W chromosome to disentangle the contribution of cis- and trans-factors to expression divergence. We found that female-specific selection does not have a significant effect on W chromosome gene-expression patterns. For different tissues, there were most parental divergence expression genes in muscle, and also more heterosis compared with two other tissues. Notably, a broader pattern of trans regulation in the W chromosome was observed, which is consistent with autosomes. Taken together, this work describes the regulatory divergence of W-linked genes between two contrasting breeds and indicates sex chromosomes have a unique regulation and expression mechanism.

The developmentally dynamic microRNA transcriptome of Glossina pallidipes tsetse flies, vectors of animal trypanosomiasis

Summary

MicroRNAs (miRNAs) are single stranded gene regulators of 18-25 bp in length. They play a crucial role in regulating several biological processes in insects. However, the functions of miRNA in Glossina pallidipes, one of the biological vectors of African animal trypanosomosis in sub-Saharan Africa, remain poorly characterized. We used a combination of both molecular biology and bioinformatics techniques to identify miRNA genes at different developmental stages (larvae, pupae, teneral and reproductive unmated adults, gravid females) and sexes of G. pallidipes. We identified 157 mature miRNA genes, including 12 novel miRNAs unique to G. pallidipes. Moreover, we identified 93 miRNA genes that were differentially expressed by sex and/or in specific developmental stages. By combining both miRanda and RNAhybrid algorithms, we identified 5550 of their target genes. Further analyses with the Gene Ontology term and KEGG pathways for these predicted target genes suggested that the miRNAs may be involved in key developmental biological processes. Our results provide the first repository of G. pallidipes miRNAs across developmental stages, some of which appear to play crucial roles in tsetse fly development. Hence, our findings provide a better understanding of tsetse biology and a baseline for exploring miRNA genes in tsetse flies.

Availability and implementation

Raw sequence data are available from NCBI Sequence Read Archives (SRA) under Bioproject accession number PRJNA590626.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Sex differences in deleterious mutational effects in Drosophila melanogaster: combining quantitative and population genetic insights

Fitness effects of deleterious mutations can differ between females and males due to: (i) sex differences in the strength of purifying selection; and (ii) sex differences in ploidy. Although sex differences in fitness effects have important broader implications (e.g., for the evolution of sex and lifespan), few studies have quantified their scope. Those that have belong to one of two distinct empirical traditions: (i) quantitative genetics, which focusses on multi-locus genetic variances in each sex, but is largely agnostic about their genetic basis; and (ii) molecular population genetics, which focusses on comparing autosomal and X-linked polymorphism, but is poorly suited for inferring contemporary sex differences. Here, we combine both traditions to present a comprehensive analysis of female and male adult reproductive fitness among 202 outbred, laboratory-adapted, hemiclonal genomes of Drosophila melanogaster. While we find no clear evidence for sex differences in the strength of purifying selection, sex differences in ploidy generate multiple signals of enhanced purifying selection for X-linked loci. These signals are present in quantitative genetic metrics-i.e., a disproportionate contribution of the X to male (but not female) fitness variation-and population genetic metrics-i.e., steeper regressions of an allele's average fitness effect on its frequency, and proportionally less nonsynonymous polymorphism on the X than autosomes. Fitting our data to models for both sets of metrics, we infer that deleterious alleles are partially recessive. Given the often-large gap between quantitative and population genetic estimates of evolutionary parameters, our study showcases the benefits of combining genomic and fitness data when estimating such parameters.

Developmental mechanisms of sex differences: from cells to organisms

Male-female differences in many developmental mechanisms lead to the formation of two morphologically and physiologically distinct sexes. Although this is expected for traits with prominent differences between the sexes, such as the gonads, sex-specific processes also contribute to traits without obvious male-female differences, such as the intestine. Here, we review sex differences in developmental mechanisms that operate at several levels of biological complexity - molecular, cellular, organ and organismal - and discuss how these differences influence organ formation, function and whole-body physiology. Together, the examples we highlight show that one simple way to gain a more accurate and comprehensive understanding of animal development is to include both sexes.

Comparative Analysis of PacBio and Oxford Nanopore Sequencing Technologies for Transcriptomic Landscape Identification of Penaeus monodon

With the advantages that long-read sequencing platforms such as Pacific Biosciences (Menlo Park, CA, USA) (PacBio) and Oxford Nanopore Technologies (Oxford, UK) (ONT) can offer, various research fields such as genomics and transcriptomics can exploit their benefits. Selecting an appropriate sequencing platform is undoubtedly crucial for the success of the research outcome, thus there is a need to compare these long-read sequencing platforms and evaluate them for specific research questions. This study aims to compare the performance of PacBio and ONT platforms for transcriptomic analysis by utilizing transcriptome data from three different tissues (hepatopancreas, intestine, and gonads) of the juvenile black tiger shrimp, Penaeus monodon. We compared three important features: (i) main characteristics of the sequencing libraries and their alignment with the reference genome, (ii) transcript assembly features and isoform identification, and (iii) correlation of the quantification of gene expression levels for both platforms. Our analyses suggest that read-length bias and differences in sequencing throughput are highly influential factors when using long reads in transcriptome studies. These comparisons can provide a guideline when designing a transcriptome study utilizing these two long-read sequencing technologies.

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.

Differential gene expression patterns during gametophyte development provide insights into sex differentiation in the dioicous kelp Saccharina japonica

BACKGROUND

In brown algae, dioicy is the prevalent sexual system, and phenotypic differences between male and female gametophytes have been found in many dioicous species. Saccharina japonica show remarkable sexual dimorphism in gametophytes before gametogenesis. A higher level of phenotypic differentiation was also found in female and male gametes after gametogenesis. However, the patterns of differential gene expression throughout gametophyte development and how these changes might relate to sex-specific fitness at the gamete stage in S. japonica are not well known.

RESULTS

In this study, differences in gene expression between male and female gametophytes in different developmental stages were investigated using comparative transcriptome analysis. Among the 20,151 genes expressed in the haploid gametophyte generation, 37.53% were sex-biased. The abundance of sex-biased genes in mature gametophytes was much higher than that in immature gametophytes, and more male-biased than female-biased genes were observed in the mature stage. The predicted functions of most sex-biased genes were closely related to the sex-specific characteristics of gametes, including cell wall biosynthesis, sperm motility, and sperm and egg recognition. In addition, 51 genes were specifically expressed in males in both stages, showing great potential as candidate male sex-determining region (SDR) genes.

CONCLUSIONS

This study describes a thorough investigation into differential gene expression between male and female gametophytes in the dioicous kelp S. japonica. A large number of sex-biased genes in mature gametophytes may be associated with the divergence of phenotypic traits and physiological functions between female gametes (eggs) and male gametes (sperm) during sexual differentiation. These genes may mainly come from new sex-biased genes that have recently evolved in the S. japonica lineage. The duplication of sex-biased genes was detected, which may increase the number of sex-biased genes after gametogenesis in S. japonica to some extent. The excess of male-biased genes over female-biased genes in the mature stage may reflect the different levels of sexual selection across sexes. This study deepens our understanding of the regulation of sex development and differentiation in the dioicous kelp S. japonica.

Support for the adaptive decoupling hypothesis from whole-transcriptome profiles of a hypermetamorphic and sexually dimorphic insect, Neodiprion lecontei

Though seemingly bizarre, the dramatic morphological and ecological transformation that occurs when immature life stages metamorphose into reproductive adults is one of the most successful developmental strategies on the planet. The adaptive decoupling hypothesis (ADH) proposes that metamorphosis is an adaptation for breaking developmental links between traits expressed in different life stages, thereby facilitating their independent evolution when exposed to opposing selection pressures. Here, we draw inspiration from the ADH to develop a conceptual framework for understanding changes in gene expression across ontogeny. We hypothesized that patterns of stage-biased and sex-biased gene expression are the product of both decoupling mechanisms and selection history. To test this hypothesis, we characterized transcriptome-wide patterns of gene-expression traits for three ecologically distinct larval stages (all male) and adult males and females of a hypermetamorphic insect (Neodiprion lecontei). We found that stage-biased gene expression was most pronounced between larval and adult males, which is consistent with the ADH. However, even in the absence of a metamorphic transition, considerable stage-biased expression was observed among morphologically and behaviourally distinct larval stages. Stage-biased expression was also observed across ecologically relevant Gene Ontology categories and genes, highlighting the role of ecology in shaping patterns of gene expression. We also found that the magnitude and prevalence of stage-biased expression far exceeded adult sex-biased expression. Overall, our results highlight how the ADH can shed light on transcriptome-wide patterns of gene expression in organisms with complex life cycles. For maximal insight, detailed knowledge of organismal ecology is also essential.

In Silico Prediction of Transcription Factor Collaborations Underlying Phenotypic Sexual Dimorphism in Zebrafish (Danio rerio)

The transcriptional regulation of gene expression in higher organisms is essential for different cellular and biological processes. These processes are controlled by transcription factors and their combinatorial interplay, which are crucial for complex genetic programs and transcriptional machinery. The regulation of sex-biased gene expression plays a major role in phenotypic sexual dimorphism in many species, causing dimorphic gene expression patterns between two different sexes. The role of transcription factor (TF) in gene regulatory mechanisms so far has not been studied for sex determination and sex-associated colour patterning in zebrafish with respect to phenotypic sexual dimorphism. To address this open biological issue, we applied bioinformatics approaches for identifying the predicted TF pairs based on their binding sites for sex and colour genes in zebrafish. In this study, we identified 25 (e.g., STAT6-GATA4; JUN-GATA4; SOX9-JUN) and 14 (e.g., IRF-STAT6; SOX9-JUN; STAT6-GATA4) potentially cooperating TFs based on their binding patterns in promoter regions for sex determination and colour pattern genes in zebrafish, respectively. The comparison between identified TFs for sex and colour genes revealed several predicted TF pairs (e.g., STAT6-GATA4; JUN-SOX9) are common for both phenotypes, which may play a pivotal role in phenotypic sexual dimorphism in zebrafish.

Ecdysone Quantification from Whole Body Samples of Drosophila melanogaster Larvae

Steroid hormones strictly control the timing of sexual maturation and final body size both in vertebrates and invertebrates. In insects, the steroid hormone ecdysone controls the timing of the molts between larval instars as well as the transition to metamorphosis. Growth during the final instar accounts for over 80% of the increase in final mass in insects, and the duration of this growth period is driven by a sequence of small ecdysone pulses that ultimately induce metamorphosis. Historically the biologically active form of ecdysone, 20-hydroxyecdysone (20E), was quantified using radio-immunoassays, bioassays, or chromatography assays. However, these assays are methodologically complicated and often time consuming. Furthermore, collecting samples for precise measurements of ecdysone concentrations using these assays is limited in small insects like Drosophila melanogaster. Here, we describe an accurate and sensitive method to collect carefully-staged third instar larvae suitable for preparing samples for ecdysone quantification using a commercially-available 20E enzyme immunoassay (EIA). Because we resynchronize larval development at the molt to the final instar, collect large samples, and weigh each sample, we are able to detect a small ecdysone peak early in the final instar known as the critical weight ecdysone peak. This method detects peaks as low as 6 pg 20E/mg larval sample, allowing us to quantify other small ecdysone peaks in flies - the necessary prerequisite for eventually determining their regulation and function.

Signals of selection beyond bottlenecks between exotic populations of the bull-headed dung beetle, Onthophagus taurus

Colonization of new environments can lead to population bottlenecks and rapid phenotypic evolution that could be due to neutral and selective processes. Exotic populations of the bull-headed dung beetle (Onthophagus taurus) have differentiated in opposite directions from native beetles in male horn-to-body size allometry and female fecundity. Here we test for genetic and transcriptional differences among two exotic and one native O. taurus populations after three generations in common garden conditions. We sequenced RNA from 24 individuals for each of the three populations including both sexes, and spanning four developmental stages for the two exotic, differentiated populations. Identifying 270,400 high-quality single nucleotide polymorphisms, we revealed a strong signal of genetic differentiation between the three populations, and evidence of recent bottlenecks within and an excess of outlier loci between exotic populations. Differences in gene expression between populations were greatest in prepupae and early adult life stages, stages during which differences in male horn development and female fecundity manifest. Finally, genes differentially expressed between exotic populations also had greater genetic differentiation and performed functions related to chitin biosynthesis and nutrient sensing, possibly underlying allometry and fecundity trait divergences. Our results suggest that beyond bottlenecks, recent introductions have led to genetic and transcriptional differences in genes correlated with observed phenotypic differences.

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

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

Differences in the contributions of sex linkage and androgen regulation to sex-biased gene expression in juvenile and adult sticklebacks

Different evolutionary interests between males and females can lead to the evolution of sexual dimorphism. However, intersex genetic correlations due to the shared genome can constrain the evolution of sexual dimorphism, resulting in intra-locus sexual conflict. One of the mechanisms resolving this conflict is sex linkage, which allows males and females to carry different alleles on sex chromosomes. Another is a regulatory mutation causing sex-biased gene expression, which is often mediated by gonadal steroids in vertebrates. How do these two mechanisms differ in the contributions to the resolution of intra-locus sexual conflict? The magnitude of sexual conflict often varies between the juvenile and adult stages. Because gonadal steroids change in titre during development, we hypothesized that gonadal steroids play a role in sexual dimorphism expression only at certain developmental stages, whereas sex linkage is more important for sexual dimorphism expressed throughout life. Our brain transcriptome analysis of juvenile and adult threespine sticklebacks showed that the majority of genes that were sex-biased in both stages were sex-linked. The relative contribution of androgen-dependent regulation to the sex-biased transcriptome increased and that of sex linkage declined in adults compared to juveniles. The magnitude of the sex differences was greater in sex-linked genes than androgen-responsive genes, suggesting that sex linkage is more effective than androgen regulation in the production of large sex differences in gene expression. Overall, our data are consistent with the hypothesis that sex linkage is effective in resolving sexual conflict throughout life, whereas androgen-dependent regulation can contribute to temporary resolution of sexual conflict.

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

Background

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

Results

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

Conclusions

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

Differential Gene Expression between Fungal Mating Types Is Associated with Sequence Degeneration

Degenerative mutations in non-recombining regions, such as in sex chromosomes, may lead to differential expression between alleles if mutations occur stochastically in one or the other allele. Reduced allelic expression due to degeneration has indeed been suggested to occur in various sex-chromosome systems. However, whether an association occurs between specific signatures of degeneration and differential expression between alleles has not been extensively tested, and sexual antagonism can also cause differential expression on sex chromosomes. The anther-smut fungus Microbotryum lychnidis-dioicae is ideal for testing associations between specific degenerative signatures and differential expression because 1) there are multiple evolutionary strata on the mating-type chromosomes, reflecting successive recombination suppression linked to mating-type loci; 2) separate haploid cultures of opposite mating types help identify differential expression between alleles; and 3) there is no sexual antagonism as a confounding factor accounting for differential expression. We found that differentially expressed genes were enriched in the four oldest evolutionary strata compared with other genomic compartments, and that, within compartments, several signatures of sequence degeneration were greater for differentially expressed than non-differentially expressed genes. Two particular degenerative signatures were significantly associated with lower expression levels within differentially expressed allele pairs: upstream insertion of transposable elements and mutations truncating the protein length. Other degenerative mutations associated with differential expression included nonsynonymous substitutions and altered intron or GC content. The association between differential expression and allele degeneration is relevant for a broad range of taxa where mating compatibility or sex is determined by genes located in large regions where recombination is suppressed.

Release from intralocus sexual conflict? Evolved loss of a male sexual trait demasculinizes female gene expression

The loss of sexual ornaments is observed across taxa, and pleiotropic effects of such losses provide an opportunity to gain insight into underlying dynamics of sex-biased gene expression and intralocus sexual conflict (IASC). We investigated this in a Hawaiian field cricket, Teleogryllus oceanicus, in which an X-linked genotype ( flatwing) feminizes males' wings and eliminates their ability to produce sexually selected songs. We profiled adult gene expression across somatic and reproductive tissues of both sexes. Despite the feminizing effect of flatwing on male wings, we found no evidence of feminized gene expression in males. Instead, female transcriptomes were more strongly affected by flatwing than males', and exhibited demasculinized gene expression. These findings are consistent with a relaxation of IASC constraining female gene expression through loss of a male sexual ornament. In a follow-up experiment, we found reduced testes mass in flatwing males, whereas female carriers showed no reduction in egg production. By contrast, female carriers exhibited greater measures of body condition. Our results suggest sex-limited phenotypic expression offers only partial resolution to IASC, owing to pleiotropic effects of the loci involved. Benefits conferred by release from intralocus conflict could help explain widespread loss of sexual ornaments across taxa.

Absence of a Faster-X Effect in Beetles (Tribolium, Coleoptera)

The faster-X effect, namely the rapid evolution of protein-coding genes on the X chromosome, has been widely reported in metazoans. However, the prevalence of this phenomenon across diverse systems and its potential causes remain largely unresolved. Analysis of sex-biased genes may elucidate its possible mechanisms: for example, in systems with X/Y males a more pronounced faster-X effect in male-biased genes than in female-biased or unbiased genes may suggest fixation of recessive beneficial mutations rather than genetic drift. Further, theory predicts that the faster-X effect should be promoted by X chromosome dosage compensation. Here, we asked whether we could detect a faster-X effect in genes of the beetle Tribolium castaneum (and T. freemani orthologs), which has X/Y sex-determination and heterogametic males. Our comparison of protein sequence divergence (dN/dS) on the X chromosome vs. autosomes indicated a rarely observed absence of a faster-X effect in this organism. Further, analyses of sex-biased gene expression revealed that the X chromosome was particularly highly enriched for ovary-biased genes, which evolved slowly. In addition, an evaluation of male X chromosome dosage compensation in the gonads and in non-gonadal somatic tissues indicated a striking lack of compensation in the testis. This under-expression in testis may limit fixation of recessive beneficial X-linked mutations in genes transcribed in these male sex organs. Taken together, these beetles provide an example of the absence of a faster-X effect on protein evolution in a metazoan, that may result from two plausible factors, strong constraint on abundant X-linked ovary-biased genes and a lack of gonadal dosage compensation.

Sex biased expression and co-expression networks in development, using the hymenopteran Nasonia vitripennis

Sexual dimorphism requires regulation of gene expression in developing organisms. These developmental differences are caused by differential expression of genes and isoforms. The effect of expressing a gene is also influenced by which other genes are simultaneously expressed (functional interactions). However, few studies have described how these processes change across development. We compare the dynamics of differential expression, isoform switching and functional interactions in the sexual development of the model parasitoid wasp Nasonia vitripennis, a system that permits genome wide analysis of sex bias from early embryos to adults. We find relatively little sex-bias in embryos and larvae at the gene level, but several sub-networks show sex-biased functional interactions in early developmental stages. These networks provide new candidates for hymenopteran sex determination, including histone modification. In contrast, sex-bias in pupae and adults is driven by the differential expression of genes. We observe sex-biased isoform switching consistently across development, but mostly in genes that are already differentially expressed. Finally, we discover that sex-biased networks are enriched by genes specific to the Nasonia clade, and that those genes possess the topological properties of key regulators. These findings suggest that regulators in sex-biased networks evolve more rapidly than regulators of other developmental networks.

Comparative Genomics Reveals Accelerated Evolution of Fright Reaction Genes in Ostariophysan Fishes

The ostariophysian fishes are the most species-rich clade in freshwaters. This diversification has been suggested to be associated with the fright reaction presented in most ostariophysians. However, the genetic forces that underlie fright reaction remains poorly understood. In the present study, through integrating behavioral, physiological, transcriptomic, and evolutionary genomic analyses, we found that the fright reaction has a broad impact on zebrafish at multiple levels, including changes in swimming behaviors, cortisol levels, and gene expression patterns. In total, 1,555 and 1,599 differentially expressed genes were identified in olfactory mucosae and brain of zebrafish, respectively, with a greater number upregulated after the fright reaction. Functional annotation showed that response to stress and signal transduction were strongly represented, which is directly associated with the fright reaction. These differentially expressed genes were shown to be evolved accelerated under the influence of positive selection, indicating that protein-coding evolution has played a major role in fright reaction. We found the basal vomeronasal type 2 receptors (v2r) gene, v2rl1, displayed significantly decrease expression after fright reaction, which suggests that v2rs may be important to detect the alarm substance and induce the fright reaction. Collectively, based on our transcriptome and evolutionary genomics analyses, we suggest that transcriptional plasticity of gene may play an important role in fright reaction in ostariophysian fishes.

Conjoint Analysis of SMRT- and Illumina-Based RNA-Sequencing Data of Fenneropenaeus chinensis Provides Insight Into Sex-Biased Expression Genes Involved in Sexual Dimorphism

Fenneropenaeus chinensis (F. chinensis) is one of the most commercially important cultured shrimps in China. The adult F. chinensis exhibit sexual dimorphism in growth and body color. In this research, we profiled the whole transcriptome of F. chinensis by using single molecule real-time-based full-length transcriptome sequencing. We further performed Illumina-based short reads RNA-seq on muscle and gonad of two sexes to detect the sex-biased expression genes. In muscle, we observed significantly more female-biased transcripts. With the differentially expressed transcripts (DETs) in muscle, some pathways related to the energy metabolism were enriched, which may be responsible for the difference of growth. We also digged out a pathway named porphyrin and chlorophyll metabolism. It was speculated to relevant to the difference of body color between the two sexes of shrimp. Interestingly, almost all DETs in these pathways were female-biased expression in muscle, which could explain the phenomenon of better growth performance and darker body color in female. In gonad, several pathways involved in reproduction were enriched. For instance, some female-biased DETs participated in the arachidonic acid metabolism, which was reported crucial in female reproduction. In conclusion, our studies identified abundant sex-biased expression transcripts and important pathways involved in sexual dimorphism by using the RNA-seq method. It provided a basis for future researches on the sexual dimorphism of F. chinensis.

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

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

The developmental origins of sex-biased expression in cardiac development

BACKGROUND

Expression patterns between males and females vary in every adult tissue, even in organs with no conspicuous dimorphisms such as the heart. While studies of male and female differences have traditionally focused on the influence of sex hormones, these do not account for all the differences at the molecular and epigenetic levels. We previously reported that a substantial number of genes were differentially expressed in male and female mouse embryonic stem (ES) cells and revealed dose-dependent enhancer activity in response to Prdm14, a key pluripotency factor expressed more highly in female ES cells. In this work, we investigated the role of Prdm14 in establishing sex-specific gene expression networks. We surveyed the sex-specific landscape in early embryogenesis with special reference to cardiac development. We generated sex-specific co-expression networks from mouse ES cells, examined the presence of sex-specific chromatin domains, and analyzed previously published datasets from different developmental time points to characterize how sex-biased gene expression waxes and wanes to evaluate whether sex-biased networks are detectable throughout heart development.

RESULTS

We performed ChIP-seq on male and female mouse ES cells to determine differences in chromatin status. Our study reveals sex-biased histone modifications, underscoring the potential for the sex chromosome complement to prime the genome differently in early development with consequences for later expression biases. Upon differentiation of ES cells to cardiac precursors, we found sex-biased expression of key transcription and epigenetic factors, some of which persisted from the undifferentiated state. Using network analyses, we also found that Prdm14 plays a prominent role in regulating a subset of dimorphic expression patterns. To determine whether sex-biased expression is present throughout cardiogenesis, we re-analyzed data from two published studies that sampled the transcriptomes of mouse hearts from 8.5 days post-coitum embryos to neonates and adults. We found sex-biased expression at every stage in heart development, and interestingly, identified a subset of genes that exhibit the same bias across multiple cardiogenic stages.

CONCLUSIONS

Overall, our results support the existence of sexually dimorphic gene expression profiles and regulatory networks at every stage of cardiac development, some of which may be established in early embryogenesis and epigenetically perpetuated.

Conservation, acquisition, and functional impact of sex-biased gene expression in mammals

Sex differences abound in human health and disease, as they do in other mammals used as models. The extent to which sex differences are conserved at the molecular level across species and tissues is unknown. We surveyed sex differences in gene expression in human, macaque, mouse, rat, and dog, across 12 tissues. In each tissue, we identified hundreds of genes with conserved sex-biased expression-findings that, combined with genomic analyses of human height, explain ~12% of the difference in height between females and males. We surmise that conserved sex biases in expression of genes otherwise operating equivalently in females and males contribute to sex differences in traits. However, most sex-biased expression arose during the mammalian radiation, which suggests that careful attention to interspecies divergence is needed when modeling human sex differences.

Sexual dimorphism and rapid turnover in gene expression in pre-reproductive seedlings of a dioecious herb

BACKGROUND AND AIMS

Sexual dimorphism in morphology, physiology or life history traits is common in dioecious plants at reproductive maturity, but it is typically inconspicuous or absent in juveniles. Although plants of different sexes probably begin to diverge in gene expression both before their reproduction commences and before dimorphism becomes readily apparent, to our knowledge transcriptome-wide differential gene expression has yet to be demonstrated for any angiosperm species.

METHODS

The present study documents differences in gene expression in both above- and below-ground tissues of early pre-reproductive individuals of the wind-pollinated dioecious annual herb, Mercurialis annua, which otherwise shows clear sexual dimorphism only at the adult stage.

KEY RESULTS

Whereas males and females differed in their gene expression at the first leaf stage, sex-biased gene expression peaked just prior to, and after, flowering, as might be expected if sexual dimorphism is partly a response to differential costs of reproduction. Sex-biased genes were over-represented among putative sex-linked genes in M. annua but showed no evidence for more rapid evolution than unbiased genes.

CONCLUSIONS

Sex-biased gene expression in M. annua occurs as early as the first whorl of leaves is produced, is highly dynamic during plant development and varies substantially between vegetative tissues.

Pleiotropy Modulates the Efficacy of Selection in Drosophila melanogaster

Pleiotropy is the well-established idea that a single mutation affects multiple phenotypes. If a mutation has opposite effects on fitness when expressed in different contexts, then genetic conflict arises. Pleiotropic conflict is expected to reduce the efficacy of selection by limiting the fixation of beneficial mutations through adaptation, and the removal of deleterious mutations through purifying selection. Although this has been widely discussed, in particular in the context of a putative "gender load," it has yet to be systematically quantified. In this work, we empirically estimate to which extent different pleiotropic regimes impede the efficacy of selection in Drosophila melanogaster. We use whole-genome polymorphism data from a single African population and divergence data from D. simulans to estimate the fraction of adaptive fixations (α), the rate of adaptation (ωA), and the direction of selection (DoS). After controlling for confounding covariates, we find that the different pleiotropic regimes have a relatively small, but significant, effect on selection efficacy. Specifically, our results suggest that pleiotropic sexual antagonism may restrict the efficacy of selection, but that this conflict can be resolved by limiting the expression of genes to the sex where they are beneficial. Intermediate levels of pleiotropy across tissues and life stages can also lead to maladaptation in D. melanogaster, due to inefficient purifying selection combined with low frequency of mutations that confer a selective advantage. Thus, our study highlights the need to consider the efficacy of selection in the context of antagonistic pleiotropy, and of genetic conflict in general.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Evolutionary and developmental dynamics of sex-biased gene expression in common frogs with proto-Y chromosomes

BACKGROUND

The patterns of gene expression on highly differentiated sex chromosomes differ drastically from those on autosomes, due to sex-specific patterns of selection and inheritance. As a result, X chromosomes are often enriched in female-biased genes (feminization) and Z chromosomes in male-biased genes (masculinization). However, it is not known how quickly sexualization of gene expression and transcriptional degeneration evolve after sex-chromosome formation. Furthermore, little is known about how sex-biased gene expression varies throughout development.

RESULTS

We sample a population of common frogs (Rana temporaria) with limited sex-chromosome differentiation (proto-sex chromosome), leaky genetic sex determination evidenced by the occurrence of XX males, and delayed gonadal development, meaning that XY individuals may first develop ovaries before switching to testes. Using high-throughput RNA sequencing, we investigate the dynamics of gene expression throughout development, spanning from early embryo to froglet stages. Our results show that sex-biased expression affects different genes at different developmental stages and increases during development, reaching highest levels in XX female froglets. Additionally, sex-biased gene expression depends on phenotypic, rather than genotypic sex, with similar expression in XX and XY males; correlates with gene evolutionary rates; and is not localized to the proto-sex chromosome nor near the candidate sex-determining gene Dmrt1.

CONCLUSIONS

The proto-sex chromosome of common frogs does not show evidence of sexualization of gene expression, nor evidence for a faster rate of evolution. This challenges the notion that sexually antagonistic genes play a central role in the initial stages of sex-chromosome evolution.

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

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

Evolution of sex-specific pace-of-life syndromes: genetic architecture and physiological mechanisms

Sex differences in life history, physiology, and behavior are nearly ubiquitous across taxa, owing to sex-specific selection that arises from different reproductive strategies of the sexes. The pace-of-life syndrome (POLS) hypothesis predicts that most variation in such traits among individuals, populations, and species falls along a slow-fast pace-of-life continuum. As a result of their different reproductive roles and environment, the sexes also commonly differ in pace-of-life, with important consequences for the evolution of POLS. Here, we outline mechanisms for how males and females can evolve differences in POLS traits and in how such traits can covary differently despite constraints resulting from a shared genome. We review the current knowledge of the genetic basis of POLS traits and suggest candidate genes and pathways for future studies. Pleiotropic effects may govern many of the genetic correlations, but little is still known about the mechanisms involved in trade-offs between current and future reproduction and their integration with behavioral variation. We highlight the importance of metabolic and hormonal pathways in mediating sex differences in POLS traits; however, there is still a shortage of studies that test for sex specificity in molecular effects and their evolutionary causes. Considering whether and how sexual dimorphism evolves in POLS traits provides a more holistic framework to understand how behavioral variation is integrated with life histories and physiology, and we call for studies that focus on examining the sex-specific genetic architecture of this integration.

Faster-X evolution: Theory and evidence from Drosophila

A faster rate of adaptive evolution of X-linked genes compared with autosomal genes can be caused by the fixation of recessive or partially recessive advantageous mutations, due to the full expression of X-linked mutations in hemizygous males. Other processes, including recombination rate and mutation rate differences between X chromosomes and autosomes, may also cause faster evolution of X-linked genes. We review population genetics theory concerning the expected relative values of variability and rates of evolution of X-linked and autosomal DNA sequences. The theoretical predictions are compared with data from population genomic studies of several species of Drosophila. We conclude that there is evidence for adaptive faster-X evolution of several classes of functionally significant nucleotides. We also find evidence for potential differences in mutation rates between X-linked and autosomal genes, due to differences in mutational bias towards GC to AT mutations. Many aspects of the data are consistent with the male hemizygosity model, although not all possible confounding factors can be excluded.

Evidence of sex-bias in gene expression in the brain transcriptome of two populations of rainbow trout (Oncorhynchus mykiss) with divergent life histories

Sex-bias in gene expression is a mechanism that can generate phenotypic variance between the sexes, however, relatively little is known about how patterns of sex-bias vary during development, and how variable sex-bias is between different populations. To that end, we measured sex-bias in gene expression in the brain transcriptome of rainbow trout (Oncorhynchus mykiss) during the first two years of development. Our sampling included from the fry stage through to when O. mykiss either migrate to the ocean or remain resident and undergo sexual maturation. Samples came from two F1 lines: One from migratory steelhead trout and one from resident rainbow trout. All samples were reared in a common garden environment and RNA sequencing (RNA-seq) was used to estimate patterns of gene expression. A total of 1,716 (4.6% of total) genes showed evidence of sex-bias in gene expression in at least one time point. The majority (96.7%) of sex-biased genes were differentially expressed during the second year of development, indicating that patterns of sex-bias in expression are tied to key developmental events, such as migration and sexual maturation. Mapping of differentially expressed genes to the O. mykiss genome revealed that the X chromosome is enriched for female upregulated genes, and this may indicate a lack of dosage compensation in rainbow trout. There were many more sex-biased genes in the migratory line than the resident line suggesting differences in patterns of gene expression in the brain between populations subjected to different forces of selection. Overall, our results suggest that there is considerable variation in the extent and identity of genes exhibiting sex-bias during the first two years of life. These differentially expressed genes may be connected to developmental differences between the sexes, and/or between adopting a resident or migratory life history.

Slow evolution of sex-biased genes in the reproductive tissue of the dioecious plant Salix viminalis

The relative rate of evolution for sex‐biased genes has often been used as a measure of the strength of sex‐specific selection. In contrast to studies in a wide variety of animals, far less is known about the molecular evolution of sex‐biased genes in plants, particularly in dioecious angiosperms. Here, we investigate the gene expression patterns and evolution of sex‐biased genes in the dioecious plant Salix viminalis. We observe lower rates of sequence evolution for male‐biased genes expressed in the reproductive tissue compared to unbiased and female‐biased genes. These results could be partially explained by the lower codon usage bias for male‐biased genes leading to elevated rates of synonymous substitutions compared to unbiased genes. However, the stronger haploid selection in the reproductive tissue of plants, together with pollen competition, would also lead to higher levels of purifying selection acting to remove deleterious variation. Future work should focus on the differential evolution of haploid‐ and diploid‐specific genes to understand the selective dynamics acting on these loci.

Evidence of Adaptive Evolution and Relaxed Constraints in Sex-Biased Genes of South American and West Indies Fruit Flies (Diptera: Tephritidae)

Several studies have demonstrated that genes differentially expressed between sexes (sex-biased genes) tend to evolve faster than unbiased genes, particularly in males. The reason for this accelerated evolution is not clear, but several explanations have involved adaptive and nonadaptive mechanisms. Furthermore, the differences of sex-biased expression patterns of closely related species are also little explored out of Drosophila. To address the evolutionary processes involved with sex-biased expression in species with incipient differentiation, we analyzed male and female transcriptomes of Anastrepha fraterculus and Anastrepha obliqua, a pair of species that have diverged recently, likely in the presence of gene flow. Using these data, we inferred differentiation indexes and evolutionary rates and tested for signals of selection in thousands of genes expressed in head and reproductive transcriptomes from both species. Our results indicate that sex-biased and reproductive-biased genes evolve faster than unbiased genes in both species, which is due to both adaptive pressure and relaxed constraints. Furthermore, among male-biased genes evolving under positive selection, we identified some related to sexual functions such as courtship behavior and fertility. These findings suggest that sex-biased genes may have played important roles in the establishment of reproductive isolation between these species, due to a combination of selection and drift, and unveil a plethora of genetic markers useful for more studies in these species and their differentiation.

The Y Chromosome Modulates Splicing and Sex-Biased Intron Retention Rates in Drosophila

The Drosophila Y chromosome is a 40-Mb segment of mostly repetitive DNA; it harbors a handful of protein-coding genes and a disproportionate amount of satellite repeats, transposable elements, and multicopy DNA arrays. Intron retention (IR) is a type of alternative splicing (AS) event by which one or more introns remain within the mature transcript. IR recently emerged as a deliberate cellular mechanism to modulate gene expression levels and has been implicated in multiple biological processes. However, the extent of sex differences in IR and the contribution of the Y chromosome to the modulation of AS and IR rates has not been addressed. Here we showed pervasive IR in the fruit fly Drosophila melanogaster with thousands of novel IR events, hundreds of which displayed extensive sex bias. The data also revealed an unsuspected role for the Y chromosome in the modulation of AS and IR. The majority of sex-biased IR events introduced premature termination codons and the magnitude of sex bias was associated with gene expression differences between the sexes. Surprisingly, an extra Y chromosome in males (X^YY genotype) or the presence of a Y chromosome in females (X^XY genotype) significantly modulated IR and recapitulated natural differences in IR between the sexes. Our results highlight the significance of sex-biased IR in tuning sex differences and the role of the Y chromosome as a source of variable IR rates between the sexes. Modulation of splicing and IR rates across the genome represent new and unexpected outcomes of the Drosophila Y chromosome.

Mating Changes Sexually Dimorphic Gene Expression in the Seed Beetle Callosobruchus maculatus

Sexually dimorphic phenotypes arise largely from sex-specific gene expression, which has mainly been characterized in sexually naïve adults. However, we expect sexual dimorphism in transcription to be dynamic and dependent on factors such as reproductive status. Mating induces many behavioral and physiological changes distinct to each sex and is therefore expected to activate regulatory changes in many sex-biased genes. Here, we first characterized sexual dimorphism in gene expression in Callosobruchus maculatus seed beetles. We then examined how females and males respond to mating and how it affects sex-biased expression, both in sex-limited (abdomen) and sex-shared (head and thorax) tissues. Mating responses were largely sex-specific and, as expected, females showed more genes responding compared with males (∼2,000 vs. ∼300 genes in the abdomen, ∼500 vs. ∼400 in the head and thorax, respectively). Of the sex-biased genes present in virgins, 16% (1,041 genes) in the abdomen and 17% (243 genes) in the head and thorax altered their relative expression between the sexes as a result of mating. Sex-bias status changed in 2% of the genes in the abdomen and 4% in the head and thorax following mating. Mating responses involved de-feminization of females and, to a lesser extent, de-masculinization of males relative to their virgin state: mating decreased rather than increased dimorphic expression of sex-biased genes. The fact that regulatory changes of both types of sex-biased genes occurred in both sexes suggests that male- and female-specific selection is not restricted to male- and female-biased genes, respectively, as is sometimes assumed.

Comparative transcriptomics reveal developmental turning points during embryogenesis of a hemimetabolous insect, the damselfly Ischnura elegans

Identifying transcriptional changes during embryogenesis is of crucial importance for unravelling evolutionary, molecular and cellular mechanisms that underpin patterning and morphogenesis. However, comparative studies focusing on early/embryonic stages during insect development are limited to a few taxa. Drosophila melanogaster is the paradigm for insect development, whereas comparative transcriptomic studies of embryonic stages of hemimetabolous insects are completely lacking. We reconstructed the first comparative transcriptome covering the daily embryonic developmental progression of the blue-tailed damselfly Ischnura elegans (Odonata), an ancient hemimetabolous representative. We identified a "core" set of 6,794 transcripts - shared by all embryonic stages - which are mainly involved in anatomical structure development and cellular nitrogen compound metabolic processes. We further used weighted gene co-expression network analysis to identify transcriptional changes during Odonata embryogenesis. Based on these analyses distinct clusters of transcriptional active sequences could be revealed, indicating that embryos at different development stages have their own transcriptomic profile according to the developmental events and leading to sequential reprogramming of metabolic and developmental genes. Interestingly, a major change in transcriptionally active sequences is correlated with katatrepsis (revolution) during mid-embryogenesis, a 180° rotation of the embryo within the egg and specific to hemimetabolous insects.

Evo-Devo-EpiR: a genome-wide search platform for epistatic control on the evolution of development

Evo-devo is a theory proposed to study how phenotypes evolve by comparing the developmental processes of different organisms or the same organism experiencing changing environments. It has been recognized that nonallelic interactions at different genes or quantitative trait loci, known as epistasis, may play a pivotal role in the evolution of development, but it has proven difficult to quantify and elucidate this role into a coherent picture. We implement a high-dimensional genome-wide association study model into the evo-devo paradigm and pack it into the R-based Evo-Devo-EpiR, aimed at facilitating the genome-wide landscaping of epistasis for the diversification of phenotypic development. By analyzing a high-throughput assay of DNA markers and their pairs simultaneously, Evo-Devo-EpiR is equipped with a capacity to systematically characterize various epistatic interactions that impact on the pattern and timing of development and its evolution. Enabling a global search for all possible genetic interactions for developmental processes throughout the whole genome, Evo-Devo-EpiR provides a computational tool to illustrate a precise genotype-phenotype map at interface between epistasis, development and evolution.