Sequence divergence, functional constraint, and selection in protein evolution

Evaluating the role of natural selection in the evolution of gene regulation

Surveys of gene expression reveal extensive variability both within and between a wide range of species. Compelling cases have been made for adaptive changes in gene regulation, but the proportion of expression divergence attributable to natural selection remains unclear. Distinguishing adaptive changes driven by positive selection from neutral divergence resulting from mutation and genetic drift is critical for understanding the evolution of gene expression. Here, we review the various methods that have been used to test for signs of selection in genomic expression data. We also discuss properties of regulatory systems relevant to neutral models of gene expression. Despite some potential caveats, published studies provide considerable evidence for adaptive changes in gene expression. Future challenges for studies of regulatory evolution will be to quantify the frequency of adaptive changes, identify the genetic basis of expression divergence and associate changes in gene expression with specific organismal phenotypes.

Protein polymorphism is negatively correlated with conservation of intronic sequences and complexity of expression patterns in Drosophila melanogaster

We report a significant negative correlation between nonsynonymous polymorphism and intron length in Drosophila melanogaster. This correlation is similar to that between protein divergence and intron length previously reported in Drosophila. We show that the relationship can be explained by the content of conserved noncoding sequences (CNS) within introns. In addition, genes with a high regulatory complexity and many genetic interactions also exhibit larger amounts of CNS within their introns and lower values of nonsynonymous polymorphism. The present study provides relevant evidence on the importance of intron content and expression patterns on the levels of coding polymorphism.

Functional nonsynonymous single nucleotide polymorphisms from the TGF-β protein interaction network

Protein complexes mediated by protein-protein interactions are essential for many cellular functions. Transforming growth factor (TGF)-β signaling involves a cascade of protein-protein interactions and malfunctioning of this pathway has been implicated in human diseases. Using an in silico approach, we analyzed the naturally occurring human genetic variations from the proteins involved in the TGF-β signaling (10 TGF-β proteins and 242 other proteins interacting with them) to identify the ones that have potential biological consequences. All proteins were searched in the dbSNP database for the presence of nonsynonymous single nucleotide polymorphisms (nsSNPs). A total of 118 validated nsSNPs from 63 proteins were retrieved and analyzed in terms of 1) evolutionary conservation status, 2) being located in a functional protein domain or motif, and 3) altering putative protein motif or phosphorylation sites. Our results indicated the presence of 31 nsSNPs that occurred at evolutionarily conserved residues, 37 nsSNPs were located in protein domains, motifs, or repeats, and 46 nsSNPs were predicted to either create or abolish putative protein motifs or phosphorylation sites. We undertook this study to analyze the human genetic variations that can affect the protein function and the TGF-β signaling. The nsSNPs reported in here can be characterized by experimental approaches to elucidate their exact biological roles and whether they are related to human disease.

Long-wavelength sensitive visual pigments of the guppy (Poecilia reticulata): six opsins expressed in a single individual

BACKGROUND

The diversity of visual systems in fish has long been of interest for evolutionary biologists and neurophysiologists, and has recently begun to attract the attention of molecular evolutionary geneticists. Several recent studies on the copy number and genomic organization of visual pigment proteins, the opsins, have revealed an increased opsin diversity in fish relative to most vertebrates, brought about through recent instances of opsin duplication and divergence. However, for the subfamily of opsin genes that mediate vision at the long-wavelength end of the spectrum, the LWS opsins, it appears that most fishes possess only one or two loci, a value comparable to most other vertebrates. Here, we characterize the LWS opsins from cDNA of an individual guppy, Poecilia reticulata, a fish that is known exhibit variation in its long-wavelength sensitive visual system, mate preferences and colour patterns.

RESULTS

We identified six LWS opsins expressed within a single individual. Phylogenetic analysis revealed that these opsins descend from duplication events both pre-dating and following the divergence of the guppy lineage from that of the bluefin killifish, Lucania goodei, the closest species for which comparable data exists. Numerous amino acid substitutions exist among these different LWS opsins, many at sites known to be important for visual pigment function, including spectral sensitivity and G-protein activation. Likelihood analyses using codon-based models of evolution reveal significant changes in selective constraint along two of the guppy LWS opsin lineages.

CONCLUSION

