Peter W. Harrison

Sexual selection drives evolution and rapid turnover of male gene expression

Significance Genes with different expression between males and females (sex-biased genes) show rapid rates of sequence and expression divergence in a range of taxa. These characteristics have led many to assume that sex-biased genes are the product of sexual selection and sexual conflict, but this assumption remains to be rigorously tested. Using a phylogenetically controlled analysis of birds that exhibit diverse levels of sexual selection, we show a rapid turnover in sex-biased gene expression primarily through evolution of male expression levels and that the degree of sexual selection predicts the proportion of male-biased genes but does not account for rates of coding sequence evolution. We also discuss the impact of allometry on gene expression studies, an issue rarely discussed in the literature. The profound and pervasive differences in gene expression observed between males and females, and the unique evolutionary properties of these genes in many species, have led to the widespread assumption that they are the product of sexual selection and sexual conflict. However, we still lack a clear understanding of the connection between sexual selection and transcriptional dimorphism, often termed sex-biased gene expression. Moreover, the relative contribution of sexual selection vs. drift in shaping broad patterns of expression, divergence, and polymorphism remains unknown. To assess the role of sexual selection in shaping these patterns, we assembled transcriptomes from an avian clade representing the full range of sexual dimorphism and sexual selection. We use these species to test the links between sexual selection and sex-biased gene expression evolution in a comparative framework. Through ancestral reconstruction of sex bias, we demonstrate a rapid turnover of sex bias across this clade driven by sexual selection and show it to be primarily the result of expression changes in males. We use phylogenetically controlled comparative methods to demonstrate that phenotypic measures of sexual selection predict the proportion of male-biased but not female-biased gene expression. Although male-biased genes show elevated rates of coding sequence evolution, consistent with previous reports in a range of taxa, there is no association between sexual selection and rates of coding sequence evolution, suggesting that expression changes may be more important than coding sequence in sexual selection. Taken together, our results highlight the power of sexual selection to act on gene expression differences and shape genome evolution.

Masculinization of gene expression is associated with exaggeration of male sexual dimorphism.

Gene expression differences between the sexes account for the majority of sexually dimorphic phenotypes, and the study of sex-biased gene expression is important for understanding the genetic basis of complex sexual dimorphisms. However, it has been difficult to test the nature of this relationship due to the fact that sexual dimorphism has traditionally been conceptualized as a dichotomy between males and females, rather than an axis with individuals distributed at intermediate points. The wild turkey (Meleagris gallopavo) exhibits just this sort of continuum, with dominant and subordinate males forming a gradient in male secondary sexual characteristics. This makes it possible for the first time to test the correlation between sex-biased gene expression and sexually dimorphic phenotypes, a relationship crucial to molecular studies of sexual selection and sexual conflict. Here, we show that subordinate male transcriptomes show striking multiple concordances with their relative phenotypic sexual dimorphism. Subordinate males were clearly male rather than intersex, and when compared to dominant males, their transcriptomes were simultaneously demasculinized for male-biased genes and feminized for female-biased genes across the majority of the transcriptome. These results provide the first evidence linking sexually dimorphic transcription and sexually dimorphic phenotypes. More importantly, they indicate that evolutionary changes in sexual dimorphism can be achieved by varying the magnitude of sex-bias in expression across a large proportion of the coding content of a genome.

The Ontogeny and Evolution of Sex-Biased Gene Expression in Drosophila melanogaster

Sexually dimorphic phenotypes are thought to largely result from sex differences in gene expression, and genes with sex-biased expression have been well characterized in adults of many species. Although most sexual dimorphisms manifest in adults, many result from sex-specific developmental trajectories, implying that juveniles may exhibit significant levels of sex-biased expression. However, it is unclear how much sex-biased expression occurs before reproductive maturity and whether preadult sex-biased genes should exhibit the same evolutionary dynamics observed for adult sex-biased genes. In order to understand the continuity, or lack thereof, and evolutionary dynamics of sex-biased expression throughout the life cycle, we examined sex-biased genes in pre-gonad tissue of two preadult stages and compared them with the adult gonad of Drosophila melanogaster. We found that the majority of the genome is sex-biased at some point in the life cycle, with some genes exhibiting conserved sex-biased expression and others displaying stage-specific sex bias. Our results also reveal a far more complex pattern of evolution for sex-biased genes throughout development. The most rapid evolutionary divergence occurred in genes expressed only in larvae within each sex, compared with continuously expressed genes. In females—but not males—this pattern appeared to be due to relaxed purifying selection in larva-limited genes. Furthermore, genes that retained male bias throughout life evolved more rapidly than stage-specific male-biased genes, due to stronger purifying selection in stage-specific genes. However, female-biased genes that were specific to larvae evolved most rapidly, a pattern that could not be definitively attributed to differences in adaptive evolution or purifying selection, suggesting that pleiotropic constraints on protein-coding sequences can arise when genes are broadly expressed across developmental stages. These results indicate that the signature of sex-specific selection can be detected well before reproductive maturity and is strongest during development.

