Detecting selection in noncoding regions of nucleotide sequences

KaKs_Calculator 3.0: Calculating Selective Pressure on Coding and Non-coding Sequences

KaKs_Calculator 3.0 is an updated toolkit that is capable of calculating selective pressure on both coding and non-coding sequences. Similar to the nonsynonymous/synonymous substitution rate ratio for coding sequences, selection on non-coding sequences can be quantified as the ratio of non-coding nucleotide substitution rate to synonymous substitution rate of adjacent coding sequences. As testified on empirical data, KaKs_Calculator 3.0 shows effectiveness to detect the strength and mode of selection operated on molecular sequences, accordingly demonstrating its great potential to achieve genome-wide scan of natural selection on diverse sequences and identification of potentially functional elements at a whole-genome scale. The package of KaKs_Calculator 3.0 is freely available for academic use only at https://ngdc.cncb.ac.cn/biocode/tools/BT000001.

Positive selection within the genomes of SARS-CoV-2 and other Coronaviruses independent of impact on protein function

BACKGROUND

The emergence of a novel coronavirus (SARS-CoV-2) associated with severe acute respiratory disease (COVID-19) has prompted efforts to understand the genetic basis for its unique characteristics and its jump from non-primate hosts to humans. Tests for positive selection can identify apparently nonrandom patterns of mutation accumulation within genomes, highlighting regions where molecular function may have changed during the origin of a species. Several recent studies of the SARS-CoV-2 genome have identified signals of conservation and positive selection within the gene encoding Spike protein based on the ratio of synonymous to nonsynonymous substitution. Such tests cannot, however, detect changes in the function of RNA molecules.

METHODS

Here we apply a test for branch-specific oversubstitution of mutations within narrow windows of the genome without reference to the genetic code.

RESULTS

We recapitulate the finding that the gene encoding Spike protein has been a target of both purifying and positive selection. In addition, we find other likely targets of positive selection within the genome of SARS-CoV-2, specifically within the genes encoding Nsp4 and Nsp16. Homology-directed modeling indicates no change in either Nsp4 or Nsp16 protein structure relative to the most recent common ancestor. These SARS-CoV-2-specific mutations may affect molecular processes mediated by the positive or negative RNA molecules, including transcription, translation, RNA stability, and evasion of the host innate immune system. Our results highlight the importance of considering mutations in viral genomes not only from the perspective of their impact on protein structure, but also how they may impact other molecular processes critical to the viral life cycle.

Selection against archaic hominin genetic variation in regulatory regions

Traces of Neandertal and Denisovan DNA persist in the modern human gene pool, but have been systematically purged by natural selection from genes and other functionally important regions. This implies that many archaic alleles harmed the fitness of hybrid individuals, but the nature of this harm is poorly understood. Here, we show that enhancers contain less Neandertal and Denisovan variation than expected given the background selection they experience, suggesting that selection acted to purge these regions of archaic alleles that disrupted their gene regulatory functions. We infer that selection acted mainly on young archaic variation that arose in Neandertals or Denisovans shortly before their contact with humans; enhancers are not depleted of older variants found in both archaic species. Some types of enhancer appear to have tolerated introgression better than others; compared with tissue-specific enhancers, pleiotropic enhancers show stronger depletion of archaic single-nucleotide polymorphisms. To some extent, evolutionary constraint is predictive of introgression depletion, but certain tissues' enhancers are more depleted of Neandertal and Denisovan alleles than expected given their comparative tolerance to new mutations. Foetal brain and muscle are the tissues whose enhancers show the strongest depletion of archaic alleles, but only brain enhancers show evidence of unusually stringent purifying selection. We conclude that epistatic incompatibilities between human and archaic alleles are needed to explain the degree of archaic variant depletion from foetal muscle enhancers, perhaps due to divergent selection for higher muscle mass in archaic hominins compared with humans.

Identifying branch-specific positive selection throughout the regulatory genome using an appropriate proxy neutral

BACKGROUND

Adaptive changes in cis-regulatory elements are an essential component of evolution by natural selection. Identifying adaptive and functional noncoding DNA elements throughout the genome is therefore crucial for understanding the relationship between phenotype and genotype.

RESULTS

We used ENCODE annotations to identify appropriate proxy neutral sequences and demonstrate that the conservativeness of the test can be modulated during the filtration of reference alignments. We applied the method to noncoding Human Accelerated Elements as well as open chromatin elements previously identified in 125 human tissues and cell lines to demonstrate its utility. Then, we evaluated the impact of query region length, proxy neutral sequence length, and branch count on test sensitivity and specificity. We found that the length of the query alignment can vary between 150 bp and 1 kb without affecting the estimation of selection, while for the reference alignment, we found that a length of 3 kb is adequate for proper testing. We also simulated sequence alignments under different classes of evolution and validated our ability to distinguish positive selection from relaxation of constraint and neutral evolution. Finally, we re-confirmed that a quarter of all non-coding Human Accelerated Elements are evolving by positive selection.

CONCLUSION

Here, we introduce a method we called adaptiPhy, which adds significant improvements to our earlier method that tests for branch-specific directional selection in noncoding sequences. The motivation for these improvements is to provide a more sensitive and better targeted characterization of directional selection and neutral evolution across the genome.

Evaluating Chromatin Accessibility Differences Across Multiple Primate Species Using a Joint Modeling Approach

Changes in transcriptional regulation are thought to be a major contributor to the evolution of phenotypic traits, but the contribution of changes in chromatin accessibility to the evolution of gene expression remains almost entirely unknown. To address this important gap in knowledge, we developed a new method to identify DNase I Hypersensitive (DHS) sites with differential chromatin accessibility between species using a joint modeling approach. Our method overcomes several limitations inherent to conventional threshold-based pairwise comparisons that become increasingly apparent as the number of species analyzed rises. Our approach employs a single quantitative test which is more sensitive than existing pairwise methods. To illustrate, we applied our joint approach to DHS sites in fibroblast cells from five primates (human, chimpanzee, gorilla, orangutan, and rhesus macaque). We identified 89,744 DHS sites, of which 41% are identified as differential between species using the joint model compared with 33% using the conventional pairwise approach. The joint model provides a principled approach to distinguishing single from multiple chromatin accessibility changes among species. We found that nondifferential DHS sites are enriched for nucleotide conservation. Differential DHS sites with decreased chromatin accessibility relative to rhesus macaque occur more commonly near transcription start sites (TSS), while those with increased chromatin accessibility occur more commonly distal to TSS. Further, differential DHS sites near TSS are less cell type-specific than more distal regulatory elements. Taken together, these results point to distinct classes of DHS sites, each with distinct characteristics of selection, genomic location, and cell type specificity.

Evolution of Subgenomic RNA Shapes Dengue Virus Adaptation and Epidemiological Fitness

Changes in dengue virus (DENV) genome affect viral fitness both clinically and epidemiologically. Even in the 3′ untranslated region (3′ UTR), mutations could affect subgenomic flaviviral RNA (sfRNA) production and its affinity for host proteins, which are necessary for successful viral replication. Indeed, we recently showed that mutations in DENV2 3′ UTR of epidemic strains increased sfRNA ability to bind host proteins and reduce interferon expression. However, whether 3′ UTR differences shape the overall DENV evolution remains incompletely understood. Herein, we combined RNA phylogeny with phylogenetics to gain insights on sfRNA evolution. We found that sfRNA structures are under purifying selection and highly conserved despite sequence divergence. Only the second flaviviral nuclease-resistant RNA (fNR2) structure of DENV2 sfRNA has undergone strong positive selection. Epidemiological reports suggest that substitutions in fNR2 may drive DENV2 epidemiological fitness, possibly through sfRNA-protein interactions. Collectively, our findings indicate that 3′ UTRs are important determinants of DENV fitness in human-mosquito cycles.