The guppy displays an unusually high number of LWS opsins compared to other fish, and to vertebrates in general. Observing both substitutions at functionally important sites and the persistence of lineages across species boundaries suggests that these opsins might have functionally different roles, especially with regard to G-protein activation. The reasons why are currently unknown, but may relate to aspects of the guppy's behavioural ecology, in which both male colour patterns and the female mate preferences for these colour patterns experience strong, highly variable selection pressures.

FUNC: a package for detecting significant associations between gene sets and ontological annotations

Background

Genome-wide expression, sequence and association studies typically yield large sets of gene candidates, which must then be further analysed and interpreted. Information about these genes is increasingly being captured and organized in ontologies, such as the Gene Ontology. Relationships between the gene sets identified by experimental methods and biological knowledge can be made explicit and used in the interpretation of results. However, it is often difficult to assess the statistical significance of such analyses since many inter-dependent categories are tested simultaneously.

Results

We developed the program package FUNC that includes and expands on currently available methods to identify significant associations between gene sets and ontological annotations. Implemented are several tests in particular well suited for genome wide sequence comparisons, estimates of the family-wise error rate, the false discovery rate, a sensitive estimator of the global significance of the results and an algorithm to reduce the complexity of the results.

Conclusion

FUNC is a versatile and useful tool for the analysis of genome-wide data. It is freely available under the GPL license and also accessible via a web service.

Isolation and characterization of LCHN: a novel factor induced by transient global ischemia in the adult rat hippocampus

Using mRNA differential display to identify cerebral ischemia-responsive mRNAs, we isolated and cloned a cDNA derived from a novel gene, that has been designated LCHN. Antisense mRNA in situ hybridization and immunoblotting confirmed LCHN expression to be induced in the rat hippocampus following transient forebrain ischemia. The deduced amino acid sequence of the novel LCHN cDNA contains an open reading frame of 455 amino acids, encoding a protein with a predicted molecular mass of approximately 51 kDa. Although LCHN is highly conserved between rat, mouse, and human, the deduced amino acid sequence of LCHN does not possess significant homology to other known genes. LCHN immunoreactivity is detected within the somatodendritic compartment of neurons, is also present on dendritic growth cones, but is not detected on astrocytes. The induction of LCHN in the hippocampus following ischemic injury may have functional consequences, as the ectopic over-expression of LCHN generated neurons with longer and more branched axons and dendrites. Taken together, these data suggest that LCHN could play a role in neuritogenesis, as well as in neuronal recovery and/or restructuring in the hippocampus following transient cerebral ischemia.

Genomics and conservation genetics

In large part, the relevance of genetics to conservation rests on the premise that neutral marker variation in populations reflects levels of detrimental and adaptive genetic variation. Despite its prominence, this tenet has been difficult to evaluate, until now. As we discuss here, genome sequence information and new technological and bioinformatics platforms now enable comprehensive surveys of neutral variation and more direct inferences of detrimental and adaptive variation in species with sequenced genomes and in 'genome-enabled' endangered taxa. Moreover, conservation schemes could begin to consider specific pathological genetic variants. A new conservation genetic agenda would utilize data from enhanced surveys of genomic variation in endangered species to better manage functional genetic variation.

Evolution of primate gene expression

It has been suggested that evolutionary changes in gene expression account for most phenotypic differences between species, in particular between humans and apes. What general rules can be described governing expression evolution? We find that a neutral model where negative selection and divergence time are the major factors is a useful null hypothesis for both transcriptome and genome evolution. Two tissues that stand out with regard to gene expression are the testes, where positive selection has exerted a substantial influence in both humans and chimpanzees, and the brain, where gene expression has changed less than in other organs but acceleration might have occurred in human ancestors.

Functional nsSNPs from carcinogenesis-related genes expressed in breast tissue: potential breast cancer risk alleles and their distribution across human populations

