Widespread Positive Selection Research Articles

Adaptive evolution in virulence effectors of the rice blast fungus Pyricularia oryzae.

Plant pathogens secrete proteins called effectors that target host cellular processes to promote disease. Recently, structural genomics has identified several families of fungal effectors that share a similar three-dimensional structure despite remarkably variable amino-acid sequences and surface properties. To explore the selective forces that underlie the sequence variability of structurally-analogous effectors, we focused on MAX effectors, a structural family of effectors that are major determinants of virulence in the rice blast fungus Pyricularia oryzae. Using structure-informed gene annotation, we identified 58 to 78 MAX effector genes per genome in a set of 120 isolates representing seven host-associated lineages. The expression of MAX effector genes was primarily restricted to the early biotrophic phase of infection and strongly influenced by the host plant. Pangenome analyses of MAX effectors demonstrated extensive presence/absence polymorphism and identified gene loss events possibly involved in host range adaptation. However, gene knock-in experiments did not reveal a strong effect on virulence phenotypes suggesting that other evolutionary mechanisms are the main drivers of MAX effector losses. MAX effectors displayed high levels of standing variation and high rates of non-synonymous substitutions, pointing to widespread positive selection shaping the molecular diversity of MAX effectors. The combination of these analyses with structural data revealed that positive selection acts mostly on residues located in particular structural elements and at specific positions. By providing a comprehensive catalog of amino acid polymorphism, and by identifying the structural determinants of the sequence diversity, our work will inform future studies aimed at elucidating the function and mode of action of MAX effectors.

Genomic analysis of hypoxia inducible factor alpha in ray-finned fishes reveals missing Ohnologs and evidence of widespread positive selection

As aquatic hypoxia worsens on a global scale, fishes will become increasingly challenged by low oxygen, and understanding the molecular basis of their response to hypoxia may help to better define the capacity of fishes to cope with this challenge. The hypoxia inducible factor (HIF) plays a critical role in the molecular response to hypoxia by activating the transcription of genes that serve to improve oxygen delivery to the tissues or enhance the capacity of tissues to function at low oxygen. The current study examines the molecular evolution of genes encoding the oxygen-dependent HIFα subunit (HIFA) in the ray-finned fishes (Actinopterygii). Genomic analyses demonstrate that several lineages retain four paralogs of HIFA predicted from two rounds of genome duplication at the base of vertebrate evolution, broaden the known distribution of teleost-specific HIFA paralogs, and provide evidence for salmonid-specific HIFA duplicates. Evolution of the HIFA gene family is characterized by widespread episodic positive selection at amino acid sites that potentially mediate protein stability, protein–protein interactions, and transcriptional regulation. HIFA transcript abundance depends upon paralog, tissue, and fish lineage. A phylogenetically-informed gene nomenclature is proposed along with avenues for future research on this critical family of transcription factors.

No evidence for widespread positive selection on double substitutions within codons in primates and yeasts.

Nucleotide substitutions in protein-coding genes can be divided into synonymous (S) and non-synonymous (N) ones that alter amino acids (including nonsense mutations causing stop codons). The S substitutions are expected to have little effect on function. The N substitutions almost always are affected by strong purifying selection that eliminates them from evolving populations. However, additional mutations of nearby bases can modulate the deleterious effect of single N substitutions and, thus, could be subjected to the positive selection. This effect has been demonstrated for mutations in the serine codons, stop codons and double N substitutions in prokaryotes. In all abovementioned cases, a novel technique was applied that allows elucidating the effects of selection on double substitutions considering mutational biases. Here, we applied the same technique to study double N substitutions in eukaryotic lineages of primates and yeast. We identified markedly fewer cases of purifying selection relative to prokaryotes and no evidence of codon double substitutions under positive selection. This is consistent with previous studies of serine codons in primates and yeast. In general, the obtained results strongly suggest that there are major differences between studied pro- and eukaryotes; double substitutions in primates and yeasts largely reflect mutational biases and are not hallmarks of selection. This is especially important in the context of detection of positive selection in codons because it has been suggested that multiple mutations in codons cause false inferences of lineage-specific site positive selection. It is likely that this concern is applicable to previously studied prokaryotes but not to primates and yeasts where markedly fewer double substitutions are affected by positive selection.