•

Dengue viruses (DENVs) preserve RNA elements in their 3′ untranslated region (UTR).

•

Quantification of natural selection revealed positive selection on DENV2 sfRNA

•

Flaviviral nuclease-resistant RNAs (fNR) in the 3′ UTRs contribute to DENV speciation

•

A highly evolving fNR structure appears to increase DENV2 epidemiological fitness

Parallel molecular evolution and adaptation in viruses

Parallel molecular evolution is the independent evolution of the same genotype or phenotype from distinct ancestors. The simple genomes and rapid evolution of many viruses mean they are useful model systems for studying parallel evolution by natural selection. Parallel adaptation occurs in the context of several viral behaviours, including cross-species transmission, drug resistance, and host immune escape, and its existence suggests that at least some aspects of virus evolution and emergence are repeatable and predictable. We introduce examples of virus parallel evolution and summarise key concepts. We outline the difficulties in detecting parallel adaptation using virus genomes, with a particular focus on phylogenetic and structural approaches, and we discuss future approaches that may improve our understanding of the phenomenon.

Selection Acting on Genomes

Populations evolve as mutations arise in individual organisms and, through hereditary transmission, may become "fixed" (shared by all individuals) in the population. Most mutations are lethal or have negative fitness consequences for the organism. Others have essentially no effect on organismal fitness and can become fixed through the neutral stochastic process known as random drift. However, mutations may also produce a selective advantage that boosts their chances of reaching fixation. Regions of genomes where new mutations are beneficial, rather than neutral or deleterious, tend to evolve more rapidly due to positive selection. Genes involved in immunity and defense are a well-known example; rapid evolution in these genes presumably occurs because new mutations help organisms to prevail in evolutionary "arms races" with pathogens. In recent years genome-wide scans for selection have enlarged our understanding of the genome evolution of various species. In this chapter, we will focus on methods to detect selection on the genome. In particular, we will discuss probabilistic models and how they have changed with the advent of new genome-wide data now available.

The diversity, evolution and epidemiology of plant viruses: A phylogenetic view

During the past four decades, the scientific community has seen an exponential advance in the number, sophistication, and quality of molecular techniques and bioinformatics tools for the genetic characterization of plant virus populations. Predating these advances, the field of Phylogenetics has significantly contributed to understand important aspects of plant virus evolution. This review aims at summarizing the impact of Phylogenetics in the current knowledge on three major aspects of plant virus evolution that have benefited from the development of phylogenetic inference: (1) The identification and classification of plant virus diversity. (2) The mechanisms and forces shaping the evolution of plant virus populations. (3) The understanding of the interaction between plant virus evolution, epidemiology and ecology. The work discussed here highlights the important role of phylogenetic approaches in the study of the dynamics of plant virus populations.

Correlated Selection on Amino Acid Deletion and Replacement in Mammalian Protein Sequences

A low ratio of nonsynonymous and synonymous substitution rates (dN/dS) at a codon is an indicator of functional constraint caused by purifying selection. Intuitively, the functional constraint would also be expected to prevent such a codon from being deleted. However, to the best of our knowledge, the correlation between the rates of deletion and substitution has never actually been estimated. Here, we use 8595 protein-coding region sequences from nine mammalian species to examine the relationship between deletion rate and dN/dS. We find significant positive correlations at the levels of both sites and genes. We compared our data against controls consisting of simulated coding sequences evolving along identical phylogenetic trees, where deletions occur independently of substitutions. A much weaker correlation was found in the corresponding simulated sequences, probably caused by alignment errors. In the real data, the correlations cannot be explained by alignment errors. Separate investigations on nonsynonymous (dN) and synonymous (dS) substitution rates indicate that the correlation is most likely due to a similarity in patterns of selection rather than in mutation rates.

Coupled enhancer and coding sequence evolution of a homeobox gene shaped leaf diversity

In this study, Vuolo et al. investigate the mechanisms underlying the genetic basis for morphological diversity in leaf shape. They show that evolution of an enhancer element in the homeobox gene REDUCED COMPLEXITY (RCO) altered leaf shape by changing gene expression from the distal leaf blade to its base.

Here we investigate mechanisms underlying the diversification of biological forms using crucifer leaf shape as an example. We show that evolution of an enhancer element in the homeobox gene REDUCED COMPLEXITY (RCO) altered leaf shape by changing gene expression from the distal leaf blade to its base. A single amino acid substitution evolved together with this regulatory change, which reduced RCO protein stability, preventing pleiotropic effects caused by its altered gene expression. We detected hallmarks of positive selection in these evolved regulatory and coding sequence variants and showed that modulating RCO activity can improve plant physiological performance. Therefore, interplay between enhancer and coding sequence evolution created a potentially adaptive path for morphological evolution.

An Ancient Fecundability-Associated Polymorphism Switches a Repressor into an Enhancer of Endometrial TAP2 Expression

Variation in female reproductive traits, such as fertility, fecundity, and fecundability, is heritable in humans, but identifying and functionally characterizing genetic variants associated with these traits has been challenging. Here, we explore the functional significance and evolutionary history of a T/C polymorphism of SNP rs2071473, which we have previously shown is an eQTL for TAP2 and significantly associated with fecundability (time to pregnancy). We replicated the association between the rs2071473 genotype and TAP2 expression by using GTEx data and demonstrated that TAP2 is expressed by decidual stromal cells at the maternal-fetal interface. Next, we showed that rs2071473 is located within a progesterone-responsive cis-regulatory element that functions as a repressor with the T allele and an enhancer with the C allele. Remarkably, we found that this polymorphism arose before the divergence of modern and archaic humans, segregates at intermediate to high frequencies across human populations, and has genetic signatures of long-term balancing selection. This variant has also previously been identified in genome-wide association studies of immune-related disease, suggesting that both alleles are maintained as a result of antagonistic pleiotropy.

DNA polymorphism and selection at the bindin locus in three Strongylocentrotus sp. (Echinoidea)

BACKGROUND

The sperm gene bindin encodes a gamete recognition protein, which plays an important role in conspecific fertilization and reproductive isolation of sea urchins. Molecular evolution of the gene has been extensively investigated with the attention focused on the protein coding regions. Intron evolution has been investigated to a much lesser extent. We have studied nucleotide variability in the complete bindin locus, including two exons and one intron, in the sea urchin Strongylocentrotus intermedius represented by two morphological forms. We have also analyzed all available bindin sequences for two other sea urchin species, S. pallidus and S. droebachiensis.

RESULTS

The results show that the bindin sequences from the two forms of S. intermedius are intermingled with no evidence of genetic divergence; however, the forms exhibit slightly different patterns in bindin variability. The level of the bindin nucleotide diversity is close for S. intermedius and S. droebachiensis, but noticeably higher for S. pallidus. The distribution of variability is non-uniform along the gene; however there are striking similarities among the species, indicating similar evolutionary trends in this gene engaged in reproductive function. The patterns of nucleotide variability and divergence are radically different in the bindin coding and intron regions. Positive selection is detected in the bindin coding region. The neutrality tests as well as the maximum likelihood approaches suggest the action of diversifying selection in the bindin intron.

CONCLUSIONS

Significant deviation from neutrality has been detected in the bindin coding region and suggested in the intron, indicating the possible functional importance of the bindin intron variability. To clarify the question concerning possible involvement of diversifying selection in the bindin intron evolution more data combining population genetic and functional approaches are necessary.

A novel strain of cynomolgus macaque cytomegalovirus: implications for host-virus co-evolution

Background

Cytomegaloviruses belong to a large, ancient, genus of DNA viruses comprised of a wide array of species-specific strains that occur in diverse array of hosts.