Incomplete sex chromosome dosage compensation in the Indian meal moth, Plodia interpunctella, based on de novo transcriptome assembly.

Males and females experience differences in gene dose for loci in the nonrecombining region of heteromorphic sex chromosomes. If not compensated, this leads to expression imbalances, with the homogametic sex on average exhibiting greater expression due to the doubled gene dose. Many organisms with heteromorphic sex chromosomes display global dosage compensation mechanisms, which equalize gene expression levels between the sexes. However, birds and Schistosoma have been previously shown to lack chromosome-wide dosage compensation mechanisms, and the status in other female heterogametic taxa including Lepidoptera remains unresolved. To further our understanding of dosage compensation in female heterogametic taxa and to resolve its status in the lepidopterans, we assessed the Indian meal moth, Plodia interpunctella. As P. interpunctella lacks a complete reference genome, we conducted de novo transcriptome assembly combined with orthologous genomic location prediction from the related silkworm genome, Bombyx mori, to compare Z-linked and autosomal gene expression levels for each sex. We demonstrate that P. interpunctella lacks complete Z chromosome dosage compensation, female Z-linked genes having just over half the expression level of males and autosomal genes. This finding suggests that the Lepidoptera and possibly all female heterogametic taxa lack global dosage compensation, although more species will need to be sampled to confirm this assertion.

The evolution of gene expression and the transcriptome–phenotype relationship

Changes in gene expression underlie the adaptive evolution in many complex phenotypes, and the recent increase in the availability of multi-species comparative transcriptome data has made it possible to scan whole transcriptomes for loci that have experienced adaptive changes in expression. However, despite the increase in data availability, current models of gene expression evolution often do not account for the complexities and inherent noise associated with transcriptome data. Additionally, in contrast to current models of gene sequence evolution, models of transcriptome evolution often lack the sophistication to effectively determine whether transcriptional differences between species or within a clade are the result of neutral or adaptive processes. In this review, we discuss the tools, methods and models that define our current understanding of the relationship between gene expression and complex phenotype evolution. Our goal is to summarize what we know about the evolution of global gene expression patterns underlying complex traits, as well to identify some of the questions that remain to be answered.

RNA-seq analysis reveals significant transcriptome changes in turbot (Scophthalmus maximus) suffering severe enteromyxosis

BackgroundEnteromyxosis caused by the intestinal myxozoan parasite Enteromyxum scophthalmi is a serious threat for turbot (Scophthalmus maximus, L.) aquaculture, causing severe catarrhal enteritis leading to a cachectic syndrome, with no therapeutic options available. There are still many aspects of host-parasite interaction and disease pathogenesis that are yet to be elucidated, and to date, no analysis of the transcriptomic changes induced by E. scophthalmi in turbot organs has been conducted. In this study, RNA-seq technology was applied to head kidney, spleen and pyloric caeca of severely infected turbot with the aim of furthering our understanding of the pathogenetic mechanisms and turbot immune response against enteromyxosis.ResultsA huge amount of information was generated with more than 23,000 identified genes in the three organs, amongst which 4,762 were differently expressed (DE) between infected and control fish. Associate gene functions were studied based on gene ontology terms and available literature, and the most interesting DE genes were classified into five categories: 1) immune and defence response; 2) apoptosis and cell proliferation; 3) iron metabolism and erythropoiesis; 4) cytoskeleton and extracellular matrix and 5) metabolism and digestive function. The analysis of down-regulated genes of the first category revealed evidences of a connexion failure between innate and adaptive immune response, especially represented by a high number of DE interferon-related genes in the three organs. Furthermore, we found an intense activation of local immune response at intestinal level that appeared exacerbated, whereas in kidney and spleen genes involved in adaptive immune response were mainly down-regulated. The apoptotic machinery was only clearly activated in pyloric caeca, while kidney and spleen showed a marked depression of genes related to erythropoiesis, probably related to disorders in iron homeostasis. The genetic signature of the causes and consequences of cachexia was also demonstrated by the down-regulation of the genes encoding structural proteins and those involved in the digestive metabolism.ConclusionsThis transcriptomic study has enabled us to gain a better understanding of the pathogenesis of enteromyxosis and identify a large number of DE target genes that bring us closer to the development of strategies designed to effectively combat this pathogen.

Variation in promiscuity and sexual selection drives avian rate of Faster‐Z evolution