Vulture Genomes Reveal Molecular Adaptations Underlying Obligate Scavenging and Low Levels of Genetic Diversity.

Obligate scavenging on the dead and decaying animal matter is a rare dietary specialization that in extant vertebrates is restricted to vultures. These birds perform essential ecological services, yet many vulture species have undergone recent steep population declines and are now endangered. To test for molecular adaptations underlying obligate scavenging in vultures, and to assess whether genomic features might have contributed to their population declines, we generated high-quality genomes of the Himalayan and bearded vultures, representing both independent origins of scavenging within the Accipitridae, alongside a sister taxon, the upland buzzard. By comparing our data to published sequences from other birds, we show that the evolution of obligate scavenging in vultures has been accompanied by widespread positive selection acting on genes underlying gastric acid production, and immunity. Moreover, we find evidence of parallel molecular evolution, with amino acid replacements shared among divergent lineages of these scavengers. Our genome-wide screens also reveal that both the Himalayan and bearded vultures exhibit low levels of genetic diversity, equating to around a half of the mean genetic diversity of other bird genomes examined. However, demographic reconstructions indicate that population declines began at around the Last Glacial Maximum, predating the well-documented dramatic declines of the past three decades. Taken together, our genomic analyses imply that vultures harbor unique adaptations for processing carrion, but that modern populations are genetically depauperate and thus especially vulnerable to further genetic erosion through anthropogenic activities.

No Evidence for Widespread Positive Selection Signatures in Common Risk Alleles Associated with Schizophrenia.

Schizophrenia poses an evolutionary-genetic paradox as it exhibits strongly negative fitness effects (early mortality and decreased fecundity), yet it persists at a prevalence of approximately 1% worldwide. Evidence from several studies have suggested that schizophrenia is evolved and maintained in part as a maladaptive byproduct of recent positive selection and adaptive evolution in human beings. However, inconsistent results have been also proposed, challenging the recent positive selection theory to explain the high population frequency of schizophrenia-associated alleles. Here, we used public domain data to locate signatures of positive selection based on genetic diversity, derived allele frequency, differentiation between populations, and long haplotypes at schizophrenia-associated single nucleotide polymorphisms (SNPs) and randomly selected SNPs (as negative controls). We found evidence for positive selection at 10 out of the 105 schizophrenia-associated SNPs, while 5 of these SNPs involved positive selection for the protective allele. Taken together, the absence of widespread positive selection signals at the schizophrenia-associated SNPs, along with the fact that half of the positive selection favored the protective allele, provide little evidence supporting the positive selection theory in schizophrenia.

Gene Expansion and Positive Selection as Bacterial Adaptations to Oligotrophic Conditions

We examined the genomic adaptations of prevalent bacterial taxa in a highly nutrient- and ion-depleted freshwater environment located in the secondary cooling water system of a nuclear research reactor. Using genome-centric metagenomics, we found that none of the prevalent bacterial taxa were related to typical freshwater bacterial lineages. We also did not identify strong signatures of genome streamlining, which has been shown to be one of the ecoevolutionary forces shaping the genome characteristics of bacterial taxa in nutrient-depleted environments. Instead, focusing on the dominant taxon, a novel Ramlibacter sp. which we propose to name Ramlibacter aquaticus, we detected extensive positive selection on genes involved in phosphorus and carbon scavenging pathways. These genes were involved in the high-affinity phosphate uptake and storage into polyphosphate granules, metabolism of nitrogen-rich organic matter, and carbon/energy storage into polyhydroxyalkanoate. In parallel, comparative genomics revealed a high number of paralogs and an accessory genome significantly enriched in environmental sensing pathways (i.e., chemotaxis and motility), suggesting extensive gene expansions in R. aquaticus The type strain of R. aquaticus (LMG 30558T) displayed optimal growth kinetics and productivity at low nutrient concentrations, as well as substantial cell size plasticity. Our findings with R. aquaticus LMG 30558T demonstrate that positive selection and gene expansions may represent successful adaptive strategies to oligotrophic environments that preserve high growth rates and cellular productivity.IMPORTANCE By combining a genome-centric metagenomic approach with a culture-based approach, we investigated the genomic adaptations of prevalent populations in an engineered oligotrophic freshwater system. We found evidence for widespread positive selection on genes involved in phosphorus and carbon scavenging pathways and for gene expansions in motility and environmental sensing to be important genomic adaptations of the abundant taxon in this system. In addition, microscopic and flow cytometric analysis of the first freshwater representative of this population (Ramlibacter aquaticus LMG 30558T) demonstrated phenotypic plasticity, possibly due to the metabolic versatility granted by its larger genome, to be a strategy to cope with nutrient limitation. Our study clearly demonstrates the need for the use of a broad set of genomic tools combined with culture-based physiological characterization assays to investigate and validate genomic adaptations.

