Towards a more nuanced understanding of the relationship between sex-biased gene expression and rates of protein-coding sequence evolution

Genes that are differentially expressed between the sexes (sex-biased genes) are among the fastest evolving genes in animal genomes. The majority of sex-biased expression is attributable to genes that are primarily expressed in sex-limited reproductive tissues, and these reproductive genes are often rapidly evolving because of intra- and intersexual selection pressures. Additionally, studies of multiple taxa have revealed that genes with sex-biased expression are also expressed in a limited number of tissues. This is worth noting because narrowly expressed genes are known to evolve faster than broadly expressed genes. Therefore, it is not clear whether sex-biased genes are rapidly evolving because they have sexually dimorphic expression, because they are expressed in sex-limited reproductive tissues, or because they are narrowly expressed. To determine the extend to which other confounding variables can explain the rapid evolution of sex-biased genes, I analyzed the rates of evolution of sex-biased genes in Drosophila melanogaster and Mus musculus in light of tissue-specific measures of expression. I find that genes with sex-biased expression in somatic tissues shared by both sexes are often evolving faster than non-sex-biased genes, but this is best explained by the narrow expression profiles of sex-biased genes. Sex-biased genes in sex-limited tissues in D. melanogaster, however, evolve faster than other narrowly expressed genes. Therefore, the rapid evolution of sex-biased genes is limited only to those genes primarily expressed in sex-limited reproductive tissues.     

Comparative analysis of testis protein evolution in rodents

Genes expressed in testes are critical to male reproductive success, affecting spermatogenesis, sperm competition, and sperm-egg interaction. Comparing the evolution of testis proteins at different taxonomic levels can reveal which genes and functional classes are targets of natural and sexual selection and whether the same genes are targets among taxa. Here we examine the evolution of testis-expressed proteins at different levels of divergence among three rodents, mouse (Mus musculus), rat (Rattus norvegicus), and deer mouse (Peromyscus maniculatus), to identify rapidly evolving genes. Comparison of expressed sequence tags (ESTs) from testes suggests that proteins with testis-specific expression evolve more rapidly on average than proteins with maximal expression in other tissues. Genes with the highest rates of evolution have a variety of functional roles including signal transduction, DNA binding, and egg-sperm interaction. Most of these rapidly evolving genes have not been identified previously as targets of selection in comparisons among more divergent mammals. To determine if these genes are evolving rapidly among closely related species, we sequenced 11 of these genes in six Peromyscus species and found evidence for positive selection in five of them. Together, these results demonstrate rapid evolution of functionally diverse testis-expressed proteins in rodents, including the identification of amino acids under lineage-specific selection in Peromyscus. Evidence for positive selection among closely related species suggests that changes in these proteins may have consequences for reproductive isolation.     

Rates of protein evolution are positively correlated with developmental timing of expression during mouse spermatogenesis

Male reproductive genes often evolve very rapidly, and sexual selection is thought to be a primary force driving this divergence. We investigated the molecular evolution of 987 genes expressed at different times during mouse spermatogenesis to determine if the rate of evolution and the intensity of positive selection vary across stages of male gamete development. Using mouse-rat orthologs, we found that rates of protein evolution were positively correlated with the developmental timing of expression. Genes expressed early in spermatogenesis had rates of divergence similar to the genome median, while genes expressed after the onset of meiosis were found to evolve much more quickly. Rates of protein evolution were fastest for genes expressed during the dramatic morphogenesis of round spermatids into spermatozoa. Late-expressed genes were also more likely to be specific to the male germline. To test for evidence of positive selection, we analyzed the ratio of nonsynonymous to synonymous changes using a maximum likelihood framework in comparisons among mouse, rat, and human. Many genes showed evidence of positive selection, and most of these genes were expressed late in spermatogenesis and were testis specific. Overall, these data suggest that the intensity of positive selection associated with the evolution of male gametes varies considerably across development and acts primarily on phenotypes that develop late in spermatogenesis.     

