Disentangling the relationship between sex-biased gene expression and X-linkage

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.

CAGEE: Computational Analysis of Gene Expression Evolution

Despite the increasing abundance of whole transcriptome data, few methods are available to analyze global gene expression across phylogenies. Here, we present a new software package (Computational Analysis of Gene Expression Evolution [CAGEE]) for inferring patterns of increases and decreases in gene expression across a phylogenetic tree, as well as the rate at which these changes occur. In contrast to previous methods that treat each gene independently, CAGEE can calculate genome-wide rates of gene expression, along with ancestral states for each gene. The statistical approach developed here makes it possible to infer lineage-specific shifts in rates of evolution across the genome, in addition to possible differences in rates among multiple tissues sampled from the same species. We demonstrate the accuracy and robustness of our method on simulated data and apply it to a data set of ovule gene expression collected from multiple self-compatible and self-incompatible species in the genus Solanum to test hypotheses about the evolutionary forces acting during mating system shifts. These comparisons allow us to highlight the power of CAGEE, demonstrating its utility for use in any empirical system and for the analysis of most morphological traits. Our software is available at https://github.com/hahnlab/CAGEE/.

Sex-biased gene expression in nutrient-sensing pathways

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

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.

Functional Diversity and Evolution of the Drosophila Sperm Proteome

Spermatozoa are central to fertilization and the evolutionary fitness of sexually reproducing organisms. As such, a deeper understanding of sperm proteomes (and associated reproductive tissues) has proven critical to the advancement of the fields of sexual selection and reproductive biology. Due to their extraordinary complexity, proteome depth-of-coverage is dependent on advancements in technology and related bioinformatics, both of which have made significant advancements in the decade since the last Drosophila sperm proteome was published. Here, we provide an updated version of the Drosophila melanogaster sperm proteome (DmSP3) using improved separation and detection methods and an updated genome annotation. Combined with previous versions of the sperm proteome, the DmSP3 contains a total of 3176 proteins, and we provide the first label-free quantitation of the sperm proteome for 2125 proteins. The top 20 most abundant proteins included the structural elements α- and β-tubulins and sperm leucyl-aminopeptidases. Both gene content and protein abundance were significantly reduced on the X chromosome, consistent with prior genomic studies of X chromosome evolution. We identified 9 of the 16 Y-linked proteins, including known testis-specific male fertility factors. We also identified almost one-half of known Drosophila ribosomal proteins in the DmSP3. The role of this subset of ribosomal proteins in sperm is unknown. Surprisingly, our expanded sperm proteome also identified 122 seminal fluid proteins (Sfps), proteins originally identified in the accessory glands. We show that a significant fraction of 'sperm-associated Sfps' are recalcitrant to concentrated salt and detergent treatments, suggesting this subclass of Sfps are expressed in testes and may have additional functions in sperm, per se. Overall, our results add to a growing landscape of both sperm and seminal fluid protein biology and in particular provides quantitative evidence at the protein level for prior findings supporting the meiotic sex-chromosome inactivation model for male-specific gene and X chromosome evolution.

Induction and inhibition of Drosophila X chromosome gene expression are both impeded by the dosage compensation complex

Sex chromosomes frequently differ from the autosomes in the frequencies of genes with sexually dimorphic or tissue-specific expression. Multiple hypotheses have been put forth to explain the unique gene content of the X chromosome, including selection against male-beneficial X-linked alleles, expression limits imposed by the haploid dosage of the X in males, and interference by the dosage compensation complex on expression in males. Here, we investigate these hypotheses by examining differential gene expression in Drosophila melanogaster following several treatments that have widespread transcriptomic effects: bacterial infection, viral infection, and abiotic stress. We found that genes that are induced (upregulated) by these biotic and abiotic treatments are frequently under-represented on the X chromosome, but so are those that are repressed (downregulated) following treatment. We further show that whether a gene is bound by the dosage compensation complex in males can largely explain the paucity of both up- and downregulated genes on the X chromosome. Specifically, genes that are bound by the dosage compensation complex, or close to a dosage compensation complex high-affinity site, are unlikely to be up- or downregulated after treatment. This relationship, however, could partially be explained by a correlation between differential expression and breadth of expression across tissues. Nonetheless, our results suggest that dosage compensation complex binding, or the associated chromatin modifications, inhibit both up- and downregulation of X chromosome gene expression within specific contexts, including tissue-specific expression. We propose multiple possible mechanisms of action for the effect, including a role of Males absent on the first, a component of the dosage compensation complex, as a dampener of gene expression variance in both males and females. This effect could explain why the Drosophila X chromosome is depauperate in genes with tissue-specific or induced expression, while the mammalian X has an excess of genes with tissue-specific expression.

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

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

Population biology of accessory gland-expressed de novo genes in Drosophila melanogaster

Early work on de novo gene discovery in Drosophila was consistent with the idea that many such genes have male-biased patterns of expression, including a large number expressed in the testis. However, there has been little formal analysis of variation in the abundance and properties of de novo genes expressed in different tissues. Here, we investigate the population biology of recently evolved de novo genes expressed in the Drosophila melanogaster accessory gland, a somatic male tissue that plays an important role in male and female fertility and the post mating response of females, using the same collection of inbred lines used previously to identify testis-expressed de novo genes, thus allowing for direct cross tissue comparisons of these genes in two tissues of male reproduction. Using RNA-seq data, we identify candidate de novo genes located in annotated intergenic and intronic sequence and determine the properties of these genes including chromosomal location, expression, abundance, and coding capacity. Generally, we find major differences between the tissues in terms of gene abundance and expression, though other properties such as transcript length and chromosomal distribution are more similar. We also explore differences between regulatory mechanisms of de novo genes in the two tissues and how such differences may interact with selection to produce differences in D. melanogaster de novo genes expressed in the two tissues.

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

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

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.

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.