Widespread positive selection on cetacean TLR extracellular domain

Toll like receptors (TLRs), key members of innate immune system, can recognize a wide diversity of pathogens and initiate both innate and adaptive immune responses in vertebrate. Cetaceans must have faced new challenges of pathogens when their terrestrial relatives transitioned from the terrestrial to aquatic environment. Here, we sequenced the extracellular domain (ECD) of 10 TLRs in cetacean lineages because this region involved in the recognition of pathogens. A total of 148 sites ranging between 5–26 codons (0.01%–4.83%) were identified to be robust candidates of positive selection at the ECD of 10 TLRs. In addition, the majority (90.54%) of these positively selected codons were found to have radical amino acid changes, which strengthen the evidence of positive selection. Importantly, more radical amino acid changes in selected sites were enriched in the period of early evolutionary transition from land to semi-aquatic and from semi-aquatic to full-aquatic habitat, which might endow cetaceans with a faster adaptation to new pathogens as they transitioned into novel habitat. Interestingly, similar selective intensity was detected in both viral and non-viral TLRs in cetaceans, which was not in line with previous studies on primates and birds that reported stronger positive selection in non-viral TLRs than in viral TLRs. This result may be explained by the fact that cetaceans might have faced diversity of bacteria and viruses during its transitions from terrestrial to aquatic environment whereas both primates and birds probably being affected by only a restricted number of related viruses due to their homogeneous habitat.

Functional and evolutionary genomic inferences in Populus through genome and population sequencing of American and European aspen

The Populus genus is one of the major plant model systems, but genomic resources have thus far primarily been available for poplar species, and primarily Populus trichocarpa (Torr. & Gray), which was the first tree with a whole-genome assembly. To further advance evolutionary and functional genomic analyses in Populus, we produced genome assemblies and population genetics resources of two aspen species, Populus tremula L. and Populus tremuloides Michx. The two aspen species have distributions spanning the Northern Hemisphere, where they are keystone species supporting a wide variety of dependent communities and produce a diverse array of secondary metabolites. Our analyses show that the two aspens share a similar genome structure and a highly conserved gene content with P. trichocarpa but display substantially higher levels of heterozygosity. Based on population resequencing data, we observed widespread positive and negative selection acting on both coding and noncoding regions. Furthermore, patterns of genetic diversity and molecular evolution in aspen are influenced by a number of features, such as expression level, coexpression network connectivity, and regulatory variation. To maximize the community utility of these resources, we have integrated all presented data within the PopGenIE web resource (PopGenIE.org).

PO-329 Rejecting neutral tumour evolution across human cancer types

IntroductionClusters of subclonal mutations attributed to subclonal expansions have been observed across most cancer types. To date it is not clear what fraction of these expansions can be ascribed to selective sweeps or genetic drifts. The widely held theory states that tumour expansion arises from a series of selective sweeps. However, Williams et al. (Nat. Genet. 48:238–224, 2016) recently claimed to have identified neutral tumour evolution in one third of 904 samples from The Cancer Genome Atlas (TCGA).Material and methodsFirst, using the same equations we simulated data of tumour undergoing selective sweeps. Second, we simulated data using branching processes, i.e. more realistic models. Third, we analysed the mutations from ~1600 exomes and ~1500 whole genomes of real tumours from TCGA and the International Cancer Genome Consortium, respectively. We grouped them into neutral and non-neutral tumours according to the rationale described by Williams et al. We then applied dN/dS, an orthogonal widely used approach to detect selection in non-synonymous mutations of coding regions. Finally, we ran dN/dS on subclonal mutations within individual cancer types.Results and discussionsApplication of the method described by Williams et al. led to over-calling of neutrality and led to random predictions (area under the curve~=50%) in simulations based on the deterministic or stochastic models, respectively. dN/dS analyses revealed significant and widespread positive selection in clonal and subclonal mutations within both neutral and non-neutral groups. Despite power limits, the subclonal mutations of many individual cancer types displayed dN/dS>1, evidence for positive selection.ConclusionAltogether, our results point to the lack of evidence for neutral tumour evolution and uncovered strong positive selection within subclonal mutations across cancers.