Rates of Evolution of Hominoid Seminal Proteins are Correlated with Function and Expression, Rather than Mating System

In species where females mate with multiple males during a single ovulatory cycle, sperm competition is hypothesized to increase the rate of adaptive evolution of proteins expressed in male reproductive tissues through recurrent selective sweeps (positive selection). The hominoids, comprising apes and humans, are a group of closely related primates with extensive variation in mating behaviors and predicted levels of sperm competition. Since previous studies of individual male reproductive genes have shown evidence of positive selection, we estimated rates of evolution of a comprehensive set of proteins expressed in ejaculated semen. Our results show that these proteins in hominoids do not have elevated rates of nonsynonymous substitutions (Ka) compared with a control dataset of nonreproductive genes. Species with greater sperm competition do not have faster rates of seminal protein evolution. Although at these broad levels our hypotheses were not confirmed, further analyses indicate specific patterns of molecular evolution. Namely, the Ka of seminal genes is more strongly correlated with measures of tissue specificity than nonreproductive genes, suggesting that the former may more readily adapt to tissue-specific functions. Proteins expressed from the seminal vesicles evolve more rapidly than those from other male reproductive tissues. Also, several gene ontology categories show elevated rates of protein evolution, not seen in the control data set. While the generalization that male reproductive genes evolve rapidly in hominoids is an oversimplification, a subset of proteins can be identified that are likely targets for adaptive evolution driven by sexual selection.     

Molecular evolutionary analysis of seminal receptacle sperm storage organ genes of Drosophila melanogaster

Sperm storage organs are common and broadly distributed among animal taxa. However, little is known about how these organs function at the molecular level. Additionally, there is a paucity of knowledge about the evolution of genes expressed in these organs. This investigation is an evolutionary expressed sequence tag (EST) study of genes expressed in the seminal receptacle, one of the sperm storage organs in Drosophila. The incidence of positive selection is higher for the seminal receptacle genes than Drosophila reproductive genes as a whole, but lower than genes associated with the spermatheca, a second type of Drosophila sperm storage organ. By identifying overrepresented classes of proteins and classes for which sperm storage function is suggested by the nature of the proteins, candidate genes were discovered. These candidates belong to protein classes such as muscle contraction, odorant binding and odorant receptor, protease inhibitor and immunity.     

Sexual selection drives weak positive selection in protamine genes and high promoter divergence,enhancing sperm competitiveness

Phenotypic adaptations may be the result of changes in gene structure or gene regulation, but little is known about the evolution of gene expression. In addition, it is unclear whether the same selective forces may operate at both levels simultaneously. Reproductive proteins evolve rapidly, but the underlying selective forces promoting such rapid changes are still a matter of debate. In particular, the role of sexual selection in driving positive selection among reproductive proteins remains controversial, whereas its potential influence on changes in promoter regions has not been explored. Protamines are responsible for maintaining DNA in a compacted form in chromosomes in sperm and the available evidence suggests that they evolve rapidly. Because protamines condense DNA within the sperm nucleus, they influence sperm head shape. Here, we examine the influence of sperm competition upon protamine 1 and protamine 2 genes and their promoters, by comparing closely related species of Mus that differ in relative testes size, a reliable indicator of levels of sperm competition. We find evidence of positive selection in the protamine 2 gene in the species with the highest inferred levels of sperm competition. In addition, sperm competition levels across all species are strongly associated with high divergence in protamine 2 promoters that, in turn, are associated with sperm swimming speed. We suggest that changes in protamine 2 promoters are likely to enhance sperm swimming speed by making sperm heads more hydrodynamic. Such phenotypic changes are adaptive because sperm swimming speed may be a major determinant of fertilization success under sperm competition. Thus, when species have diverged recently, few changes in gene-coding sequences are found, while high divergence in promoters seems to be associated with the intensity of sexual selection.     

Olfactory receptor gene evolution is unusually rapid across Tetrapoda and outpaces chemosensory phenotypic change