Chromosome-scale assembly of the genome of Salix dunnii reveals a male-heterogametic sex determination system on chromosome 7

Sex determination systems in plants can involve either female or male heterogamety (ZW or XY, respectively). Here we used Illumina short reads, Oxford Nanopore Technologies (ONT) long reads and Hi-C reads to assemble the first chromosome-scale genome of a female willow tree (Salix dunnii), and to predict genes using transcriptome sequences and available databases. The final genome sequence of 328 Mb in total was assembled in 29 scaffolds, and includes 31,501 predicted genes. Analyses of short-read sequence data that included female and male plants suggested a male heterogametic sex-determining factor on chromosome 7, implying that, unlike the female heterogamety of most species in the genus Salix, male heterogamety evolved in the subgenus Salix. The S. dunnii sex-linked region occupies about 3.21 Mb of chromosome 7 in females (representing its position in the X chromosome), probably within a pericentromeric region. Our data suggest that this region is enriched for transposable element insertions, and about one-third of its 124 protein-coding genes were gained via duplications from other genome regions. We detect purifying selection on the genes that were ancestrally present in the region, though some have been lost. Transcriptome data from female and male individuals show more male- than female-biased genes in catkin and leaf tissues, and indicate enrichment for male-biased genes in the pseudo-autosomal regions. Our study provides valuable genomic resources for further studies of sex-determining regions in the family Salicaceae, and sex chromosome evolution.

A-to-I RNA editing in honeybees shows signals of adaptation and convergent evolution

Social insects exhibit extensive phenotypic diversities among the genetically similar individuals, suggesting a role for the epigenetic regulations beyond the genome level. The ADAR-mediated adenosine-to-inosine (A-to-I) RNA editing, an evolutionarily conserved mechanism, facilitates adaptive evolution by expanding proteomic diversities. Here, we characterize the A-to-I RNA editome of honeybees (Apis mellifera), identifying 407 high-confidence A-to-I editing sites. Editing is most abundant in the heads and shows signatures for positive selection. Editing behavior differs between foragers and nurses, suggesting a role for editing in caste differentiation. Although only five sites are conserved between bees and flies, an unexpectedly large number of genes exhibit editing in both species, albeit at different locations, including the nonsynonymous auto-editing of Adar. This convergent evolution, where the same target genes independently acquire recoding events in distant diverged clades, together with the signals of adaptation observed in honeybees alone, further supports the notion of recoding being adaptive.

•

Nonsynonymous editing sites in honeybees were under positive selection

•

Differential editing may contribute to the phenotypic diversity between sub-castes

•

Target genes acquire editing in different clades, suggesting convergent evolution

The Influence of Chromosomal Environment on X-Linked Gene Expression in Drosophila melanogaster

Sex chromosomes often differ from autosomes with respect to their gene expression and regulation. In Drosophila melanogaster, X-linked genes are dosage compensated by having their expression upregulated in the male soma, a process mediated by the X-chromosome-specific binding of the dosage compensation complex (DCC). Previous studies of X-linked gene expression found a negative correlation between a gene’s male-to-female expression ratio and its distance to the nearest DCC binding site in somatic tissues, including head and brain, which suggests that dosage compensation influences sex-biased gene expression. A limitation of the previous studies, however, was that they focused on endogenous X-linked genes and, thus, could not disentangle the effects of chromosomal position from those of gene-specific regulation. To overcome this limitation, we examined the expression of an exogenous reporter gene inserted at many locations spanning the X chromosome. We observed a negative correlation between the male-to-female expression ratio of the reporter gene and its distance to the nearest DCC binding site in somatic tissues, but not in gonads. A reporter gene’s location relative to a DCC binding site had greater influence on its expression than the local regulatory elements of neighboring endogenous genes, suggesting that intra-chromosomal variation in the strength of dosage compensation is a major determinant of sex-biased gene expression. Average levels of sex-biased expression did not differ between head and brain, but there was greater positional effect variation in the brain, which may explain the observed excess of endogenous sex-biased genes located on the X chromosome in this tissue.

SAGD: a comprehensive sex-associated gene database from transcriptomes