Methods

In this study we sequenced the ~217 Kb genome of a cytomegalovirus isolated from a Mauritius cynomolgus macaque, CyCMV Mauritius, and compared it to previously sequenced cytomegaloviruses from a cynomolgus macaque of Filipino origin (CyCMV Ottawa) and two from Indian rhesus macaques (RhCMV 180.92 and RhCMV 68–1).

Results

Though more closely related to CyCMV Ottawa, CyCMV Mauritius is less genetically distant from both RhCMV strains than is CyCMV Ottawa. Several individual genes, including homologues of CMV genes RL11B, UL123, UL83b, UL84 and a homologue of mammalian COX-2, show a closer relationship between homologues of CyCMV Mauritius and the RhCMVs than between homologues of CyCMV Mauritius and CyCMV Ottawa. A broader phylogenetic analysis of 12 CMV strains from eight species recovers evolutionary relationships among viral strains that mirror those amongst the host species, further demonstrating co-evolution of host and virus.

Conclusions

Phylogenetic analyses of rhesus and cynomolgus macaque CMV genome sequences demonstrate co-speciation of the virus and host.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2588-3) contains supplementary material, which is available to authorized users.

Estimate of within population incremental selection through branch imbalance in lineage trees

Incremental selection within a population, defined as limited fitness changes following mutation, is an important aspect of many evolutionary processes. Strongly advantageous or deleterious mutations are detected using the synonymous to non-synonymous mutations ratio. However, there are currently no precise methods to estimate incremental selection. We here provide for the first time such a detailed method and show its precision in multiple cases of micro-evolution. The proposed method is a novel mixed lineage tree/sequence based method to detect within population selection as defined by the effect of mutations on the average number of offspring. Specifically, we propose to measure the log of the ratio between the number of leaves in lineage trees branches following synonymous and non-synonymous mutations. The method requires a high enough number of sequences, and a large enough number of independent mutations. It assumes that all mutations are independent events. It does not require of a baseline model and is practically not affected by sampling biases. We show the method's wide applicability by testing it on multiple cases of micro-evolution. We show that it can detect genes and inter-genic regions using the selection rate and detect selection pressures in viral proteins and in the immune response to pathogens.

Evolutionary suppression of erythropoiesis via the modulation of TGF-β signalling in an Antarctic icefish

The Antarctic icefish, a family (Channichthyidae) of teleosts within the perciform suborder Notothenioidei, are the only known vertebrates without oxygen-transporting haemoglobins and that are largely devoid of circulating erythrocytes. To elucidate the evo-devo mechanisms underpinning the suppressed erythropoiesis in the icefish, we conducted comparative studies on the transcriptomes and microRNAomes of the primary haematopoietic tissues between an icefish (Chionodraco hamatus) and two red-blooded notothenioids (Trematomus bernacchii and Gymnodraco acuticeps). We identified substantial remodelling of the haematopoietic programs in the icefish through which erythropoiesis is selectively suppressed. Experimental verification showed that erythropoietic suppression in the icefish may be attributable to the upregulation of TGF-β signalling, which coincides with reductions in multiple transcription factors essential for erythropoiesis and the upregulation of hundreds of microRNAs, the majority (> 80%) of which potentially target erythropoiesis regulating factors. Of the six microRNAs selected for verification, three miRNAs (miR-152, miR-1388 and miR-16b) demonstrated suppressive functions on GATA1 and ALAS2, which are two factors important for erythroid differentiation, resulting in reduced numbers of erythroids in microinjected zebra fish embryos. Codon substitution analyses of the genes of the TGF-β superfamily revealed signs of positive selection in TGF-β1 and endoglin in the lineages leading to Antarctic notothenioids. Both genes are previously known to function in erythropoietic suppression. These findings implied a general trend of erythropoietic suppression in the cold-adapted notothenioid lineages through evolutionary modulation of the multi-functional TGF-β signalling pathway. This trend is more pronounced in the haemoglobin-less icefish, which may pre-emptively hinder the otherwise defective erythroids from production.

Maintenance of sex-related genes and the co-occurrence of both mating types in Verticillium dahliae

Verticillium dahliae is a cosmopolitan, soilborne fungus that causes a significant wilt disease on a wide variety of plant hosts including economically important crops, ornamentals, and timber species. Clonal expansion through asexual reproduction plays a vital role in recurring plant epidemics caused by this pathogen. The recent discovery of recombination between clonal lineages and preliminary investigations of the meiotic gene inventory of V. dahliae suggest that cryptic sex appears to be rare in this species. Here we expanded on previous findings on the sexual nature of V. dahliae. Only 1% of isolates in a global collection of 1120 phytopathogenic V. dahliae isolates contained the MAT1-1 idiomorph, whereas 99% contained MAT1-2. Nine unique multilocus microsatellite types comprised isolates of both mating types, eight of which were collected from the same substrate at the same time. Orthologs of 88 previously characterized sex-related genes from fungal model systems in the Ascoymycota were identified in the genome of V. dahliae, out of 93 genes investigated. Results of RT-PCR experiments using both mating types revealed that 10 arbitrarily chosen sex-related genes, including MAT1-1-1 and MAT1-2-1, were constitutively expressed in V. dahliae cultures grown under laboratory conditions. Ratios of non-synonymous (amino-acid altering) to synonymous (silent) substitutions in V. dahliae MAT1-1-1 and MAT1-2-1 sequences were indistinguishable from the ratios observed in the MAT genes of sexual fungi in the Pezizomycotina. Patterns consistent with strong purifying selection were also observed in 18 other arbitrarily chosen V. dahliae sex-related genes, relative to the patterns in orthologs from fungi with known sexual stages. This study builds upon recent findings from other laboratories and mounts further evidence for an ancestral or cryptic sexual stage in V. dahliae.

Natural selection on coding and noncoding DNA sequences is associated with virulence genes in a plant pathogenic fungus

Natural selection leaves imprints on DNA, offering the opportunity to identify functionally important regions of the genome. Identifying the genomic regions affected by natural selection within pathogens can aid in the pursuit of effective strategies to control diseases. In this study, we analyzed genome-wide patterns of selection acting on different classes of sequences in a worldwide sample of eight strains of the model plant-pathogenic fungus Colletotrichum graminicola. We found evidence of selective sweeps, balancing selection, and positive selection affecting both protein-coding and noncoding DNA of pathogenicity-related sequences. Genes encoding putative effector proteins and secondary metabolite biosynthetic enzymes show evidence of positive selection acting on the coding sequence, consistent with an Arms Race model of evolution. The 5' untranslated regions (UTRs) of genes coding for effector proteins and genes upregulated during infection show an excess of high-frequency polymorphisms likely the consequence of balancing selection and consistent with the Red Queen hypothesis of evolution acting on these putative regulatory sequences. Based on the findings of this work, we propose that even though adaptive substitutions on coding sequences are important for proteins that interact directly with the host, polymorphisms in the regulatory sequences may confer flexibility of gene expression in the virulence processes of this important plant pathogen.

The evolutionary jigsaw puzzle of the surviving trout (Salmo trutta L. complex) diversity in the Italian region. A multilocus Bayesian approach