Although highly penetrant alleles of BRCA1 and BRCA2 have been shown to predispose to breast cancer, the majority of breast cancer cases are assumed to result from the presence of low-moderate penetrant alleles and environmental carcinogens. Non-synonymous single nucleotide polymorphisms (nsSNPs) are hypothesised to contribute to disease susceptibility and approximately 30 per cent of them are predicted to have a biological significance. In this study, we have applied a bioinformatics-based strategy to identify breast cancer-related nsSNPs from 981 carcinogenesis-related genes expressed in breast tissue. Our results revealed a total of 367 validated nsSNPs, 109 (29.7 per cent) of which are predicted to affect the protein function (functional nsSNPs), suggesting that these nsSNPs are likely to influence the development and homeostasis of breast tissue and hence contribute to breast cancer susceptibility. Sixty-seven of the functional nsSNPs presented as commonly occurring nsSNPs (minor allele frequencies ≥ 5 per cent), representing excellent candidates for breast cancer susceptibility. Additionally, a non-uniform distribution of the common functional nsSNPs among different human populations was observed: 15 nsSNPs were reported to be present in all populations analysed, whereas another set of 15 nsSNPs was specific to particular population(s). We propose that the nsSNPs analysed in this study constitute a unique resource of potential genetic factors for breast cancer susceptibility. Furthermore, the variations in functional nsSNP allele frequencies across major population backgrounds may point to the potential variability of the molecular basis of breast cancer predisposition and treatment response among different human populations.

Positive selection in the evolution of cancer

We hypothesize that forms of antagonistic coevolution have forged strong links between positive selection at the molecular level and increased cancer risk. By this hypothesis, evolutionary conflict between males and females, mothers and foetuses, hosts and parasites, and other parties with divergent fitness interests has led to rapid evolution of genetic systems involved in control over fertilization and cellular resources. The genes involved in such systems promote cancer risk as a secondary effect of their roles in antagonistic coevolution, which generates evolutionary disequilibrium and maladaptation. Evidence from two sources: (1) studies on specific genes, including SPANX cancer/testis antigen genes, several Y-linked genes, the pem homebox gene, centromeric histone genes, the breast cancer gene BRCA1, the angiogenesis gene ANG, cadherin genes, cytochrome P450 genes, and viral oncogenes; and (2) large-scale database studies of selection on different functional categories of genes, supports our hypothesis. These results have important implications for understanding the evolutionary underpinnings of cancer and the dynamics of antagonistically-coevolving molecular systems.

Multiple mechanisms promote the retained expression of gene duplicates in the tetraploid frog Xenopus laevis

Gene duplication provides a window of opportunity for biological variants to persist under the protection of a co-expressed copy with similar or redundant function. Duplication catalyzes innovation (neofunctionalization), subfunction degeneration (subfunctionalization), and genetic buffering (redundancy), and the genetic survival of each paralog is triggered by mechanisms that add, compromise, or do not alter protein function. We tested the applicability of three types of mechanisms for promoting the retained expression of duplicated genes in 290 expressed paralogs of the tetraploid clawed frog, Xenopus laevis. Tests were based on explicit expectations concerning the ka/ks ratio, and the number and location of nonsynonymous substitutions after duplication. Functional constraints on the majority of paralogs are not significantly different from a singleton ortholog. However, we recover strong support that some of them have an asymmetric rate of nonsynonymous substitution: 6% match predictions of the neofunctionalization hypothesis in that (1) each paralog accumulated nonsynonymous substitutions at a significantly different rate and (2) the one that evolves faster has a higher ka/ks ratio than the other paralog and than a singleton ortholog. Fewer paralogs (3%) exhibit a complementary pattern of substitution at the protein level that is predicted by enhancement or degradation of different functional domains, and the remaining 13% have a higher average ka/ks ratio in both paralogs that is consistent with altered functional constraints, diversifying selection, or activity-reducing mutations after duplication. We estimate that these paralogs have been retained since they originated by genome duplication between 21 and 41 million years ago. Multiple mechanisms operate to promote the retained expression of duplicates in the same genome, in genes in the same functional class, over the same period of time following duplication, and sometimes in the same pair of paralogs. None of these paralogs are superfluous; degradation or enhancement of different protein subfunctions and neofunctionalization are plausible hypotheses for the retained expression of some of them. Evolution of most X. laevis paralogs, however, is consistent with retained expression via mechanisms that do not radically alter functional constraints, such as selection to preserve post-duplication stoichiometry or temporal, quantitative, or spatial subfunctionalization.