Many animal species present sex differences. Sex-associated genes (SAGs), which have female-biased or male-biased expression, have major influences on the remarkable sex differences in important traits such as growth, reproduction, disease resistance and behaviors. However, the SAGs resulting in the vast majority of phenotypic sex differences are still unknown. To provide a useful resource for the functional study of SAGs, we manually curated public RNA-seq datasets with paired female and male biological replicates from the same condition and systematically re-analyzed the datasets using standardized methods. We identified 27,793 female-biased SAGs and 64,043 male-biased SAGs from 2,828 samples of 21 species, including human, chimpanzee, macaque, mouse, rat, cow, horse, chicken, zebrafish, seven fly species and five worm species. All these data were cataloged into SAGD, a user-friendly database of SAGs (http://bioinfo.life.hust.edu.cn/SAGD) where users can browse SAGs by gene, species, drug and dataset. In SAGD, the expression, annotation, targeting drugs, homologs, ontology and related RNA-seq datasets of SAGs are provided to help researchers to explore their functions and potential applications in agriculture and human health.

Sex Chromosome Evolution in Muscid Flies

Sex chromosomes and sex determining genes can evolve fast, with the sex-linked chromosomes often differing between closely related species. Population genetics theory has been developed and tested to explain the rapid evolution of sex chromosomes and sex determination. However, we do not know why the sex chromosomes are divergent in some taxa and conserved in others. Addressing this question requires comparing closely related taxa with conserved and divergent sex chromosomes to identify biological features that could explain these differences. Cytological karyotypes suggest that muscid flies (e.g., house fly) and blow flies are such a taxonomic pair. The sex chromosomes appear to differ across muscid species, whereas they are conserved across blow flies. Despite the cytological evidence, we do not know the extent to which muscid sex chromosomes are independently derived along different evolutionary lineages. To address that question, we used genomic and transcriptomic sequence data to identify young sex chromosomes in two closely related muscid species, horn fly (Haematobia irritans) and stable fly (Stomoxys calcitrans). We provide evidence that the nascent sex chromosomes of horn fly and stable fly were derived independently from each other and from the young sex chromosomes of the closely related house fly (Musca domestica). We present three different scenarios that could have given rise to the sex chromosomes of horn fly and stable fly, and we describe how the scenarios could be distinguished. Distinguishing between these scenarios in future work could identify features of muscid genomes that promote sex chromosome divergence.

Repeated sex chromosome evolution in vertebrates supported by expanded avian sex chromosomes

Sex chromosomes have evolved from the same autosomes multiple times across vertebrates, suggesting that selection for recombination suppression has acted repeatedly and independently on certain genetic backgrounds. Here, we perform comparative genomics of a bird clade (larks and their sister lineage; Alaudidae and Panuridae) where multiple autosome-sex chromosome fusions appear to have formed expanded sex chromosomes. We detected the largest known avian sex chromosome (195.3 Mbp) and show that it originates from fusions between parts of four avian chromosomes: Z, 3, 4A and 5. Within these four chromosomes, we found evidence of five evolutionary strata where recombination had been suppressed at different time points, and show that stratum age explained the divergence rate of Z-W gametologs. Next, we analysed chromosome content and found that chromosome 3 was significantly enriched for genes with predicted sex-related functions. Finally, we demonstrate extensive homology to sex chromosomes in other vertebrate lineages: chromosomes Z, 3, 4A and 5 have independently evolved into sex chromosomes in fish (Z), turtles (Z, 5), lizards (Z, 4A), mammals (Z, 4A) and frogs (Z, 3, 4A, 5). Our results provide insights into and support for repeated evolution of sex chromosomes in vertebrates.

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.

Minimal Effects of Proto-Y Chromosomes on House Fly Gene Expression in Spite of Evidence that Selection Maintains Stable Polygenic Sex Determination

Sex determination, the developmental process by which organismal sex is established, evolves fast, often due to changes in the master regulators at the top of the pathway. Additionally, in species with polygenic sex determination, multiple different master regulators segregate as polymorphisms. Understanding the forces that maintain polygenic sex determination can be informative of the factors that drive the evolution of sex determination. The house fly, Musca domestica, is a well-suited model to those ends because natural populations harbor male-determining loci on each of the six chromosomes and a biallelic female determiner. To investigate how natural selection maintains polygenic sex determination in the house fly, we assayed the phenotypic effects of proto-Y chromosomes by performing mRNA-sequencing experiments to measure gene expression in house fly males carrying different proto-Y chromosomes. We find that the proto-Y chromosomes have similar effects as a nonsex-determining autosome. In addition, we created sex-reversed males without any proto-Y chromosomes and they had nearly identical gene expression profiles as genotypic males. Therefore, the proto-Y chromosomes have a minor effect on male gene expression, consistent with previously described minimal X-Y sequence differences. Despite these minimal differences, we find evidence for a disproportionate effect of one proto-Y chromosome on male-biased expression, which could be partially responsible for fitness differences between males with different proto-Y chromosome genotypes. Therefore our results suggest that, if natural selection maintains polygenic sex determination in house fly via gene expression differences, the phenotypes under selection likely depend on a small number of genetic targets.

Genetic constraints on microevolutionary divergence of sex-biased gene expression

The evolution of sex-specific phenotypes is an important dimension of diversification and local adaptation. The sex-dependent regulation of gene expression is considered a key genomic mechanism facilitating sex-dependent adaptation. In many species, genes with male-biased expression evolve faster in DNA sequence and expression level than genes with female-biased or sexually monomorphic expression. While positive selection may be responsible for rapid DNA sequence evolution, why expression of male-biased genes also evolves rapidly remains unclear. Beyond sex differences in selection, some aspects of the genetic architecture of gene expression could contribute to the rapid evolution of male-biased gene expression. First, male-biased genes might simply have greater standing genetic variance than female-biased genes. Second, male-biased genes could be less constrained by pleiotropy, either within or between sexes. Here, we evaluate these alternative explanations on an intraspecific scale using a series of quantitative genetic experiments conducted on natural variation in male and female gene expression in the fly Drosophila serrata Male-biased genes had significantly higher genetic variance than female-biased genes and were generally more narrowly expressed across tissues, suggesting lower within-individual pleiotropy. However, consistent with stronger constraints due to between-sex pleiotropy, their between-sex genetic correlations, r_MF, were higher than for female-biased genes and more strongly negatively associated with sex bias. Using an extensive clinal dataset, we tested whether sex differences in gene expression divergence among populations have been shaped by pleiotropy. Here too, male-biased gene divergence was more strongly associated with between-sex pleiotropy than was female-biased gene divergence. Systematic differences in genetic variance and pleiotropy may be important factors influencing sex-specific adaptation arising through changes in gene expression.This article is part of the theme issue 'Linking local adaptation with the evolution of 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.

Early Sex-Chromosome Evolution in the Diploid Dioecious Plant Mercurialis annua

Suppressed recombination allows divergence between homologous sex chromosomes and the functionality of their genes. Here, we reveal patterns of the earliest stages of sex-chromosome evolution in the diploid dioecious herb Mercurialis annua on the basis of cytological analysis, de novo genome assembly and annotation, genetic mapping, exome resequencing of natural populations, and transcriptome analysis. The genome assembly contained 34,105 expressed genes, of which 10,076 were assigned to linkage groups. Genetic mapping and exome resequencing of individuals across the species range both identified the largest linkage group, LG1, as the sex chromosome. Although the sex chromosomes of M. annua are karyotypically homomorphic, we estimate that about one-third of the Y chromosome, containing 568 transcripts and spanning 22.3 cM in the corresponding female map, has ceased recombining. Nevertheless, we found limited evidence for Y-chromosome degeneration in terms of gene loss and pseudogenization, and most X- and Y-linked genes appear to have diverged in the period subsequent to speciation between M. annua and its sister species M. huetii, which shares the same sex-determining region. Taken together, our results suggest that the M. annua Y chromosome has at least two evolutionary strata: a small old stratum shared with M. huetii, and a more recent larger stratum that is probably unique to M. annua and that stopped recombining ∼1 MYA. Patterns of gene expression within the nonrecombining region are consistent with the idea that sexually antagonistic selection may have played a role in favoring suppressed recombination.

Sexual conflict through mother's curse and father's curse

In contrast with autosomes, lineages of sex chromosomes reside for different amounts of time in males and females, and this transmission asymmetry makes them hotspots for sexual conflict. Similarly, the maternal inheritance of the mitochondrial genome (mtDNA) means that mutations that are beneficial in females can spread in a population even if they are deleterious in males, a form of sexual conflict known as Mother's Curse. While both Mother's Curse and sex chromosome induced sexual conflict have been well studied on their own, the interaction between mitochondrial genes and genes on sex chromosomes is poorly understood. Here, we use analytical models and computer simulations to perform a comprehensive examination of how transmission asymmetries of nuclear, mitochondrial, and sex chromosome-linked genes may both cause and resolve sexual conflicts. For example, the accumulation of male-biased Mother's Curse mtDNA mutations will lead to selection in males for compensatory nuclear modifier loci that alleviate the effect. We show how the Y chromosome, being strictly paternally transmitted provides a particularly safe harbor for such modifiers. This analytical framework also allows us to discover a novel kind of sexual conflict, by which Y chromosome-autosome epistasis may result in the spread of male beneficial but female deleterious mutations in a population. We christen this phenomenon Father's Curse. Extending this analytical framework to ZW sex chromosome systems, where males are the heterogametic sex, we also show how W-autosome epistasis can lead to a novel kind of nuclear Mother's Curse. Overall, this study provides a comprehensive framework to understand how genetic transmission asymmetries may both cause and resolve sexual conflicts.

Genes Relocated Between Drosophila Chromosome Arms Evolve Under Relaxed Selective Constraints Relative to Non-Relocated Genes

Gene duplication creates a second copy of a gene either in tandem to the ancestral locus or dispersed to another chromosomal location. When the ancestral copy of a dispersed duplicate is lost from the genome, it creates the appearance that the gene was "relocated" from the ancestral locus to the derived location. Gene relocations may be as common as canonical dispersed duplications in which both the ancestral and derived copies are retained. Relocated genes appear to be under more selective constraints than the derived copies of canonical duplications, and they are possibly as conserved as single-copy non-relocated genes. To test this hypothesis, we combined comparative genomics, population genetics, gene expression, and functional analyses to assess the selection pressures acting on relocated, duplicated, and non-relocated single-copy genes in Drosophila genomes. We find that relocated genes evolve faster than single-copy non-relocated genes, and there is no evidence that this faster evolution is driven by positive selection. In addition, relocated genes are less essential for viability and male fertility than single-copy non-relocated genes, suggesting that relocated genes evolve fast because of relaxed selective constraints. However, relocated genes evolve slower than the derived copies of canonical dispersed duplicated genes. We therefore conclude that relocated genes are under more selective constraints than canonical duplicates, but are not as conserved as single-copy non-relocated genes.

Transcriptional modulation of entire chromosomes: dosage compensation

Dosage compensation is a regulatory system designed to equalize the transcription output of the genes of the sex chromosomes that are present in different doses in the sexes (X or Z chromosome, depending on the animal species involved). Different mechanisms of dosage compensation have evolved in different animal groups. In Drosophila males, a complex (male-specific lethal) associates with the X chromosome and enhances the activity of most X-linked genes by increasing the rate of RNAPII elongation. In Caenorhabditis, a complex (dosage compensation complex) that contains a number of proteins involved in condensing chromosomes decreases the level of transcription of both X chromosomes in the XX hermaphrodite. In mammals, dosage compensation is achieved by the inactivation, early during development, of most X-linked genes on one of the two X chromosomes in females. The mechanism involves the synthesis of an RNA (Tsix) that protects one of the two Xs from inactivation, and of another RNA (Xist) that coats the other X chromosome and recruits histone and DNA modifying enzymes. This review will focus on the current progress in understanding the dosage compensation mechanisms in the three taxa where it has been best studied at the molecular level: flies, round worms and mammals.

Evolution of Sex Chromosome Dosage Compensation in Animals: A Beautiful Theory, Undermined by Facts and Bedeviled by Details

Many animals with genetic sex determination harbor heteromorphic sex chromosomes, where the heterogametic sex has half the gene dose of the homogametic sex. This imbalance, if reflected in the abundance of transcripts or proteins, has the potential to deleteriously disrupt interactions between X-linked and autosomal loci in the heterogametic sex. Classical theory predicts that molecular mechanisms will evolve to provide dosage compensation that recovers expression levels comparable to ancestral expression prior to sex chromosome divergence. Such dosage compensating mechanisms may also, secondarily, result in balanced sex-linked gene expression between males and females. However, numerous recent studies addressing sex chromosome dosage compensation (SCDC) in a diversity of animals have yielded a surprising array of patterns concerning dosage compensation in the heterogametic sex, as well as dosage balance between sexes. These results substantially contradict longstanding theory, catalyzing both novel perspectives and new approaches in dosage compensation research. In this review, we summarize the theory, analytical approaches, and recent results concerning evolutionary patterns of SCDC in animals. We also discuss methodological challenges and discrepancies encountered in this research, which often underlie conflicting results. Finally, we discuss what outstanding questions and opportunities exist for future research on SCDC.

Distinct core promoter codes drive transcription initiation at key developmental transitions in a marine chordate

Background

Development is largely driven by transitions between transcriptional programs. The initiation of transcription at appropriate sites in the genome is a key component of this and yet few rules governing selection are known. Here, we used cap analysis of gene expression (CAGE) to generate bp-resolution maps of transcription start sites (TSSs) across the genome of Oikopleura dioica, a member of the closest living relatives to vertebrates.

Results

Our TSS maps revealed promoter features in common with vertebrates, as well as striking differences, and uncovered key roles for core promoter elements in the regulation of development. During spermatogenesis there is a genome-wide shift in mode of transcription initiation characterized by a novel core promoter element. This element was associated with > 70% of male-specific transcription, including the use of cryptic internal promoters within operons. In many cases this led to the exclusion of trans-splice sites, revealing a novel mechanism for regulating which mRNAs receive the spliced leader. Binding of the cell cycle regulator, E2F1, is enriched at the TSS of maternal genes in endocycling nurse nuclei. In addition, maternal promoters lack the TATA-like element found in zebrafish and have broad, rather than sharp, architectures with ordered nucleosomes. Promoters of ribosomal protein genes lack the highly conserved TCT initiator. We also report an association between DNA methylation on transcribed gene bodies and the TATA-box.

Conclusions

Our results reveal that distinct functional promoter classes and overlapping promoter codes are present in protochordates like in vertebrates, but show extraordinary lineage-specific innovations. Furthermore, we uncover a genome-wide, developmental stage-specific shift in the mode of TSS selection. Our results provide a rich resource for the study of promoter structure and evolution in Metazoa.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4504-5) contains supplementary material, which is available to authorized users.