Mediterranean trout populations display a diversity of phenotypes, representing a valuable model for the study of adaptation and a puzzling dilemma for taxonomists and biogeographers, which is further entangled by the widespread introgression of allochthonous genes. In this paper we analysed DNA polymorphism at multiple loci (sequence variation of the mitochondrial control region and eight nuclear fragments and length variation at eleven nuclear microsatellite loci) in representative samples of the autochthonous taxonomic diversity described in Italian trout populations (Salmo marmoratus, S. carpio, S. cenerinus, S. cettii and S. fibreni) and in samples from hatchery-originated strains of Atlantic S. trutta. We employed model-based clustering and Approximate Bayesian Computation in order to: (i) describe the phylogeographic structure of Italian autochthonous trout populations; (ii) evaluate a set of evolutionary/biogeographic models. The inclusion of hatchery-originated strains allowed to account for man-mediated allochthonous introgression in Italian populations. Our results (i) showed that the analysed sample consists of two main autochthonous evolutionary lineages, including the marble trout populations on one side ('marble' lineage) and the three peninsular populations of S. cettii, S. cenerinus and S. fibreni on the other side ('peninsular' lineage); (ii) indicated that S. carpio originated from a 'peninsular' population, with a possible, limited contribution from the 'marble' lineage; (iii) pointed out that the 'marble' lineage started diverging before the separation of the 'peninsular' lineage from Atlantic S. trutta; (iv) suggested that a model of divergence involving gene flow from the 'peninsular' population into the ancestral gene pool of 'marble' trout is most consistent with the genetic data; (v) provided evidence that the autochthonous trout gene pools in the Tyrrhenian and Adriatic basins of the Italian peninsula started diverging very recently (most likely after the last glacial maximum).

Genetic comparisons yield insight into the evolution of enamel thickness during human evolution

Enamel thickness varies substantially among extant hominoids and is a key trait with significance for interpreting dietary adaptation, life history trajectory, and phylogenetic relationships. There is a strong link in humans between enamel formation and mutations in the exons of the four genes that code for the enamel matrix proteins and the associated protease. The evolution of thick enamel in humans may have included changes in the regulation of these genes during tooth development. The cis-regulatory region in the 5' flank (upstream non-coding region) of MMP20, which codes for enamelysin, the predominant protease active during enamel secretion, has previously been shown to be under strong positive selection in the lineages leading to both humans and chimpanzees. Here we examine evidence for positive selection in the 5' flank and 3' flank of AMELX, AMBN, ENAM, and MMP20. We contrast the human sequence changes with other hominoids (chimpanzees, gorillas, orangutans, gibbons) and rhesus macaques (outgroup), a sample comprising a range of enamel thickness. We find no evidence for positive selection in the protein-coding regions of any of these genes. In contrast, we find strong evidence for positive selection in the 5' flank region of MMP20 and ENAM along the lineage leading to humans, and in both the 5' flank and 3' flank regions of MMP20 along the lineage leading to chimpanzees. We also identify putative transcription factor binding sites overlapping some of the species-specific nucleotide sites and we refine which sections of the up- and downstream putative regulatory regions are most likely to harbor important changes. These non-coding changes and their potential for differential regulation by transcription factors known to regulate tooth development may offer insight into the mechanisms that allow for rapid evolutionary changes in enamel thickness across closely-related species, and contribute to our understanding of the enamel phenotype in hominoids.

Cell cycle, oncogenic and tumor suppressor pathways regulate numerous long and macro non-protein-coding RNAs

Background

The genome is pervasively transcribed but most transcripts do not code for proteins, constituting non-protein-coding RNAs. Despite increasing numbers of functional reports of individual long non-coding RNAs (lncRNAs), assessing the extent of functionality among the non-coding transcriptional output of mammalian cells remains intricate. In the protein-coding world, transcripts differentially expressed in the context of processes essential for the survival of multicellular organisms have been instrumental in the discovery of functionally relevant proteins and their deregulation is frequently associated with diseases. We therefore systematically identified lncRNAs expressed differentially in response to oncologically relevant processes and cell-cycle, p53 and STAT3 pathways, using tiling arrays.

Results

We found that up to 80% of the pathway-triggered transcriptional responses are non-coding. Among these we identified very large macroRNAs with pathway-specific expression patterns and demonstrated that these are likely continuous transcripts. MacroRNAs contain elements conserved in mammals and sauropsids, which in part exhibit conserved RNA secondary structure. Comparing evolutionary rates of a macroRNA to adjacent protein-coding genes suggests a local action of the transcript. Finally, in different grades of astrocytoma, a tumor disease unrelated to the initially used cell lines, macroRNAs are differentially expressed.

Conclusions

It has been shown previously that the majority of expressed non-ribosomal transcripts are non-coding. We now conclude that differential expression triggered by signaling pathways gives rise to a similar abundance of non-coding content. It is thus unlikely that the prevalence of non-coding transcripts in the cell is a trivial consequence of leaky or random transcription events.

Multi Step Selection in Ig H Chains is Initially Focused on CDR3 and Then on Other CDR Regions

AFFINITY MATURATION OCCURS THROUGH TWO SELECTION PROCESSES

the choice of appropriate clones (clonal selection), and the internal evolution within clones, induced by somatic hyper-mutations, where high affinity mutants are selected for. When a final population of immunoglobulin sequences is observed, the genetic composition of this population is affected by a combination of these two processes. Different immune induced diseases can result from the failure of regulation of clonal selection or of the regulation of the within clone affinity maturation. In order to understand each of these processes separately, we propose a mixed lineage tree/sequence based method to detect within clone selection as defined by the effect of mutations on the average number of offspring. Specifically, we measure the imbalance in the number of leaves in lineage trees branches following synonymous and non-synonymous (NS) mutations. If a mutation is positively selected, we expect the number of leaves in the sub-tree below this mutation to be larger than in the parallel sub-tree without the mutation. The ratio between the number of leaves in such branches following NS mutations can be used to measure selection within a clone. We apply this method to the sampled Ig repertoire from multiple healthy volunteers and show that within clone selection is positive in the CDR2 region and either positive or negative in the CDR3 and FWR3 regions. Selection occurs already at the IgM isotype level mainly in the DH gene region, with a strong negative selection in the join region. This is followed in the later memory stages in the CDR2 region. We have not studied here the FWR1 and CDR1 regions. An important advantage of this method is that it is very weakly affected by the baseline mutation model or by sampling biases, as are most synonymous to NS mutations ratio based methods.

Selection and adaptation in the human genome

An enduring goal of evolutionary biology is to understand how natural selection has shaped patterns of polymorphism and divergence within and between species and to map the genetic basis of adaptations. The rapid maturation of next-generation sequencing technology has generated a deluge of genomics data from nonhuman primates, extinct hominins, and diverse human populations. These emerging genome data sets have simultaneously broadened our understanding of human evolution and sharply defined existing gaps in knowledge about the mechanistic basis of evolutionary change. In this review, we summarize recent insights into how natural selection has influenced the human genome across different timescales. Although the path to a more comprehensive understanding of selection and adaptation in humans remains arduous, some general insights are beginning to emerge, such as the importance of adaptive regulatory evolution, the absence of pervasive classic selective sweeps, and the potential roles that selection from standing variation and polygenic adaptation have likely played in recent human evolutionary history.

Putting pleiotropy and selection into context defines a new paradigm for interpreting genetic data

BACKGROUND

Natural selection shapes many human genes, including some related to complex diseases. Understanding how selection affects genes, especially pleiotropic ones, may be important in evaluating disease associations and the role played by environmental variation. This may be of particular interest for genes with antagonistic roles that cause divergent patterns of selection. The lectin-like low-density lipoprotein 1 receptor, encoded by OLR1, is exemplary. It has antagonistic functions in the cardiovascular and immune systems because the same protein domain binds oxidized low-density lipoprotein and bacterial cell wall proteins, the former contributing to atherosclerosis and the latter presumably protecting from infection. We studied patterns of selection in this gene, in humans and nonhuman primates, to determine whether variable selection can lead to conflicting results in cardiovascular disease association studies.

METHODS AND RESULTS