Chemosensation is the most ubiquitous sense in animals, enacted by the products of complex gene families that detect environmental chemical cues and larger-scale sensory structures that process these cues. While there is a general conception that olfactory receptor (OR) genes evolve rapidly, the universality of this phenomenon across vertebrates, and its magnitude, are unclear. The supposed correlation between molecular rates of chemosensory evolution and phenotypic diversity of chemosensory systems is largely untested. We combine comparative genomics and sensory morphology to test whether OR genes and olfactory phenotypic traits evolve at faster rates than other genes or traits. Using published genomes, we identified ORs in 21 tetrapods, including amphibians, reptiles, birds, and mammals and compared their rates of evolution to those of orthologous non-OR protein-coding genes. We found that, for all clades investigated, most OR genes evolve nearly an order of magnitude faster than other protein-coding genes, with many OR genes showing signatures of diversifying selection across nearly all taxa in this study. This rapid rate of evolution suggests that chemoreceptor genes are in “evolutionary overdrive,” perhaps evolving in response to the ever-changing chemical space of the environment. To obtain complementary morphological data, we stained whole fixed specimens with iodine, µCT-scanned the specimens, and digitally segmented chemosensory and nonchemosensory brain regions. We then estimated phenotypic variation within traits and among tetrapods. While we found considerable variation in chemosensory structures, they were no more diverse than nonchemosensory regions. We suggest chemoreceptor genes evolve quickly in reflection of an ever-changing chemical space, whereas chemosensory phenotypes and processing regions are more conserved because they use a standardized or constrained architecture to receive and process a range of chemical cues.     

Conflict theory of genomic imprinting in mammals

In mammals, some embryonic genes are expressed differently depending on whether they are inherited from the sperm or egg, a phenomenon known as genomic imprinting. The information on the parental origin is transmitted by an epigenetic mark. Both the molecular mechanisms and evolutionary processes of genomic imprinting have been studied extensively. Here, I illustrate the simplest evolutionary dynamics of imprinting evolution based on the “conflict theory,” by considering the evolution of a gene encoding an embryonic growth factor controlling the maternal resource supply. It demonstrates that (a) the autosomal genes controlling placenta development to modify maternal resource acquisition may evolve a strong asymmetry of gene expression, provided the mother has some chance of accepting multiple males. (b) The genomic imprinting may not evolve if there is a small fraction of recessive deleterious mutations on the gene. (c) The growth-enhancing genes should evolve to paternally expressed, while the growth-suppressing genes should evolve to maternally expressed. (d) The X-linked genes also evolve genomic imprinting, but the main evolutionary force is the sex difference in the optimal embryonic size. I discuss other aberrations that can be explained by the modified versions of the basic model.     

Genomic and transcriptomic analyses reveal selection of genes for puberty in Bama Xiang pigs

The Bama Xiang pig(BMX) is a famous early-maturing Chinese indigenous breed with a two-end black coat.To uncover the genetic basis of the BMX phenotype,we conducted comparative genomic analyses between BMX and East Asian wild boars and Laiwu pigs,respectively. Genes under positive selection were enriched in pathways associated with gonadal hormone and melanin synthesis, consistent with the phenotypic changes observed during development in BMX pigs. We also performed differentially expressed gene analysis based on RNA-seq data from pituitary tissues of BMX and Large White pigs. The CTTNBP2 NL, FRS2,KANK4, and KATNAL1 genes were under selection and exhibited expressional changes in the pituitary tissue, which may affect BMX pig puberty. Our study demonstrated the positive selection of early maturity in the development of BMX pigs and advances our knowledge on the role of regulatory elements in puberty evolution in pigs.     

Mammalian sperm proteins are rapidly evolving: evidence of positive selection in functionally diverse genes