Genomic Loss and Silencing on the Y Chromosomes of Rumex

Across many unrelated lineages of plants and animals, Y chromosomes show a recurrent pattern of gene degeneration and loss, but the relative importance of inefficient selection, adaptive gene silencing, and neutral genetic drift in causing degeneration remain poorly understood. Here, we use next-generation genome and transcriptome sequencing to investigate patterns of ongoing Y chromosome degeneration in two annual plant species of Rumex (Polygonaceae) differing in their degree of degeneration and sex chromosome heteromorphism. We find evidence for both gene loss as well as silencing in these young plant sex chromosomes. Our analyses revealed significantly more gene deletion relative to silencing in R. rothschildianus, which has had a larger nonrecombining region for a longer period than R. hastatulus, consistent with this system being at a more advanced stage of degeneration. Intra- and interspecific comparisons of genomic coverage and heterozygosity indicated that loss of expression precedes gene deletion, implying that the final stages of mutation accumulation and gene loss may often occur neutrally. We found no evidence for adaptive silencing of genes that have lost expression. Our results suggest that the initial spread of deleterious regulatory variants and/or epigenetic silencing may be an important driver of early degeneration of Y chromosomes.

Rapidly Evolving Toll-3/4 Genes Encode Male-Specific Toll-Like Receptors in Drosophila