Evolutionary rates of mammalian telomere-stability genes correlate with karyotype features and female germline expression.

Telomeres protect the ends of eukaryotic chromosomes and are essential for cell viability. In mammals, telomere dynamics vary with life history traits (e.g. body mass and longevity), suggesting differential selection depending on physiological characteristics. Telomeres, in analogy to centromeric regions, also represent candidate meiotic drivers and subtelomeric DNA evolves rapidly. We analyzed the evolutionary history of mammalian genes implicated in telomere homeostasis (TEL genes). We detected widespread positive selection and we tested two alternative hypotheses: (i) fast evolution is driven by changes in life history traits; (ii) a conflict with selfish DNA elements at the female meiosis represents the underlying selective pressure. By accounting for the phylogenetic relationships among mammalian species, we show that life history traits do not contribute to shape diversity of TEL genes. Conversely, the evolutionary rate of TEL genes correlates with expression levels during meiosis and episodes of positive selection across mammalian species are associated with karyotype features (number of chromosome arms). We thus propose a telomere drive hypothesis, whereby (sub)telomeres and telomere-binding proteins are engaged in an intra-genomic conflict similar to the one described for centromeres.

Extensive Positive Selection Drives the Evolution of Nonstructural Proteins in Lineage C Betacoronaviruses.

Middle East respiratory syndrome-related coronavirus (MERS-CoV) spreads to humans via zoonotic transmission from camels. MERS-CoV belongs to lineage C of betacoronaviruses (betaCoVs), which also includes viruses isolated from bats and hedgehogs. A large portion of the betaCoV genome consists of two open reading frames (ORF1a and ORF1b) that are translated into polyproteins. These are cleaved by viral proteases to generate 16 nonstructural proteins (nsp1 to nsp16) which compose the viral replication-transcription complex. We investigated the evolution of ORF1a and ORF1b in lineage C betaCoVs. Results indicated widespread positive selection, acting mostly on ORF1a. The proportion of positively selected sites in ORF1a was much higher than that previously reported for the surface-exposed spike protein. Selected sites were unevenly distributed, with nsp3 representing the preferential target. Several pairs of coevolving sites were also detected, possibly indicating epistatic interactions; most of these were located in nsp3. Adaptive evolution at nsp3 is ongoing in MERS-CoV strains, and two selected sites (G720 and R911) were detected in the protease domain. While position 720 is variable in camel-derived viruses, suggesting that the selective event does not represent a specific adaptation to humans, the R911C substitution was observed only in human-derived MERS-CoV isolates, including the viral strain responsible for the recent South Korean outbreak. It will be extremely important to assess whether these changes affect host range or other viral phenotypes. More generally, data herein indicate that CoV nsp3 represents a major selection target and that nsp3 sequencing should be envisaged in monitoring programs and field surveys. Both severe acute respiratory syndrome coronavirus (SARS-CoV) and MERS-CoV originated in bats and spread to humans via an intermediate host. This clearly highlights the potential for coronavirus host shifting and the relevance of understanding the molecular events underlying the adaptation to new host species. We investigated the evolution of ORF1a and ORF1b in lineage C betaCoVs and in 87 sequenced MERS-CoV isolates. Results indicated widespread positive selection, stronger in ORF1a than in ORF1b. Several selected sites were found to be located in functionally relevant protein regions, and some of them corresponded to functional mutations in other coronaviruses. The proportion of selected sites we identified in ORF1a is much higher than that for the surface-exposed spike protein. This observation suggests that adaptive evolution in ORF1a might contribute to host shifts or immune evasion. Data herein also indicate that genetic diversity at nonstructural proteins should be taken into account when antiviral compounds are developed.