Gene duplication plays a fundamental role in biological innovation but it is not clear how both copies of a duplicated gene manage to circumvent degradation by mutation if neither is unique. This study explores genetic mechanisms that could make each copy of a duplicate gene different, and therefore distinguishable and potentially preserved by natural selection. It is based on DNA sequences of the protein-coding region of 290 expressed duplicated genes in a frog, Xenopus laevis, that underwent complete duplication of its entire genome. Results provide evidence for multiple mechanisms acting within the same genome, within the same functional classes of genes, within the same period of time following duplication, and even on the same set of duplicated genes. Each copy of a duplicate gene may be subject to distinct evolutionary constraints, and this could be associated with degradation or enhancement of function. Functional constraints of most of these duplicates, however, are not substantially different from a single copy gene; their persistence in the first dozens of millions of years after duplication may more frequently be explained by mechanisms acting on their expression rather than their function.

Identification of genes subject to positive selection in uropathogenic strains of Escherichia coli: a comparative genomics approach

Escherichia coli is a model laboratory bacterium, a species that is widely distributed in the environment, as well as a mutualist and pathogen in its human hosts. As such, E. coli represents an attractive organism to study how environment impacts microbial genome structure and function. Uropathogenic E. coli (UPEC) must adapt to life in several microbial communities in the human body, and has a complex life cycle in the bladder when it causes acute or recurrent urinary tract infection (UTI). Several studies designed to identify virulence factors have focused on genes that are uniquely represented in UPEC strains, whereas the role of genes that are common to all E. coli has received much less attention. Here we describe the complete 5,065,741-bp genome sequence of a UPEC strain recovered from a patient with an acute bladder infection and compare it with six other finished E. coli genome sequences. We searched 3,470 ortholog sets for genes that are under positive selection only in UPEC strains. Our maximum likelihood-based analysis yielded 29 genes involved in various aspects of cell surface structure, DNA metabolism, nutrient acquisition, and UTI. These results were validated by resequencing a subset of the 29 genes in a panel of 50 urinary, periurethral, and rectal E. coli isolates from patients with UTI. These studies outline a computational approach that may be broadly applicable for studying strain-specific adaptation and pathogenesis in other bacteria.

UVPAR: fast detection of functional shifts in duplicate genes

BACKGROUND

The imprint of natural selection on gene sequences is often difficult to detect. A plethora of methods have been devised to detect genetic changes due to selective processes. However, many of those methods depend heavily on underlying assumptions regarding the mode of change of DNA sequences and often require sophisticated mathematical treatments that made them computationally slow. The development of fast and effective methods to detect modifications in the selective constraints of genes is therefore of great interest.

RESULTS

We describe UVPAR, a program designed to quickly test for changes in the functional constraints of duplicate genes. Starting with alignments of the proteins encoded by couples of duplicate genes in two different species, UVPAR detects the regions in which modifications of the functional constraints in the paralogs occurred since both species diverged. Sequences can be analyzed with UVPAR in just a few minutes on a standard PC computer. To demonstrate the power of the program, we first show how the results obtained with UVPAR compare to those based on other approaches, using data for vertebrate Hox genes. We then describe a comprehensive study of the RBR family of ubiquitin ligases in which we have performed 529 analyses involving 14 duplicate genes in seven model species. A significant increase in the number of functional shifts was observed for the species Danio rerio and for the gene Ariadne-2.

CONCLUSION

These results show that UVPAR can be used to generate sensitive analyses to detect changes in the selection constraints acting on paralogs. The high speed of the program allows its application to genome-scale analyses.

The accumulation of deleterious mutations in rice genomes: a hypothesis on the cost of domestication

The extent of molecular differentiation between domesticated animals or plants and their wild relatives is postulated to be small. The availability of the complete genome sequences of two subspecies of the Asian rice, Oryza sativa (indica and japonica) and their wild relatives have provided an unprecedented opportunity to study divergence following domestication. We observed significantly more amino acid substitutions during rice domestication than can be expected from a comparison among wild species. This excess is disproportionately larger for the more radical kinds of amino acid changes (e.g. Cys<-->Tyr). We estimate that approximately a quarter of the amino acid differences between rice cultivars are deleterious, not accountable by the relaxation of selective constraints. This excess is negatively correlated with the rate of recombination, suggesting that 'hitchhiking' has occurred. We hypothesize that during domestication artificial selection increased the frequency of many deleterious mutations.

Human SNPs resulting in premature stop codons and protein truncation