A growing number of genes involved in sex and reproduction have been demonstrated to be rapidly evolving. Here, we show that genes expressed solely in spermatozoa represent a highly diverged subset among mouse and human tissue-specific orthologs. The average rate of nonsynonymous substitutions per site (K(a)) is significantly higher in sperm proteins (mean K(a) = 0.18; N = 35) than in proteins expressed specifically in all other tissues (mean K(a) = 0.074; N = 473). No differences, however, are found in the synonymous substitution rate (K(s)) between tissues, suggesting that selective forces, and not mutation rate, explain the high rate of replacement substitutions in sperm proteins. Four out of 19 sperm-specific genes with characterized function demonstrated evidence of strong positive Darwinian selection, including a protein involved in gene regulation, Protamine-1 (PRM1), a protein involved in glycolysis, GAPDS, and two egg-binding proteins, Adam-2 precursor (ADAM2) and sperm-adhesion molecule-1 (SAM1). These results demonstrate the rapid evolution of sperm-specific genes and highlight the molecular action of sexual selection on a variety of characters involved in mammalian sperm function.     

The evolution of sperm competition genes: The effect of mating system on levels of genetic variation within and between species

It is widely established that proteins involved in reproduction diverge between species more quickly than other proteins. For male sperm proteins, rapid divergence is believed to be caused by postcopulatory sexual selection and/or sexual conflict. Here, we derive the expected levels of gene diversity within populations and divergence between them for male sperm protein genes evolving by postcopulatory, prezygotic fertility competition, i.e. the function imputed for some sperm and seminal fluid genes. We find that, at the mutation‐selection equilibrium, both gene diversity within species and divergence between them are elevated relative to genes with similar selection coefficients expressed by both sexes. We show that their expected level of diversity is a function of the harmonic mean number of mates per female, which affects the strength of fertility selection stemming from male–male sperm competition. Our predictions provide a null hypothesis for distinguishing between other selective hypotheses accounting for the rapid evolution of male reproductive genes.     

Maternal effect genes and the evolution of sociality in haplo-diploid organisms

Maternal care and female-biased sex ratios are considered by many to be essential prerequisites for the evolution of eusocial behaviors among the hymenoptera. Using population genetic models, I investigate the evolution of genes that have positive maternal effects but negative, direct effects on offspring fitness. I find that, under many conditions, such genes evolve more easily in haplo-diploids than in diplo-diploids. In fact, the conditions are less restrictive than those of kin selection theory, which postulate genes with negative direct effects but positive sib-social effects. For example, the conditions permitting the evolution of maternal effect genes are not affected if females mate multiply, whereas multiple mating reduces the efficacy of kin selection by reducing genetic relatedness within colonies. Inbreeding also differentially facilitates evolution of maternal effect genes in haplo-diploids relative to diplo-diploids, although it does not differentially affect the evolution of sib-altruism genes. Furthermore, when the direct, deleterious pleiotropic effect is restricted to sons, a maternal effect gene can evolve when the beneficial maternal effect is less than half (with inbreeding, much less) of the deleterious effect on sons. For kin selection, however, the sib-social benefits must always exceed the direct costs because genetic relatedness is always less than or equal to 1.0. The results suggest that haplo-diploidy facilitates (1) the evolution of maternal care, and (2) the evolution of maternal effect genes with antagonistic pleiotropic effects on sons. The latter effect may help explain the tendency toward female-biased sex ratios in haplo-diploids, especially those with inbreeding. I conclude that haplo-diploidy not only facilitates the evolution of sister-sister altruism by kin selection but also facilitates the evolution of maternal care and female-biased sex ratios, two prerequisites for eusociality.     

Relaxed genetic constraint is ancestral to the evolution of phenotypic plasticity