Widespread Positive Selection Drives Differentiation of Centromeric Proteins in the Drosophila melanogaster subgroup.

Rapid evolution of centromeric satellite repeats is thought to cause compensatory amino acid evolution in interacting centromere-associated kinetochore proteins. Cid, a protein that mediates kinetochore/centromere interactions, displays particularly high amino acid turnover. Rapid evolution of both Cid and centromeric satellite repeats led us to hypothesize that the apparent compensatory evolution may extend to interacting partners in the Condensin I complex (i.e., SMC2, SMC4, Cap-H, Cap-D2, and Cap-G) and HP1s. Missense mutations in these proteins often result in improper centromere formation and aberrant chromosome segregation, thus selection for maintained function and coevolution among proteins of the complex is likely strong. Here, we report evidence of rapid evolution and recurrent positive selection in seven centromere-associated proteins in species of the Drosophila melanogaster subgroup, and further postulate that positive selection on these proteins could be a result of centromere drive and compensatory changes, with kinetochore proteins competing for optimal spindle attachment.

OASes and STING: adaptive evolution in concert.

OAS (2′–5′-oligoadenylate synthases) proteins and cyclic GMP–AMP synthase (cGAS, gene symbol: MB21D1) patrol the cytoplasm for the presence of foreign nucleic acids. Upon binding to double-stranded RNA or double-stranded DNA, OAS proteins and cGAS produce nucleotide second messengers to activate RNase L and STING (stimulator of interferon genes, gene symbol: TMEM173), respectively; this leads to the initiation of antiviral responses. We analyzed the evolutionary history of the MB21D1–TMEM173 and OAS–RNASEL axes in primates and bats and found evidence of widespread positive selection in both orders. In TMEM173, residue 230, a major determinant of response to natural ligands and to mimetic drugs (e.g., DMXAA), was positively selected in Primates and Chiroptera. In both orders, selection also targeted an α-helix/loop element in RNase L that modulates the enzyme preference for single-stranded RNA versus stem loops. Analysis of positively selected sites in OAS1, OAS2, and MB21D1 revealed parallel evolution, with the corresponding residues being selected in different genes. As this cannot result from gene conversion, these data suggest that selective pressure acting on OAS and MB21D1 genes is related to nucleic acid recognition and to the specific mechanism of enzyme activation, which requires a conformational change. Finally, a population genetics-phylogenetics analysis in humans, chimpanzees, and gorillas detected several positively selected sites in most genes. Data herein shed light into species-specific differences in infection susceptibility and in response to synthetic compounds, with relevance for the design of synthetic compounds as vaccine adjuvants.

Evidence for widespread positive and negative selection in coding and conserved noncoding regions of Capsella grandiflora.

The extent that both positive and negative selection vary across different portions of plant genomes remains poorly understood. Here, we sequence whole genomes of 13 Capsella grandiflora individuals and quantify the amount of selection across the genome. Using an estimate of the distribution of fitness effects, we show that selection is strong in coding regions, but weak in most noncoding regions, with the exception of 5′ and 3′ untranslated regions (UTRs). However, estimates of selection on noncoding regions conserved across the Brassicaceae family show strong signals of selection. Additionally, we see reductions in neutral diversity around functional substitutions in both coding and conserved noncoding regions, indicating recent selective sweeps at these sites. Finally, using expression data from leaf tissue we show that genes that are more highly expressed experience stronger negative selection but comparable levels of positive selection to lowly expressed genes. Overall, we observe widespread positive and negative selection in coding and regulatory regions, but our results also suggest that both positive and negative selection on plant noncoding sequence are considerably rarer than in animal genomes.

Molecular adaptation during a rapid adaptive radiation.

“Explosive” adaptive radiations on islands remain one of the most puzzling evolutionary phenomena and the evolutionary genetic processes behind such radiations remain unclear. Rapid morphological and ecological evolution during island radiations suggests that many genes may be under fairly strong selection, although this remains untested. Here, we report that during a rapid recent diversification in the Hawaiian endemic plant genus Schiedea (Caryophyllaceae), 5 in 36 studied genes evolved under positive selection. Positively selected genes are involved in defence mechanisms, photosynthesis, and reproduction. Comparison with eight mainland plant groups demonstrates both the relaxation of purifying selection and more widespread positive selection in Hawaiian Schiedea. This provides compelling evidence that adaptive evolution of protein-coding genes may play a significant role during island adaptive radiations.