We analyzed sequences from 11 nonhuman primate species, as well as single-nucleotide polymorphisms and sequence data from multiple human populations. Results indicate that the derived allele is favored across primate lineages (probably because of recent positive selection). However, both the derived and ancestral alleles were maintained in human populations, especially European ones (possibly because of balancing selection derived from dual roles of LOX-1). Balancing selection likely reflects response to diverse environmental pressures among humans.

CONCLUSIONS

These data indicate that differential selection patterns, within and between species, in OLR1 render association studies difficult to replicate even if the gene is etiologically connected to cardiovascular disease. Selection analyses can identify genes exhibiting gene-environment interactions critical for unraveling disease association.

Rapid adversarial co-evolution of viruses and cellular restriction factors

Since the discovery of viruses over a century ago, virologists have recognized that host genetics plays a major role in viral tropism and the distribution of viruses in nature. Traditionally, studies of tropism have centered on identification of cellular factors required for viral replication, such as cell-surface entry receptors. However, over the past 20 years, there has been a steady increase in the identification and characterization of restriction factors (RFs), here defined as dominant cellular factors that have evolved specifically to interfere with viral replication. Genetic studies suggest that restriction factors impose significant barriers to interspecies movement of viruses and are therefore critical determinants of viral tropism. Furthermore, the scope of the ever-expanding list of restriction factors, and the variety of antiviral mechanisms they represent, testifies to the extraordinary impact viruses have had on organismal evolution-an impact hitherto underappreciated by evolutionary biologists and virologists alike. Recent studies of RF-encoding genes that combine molecular evolutionary analysis with functional assays illustrate the potential for asking questions about virus-host interactions as they play out in natural populations and across evolutionary timescales. Most notably, it has become common to apply tests of positive selection to RF genes and couple these analyses with virological assays, to reveal evidence for antagonistic virus-host co-evolution. Herein, I summarize recent work on the evolutionary genetics of mammalian RFs, particularly those of humans, non-human primates, and model organisms, and how RFs can reveal the influence of virus-host interactions on organismal evolution. Because intensive investigation of RF evolution is fairly new (and because there is still much to learn), the discussion is organized around five broad, outstanding questions that will need to be answered before we can fully appreciate the evolutionary biology of restriction.

Sex in cheese: evidence for sexuality in the fungus Penicillium roqueforti

Although most eukaryotes reproduce sexually at some moment of their life cycle, as much as a fifth of fungal species were thought to reproduce exclusively asexually. Nevertheless, recent studies have revealed the occurrence of sex in some of these supposedly asexual species. For industrially relevant fungi, for which inoculums are produced by clonal-subcultures since decades, the potentiality for sex is of great interest for strain improvement strategies. Here, we investigated the sexual capability of the fungus Penicillium roqueforti, used as starter for blue cheese production. We present indirect evidence suggesting that recombination could be occurring in this species. The screening of a large sample of strains isolated from diverse substrates throughout the world revealed the existence of individuals of both mating types, even in the very same cheese. The MAT genes, involved in fungal sexual compatibility, appeared to evolve under purifying selection, suggesting that they are still functional. The examination of the recently sequenced genome of the FM 164 cheese strain enabled the identification of the most important genes known to be involved in meiosis, which were found to be highly conserved. Linkage disequilibria were not significant among three of the six marker pairs and 11 out of the 16 possible allelic combinations were found in the dataset. Finally, the detection of signatures of repeat induced point mutations (RIP) in repeated sequences and transposable elements reinforces the conclusion that P. roqueforti underwent more or less recent sex events. In this species of high industrial importance, the induction of a sexual cycle would open the possibility of generating new genotypes that would be extremely useful to diversify cheese products.

Extensive evolutionary changes in regulatory element activity during human origins are associated with altered gene expression and positive selection

Understanding the molecular basis for phenotypic differences between humans and other primates remains an outstanding challenge. Mutations in non-coding regulatory DNA that alter gene expression have been hypothesized as a key driver of these phenotypic differences. This has been supported by differential gene expression analyses in general, but not by the identification of specific regulatory elements responsible for changes in transcription and phenotype. To identify the genetic source of regulatory differences, we mapped DNaseI hypersensitive (DHS) sites, which mark all types of active gene regulatory elements, genome-wide in the same cell type isolated from human, chimpanzee, and macaque. Most DHS sites were conserved among all three species, as expected based on their central role in regulating transcription. However, we found evidence that several hundred DHS sites were gained or lost on the lineages leading to modern human and chimpanzee. Species-specific DHS site gains are enriched near differentially expressed genes, are positively correlated with increased transcription, show evidence of branch-specific positive selection, and overlap with active chromatin marks. Species-specific sequence differences in transcription factor motifs found within these DHS sites are linked with species-specific changes in chromatin accessibility. Together, these indicate that the regulatory elements identified here are genetic contributors to transcriptional and phenotypic differences among primate species.

The human genome shares a remarkable amount of genomic sequence with our closest living primate relatives. Researchers have long sought to understand what regions of the genome are responsible for unique species-specific traits. Previous studies have shown that many genes are differentially expressed between species, but the regulatory elements contributing to these differences are largely unknown. Here we report a genome-wide comparison of active gene regulatory elements in human, chimpanzee, and macaque, and we identify hundreds of regulatory elements that have been gained or lost in the human or chimpanzee genomes since their evolutionary divergence. These elements contain evidence of natural selection and correlate with species-specific changes in gene expression. Polymorphic DNA bases in transcription factor motifs that we found in these regulatory elements may be responsible for the varied biological functions across species. This study directly links phenotypic and transcriptional differences between species with changes in chromatin structure.

Transformation of a transposon into a derived prolactin promoter with function during human pregnancy

Transposable elements (TEs) are known to provide DNA for host regulatory functions, but the mechanisms underlying the transformation of TEs into cis-regulatory elements are unclear. In humans two TEs--MER20 and MER39--contribute the enhancer/promoter for decidual prolactin (dPRL), which is dramatically induced during pregnancy. We show that evolution of the strong human dPRL promoter was a multistep process that took millions of years. First, MER39 inserted near MER20 in the primate/rodent ancestor, and then there were two phases of activity enhancement in primates. Through the mapping of causal nucleotide substitutions, we demonstrate that strong promoter activity in apes involves epistasis between transcription factor binding sites (TFBSs) ancestral to MER39 and derived sites. We propose a mode of molecular evolution that describes the process by which MER20/MER39 was transformed into a strong promoter, called "epistatic capture." Epistatic capture is the stabilization of a TFBS that is ancestral but variable in outgroup lineages, and is fixed in the ingroup because of epistatic interactions with derived TFBSs. Finally, we note that evolution of human promoter activity coincides with the emergence of a unique reproductive character in apes, highly invasive placentation. Because prolactin communicates with immune cells during pregnancy, which regulate fetal invasion into maternal tissues, we speculate that ape dPRL promoter activity evolved in response to increased invasiveness of ape fetal tissue.

Methods to detect selection on noncoding DNA

Vast tracts of noncoding DNA contain elements that regulate gene expression in higher eukaryotes. Describing these regulatory elements and understanding how they evolve represent major challenges for biologists. Advances in the ability to survey genome-scale DNA sequence data are providing unprecedented opportunities to use evolutionary models and computational tools to identify functionally important elements and the mode of selection acting on them in multiple species. This chapter reviews some of the current methods that have been developed and applied on noncoding DNA, what they have shown us, and how they are limited. Results of several recent studies reveal that a significantly larger fraction of noncoding DNA in eukaryotic organisms is likely to be functional than previously believed, implying that the functional annotation of most noncoding DNA in these organisms is largely incomplete. In Drosophila, recent studies have further suggested that a large fraction of noncoding DNA divergence observed between species may be the product of recurrent adaptive substitution. Similar studies in humans have revealed a more complex pattern, with signatures of recurrent positive selection being largely concentrated in conserved noncoding DNA elements. Understanding these patterns and the extent to which they generalize to other organisms awaits the analysis of forthcoming genome-scale polymorphism and divergence data from more species.