Phenotypic plasticity--the capacity of a single genotype to produce different phenotypes in response to varying environmental conditions--is widespread. Yet, whether, and how, plasticity impacts evolutionary diversification is unclear. According to a widely discussed hypothesis, plasticity promotes rapid evolution because genes expressed differentially across different environments (i.e., genes with "biased" expression) experience relaxed genetic constraint and thereby accumulate variation faster than do genes with unbiased expression. Indeed, empirical studies confirm that biased genes evolve faster than unbiased genes in the same genome. An alternative hypothesis holds, however, that the relaxed constraint and faster evolutionary rates of biased genes may be a precondition for, rather than a consequence of, plasticity's evolution. Here, we evaluated these alternative hypotheses by characterizing evolutionary rates of biased and unbiased genes in two species of frogs that exhibit a striking form of phenotypic plasticity. We also characterized orthologs of these genes in four species of frogs that had diverged from the two plastic species before the plasticity evolved. We found that the faster evolutionary rates of biased genes predated the evolution of the plasticity. Furthermore, biased genes showed greater expression variance than did unbiased genes, suggesting that they may be more dispensable. Phenotypic plasticity may therefore evolve when dispensable genes are co-opted for novel function in environmentally induced phenotypes. Thus, relaxed genetic constraint may be a cause--not a consequence--of the evolution of phenotypic plasticity, and thereby contribute to the evolution of novel traits.     

Sexual selection and the adaptive evolution of mammalian ejaculate proteins

An elevated rate of substitution characterizes the molecular evolution of reproductive proteins from a wide range of taxa. Although the selective pressures explaining this rapid evolution are yet to be resolved, recent evidence implicates sexual selection as a potentially important explanatory factor. To investigate this hypothesis, we sought evidence of a high rate of adaptive gene evolution linked to postcopulatory sexual selection in muroid rodents, a model vertebrate group displaying a broad range of mating systems. Specifically, we sequenced 7 genes from diverse rodents that are expressed in the testes, prostate, or seminal vesicles, products of which have the potential to act in sperm competition. We inferred positive Darwinian selection in these genes by estimation of the ratio of nonsynonymous (d(N), amino acid changing) to synonymous (d(S), amino acid retaining) substitution rates (omega = d(N)/d(S)). Next, we tested whether variation in this ratio among lineages could be attributed to interspecific variation in mating systems, as inferred from the variation in these rodents' relative testis sizes (RTS). Four of the 7 genes examined (Prm1, Sva, Acrv1, and Svs2, but not Svp2, Msmb, or Spink3) exhibit unambiguous evidence of positive selection. One of these, the seminal vesicle-derived protein Svs2, also shows some evidence for a concentration of positive selection in those lineages in which sperm competition is common. However, this was not a general trend among all the rodent genes we examined. Using the same methods, we then reanalyzed previously published data on 2 primate genes, SEMG1 and SEMG2. Although SEMG2 also shows evidence of positive selection concentrated in lineages subject to high levels of sperm competition, no such trend was found for SEMG1. Overall, despite a high rate of positive selection being a feature of many ejaculate proteins, these results indicate that the action of sexual selection potentially responsible for elevated evolutionary rates may be difficult to detect on a gene-by-gene basis. Although the extreme diversity of reproductive phenotypes exhibited in nature attests to the power of sexual selection, the extent to which this force predominates in driving the rapid molecular evolution of reproductive genes therefore remains to be determined.     

Genomic characterization of multiple clinical phenotypes of cancer using multivariate linear regression models

MOTIVATION: The development of gene expression microarray technology has allowed the identification of differentially expressed genes between different clinical phenotypic classes of cancer from a large pool of candidate genes. Although many class comparisons concerned only a single phenotype, simultaneous assessment of the relationship between gene expression and multiple phenotypes would be warranted to better understand the underlying biological structure. RESULTS: We develop a method to select genes related to multiple clinical phenotypes based on a set of multivariate linear regression models. For each gene, we perform model selection based on the doubly-adjusted R-square statistic and use the maximum of this statistic for gene selection. The method can substantially improve the power in gene selection, compared with a conventional method that uses a single model exclusively for gene selection. Application to a bladder cancer study to correlate pre-treatment gene expressions with pathological stage and grade is given. The methods would be useful for screening for genes related to multiple clinical phenotypes. AVAILABILITY: SAS and MATLAB codes are available from author upon request.