Evidence for Widespread Positive and Purifying Selection Across the European Rabbit (Oryctolagus cuniculus) Genome

The nearly neutral theory of molecular evolution predicts that the efficacy of both positive and purifying selection is a function of the long-term effective population size (N(e)) of a species. Under this theory, the efficacy of natural selection should increase with N(e). Here, we tested this simple prediction by surveying ~1.5 to 1.8 Mb of protein coding sequence in the two subspecies of the European rabbit (Oryctolagus cuniculus algirus and O. c. cuniculus), a mammal species characterized by high levels of nucleotide diversity and N(e) estimates for each subspecies on the order of 1 × 10(6). When the segregation of slightly deleterious mutations and demographic effects were taken into account, we inferred that >60% of amino acid substitutions on the autosomes were driven to fixation by positive selection. Moreover, we inferred that a small fraction of new amino acid mutations (<4%) are effectively neutral (defined as 0 < N(e)s < 1) and that this fraction was negatively correlated with a gene's expression level. Consistent with models of recurrent adaptive evolution, we detected a negative correlation between levels of synonymous site polymorphism and the rate of protein evolution, although the correlation was weak and nonsignificant. No systematic X chromosome-autosome difference was found in the efficacy of selection. For example, the proportion of adaptive substitutions was significantly higher on the X chromosome compared with the autosomes in O. c. algirus but not in O. c. cuniculus. Our findings support widespread positive and purifying selection in rabbits and add to a growing list of examples suggesting that differences in N(e) among taxa play a substantial role in determining rates and patterns of protein evolution.

A Hierarchical Bayesian Model for Next-Generation Population Genomics

The demography of populations and natural selection shape genetic variation across the genome and understanding the genomic consequences of these evolutionary processes is a fundamental aim of population genetics. We have developed a hierarchical Bayesian model to quantify genome-wide population structure and identify candidate genetic regions affected by selection. This model improves on existing methods by accounting for stochastic sampling of sequences inherent in next-generation sequencing (with pooled or indexed individual samples) and by incorporating genetic distances among haplotypes in measures of genetic differentiation. Using simulations we demonstrate that this model has a low false-positive rate for classifying neutral genetic regions as selected genes (i.e., Φ(ST) outliers), but can detect recent selective sweeps, particularly when genetic regions in multiple populations are affected by selection. Nonetheless, selection affecting just a single population was difficult to detect and resulted in a high false-negative rate under certain conditions. We applied the Bayesian model to two large sets of human population genetic data. We found evidence of widespread positive and balancing selection among worldwide human populations, including many genetic regions previously thought to be under selection. Additionally, we identified novel candidate genes for selection, several of which have been linked to human diseases. This model will facilitate the population genetic analysis of a wide range of organisms on the basis of next-generation sequence data.

A Newly Classified Vertebrate Calpain Protease, Directly Ancestral to CAPN1 and 2, Episodically Evolved a Restricted Physiological Function in Placental Mammals

The most studied members of the calpain protease superfamily are CAPN1 and 2, which are conserved across vertebrates. Another similar family member called mu/m-CAPN has been identified in birds alone. Here, we establish that mu/m-CAPN shares one-to-one orthology with CAPN11, previously described only in eutherians (placental mammals). We use the name CAPN11 for this family member and identify orthologues across vertebrate lineages, which form a monophyletic phylogenetic clade directly ancestral to CAPN1 and 2. In lineages branching before therians (live-bearing mammals), the CAPN11 coding region has evolved under strong purifying selection, with low nonsynonymous (d(N)) versus synonymous (d(S)) substitution rates (d(N)/d(S) = 0.076 across pretherians), and its transcripts were detected widely across different tissues. These characteristics are present in CAPN1 and 2 across vertebrate lineages and indicate that pretherian CAPN11 likewise has conserved a wide physiological function. However, an approximately 7-fold elevation in d(N)/d(S) is evident along the CAPN11 branch splitting eutherians from platypus, paralleled by a shift to "testis-specific" gene regulation. Estimates of d(N)/d(S) in eutherians were approximately 3-fold elevated compared with pretherians and coding and transcriptional-level evidence suggests that CAPN11 is functionally absent in marsupials. Many CAPN11 sites are functionally constrained in eutherians to conserve a residue with radically different biochemical properties to a fixed state shared between pretherian CAPN11 and CAPN1 and 2. Protein homology modeling demonstrated that many such eutherian-specific residue replacements modify or ablate interactions with the calpain inhibitor calpastatin that are observed in both pretherian orthologues and CAPN1/2. We propose a model akin to the Dykhuizen-Hartl effect, where inefficient purifying selection and increased genetic drift associated with a reduction in effective population size, drove the fixation of mutations in regulatory and coding regions of CAPN11 of a common marsupial-eutherian ancestor. A subset of these changes had a cumulative adaptive advantage in a eutherian ancestor because of lineage-specific aspects of sperm physiology, whereas in marsupials, no advantage was realized and the gene was disabled. This work supports that functional divergence among gene family member orthologues is possible in the absence of widespread positive selection.