Genetic correlates of the evolving primate brain

The tremendous shifts in the size, structure, and function of the brain during primate evolution are ultimately caused by changes at the genetic level. Understanding what these changes are and how they effect the phenotypic changes observed lies at the heart of understanding evolutionary change. This chapter focuses on understanding the genetic basis of primate brain evolution, considering the substrates and mechanisms through which genetic change occurs. It also discusses the implications that our current understandings and tools have for what we have already discovered and where our studies will head in the future. While genetic and genomic studies have identified many regions undergoing positive selection during primate evolution, the findings are certainly not exhaustive and functional relevance remains to be confirmed. Nevertheless, a strong foundation has been built upon which future studies will emerge.

Selection on the protein-coding genome

Populations evolve as mutations arise in individual organisms and, through hereditary transmission, may become "fixed" (shared by all individuals) in the population. Most mutations are lethal or have negative fitness consequences for the organism. Others have essentially no effect on organismal fitness and can become fixed through the neutral stochastic process known as random drift. However, mutations may also produce a selective advantage that boosts their chances of reaching fixation. Regions of genes where new mutations are beneficial, rather than neutral or deleterious, tend to evolve more rapidly due to positive selection. Genes involved in immunity and defense are a well-known example; rapid evolution in these genes presumably occurs because new mutations help organisms to prevail in evolutionary "arms races" with pathogens. In recent years, genome-wide scans for selection have enlarged our understanding of the evolution of the protein-coding regions of the various species. In this chapter, we focus on the methods to detect selection in protein-coding genes. In particular, we discuss probabilistic models and how they have changed with the advent of new genome-wide data now available.

A potential role for glucose transporters in the evolution of human brain size

Differences in cognitive abilities and the relatively large brain are among the most striking differences between humans and their closest primate relatives. The energy trade-off hypothesis predicts that a major shift in energy allocation among tissues occurred during human origins in order to support the remarkable expansion of a metabolically expensive brain. However, the molecular basis of this adaptive scenario is unknown. Two glucose transporters (SLC2A1 and SLC2A4) are promising candidates and present intriguing mutations in humans, resulting, respectively, in microcephaly and disruptions in whole-body glucose homeostasis. We compared SLC2A1 and SLC2A4 expression between humans, chimpanzees and macaques, and found compensatory and biologically significant expression changes on the human lineage within cerebral cortex and skeletal muscle, consistent with mediating an energy trade-off. We also show that these two genes are likely to have undergone adaptation and participated in the development and maintenance of a larger brain in the human lineage by modulating brain and skeletal muscle energy allocation. We found that these two genes show human-specific signatures of positive selection on known regulatory elements within their 5'-untranslated region, suggesting an adaptation of their regulation during human origins. This study represents the first case where adaptive, functional and genetic lines of evidence implicate specific genes in the evolution of human brain size.

Evidence that purifying selection acts on promoter sequences

We tested whether functionally important sites in bacterial, yeast, and animal promoters are more conserved than their neighbors. We found that substitutions are predominantly seen in less important sites and that those that occurred tended to have less impact on gene expression than possible alternatives. These results suggest that purifying selection operates on promoter sequences.

Recombination and positive selection contributed to the evolution of Listeria monocytogenes lineages III and IV, two distinct and well supported uncommon L. monocytogenes lineages

Listeriamonocytogenes lineages III and IV represent two uncommon lineages of the human and animal pathogen L. monocytogenes, characterized by occurrence of unusual phenotypic and genetic characteristics that differentiate them from the common lineages I and II. To gain further insights into the evolution of lineages III and IV, we amplified and sequenced housekeeping genes (i.e., gap, prs, purM, ribC, and sigB), internalin genes (i.e., inlA, inlB, inlC, inlG, inlC2, inlD, inlE, inlF, and inlH) and the virulence gene cluster containing prfA, plcA, hly, mpl, actA, and plcB for lineages III (n = 7) and IV (n = 4) isolates. Phylogenetic analyses of the sequences obtained along with previously reported sequence data for 40 isolates representing lineages I (n = 18), II (n = 21), and III (n = 1), showed that lineages III and IV represent divergent and monophyletic lineages. The virulence gene cluster as well as the inlAB operon were present in all isolates, with inlF absent from all lineages III and IV isolates. While all lineage IV isolates contained only inlC (in addition to inlAB), lineage III isolates showed considerable diversity with regard to internalin gene presence, including presence of (i) only inlC (n = 2), (ii) inlC and inlGC2DE (n = 3), (iii) only inlGC2DE (n = 2), and (iv) inlC and inlC2DE (n = 1). In addition to evidence for horizontal gene transfer events, among lineages III and IV isolates, in prs, actA, plcB, mpl, inlA, inlB, inlG, inlD, and inlE, we also found significant evidence for positive selection in the hly promoter region and, along the lineages III and IV branches, for actA (including in sites recognized for interactions with proteins involved in actin tail polymerization). In conclusion, lineages III and IV represent two distinct monophyletic groups with contributions of intragenic recombination to the evolution of their internalin genes as well as contributions of positive selection to evolution of the virulence genes island.

Complex interplay of evolutionary forces in the ladybird homeobox genes of Drosophila melanogaster

Tandemly arranged paralogous genes lbe and lbl are members of the Drosophila NK homeobox family. We analyzed population samples of Drosophila melanogaster from Africa, Europe, North and South America, and single strains of D. sechellia, D. simulans, and D. yakuba within two linked regions encompassing partial sequences of lbe and lbl. The evolution of lbe and lbl is highly constrained due to their important regulatory functions. Despite this, a variety of forces have shaped the patterns of variation in lb genes: recombination, intragenic gene conversion and natural selection strongly influence background variation created by linkage disequilibrium and dimorphic haplotype structure. The two genes exhibited similar levels of nucleotide diversity and positive selection was detected in the noncoding regions of both genes. However, synonymous variability was significantly higher for lbe: no nonsynonymous changes were observed in this gene. We argue that balancing selection impacts some synonymous sites of the lbe gene. Stability of mRNA secondary structure was significantly different between the lbe (but not lbl) haplotype groups and may represent a driving force of balancing selection in epistatically interacting synonymous sites. Balancing selection on synonymous sites may be the first, or one of a few such observations, in Drosophila. In contrast, recurrent positive selection on lbl at the protein level influenced evolution at three codon sites. Transcription factor binding-site profiles were different for lbe and lbl, suggesting that their developmental functions are not redundant. Combined with our previous results on nucleotide variation in esterase and other homeobox genes, these results suggest that interplay of balancing and directional selection may be a general feature of molecular evolution in Drosophila and other eukaryote genomes.

Statistical methods for detecting natural selection from genomic data

In the study of molecular and phenotypic evolution, understanding the relative importance of random genetic drift and positive selection as the mechanisms for driving divergences between populations and maintaining polymorphisms within populations has been a central issue. A variety of statistical methods has been developed for detecting natural selection operating at the amino acid and nucleotide sequence levels. These methods may be largely classified into those aimed at detecting recurrent and/or recent/ongoing natural selection by utilizing the divergence and/or polymorphism data. Using these methods, pervasive positive selection has been identified for protein-coding and non-coding sequences in the genomic analysis of some organisms. However, many of these methods have been criticized by using computer simulation and real data analysis to produce excessive false-positives and to be sensitive to various disturbing factors. Importantly, some of these methods have been invalidated experimentally. These facts indicate that many of the statistical methods for detecting natural selection are unreliable. In addition, the signals that have been believed as the evidence for fixations of advantageous mutations due to positive selection may also be interpreted as the evidence for fixations of deleterious mutations due to random genetic drift. The genomic diversity data are rapidly accumulating in various organisms, and detection of natural selection may play a critical role for clarifying the relative role of random genetic drift and positive selection in molecular and phenotypic evolution. It is therefore important to develop reliable statistical methods that are unbiased as well as robust against various disturbing factors, for inferring natural selection.