Single nucleotide polymorphisms (SNPs) constitute the most common type of genetic variation in humans. SNPs introducing premature termination codons (PTCs), herein called X-SNPs, can alter the stability and function of transcripts and proteins and thus are considered to be biologically important. Initial studies suggested a strong selection against such variations/mutations. In this study, we undertook a genome-wide systematic screening to identify human X-SNPs using the dbSNP database. Our results demonstrated the presence of 28 X-SNPs from 28 genes with known minor allele frequencies. Eight X-SNPs (28.6 per cent) were predicted to cause transcript degradation by nonsense-mediated mRNA decay. Seventeen X-SNPs (60.7 per cent) resulted in moderate to severe truncation at the C-terminus of the proteins (deletion of >50 per cent of the amino acids). The majority of the X-SNPs (78.6 per cent) represent commonly occurring SNPs, by contrast with the rarely occurring disease-causing PTC mutations. Interestingly, X-SNPs displayed a non-uniform distribution across human populations: eight X-SNPs were reported to be prevalent across three different human populations, whereas six X-SNPs were found exclusively in one or two population(s). In conclusion, we have systematically investigated human SNPs introducing PTCs with respect to their possible biological consequences, distributions across different human populations and evolutionary aspects. We believe that the SNPs reported here are likely to affect gene/protein function, although their biological and evolutionary roles need to be further investigated.

Selective sweep mapping of genes with large phenotypic effects

Many domestic dog breeds have originated through fixation of discrete mutations by intense artificial selection. As a result of this process, markers in the proximity of genes influencing breed-defining traits will have reduced variation (a selective sweep) and will show divergence in allele frequency. Consequently, low-resolution genomic scans can potentially be used to identify regions containing genes that have a major influence on breed-defining traits. We model the process of breed formation and show that the probability of two or three adjacent marker loci showing a spurious signal of selection within at least one breed (i.e., Type I error or false-positive rate) is low if highly variable and moderately spaced markers are utilized. We also use simulations with selection to demonstrate that even a moderately spaced set of highly polymorphic markers (e.g., one every 0.8 cM) has high power to detect regions targeted by strong artificial selection in dogs. Further, we show that a gene responsible for black coat color in the Large Munsterlander has a 40-Mb region surrounding the gene that is very low in heterozygosity for microsatellite markers. Similarly, we survey 302 microsatellite markers in the Dachshund and find three linked monomorphic microsatellite markers all within a 10-Mb region on chromosome 3. This region contains the FGFR3 gene, which is responsible for achondroplasia in humans, but not in dogs. Consequently, our results suggest that the causative mutation is a gene or regulatory region closely linked to FGFR3.

Repositioning-dependent fate of duplicate genes

Gene duplication is the main source of evolutionary novelties. However, the problem with duplicates is that the purifying selection overlooks deleterious mutations in the redundant sequence, which therefore, instead of gaining a new function, often degrades into a functionless pseudogene. This risk of functional loss instead of gain is much higher for small populations of higher organisms with a slow and complex development. We propose that it is the epigenetic tissue/stage-complementary silencing of duplicates that makes them exposable to the purifying selection, thus saving them from pseudogenization and opening the way towards new function(s). Our genome-wide analyses of gene duplicates in several eukaryotic species combined with the phylogenetic comparison of vertebrate alpha- and beta-globin gene clusters strongly support this epigenetic complementation (EC) model. The distinctive condition for a new duplicate to survive by the EC mechanism seems to be its repositioning to an ectopic site, which is accompanied by changes in the rate and direction of mutagenesis. The most distinguished in this respect is the human genome. In this review, we extend and discuss the data on the EC- and repositioning-dependent fate of gene duplicates with the special emphasis on the problem of detecting brief postduplication period of adaptive evolution driven by positive selection. Accordingly, we propose a new CpG-focused measure of selection that is insensitive to translocation-caused biases in mutagenesis.

Simple evolutionary pathways to complex proteins

A recent paper in this journal has challenged the idea that complex adaptive features of proteins can be explained by known molecular, genetic, and evolutionary mechanisms. It is shown here that the conclusions of this prior work are an artifact of unwarranted biological assumptions, inappropriate mathematical modeling, and faulty logic. Numerous simple pathways exist by which adaptive multi-residue functions can evolve on time scales of a million years (or much less) in populations of only moderate size. Thus, the classical evolutionary trajectory of descent with modification is adequate to explain the diversification of protein functions.