Animal Toll-like receptors (TLRs) have evolved through a pattern of duplication and divergence. Whereas mammalian TLRs directly recognize microbial ligands, Drosophila Tolls bind endogenous ligands downstream of both developmental and immune signaling cascades. Here, we find that most Toll genes in Drosophila evolve slowly with little gene turnover (gains/losses), consistent with their important roles in development and indirect roles in microbial recognition. In contrast, we find that the Toll-3/4 genes have experienced an unusually rapid rate of gene gains and losses, resulting in lineage-specific Toll-3/4s and vastly different gene repertoires among Drosophila species, from zero copies (e.g., D. mojavensis) to nineteen copies (e.g., D. willistoni). In D. willistoni, we find strong evidence for positive selection in Toll-3/4 genes, localized specifically to an extracellular region predicted to overlap with the binding site of Spätzle, the only known ligand of insect Tolls. However, because Spätzle genes are not experiencing similar selective pressures, we hypothesize that Toll-3/4s may be rapidly evolving because they bind to a different ligand, akin to TLRs outside of insects. We further find that most Drosophila Toll-3/4 genes are either weakly expressed or expressed exclusively in males, specifically in the germline. Unlike other Toll genes in D. melanogaster, Toll-3, and Toll-4 have apparently escaped from essential developmental roles, as knockdowns have no substantial effects on viability or male fertility. Based on these findings, we propose that the Toll-3/4 genes represent an exceptionally rapidly evolving lineage of Drosophila Toll genes, which play an unusual, as-yet-undiscovered role in the male germline.

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.

The Z chromosome is enriched for sperm proteins in two divergent species of Lepidoptera

Genes that promote sexual conflict, such as those with a sex-limited fitness benefit, are expected to accumulate differentially on sex chromosomes relative to autosomes. Few tests of this hypothesis exist for male homogametic (ZZ) taxa, however, and most use RNA expression data to identify such genes. Here, we employ a different identification method by using proteomic analysis of sperm cells to identify genes with a sex-limited benefit. We tested for a bias in genomic location of sperm protein genes in two species of Lepidoptera. An excess of sperm protein genes was identified on the Z chromosomes of both the Carolina sphinx moth (Manduca sexta) and the monarch butterfly (Danaus plexippus). Taking into consideration a Z-autosome fusion in monarchs, we discover that the ancestrally sex-linked portion of the genome is the source of this enrichment, while the newly sex-linked portion still appears similar to autosomes in relative abundance of sperm protein genes. Together, these results point to an enrichment of male-beneficial genes on the Z chromosome and demonstrate the usefulness of proteomic datasets in sexual conflict research.