Enrichment of HapMap recombination hotspot predictions around human nervous system genes: evidence for positive selection ?

Channels and developmental genes belong to the molecular key players in the human central nervous system (CNS). Mutations in these genes often cause monogenic neurological disease and interspecies comparisons had shown reduced divergence. On the other hand, accelerated evolution of genes with roles in neurotransmission and development had indicated widespread positive selection in hominids. In the present study, we hypothesized that recombination hotspots could be enriched at genes with particularly important role in the CNS, because at those loci beneficial mutations may occur on a highly constrained background and consequently increased recombination could promote their fixation. To test this hypothesis, we retrieved CNS genes based on keyword search, expression data and expert knowledge. Consistent with our hypothesis, we find an enrichment of hotspot predictions around genes that are retrieved by all three strategies. Moreover, when comparing human genes based on their Gene Ontology annotations, we find hotspot predictions preferentially located around channels and neurodevelopmental genes. Taken together with the distinct sequence evolution that was reported by comparative genomic studies, this finding indicates continued positive selection at many CNS gene loci. In support of this interpretation, we also find an enrichment of recombination hotspot predictions around conserved noncoding regions that were reported to display a signature of accelerated evolution in the human lineage. Widespread positive selection acting on CNS gene loci could relate to the high prevalence of human nervous system disorders with genetically complex inheritance, potentially under an ancestral susceptibility allele model.

Widespread Positive Selection in Synonymous Sites of Mammalian Genes

Evolution of protein sequences is largely governed by purifying selection, with a small fraction of proteins evolving under positive selection. The evolution at synonymous positions in protein-coding genes is not nearly as well understood, with the extent and types of selection remaining, largely, unclear. A statistical test to identify purifying and positive selection at synonymous sites in protein-coding genes was developed. The method compares the rate of evolution at synonymous sites (Ks) to that in intron sequences of the same gene after sampling the aligned intron sequences to mimic the statistical properties of coding sequences. We detected purifying selection at synonymous sites in approximately 28% of the 1,562 analyzed orthologous genes from mouse and rat, and positive selection in approximately 12% of the genes. Thus, the fraction of genes with readily detectable positive selection at synonymous sites is much greater than the fraction of genes with comparable positive selection at nonsynonymous sites, i.e., at the level of the protein sequence. Unlike other genes, the genes with positive selection at synonymous sites showed no correlation between Ks and the rate of evolution in nonsynonymous sites (Ka), indicating that evolution of synonymous sites under positive selection is decoupled from protein evolution. The genes with purifying selection at synonymous sites showed significant anticorrelation between Ks and expression level and breadth, indicating that highly expressed genes evolve slowly. The genes with positive selection at synonymous sites showed the opposite trend, i.e., highly expressed genes had, on average, higher Ks. For the genes with positive selection at synonymous sites, a significantly lower mRNA stability is predicted compared to the genes with negative selection. Thus, mRNA destabilization could be an important factor driving positive selection in nonsynonymous sites, probably, through regulation of expression at the level of mRNA degradation and, possibly, also translation rate. So, unexpectedly, we found that positive selection at synonymous sites of mammalian genes is substantially more common than positive selection at the level of protein sequences. Positive selection at synonymous sites might act through mRNA destabilization affecting mRNA levels and translation.