Helicobacter pylori adhesion to carbohydrates

Adherence of bacterial pathogens to host tissues contributes to colonization and virulence and typically involves specific interactions between bacterial proteins called adhesins and cognate oligosaccharide (glycan) or protein motifs in the host that are used as receptors. A given pathogen may have multiple adhesins, each specific for a different set of receptors and, potentially, with different roles in infection and disease. This chapter provides strategies for identifying and analyzing host glycan receptors and the bacterial adhesins that exploit them as receptors, with particular reference to adherence of the gastric pathogen Helicobacter pylori.

Genome sequence, comparative analysis and haplotype structure of the domestic dog

Here we report a high-quality draft genome sequence of the domestic dog (Canis familiaris), together with a dense map of single nucleotide polymorphisms (SNPs) across breeds. The dog is of particular interest because it provides important evolutionary information and because existing breeds show great phenotypic diversity for morphological, physiological and behavioural traits. We use sequence comparison with the primate and rodent lineages to shed light on the structure and evolution of genomes and genes. Notably, the majority of the most highly conserved non-coding sequences in mammalian genomes are clustered near a small subset of genes with important roles in development. Analysis of SNPs reveals long-range haplotypes across the entire dog genome, and defines the nature of genetic diversity within and across breeds. The current SNP map now makes it possible for genome-wide association studies to identify genes responsible for diseases and traits, with important consequences for human and companion animal health.

Synonymous SNPs Provide Evidence for Selective Constraint on Human Exonic Splicing Enhancers

The human SNP database was used to detect selection on 238 hexamers previously identified as exonic splicing enhancers (ESEs). We compared the distribution of the 238 putative ESEs in biallelic and triallelic SNPs within five different functional categories of the SNP database: synonymous, nonsynonymous, introns, UTRs, and nongenic SNPs. Since true ESEs do not function outside of exons, SNPs that disrupt ESE motifs were expected to be more common in nonexonic portions of the genome. Our results supported this expectation: ESEs were least prevalent within synonymous SNPs and most common in nongenic SNPs. There were approximately 11% fewer ESEs within synonymous biallelic SNPs than expected under no selective constraint. We also compared the frequency of neutral SNPs, those where neither allele was an ESE, with deleterious SNPs, those where one or more alleles was an ESE, across the five different functional classes of SNPs. In comparison with the other functional classes of SNPs, synonymous SNPs contained an excess of neutral variants (+1.64% and +6.04% for biallelic and triallelic SNPs, respectively) and a dearth of deleterious variants (-13.11% and -52.39% for biallelic and triallelic SNPs, respectively). The observed patterns were consistent with purifying selection on the 238 hexamers to maintain their function as ESEs. However, in contrast to previous work, we did not find evidence for selection to maintain ESE function at nonsynonymous SNPs because selection at the protein level probably obscured any difference at the level of ESE function.

Rapid nonsynonymous evolution of the iron-sulfur protein in anthropoid primates

Cytochrome c (CYC) and 9 of the 13 subunits of cytochrome c oxidase (complex IV; COX) were previously shown to have accelerated rates of nonsynonymous substitution in anthropoid primates. Cytochrome b, the mtDNA encoded subunit of ubiquinol-cytochrome c reductase (complex III), also showed an accelerated nonsynonymous substitution rate in anthropoid primates but rate information about the nuclear encoded subunits of complex III has been lacking. We now report that phylogenetic and relative rates analysis of a nuclear encoded catalytically active subunit of complex III, the iron-sulfur protein (ISP), shows an accelerated rate of amino acid replacement similar to cytochrome b. Because both ISP and subunit 9, whose function is not directly related to electron transport, are produced by cleavage into two subunits of the initial translation product of a single gene, it is probable that these two subunits of complex III have essentially identical underlying rates of mutation. Nevertheless, we find that the catalytically active ISP has an accelerated rate of amino acid replacement in anthropoid primates whereas the catalytically inactive subunit 9 does not.

Chromosomal rearrangements are associated with higher rates of molecular evolution in mammals