The house fly Y Chromosome is young and minimally differentiated from its ancient X Chromosome partner

Canonical ancient sex chromosome pairs consist of a gene rich X (or Z) Chromosome and a male-limited (or female-limited) Y (or W) Chromosome that is gene poor. In contrast to highly differentiated sex chromosomes, nascent sex chromosome pairs are homomorphic or very similar in sequence content. Nascent sex chromosomes can arise if an existing sex chromosome fuses to an autosome or an autosome acquires a new sex-determining locus/allele. Sex chromosomes often differ between closely related species and can even be polymorphic within species, suggesting that nascent sex chromosomes arise frequently over the course of evolution. Previously documented sex chromosome transitions involve changes to both members of the sex chromosome pair (X and Y, or Z and W). The house fly has sex chromosomes that resemble the ancestral fly karyotype that originated ∼100 million yr ago; therefore, the house fly is expected to have X and Y Chromosomes with different gene content. We tested this hypothesis using whole-genome sequencing and transcriptomic data, and we discovered little evidence for genetic differentiation between the X and Y in house fly. We propose that the house fly has retained the ancient X Chromosome, but the ancestral Y was replaced by an X Chromosome carrying a new male determining gene. Our proposed hypothesis provides a mechanism for how one member of a sex chromosome pair can experience evolutionary turnover while the other member remains unaffected.

X-linkage is not a general inhibitor of tissue-specific gene expression in Drosophila melanogaster

As a consequence of its difference in copy number between males and females, the X chromosome is subject to unique evolutionary forces and gene regulatory mechanisms. Previous studies of Drosophila melanogaster have shown that the expression of X-linked, testis-specific reporter genes is suppressed in the male germline. However, it is not known whether this phenomenon is restricted to testis-expressed genes or if it is a more general property of genes with tissue-specific expression, which are also underrepresented on the X chromosome. To test this, we compared the expression of three tissue-specific reporter genes (ovary, accessory gland and Malpighian tubule) inserted at various autosomal and X-chromosomal locations. In contrast to testis-specific reporter genes, we found no reduction of X-linked expression in any of the other tissues. In accessory gland and Malpighian tubule, we detected higher expression of the X-linked reporter genes, which suggests that they are at least partially dosage compensated. We found no difference in the tissue-specificity of X-linked and autosomal reporter genes. These findings indicate that, in general, the X chromosome is not a detrimental environment for tissue-specific gene expression and that the suppression of X-linked expression is limited to the male germline.

Convergent recombination suppression suggests role of sexual selection in guppy sex chromosome formation

Sex chromosomes evolve once recombination is halted between a homologous pair of chromosomes. The dominant model of sex chromosome evolution posits that recombination is suppressed between emerging X and Y chromosomes in order to resolve sexual conflict. Here we test this model using whole genome and transcriptome resequencing data in the guppy, a model for sexual selection with many Y-linked colour traits. We show that although the nascent Y chromosome encompasses nearly half of the linkage group, there has been no perceptible degradation of Y chromosome gene content or activity. Using replicate wild populations with differing levels of sexually antagonistic selection for colour, we also show that sexual selection leads to greater expansion of the non-recombining region and increased Y chromosome divergence. These results provide empirical support for longstanding models of sex chromosome catalysis, and suggest an important role for sexual selection and sexual conflict in genome evolution.

The transcriptional architecture of phenotypic dimorphism

The profound differences in gene expression between the sexes are increasingly used to study the molecular basis of sexual dimorphism, sexual selection and sexual conflict. Studies of transcriptional architecture, based on comparisons of gene expression, have also been implemented for a wide variety of other intra-specific polymorphisms. These efforts are based on key assumptions regarding the relationship between transcriptional architecture, phenotypic variation and the target of selection. Some of these assumptions are better supported by available evidence than others. In all cases, the evidence is largely circumstantial, leaving considerable gaps in our understanding of the relationship between transcriptional and phenotypic dimorphism.

Evolution of sex-biased gene expression in a dioecious plant

Separate sexes and sex-biased gene expression have repeatedly evolved in animals and plants, but the underlying changes in gene expression remain unknown. Here, we studied a pair of plant species, one in which separate sexes and sex chromosomes evolved recently and one which maintained hermaphrodite flowers resembling the ancestral state, to reconstruct expression changes associated with the evolution of dioecy. We found that sex-biased gene expression has evolved in autosomal and sex-linked genes in the dioecious species. Most expression changes relative to hermaphrodite flowers occurred in females rather than males, with higher and lower expression in females leading to female-biased and male-biased expression, respectively. Expression changes were more common in genes located on the sex chromosomes than the autosomes and led to feminization of the X chromosome and masculinization of the Y chromosome. Our results support a scenario in which sex-biased gene expression evolved during the evolution of dioecy to resolve intralocus sexual conflicts over the allocation of resources.

De Novo Transcriptome Assembly and Sex-Biased Gene Expression in the Cyclical Parthenogenetic Daphnia galeata

Daphnia species have become models for ecological genomics and exhibit interesting features, such as high phenotypic plasticity and a densely packed genome with many lineage-specific genes. They are also cyclic parthenogenetic, with alternating asexual and sexual cycles and environmental sex determination. Here, we present a de novo transcriptome assembly of over 32,000 D. galeata genes and use it to investigate gene expression in females and spontaneously produced males of two clonal lines derived from lakes in Germany and the Czech Republic. We find that only a low percentage (18%) of genes shows sex-biased expression and that there are many more female-biased gene (FBG) than male-biased gene (MBG). Furthermore, FBGs tend to be more conserved between species than MBGs in both sequence and expression. These patterns may be a consequence of cyclic parthenogenesis leading to a relaxation of purifying selection on MBGs. The two clonal lines show considerable differences in both number and identity of sex-biased genes, suggesting that they may have reproductive strategies differing in their investment in sexual reproduction. Orthologs of key genes in the sex determination and juvenile hormone pathways, which are thought to be important for the transition from asexual to sexual reproduction, are present in D. galeata and highly conserved among Daphnia species.