Detection of lineage-specific evolutionary changes among primate species

BACKGROUND

Comparison of the human genome with other primates offers the opportunity to detect evolutionary events that created the diverse phenotypes among the primate species. Because the primate genomes are highly similar to one another, methods developed for analysis of more divergent species do not always detect signs of evolutionary selection.

RESULTS

We have developed a new method, called DivE, specifically designed to find regions that have evolved either more or less rapidly than expected, for any clade within a set of very closely related species. Unlike some previous methods, DivE does not rely on rates of synonymous and nonsynonymous substitution, which enables it to detect evolutionary events in noncoding regions. We demonstrate using simulated data that DivE compares favorably to alternative methods, and we then apply DivE to the ENCODE regions in 14 primate species. We identify thousands of regions in these primates, ranging from 50 to >10000 bp in length, that appear to have experienced either constrained or accelerated rates of evolution. In particular, we detected 4942 regions that have potentially undergone positive selection in one or more primate species. Most of these regions occur outside of protein-coding genes, although we identified 20 proteins that have experienced positive selection.

CONCLUSIONS

DivE provides an easy-to-use method to predict both positive and negative selection in noncoding DNA, that is particularly well-suited to detecting lineage-specific selection in large genomes.

The sequence, structure and evolutionary features of HOTAIR in mammals

Background

An increasing number of long noncoding RNAs (lncRNAs) have been identified recently. Different from all the others that function in cis to regulate local gene expression, the newly identified HOTAIR is located between HoxC11 and HoxC12 in the human genome and regulates HoxD expression in multiple tissues. Like the well-characterised lncRNA Xist, HOTAIR binds to polycomb proteins to methylate histones at multiple HoxD loci, but unlike Xist, many details of its structure and function, as well as the trans regulation, remain unclear. Moreover, HOTAIR is involved in the aberrant regulation of gene expression in cancer.

Results

To identify conserved domains in HOTAIR and study the phylogenetic distribution of this lncRNA, we searched the genomes of 10 mammalian and 3 non-mammalian vertebrates for matches to its 6 exons and the two conserved domains within the 1800 bp exon6 using Infernal. There was just one high-scoring hit for each mammal, but many low-scoring hits were found in both mammals and non-mammalian vertebrates. These hits and their flanking genes in four placental mammals and platypus were examined to determine whether HOTAIR contained elements shared by other lncRNAs. Several of the hits were within unknown transcripts or ncRNAs, many were within introns of, or antisense to, protein-coding genes, and conservation of the flanking genes was observed only between human and chimpanzee. Phylogenetic analysis revealed discrete evolutionary dynamics for orthologous sequences of HOTAIR exons. Exon1 at the 5' end and a domain in exon6 near the 3' end, which contain domains that bind to multiple proteins, have evolved faster in primates than in other mammals. Structures were predicted for exon1, two domains of exon6 and the full HOTAIR sequence. The sequence and structure of two fragments, in exon1 and the domain B of exon6 respectively, were identified to robustly occur in predicted structures of exon1, domain B of exon6 and the full HOTAIR in mammals.

Conclusions

HOTAIR exists in mammals, has poorly conserved sequences and considerably conserved structures, and has evolved faster than nearby HoxC genes. Exons of HOTAIR show distinct evolutionary features, and a 239 bp domain in the 1804 bp exon6 is especially conserved. These features, together with the absence of some exons and sequences in mouse, rat and kangaroo, suggest ab initio generation of HOTAIR in marsupials. Structure prediction identifies two fragments in the 5' end exon1 and the 3' end domain B of exon6, with sequence and structure invariably occurring in various predicted structures of exon1, the domain B of exon6 and the full HOTAIR.

Integrating evolutionary and functional approaches to infer adaptation at specific loci

Inferences about adaptation at specific loci are often exclusively based on the static analysis of DNA sequence variation. Ideally,population-genetic evidence for positive selection serves as a stepping-off point for experimental studies to elucidate the functional significance of the putatively adaptive variation. We argue that inferences about adaptation at specific loci are best achieved by integrating the indirect, retrospective insights provided by population-genetic analyses with the more direct, mechanistic insights provided by functional experiments. Integrative studies of adaptive genetic variation may sometimes be motivated by experimental insights into molecular function, which then provide the impetus to perform population genetic tests to evaluate whether the functional variation is of adaptive significance. In other cases, studies may be initiated by genome scans of DNA variation to identify candidate loci for recent adaptation. Results of such analyses can then motivate experimental efforts to test whether the identified candidate loci do in fact contribute to functional variation in some fitness-related phenotype. Functional studies can provide corroborative evidence for positive selection at particular loci, and can potentially reveal specific molecular mechanisms of adaptation.

Nonadaptive processes in primate and human evolution

Evolutionary biology has tended to focus on adaptive evolution by positive selection as the primum mobile of evolutionary trajectories in species while underestimating the importance of nonadaptive evolutionary processes. In this review, I describe evidence that suggests that primate and human evolution has been strongly influenced by nonadaptive processes, particularly random genetic drift and mutation. This is evidenced by three fundamental effects: a relative relaxation of selective constraints (i.e., purifying selection), a relative increase in the fixation of slightly deleterious mutations, and a general reduction in the efficacy of positive selection. These effects are observed in protein-coding, regulatory regions, and in gene expression data, as well as in an augmentation of fixation of large-scale mutations, including duplicated genes, mobile genetic elements, and nuclear mitochondrial DNA. The evidence suggests a general population-level explanation such as a reduction in effective population size (N(e)). This would have tipped the balance between the evolutionary forces of natural selection and random genetic drift toward genetic drift for variants having small selective effects. After describing these proximate effects, I describe the potential consequences of these effects for primate and human evolution. For example, an increase in the fixation of slightly deleterious mutations could potentially have led to an increase in the fixation rate of compensatory mutations that act to suppress the effects of slightly deleterious substitutions. The potential consequences of compensatory evolution for the evolution of novel gene functions and in potentially confounding the detection of positively selected genes are explored. The consequences of the passive accumulation of large-scale genomic mutations by genetic drift are unclear, though evidence suggests that new gene copies as well as insertions of transposable elements into genes can potentially lead to adaptive phenotypes. Finally, because a decrease in selective constraint at the genetic level is expected to have effects at the morphological level, I review studies that compare rates of morphological change in various mammalian and island populations where N(e) is reduced. Furthermore, I discuss evidence that suggests that craniofacial morphology in the Homo lineage has shifted from an evolutionary rate constrained by purifying selection toward a neutral evolutionary rate.

Comparative genetic approaches to the evolution of human brain and behavior

With advances in genomic technologies, the amount of genetic data available to scientists today is vast. Genomes are now available or planned for 14 different primate species and complete resequencing of numerous human individuals from numerous populations is underway. Moreover, high-throughput deep sequencing is quickly making whole genome efforts within the reach of single laboratories allowing for unprecedented studies. Comparative genetic approaches to the identification of the underlying basis of human brain, behavior, and cognitive ability are moving to the forefront. Two approaches predominate: inter-species divergence comparisons and intra-species polymorphism studies. These methodological differences are useful for different time scales of evolution and necessarily focus on different evolutionary events in the history of primate and hominin evolution. Inter-species divergence is more useful in studying large scale primate, or hominoid, evolution whereas intra-species polymorphism can be more illuminating of recent hominin evolution. These differences in methodological utility also extend to studies of differing genetic substrates; current divergence studies focus primarily on protein evolution whereas polymorphism studies are substrate ambivalent. Some of the issues inherent in these studies can be ameliorated by current sequencing capabilities whereas others remain intractable. New avenues are also being opened that allow for the incorporation of novel substrates and approaches. In the post-genomic era, the study of human evolution, specifically as it relates to the brain, is becoming more complete focusing increasingly on the totality of the system and better conceptualizing the entirety of the genetic changes that have lead to the human phenotype today.