Higher rates of coding sequence evolution have been observed on the Z chromosome relative to the autosomes across a wide range of species. However, despite a considerable body of theory, we lack empirical evidence explaining variation in the strength of the Faster‐Z Effect. To assess the magnitude and drivers of Faster‐Z Evolution, we assembled six de novo transcriptomes, spanning 90 million years of avian evolution. Our analysis combines expression, sequence and polymorphism data with measures of sperm competition and promiscuity. In doing so, we present the first empirical evidence demonstrating the positive relationship between Faster‐Z Effect and measures of promiscuity, and therefore variance in male mating success. Our results from multiple lines of evidence indicate that selection is less effective on the Z chromosome, particularly in promiscuous species, and that Faster‐Z Evolution in birds is due primarily to genetic drift. Our results reveal the power of mating system and sexual selection in shaping broad patterns in genome evolution.

Positive Selection Underlies Faster-Z Evolution of Gene Expression in Birds

The elevated rate of evolution for genes on sex chromosomes compared with autosomes (Fast-X or Fast-Z evolution) can result either from positive selection in the heterogametic sex or from nonadaptive consequences of reduced relative effective population size. Recent work in birds suggests that Fast-Z of coding sequence is primarily due to relaxed purifying selection resulting from reduced relative effective population size. However, gene sequence and gene expression are often subject to distinct evolutionary pressures; therefore, we tested for Fast-Z in gene expression using next-generation RNA-sequencing data from multiple avian species. Similar to studies of Fast-Z in coding sequence, we recover clear signatures of Fast-Z in gene expression; however, in contrast to coding sequence, our data indicate that Fast-Z in expression is due to positive selection acting primarily in females. In the soma, where gene expression is highly correlated between the sexes, we detected Fast-Z in both sexes, although at a higher rate in females, suggesting that many positively selected expression changes in females are also expressed in males. In the gonad, where intersexual correlations in expression are much lower, we detected Fast-Z for female gene expression, but crucially, not males. This suggests that a large amount of expression variation is sex-specific in its effects within the gonad. Taken together, our results indicate that Fast-Z evolution of gene expression is the product of positive selection acting on recessive beneficial alleles in the heterogametic sex. More broadly, our analysis suggests that the adaptive potential of Z chromosome gene expression may be much greater than that of gene sequence, results which have important implications for the role of sex chromosomes in speciation and sexual selection.

INDEPENDENT STRATUM FORMATION ON THE AVIAN SEX CHROMOSOMES REVEALS INTER CHROMOSOMAL GENE CONVERSION AND PREDOMINANCE OF PURIFYING SELECTION ON THE W CHROMOSOME

We used a comparative approach spanning three species and 90 million years to study the evolutionary history of the avian sex chromosomes. Using whole transcriptomes, we assembled the largest cross‐species dataset of W‐linked coding content to date. Our results show that recombination suppression in large portions of the avian sex chromosomes has evolved independently, and that long‐term sex chromosome divergence is consistent with repeated and independent inversions spreading progressively to restrict recombination. In contrast, over short‐term periods we observe heterogeneous and locus‐specific divergence. We also uncover four instances of gene conversion between both highly diverged and recently evolved gametologs, suggesting a complex mosaic of recombination suppression across the sex chromosomes. Lastly, evidence from 16 gametologs reveal that the W chromosome is evolving with a significant contribution of purifying selection, consistent with previous findings that W‐linked genes play an important role in encoding sex‐specific fitness.

The plover neurotranscriptome assembly: transcriptomic analysis in an ecological model species without a reference genome

We assembled a de novo transcriptome of short‐read Illumina RNA‐Seq data generated from telencephalon and diencephalon tissue samples from the Kentish plover, Charadrius alexandrinus. This is a species of considerable interest in behavioural ecology for its highly variable mating system and parental behaviour, but it lacks genomic resources and is evolutionarily distant from the few available avian draft genome sequences. We assembled and identified over 21 000 transcript contigs with significant expression in our samples, showing high homology to exonic sequences in avian draft genomes. From these, we identified >31 000 high‐quality SNPs and > 2500 simple sequence repeats (SSRs). We also analysed expression patterns in our data to identify potential candidate genes related to differences in male and female behaviour, identifying over 200 nonoverlapping putative autosomal transcripts that show significant expression differences between males and females. Gene ontology analysis revealed that female‐biased transcripts were significantly enriched for cerebral functions related to learning, cognition and memory, and male‐biased transcripts were mostly enriched for terms related to neural function such as neuron projection and synapses. This data set provides one of the first de novo transcriptome assemblies from non‐normalized short‐read next‐generation data and outlines an effective strategy for measuring sequence and expression variability simultaneously without the aid of a reference genome.