Repeated Duplication of Argonaute2 Is Associated with Strong Selection and Testis Specialization in Drosophila

Argonaute2 (Ago2) is a rapidly evolving nuclease in the Drosophila melanogaster RNA interference (RNAi) pathway that targets viruses and transposable elements in somatic tissues. Here we reconstruct the history of Ago2 duplications across the D. obscura group and use patterns of gene expression to infer new functional specialization. We show that some duplications are old, shared by the entire species group, and that losses may be common, including previously undetected losses in the lineage leading to D. pseudoobscura We find that while the original (syntenic) gene copy has generally retained the ancestral ubiquitous expression pattern, most of the novel Ago2 paralogs have independently specialized to testis-specific expression. Using population genetic analyses, we show that most testis-specific paralogs have significantly lower genetic diversity than the genome-wide average. This suggests recent positive selection in three different species, and model-based analyses provide strong evidence of recent hard selective sweeps in or near four of the six D. pseudoobscura Ago2 paralogs. We speculate that the repeated evolution of testis specificity in obscura group Ago2 genes, combined with their dynamic turnover and strong signatures of adaptive evolution, may be associated with highly derived roles in the suppression of transposable elements or meiotic drive. Our study highlights the lability of RNAi pathways, even within well-studied groups such as Drosophila, and suggests that strong selection may act quickly after duplication in RNAi pathways, potentially giving rise to new and unknown RNAi functions in nonmodel species.

Sex Differences in Drosophila Somatic Gene Expression: Variation and Regulation by doublesex

Sex differences in gene expression have been widely studied in Drosophila melanogaster. Sex differences vary across strains, but many molecular studies focus on only a single strain, or on genes that show sexually dimorphic expression in many strains. How extensive variability is and whether this variability occurs among genes regulated by sex determination hierarchy terminal transcription factors is unknown. To address these questions, we examine differences in sexually dimorphic gene expression between two strains in Drosophila adult head tissues. We also examine gene expression in doublesex (dsx) mutant strains to determine which sex-differentially expressed genes are regulated by DSX, and the mode by which DSX regulates expression. We find substantial variation in sex-differential expression. The sets of genes with sexually dimorphic expression in each strain show little overlap. The prevalence of different DSX regulatory modes also varies between the two strains. Neither the patterns of DSX DNA occupancy, nor mode of DSX regulation explain why some genes show consistent sex-differential expression across strains. We find that the genes identified as regulated by DSX in this study are enriched with known sites of DSX DNA occupancy. Finally, we find that sex-differentially expressed genes and genes regulated by DSX are highly enriched on the fourth chromosome. These results provide insights into a more complete pool of potential DSX targets, as well as revealing the molecular flexibility of DSX regulation.

Testis-specific RNA-Seq of Calligrapha (Chrysomelidae) as a transcriptomic resource for male-biased gene inquiry in Coleoptera

We report the architecture of testis transcriptomes of four closely related species of Calligrapha (Chrysomelidae) beetles, which diverged during the last 3 million years. Five cDNA libraries were sequenced using Illumina HiSeq technology, retrieving 102 884-176 514 assembled contigs, of which ~33-45% of these longer than 499 nt were functionally annotated. Annotation and sequence similarity comparisons of these libraries revealed high homogeneity in gene composition and the presence of several functional candidates related to reproduction or reproductive processes (0.72-1.08% of annotated sequences). Stringent sequence similarity analyses of these transcriptomes against empirically demonstrated male-biased genes in Drosophila melanogaster and Tribolium castaneum allowed the identification of 77 homologues in Calligrapha, possible candidates of male-biased expression. Some of these genes - including CG9313, Tektin-A or tomboy40 - were confirmed as orthologs of these male-biased genes using phylogenetic inference and available model insect data, increasing our confidence that they represent functional homologues too. Our transcriptomes are a valuable transcriptomic resource for the analysis of male-biased genes in Calligrapha, which has the added interest of including several female-only species. But it simultaneously represents a landmark for similar studies in Coleoptera, broadening the taxonomic diversity currently represented by the model species T. castaneum, and incipient genomic data in other herbivorous lineages, including weevils, longhorn beetles and leaf beetles.

Comparative Expression Dynamics of Intergenic Long Noncoding RNAs in the Genus Drosophila

Thousands of long noncoding RNAs (lncRNAs) have been annotated in eukaryotic genomes, but comparative transcriptomic approaches are necessary to understand their biological impact and evolution. To facilitate such comparative studies in Drosophila, we identified and characterized lncRNAs in a second Drosophilid—the evolutionary model Drosophila pseudoobscura. Using RNA-Seq and computational filtering of protein-coding potential, we identified 1,589 intergenic lncRNA loci in D. pseudoobscura. We surveyed multiple sex-specific developmental stages and found, like in Drosophila melanogaster, increasingly prolific lncRNA expression through male development and an overrepresentation of lncRNAs in the testes. Other trends seen in D. melanogaster, like reduced pupal expression, were not observed. Nonrandom distributions of female-biased and non-testis-specific male-biased lncRNAs between the X chromosome and autosomes are consistent with selection-based models of gene trafficking to optimize genomic location of sex-biased genes. The numerous testis-specific lncRNAs, however, are randomly distributed between the X and autosomes, and we cannot reject the hypothesis that many of these are likely to be spurious transcripts. Finally, using annotated lncRNAs in both species, we identified 134 putative lncRNA homologs between D. pseudoobscura and D. melanogaster and find that many have conserved developmental expression dynamics, making them ideal candidates for future functional analyses.