Mito-nuclear genetic comparison in a Wolbachia infected weevil: insights on reproductive mode, infection age and evolutionary forces shaping genetic variation

BACKGROUND

Maternally inherited endosymbionts like Wolbachia pipientis are in linkage disequilibrium with the mtDNA of their hosts. Therefore, they can induce selective sweeps, decreasing genetic diversity over many generations. This sex ratio distorter, that is involved in the origin of parthenogenesis and other reproductive alterations, infects the parthenogenetic weevil Naupactus cervinus, a serious pest of ornamental and fruit plants.

RESULTS

Molecular evolution analyses of mitochondrial (COI) and nuclear (ITS1) sequences from 309 individuals of Naupactus cervinus sampled over a broad range of its geographical distribution were carried out. Our results demonstrate lack of recombination in the nuclear fragment, non-random association between nuclear and mitochondrial genomes and the consequent coevolution of both genomes, being an indirect evidence of apomixis. This weevil is infected by a single Wolbachia strain, which could have caused a moderate bottleneck in the invaded population which survived the initial infection.

CONCLUSIONS

Clonal reproduction and Wolbachia infection induce the coevolution of bacterial, mitochondrial and nuclear genomes. The time elapsed since the Wolbachia invasion would have erased the traces of the demographic crash in the mtDNA, being the nuclear genome the only one that retained the signal of the bottleneck. The amount of genetic change accumulated in the mtDNA and the high prevalence of Wolbachia in all populations of N. cervinus agree with the hypothesis of an ancient infection. Wolbachia probably had great influence in shaping the genetic diversity of N. cervinus. However, it would have not caused the extinction of males, since sexual and asexual infected lineages coexisted until recent times.

Molecular analysis of phosphomannomutase (PMM) genes reveals a unique PMM duplication event in diverse Triticeae species and the main PMM isozymes in bread wheat tissues

BACKGROUND

Phosphomannomutase (PMM) is an essential enzyme in eukaryotes. However, little is known about PMM gene and function in crop plants. Here, we report molecular evolutionary and biochemical analysis of PMM genes in bread wheat and related Triticeae species.

RESULTS

Two sets of homologous PMM genes (TaPMM-1 and 2) were found in bread wheat, and two corresponding PMM genes were identified in the diploid progenitors of bread wheat and many other diploid Triticeae species. The duplication event yielding PMM-1 and 2 occurred before the radiation of diploid Triticeae genomes. The PMM gene family in wheat and relatives may evolve largely under purifying selection. Among the six TaPMM genes, the transcript levels of PMM-1 members were comparatively high and their recombinant proteins were all enzymatically active. However, PMM-2 homologs exhibited lower transcript levels, two of which were also inactive. TaPMM-A1, B1 and D1 were probably the main active isozymes in bread wheat tissues. The three isozymes differed from their counterparts in barley and Brachypodium distachyon in being more tolerant to elevated test temperatures.

CONCLUSION

Our work identified the genes encoding PMM isozymes in bread wheat and relatives, uncovered a unique PMM duplication event in diverse Triticeae species, and revealed the main active PMM isozymes in bread wheat tissues. The knowledge obtained here improves the understanding of PMM evolution in eukaryotic organisms, and may facilitate further investigations of PMM function in the temperature adaptability of bread wheat.

Molecular evolution of the endosperm starch synthesis pathway genes in rice (Oryza sativa L.) and its wild ancestor, O. rufipogon L

The evolution of metabolic pathways is a fundamental but poorly understood aspect of evolutionary change. One approach for understanding the complexity of pathway evolution is to examine the molecular evolution of genes that together comprise an integrated metabolic pathway. The rice endosperm starch biosynthetic pathway is one of the most thoroughly characterized metabolic pathways in plants, and starch is a trait that has evolved in response to strong selection during rice domestication. In this study, we have examined six key genes (AGPL2, AGPS2b, SSIIa, SBEIIb, GBSSI, ISA1) in the rice endosperm starch biosynthesis pathway to investigate the evolution of these genes before and after rice domestication. Genome-wide sequence tagged sites data were used as a neutral reference to overcome the problems of detecting selection in species with complex demographic histories such as rice. Five variety groups of Oryza sativa (aus, indica, tropical japonica, temperate japonica, aromatic) and its wild ancestor (O. rufipogon) were sampled. Our results showed evidence of purifying selection at AGPL2 in O. rufipogon and strong evidence of positive selection at GBSSI in temperate japonica and tropical japonica varieties and at GBSSI and SBEIIb in aromatic varieties. All the other genes showed a pattern consistent with neutral evolution in both cultivated rice and its wild ancestor. These results indicate the important role of positive selection in the evolution of starch genes during rice domestication. We discuss the role of SBEIIb and GBSSI in the evolution of starch quality during rice domestication and the power and limitation of detecting selection using genome-wide data as a neutral reference.

An effective model for natural selection in promoters

We have produced an evolutionary model for promoters, analogous to the commonly used synonymous/nonsynonymous mutation models for protein-coding sequences. Although our model, called Sunflower, relies on some simple assumptions, it captures enough of the biology of transcription factor action to show clear correlation with other biological features. Sunflower predicts a binding profile of transcription factors to DNA sequences, in which different factors compete for the same potential binding sites. The parametrized model simultaneously estimates a continuous measurement of binding occupancy across the genomic sequence for each factor. We can then introduce a localized mutation, rerun the binding model, and record the difference in binding profiles. A single mutation can alter interactions both upstream and downstream of its position due to potential overlapping binding sites, and our statistic captures this domino effect. Over evolutionary time, we observe a clear excess of low-scoring mutations fixed in promoters, consistent with most changes being neutral. However, this is not consistent across all promoters, and some promoters show more rapid divergence. This divergence often occurs in the presence of relatively constant protein-coding divergence. Interestingly, different classes of promoters show different sensitivity to mutations, with phosphorylation-related genes having promoters inherently more sensitive to mutations than immune genes. Although there have previously been a number of models attempting to handle transcription factor binding, Sunflower provides a richer biological model, incorporating weak binding sites and the possibility of competition. The results show the first clear correlations between such a model and evolutionary processes.

Statistical tests for natural selection on regulatory regions based on the strength of transcription factor binding sites

Background

Although cis-regulatory changes play an important role in evolution, it remains difficult to establish the contribution of natural selection to regulatory differences between species. For protein coding regions, powerful tests of natural selection have been developed based on comparisons of synonymous and non-synonymous substitutions, and analogous tests for regulatory regions would be of great utility.

Results

Here, tests for natural selection on regulatory regions are proposed based on nucleotide substitutions that occur in characterized transcription factor binding sites (an important type functional element within regulatory regions). In the absence of selection, these substitutions will tend to reduce the strength of existing binding sites. On the other hand, purifying selection will act to preserve the binding sites in regulatory regions, while positive selection can act to create or destroy binding sites, as well as change their strength. Using standard models of binding site strength and molecular evolution in the absence of selection, this intuition can be used to develop statistical tests for natural selection. Application of these tests to two well-characterized regulatory regions in Drosophila provides evidence for purifying selection.

Conclusion

This demonstrates that it is possible to develop tests for selection on regulatory regions based on the specific functional constrains on these sequences.