Evolutionary rates are not uniformly distributed across the genome. Knowledge about the biological causes of this observation is still incomplete, but its exploration has provided valuable insight into the genomical, historical and demographical variables that influence rates of genetic divergence. Recent studies suggest a possible association between chromosomal rearrangements and regions of greater divergence, but evidence is limited and contradictory. Here, we test the hypothesis of a relationship between chromosomal rearrangements and higher rates of molecular evolution by studying the genomic distribution of divergence between 12,000 human-mouse orthologous genes. Our results clearly show that genes located in genomic regions that have been highly rearranged between the two species present higher rates of synonymous (0.7686 vs. 0.7076) and non-synonymous substitution (0.1014 vs. 0.0871), and that synonymous substitution rates are higher in genes close to the breakpoints of individual rearrangements. The many potential causes of such striking are discussed, particularly in the light of speciation models suggesting that chromosomal rearrangements may have contributed to some of the speciation processes along the human and mouse lineages. Still, there are other possible causes and further research is needed to properly explore them.

Comprehensive search for intra- and inter-specific sequence polymorphisms among coding envelope genes of retroviral origin found in the human genome: genes and pseudogenes

Background

The human genome carries a high load of proviral-like sequences, called Human Endogenous Retroviruses (HERVs), which are the genomic traces of ancient infections by active retroviruses. These elements are in most cases defective, but open reading frames can still be found for the retroviral envelope gene, with sixteen such genes identified so far. Several of them are conserved during primate evolution, having possibly been co-opted by their host for a physiological role.

Results

To characterize further their status, we presently sequenced 12 of these genes from a panel of 91 Caucasian individuals. Genomic analyses reveal strong sequence conservation (only two non synonymous Single Nucleotide Polymorphisms [SNPs]) for the two HERV-W and HERV-FRD envelope genes, i.e. for the two genes specifically expressed in the placenta and possibly involved in syncytiotrophoblast formation. We further show – using an ex vivo fusion assay for each allelic form – that none of these SNPs impairs the fusogenic function. The other envelope proteins disclose variable polymorphisms, with the occurrence of a stop codon and/or frameshift for most – but not all – of them. Moreover, the sequence conservation analysis of the orthologous genes that can be found in primates shows that three env genes have been maintained in a fully coding state throughout evolution including envW and envFRD.

Conclusion

Altogether, the present study strongly suggests that some but not all envelope encoding sequences are bona fide genes. It also provides new tools to elucidate the possible role of endogenous envelope proteins as susceptibility factors in a number of pathologies where HERVs have been suspected to be involved.

Parallel patterns of evolution in the genomes and transcriptomes of humans and chimpanzees

The determination of the chimpanzee genome sequence provides a means to study both structural and functional aspects of the evolution of the human genome. Here we compare humans and chimpanzees with respect to differences in expression levels and protein-coding sequences for genes active in brain, heart, liver, kidney, and testis. We find that the patterns of differences in gene expression and gene sequences are markedly similar. In particular, there is a gradation of selective constraints among the tissues so that the brain shows the least differences between the species whereas liver shows the most. Furthermore, expression levels as well as amino acid sequences of genes active in more tissues have diverged less between the species than have genes active in fewer tissues. In general, these patterns are consistent with a model of neutral evolution with negative selection. However, for X-chromosomal genes expressed in testis, patterns suggestive of positive selection on sequence changes as well as expression changes are seen. Furthermore, although genes expressed in the brain have changed less than have genes expressed in other tissues, in agreement with previous work we find that genes active in brain have accumulated more changes on the human than on the chimpanzee lineage.

Phosphorylation states of cell cycle and DNA repair proteins can be altered by the nsSNPs

Background

Phosphorylation is a reversible post-translational modification that affects the intrinsic properties of proteins, such as structure and function. Non-synonymous single nucleotide polymorphisms (nsSNPs) result in the substitution of the encoded amino acids and thus are likely to alter the phosphorylation motifs in the proteins.

Methods

In this study, we used the web-based NetPhos tool to predict candidate nsSNPs that either introduce or remove putative phosphorylation sites in proteins that act in DNA repair and cell cycle pathways.

Results

Our results demonstrated that a total of 15 nsSNPs (16.9%) were likely to alter the putative phosphorylation patterns of 14 proteins. Three of these SNPs (CDKN1A-S31R, OGG1-S326C, and XRCC3-T241M) have already found to be associated with altered cancer risk. We believe that this set of nsSNPs constitutes an excellent resource for further molecular and genetic analyses.