Contrasting Levels of Molecular Evolution on the Mouse X Chromosome

The mammalian X chromosome has unusual evolutionary dynamics compared to autosomes. Faster-X evolution of spermatogenic protein-coding genes is known to be most pronounced for genes expressed late in spermatogenesis, but it is unclear if these patterns extend to other forms of molecular divergence. We tested for faster-X evolution in mice spanning three different forms of molecular evolution-divergence in protein sequence, gene expression, and DNA methylation-across different developmental stages of spermatogenesis. We used FACS to isolate individual cell populations and then generated cell-specific transcriptome profiles across different stages of spermatogenesis in two subspecies of house mice (Mus musculus), thereby overcoming a fundamental limitation of previous studies on whole tissues. We found faster-X protein evolution at all stages of spermatogenesis and faster-late protein evolution for both X-linked and autosomal genes. In contrast, there was less expression divergence late in spermatogenesis (slower late) on the X chromosome and for autosomal genes expressed primarily in testis (testis-biased). We argue that slower-late expression divergence reflects strong regulatory constraints imposed during this critical stage of sperm development and that these constraints are particularly acute on the tightly regulated sex chromosomes. We also found slower-X DNA methylation divergence based on genome-wide bisulfite sequencing of sperm from two species of mice (M. musculus and M. spretus), although it is unclear whether slower-X DNA methylation reflects development constraints in sperm or other X-linked phenomena. Our study clarifies key differences in patterns of regulatory and protein evolution across spermatogenesis that are likely to have important consequences for mammalian sex chromosome evolution, male fertility, and speciation.

The potential role of sexual conflict and sexual selection in shaping the genomic distribution of Mito-nuclear genes

Mitochondrial interactions with the nuclear genome represent one of life's most important co-evolved mutualisms. In many organisms, mitochondria are maternally inherited, and in these cases, co-transmission between the mitochondrial and nuclear genes differs across different parts of the nuclear genome, with genes on the X chromosome having two-third probability of co-transmission, compared with one-half for genes on autosomes. These asymmetrical inheritance patterns of mitochondria and different parts of the nuclear genome have the potential to put certain gene combinations in inter-genomic co-adaptation or conflict. Previous work in mammals found strong evidence that the X chromosome has a dearth of genes that interact with the mitochondria (mito-nuclear genes), suggesting that inter-genomic conflict might drive genes off the X onto the autosomes for their male-beneficial effects. Here, we developed this idea to test coadaptation and conflict between mito-nuclear gene combinations across phylogenetically independent sex chromosomes on a far broader scale. We found that, in addition to therian mammals, only Caenorhabditis elegans showed an under-representation of mito-nuclear genes on the sex chromosomes. The remaining species studied showed no overall bias in their distribution of mito-nuclear genes. We discuss possible factors other than inter-genomic conflict that might drive the genomic distribution of mito-nuclear genes.

Dosage Compensation and the Distribution of Sex-Biased Gene Expression in Drosophila: Considerations and Genomic Constraints

Several studies in Drosophila have shown a paucity of male-biased genes (i.e., genes that express higher in males than in females) on the X chromosome. Dosage compensation (DC) is a regulatory mechanism of gene expression triggered in males that hypertranscribes the X-linked genes to the level of transcription in females. There are currently two different hypotheses about the effects of DC on the distribution of male-biased genes: (1) it might limit male-expression level, or (2) it might interfere with the male upregulation of gene expression. Here, we used previously published gene expression datasets to reevaluate both hypotheses and introduce a mutually exclusive prediction that helped us to reject the hypothesis that the paucity of male-biased genes in the X chromosome is due to a limit in the male-expression level. Our analysis also uncovers unanticipated details about how DC interferes with the genomic distribution of both, male-biased and female-biased genes. We suggest that DC actually interferes with female downregulation of gene expression and not male upregulation, as previously suggested.

Spermatogenesis Drives Rapid Gene Creation and Masculinization of the X Chromosome in Stalk-Eyed Flies (Diopsidae)

Throughout their evolutionary history, genomes acquire new genetic material that facilitates phenotypic innovation and diversification. Developmental processes associated with reproduction are particularly likely to involve novel genes. Abundant gene creation impacts the evolution of chromosomal gene content and general regulatory mechanisms such as dosage compensation. Numerous studies in model organisms have found complex and, at times contradictory, relationships among these genomic attributes highlighting the need to examine these patterns in other systems characterized by abundant sexual selection. Therefore, we examined the association among novel gene creation, tissue-specific gene expression, and chromosomal gene content within stalk-eyed flies. Flies in this family are characterized by strong sexual selection and the presence of a newly evolved X chromosome. We generated RNA-seq transcriptome data from the testes for three species within the family and from seven additional tissues in the highly dimorphic species, Teleopsis dalmanni. Analysis of dipteran gene orthology reveals dramatic testes-specific gene creation in stalk-eyed flies, involving numerous gene families that are highly conserved in other insect groups. Identification of X-linked genes for the three species indicates that the X chromosome arose prior to the diversification of the family. The most striking feature of this X chromosome is that it is highly masculinized, containing nearly twice as many testes-specific genes as expected based on its size. All the major processes that may drive differential sex chromosome gene content—creation of genes with male-specific expression, development of male-specific expression from pre-existing genes, and movement of genes with male-specific expression—are elevated on the X chromosome of T. dalmanni. This masculinization occurs despite evidence that testes expressed genes do not achieve the same levels of gene expression on the X chromosome as they do on the autosomes.