Conclusion

The novel systematic approach used in this study will accelerate the understanding of how naturally occurring human SNPs may alter protein function through the modification of phosphorylation mechanisms and contribute to disease susceptibility.

Hypervariable noncoding sequences in Saccharomyces cerevisiae

Compared to protein-coding sequences, the evolution of noncoding sequences and the selective constraints placed on these sequences is not well characterized. To compare the evolution of coding and noncoding sequences, we have conducted a survey for DNA polymorphism at five randomly chosen loci among a diverse collection of 81 strains of Saccharomyces cerevisiae. Average rates of both polymorphism and divergence are 40% lower at noncoding sites and 90% lower at nonsynonymous sites in comparison to synonymous sites. Although noncoding and coding sequences show substantial variability in ratios of polymorphism to divergence, two of the loci, MLS1 and PDR10, show a higher rate of polymorphism at noncoding compared to synonymous sites. The high rate of polymorphism is not accompanied by a high rate of divergence and is limited to a few small regions. These hypervariable regions include sites with three segregating bases at a single site and adjacent polymorphic sites. We show that this clustering of polymorphic sites is significantly greater than one would expect on the basis of the spacing between polymorphic fourfold degenerate sites. Although hypervariable noncoding sequences could result from selection on regulatory mutations, they could also result from transient mutational hotspots.

Autosomal recessive primary microcephaly (MCPH): a review of clinical, molecular, and evolutionary findings

Autosomal recessive primary microcephaly (MCPH) is a neurodevelopmental disorder. It is characterized by two principal features, microcephaly present at birth and nonprogressive mental retardation. The microcephaly is the consequence of a small but architecturally normal brain, and it is the cerebral cortex that shows the greatest size reduction. There are at least seven MCPH loci, and four of the genes have been identified: MCPH1, encoding Microcephalin; MCPH3, encoding CDK5RAP2; MCPH5, encoding ASPM; and MCPH6, encoding CENPJ. These findings are starting to have an impact on the clinical management of families affected with MCPH. Present data suggest that MCPH is the consequence of deficient neurogenesis within the neurogenic epithelium. Evolutionary interest in MCPH has been sparked by the suggestion that changes in the MCPH genes might also be responsible for the increase in brain size during human evolution. Indeed, evolutionary analyses of Microcephalin and ASPM reveal evidence for positive selection during human and great ape evolution. So an understanding of this rare genetic disorder may offer us significant insights into neurogenic mitosis and the evolution of the most striking differences between us and our closest living relatives: brain size and cognitive ability.

Candidate nsSNPs that can affect the functions and interactions of cell cycle proteins

Nonsynonymous single nucleotide polymorphisms (nsSNPs) alter the encoded amino acid sequence, and are thus likely to affect the function of the proteins, and represent potential disease-modifiers. There is an enormous number of nsSNPs in the human population, and the major challenge lies in distinguishing the functionally significant and potentially disease-related ones from the rest. In this study, we analyzed the genetic variations that can alter the functions and the interactions of a group of cell cycle proteins (n = 60) and the proteins interacting with them (n = 26) using computational tools. As a result, we extracted 249 nsSNPs from 77 cell cycle proteins and their interaction partners from public SNP databases. Only 31 (12.4%) of the nsSNPs were validated. The majority (64.5%) of the validated SNPs were rare (minor allele frequencies < 5%). Evolutionary conservation analysis using the SIFT tool suggested that 16.1% of the validated nsSNPs may disrupt the protein function. In addition, 58% of the validated nsSNPs were located in functional protein domains/motifs, which together with the evolutionary conservation analysis enabled us to infer possible biological consequences of the nsSNPs in our set. Our study strongly suggests the presence of naturally occurring genetic variations in the cell cycle proteins that may affect their interactions and functions with possible roles in complex human diseases, such as cancer.

Comparative primate genomics

With the completion of the human genome sequence and the advent of technologies to study functional aspects of genomes, molecular comparisons between humans and other primates have gained momentum. The comparison of the human genome to the genomes of species closely related to humans allows the identification of genomic features that set primates apart from other mammals and of features that set certain primates notably humans apart from other primates. In this article, we review recent progress in these areas with an emphasis on how comparative approaches may be used to identify functionally relevant features unique to the human genome.