Higher intensity of purifying selection on >90% of the human genes revealed by the intrinsic replacement mutation rates

For over 3 decades, the rate of replacement mutations has been assumed to be equal to, and estimated from, the rate of "strictly" neutral sequence divergence in noncoding regions and in silent-codon positions where mutations do not alter the amino acid encoded. This assumption is fundamental to estimating the fraction of harmful protein mutations and to identifying adaptive evolution at individual codons and proteins. We show that the assumption is not justifiable because a much larger fraction of codon positions is involved in hypermutable CpG dinucleotides as compared with the introns, leading to a higher expected replacement mutation rate per site in a vast majority of the genes. Consideration of this difference reveals a higher intensity of purifying natural selection than previously inferred in human genes. We also show that a much smaller number of genes are expected to be evolving with positive selection than that predicted using sequence divergence at intron and silent positions in the human genome. These patterns indicate the need for using new approaches for estimating rates of amino acid-altering mutations in order to find positively selected genes and codons in genomes that contain hypermutable CpG's.     

Long Simple Sequence Repeats in Host-Adapted Pathogens Localize Near Genes Encoding Antigens,Housekeeping Genes,and Pseudogenes

Simple sequence repeats (SSRs) in DNA sequences are tandem iterations of a single nucleotide or a short oligonucleotide. SSRs are subject to slipped-strand mutations and a common source of phase variation in bacteria and antigenic variation in pathogens. Significantly long SSRs are generally rare in prokaryotic genomes, and long SSRs composed of iterations of mono-, di-, tri-, and tetranucleotides are mostly restricted to host-adapted pathogens. We present new results concerning associations between long SSRs and genes related to different cellular functions in genomes of host-adapted pathogens. We found that in the majority of the analyzed genomes, at least some of the genes associated with SSRs encode potential antigens, which is expected if the primary function of SSRs is their contribution to antigenic variation. However, we also found a number of long SSRs associated with housekeeping genes, including rRNA and tRNA genes, genes encoding ribosomal proteins, amino acyl-tRNA synthetases, chaperones, and important metabolic enzymes. Many of these genes are probably essential and it is unlikely that they are phase-variable. Few statistically significant associations between SSRs and gene functional classifications were detected, suggesting that most long SSRs are not related to a particular cellular function or process. Long SSRs in Mycobacterium leprae are mostly associated with pseudogenes and may be contributing to gene loss following the adaptation to an obligate pathogenic lifestyle. We speculate that LSSRs may have played a similar role in genome reduction of other host-adapted pathogens.     

Analysis of distribution indicates diverse functions of simple sequence repeats in Mycoplasma genomes

Simple sequence repeats (SSRs) composed of extensive tandem iterations of a single nucleotide or a short oligonucleotide are rare in most bacterial genomes, but they are common among Mycoplasma. Some of these repeats act as contingency loci in association with families of surface antigens. By contraction or expansion during replication, these SSRs increase genetic variance of the population and facilitate avoidance of the immune response of the host. Occurrence and distribution of SSRs are analyzed in complete genomes of 11 Mycoplasma and 3 related Mollicutes in order to gain insights into functional and evolutionary diversity of the SSRs in Mycoplasma. The results revealed an unexpected variety of SSRs with respect to their distribution and composition and suggest that it is unlikely that all SSRs function as contingency loci or recombination hot spots. Various types of SSRs are most abundant in Mycoplasma hyopneumoniae, whereas Mycoplasma penetrans, Mycoplasma mobile, and Mycoplasma synoviae do not contain unusually long SSRs. Mycoplasma hyopneumoniae and Mycoplasma pulmonis feature abundant short adenine and thymine runs periodically spaced at 11 and 12 bp, respectively, which likely affect the supercoiling propensities of the DNA molecule. Physiological roles of long adenine and thymine runs in M. hyopneumoniae appear independent of location upstream or downstream of genes, unlike contingency loci that are typically located in protein-coding regions or upstream regulatory regions. Comparisons among 3 M. hyopneumoniae strains suggest that the adenine and thymine runs are rarely involved in genome rearrangements. The results indicate that the SSRs in the Mycoplasma genomes play diverse roles, including modulating gene expression as contingency loci, facilitating genome rearrangements via recombination, affecting protein structure and possibly protein-protein interactions, and contributing to the organization of the DNA molecule in the cell.     

A simple method for estimating the intensity of purifying selection in protein-coding genes.

We propose a method by which the intensity of purifying selection on a functional protein-coding gene is estimated by using three aligned homologous sequences: a processed pseudogene (psi), a functional paralog from the same species (g), and a functional ortholog from a different species (o). For each such trio, we calculate the numbers of nucleotide substitutions along the branches leading to psi and g, i.e., K psi and K(g). If we assume that the mutation rates are the same in the genes and the pseudogenes and that mutations occurring in a pseudogene do not affect the fitness of the organism, we can show that the fraction of mutations that are selectively neutral, fg, is equal to the ratio K(g)/K psi. Since advantageous mutations occur only very rarely, such that they do not contribute significantly to the rate of molecular evolution, the fraction of deleterious mutations that are subject to purifying selection is 1-fg. Therefore, the K(g)/K psi ratio can be used directly to estimate the intensity of purifying selection, thereby isolating its effects on the rate of evolution from those of mutation. We compared the selection intensities of 12 orthologous protein-coding pairs from humans and murids. As expected, the fraction of mutations that are subject to purifying selection is strongest in the second codon position and weakest in the third. Interestingly, the mean fractions of effectively neutral mutations in the third codon position were only 41% and 42% for murids and humans, respectively, indicating that many synonymous mutations are subject to selective constraint. In several orthologous genes, we found that the intensity of purifying selection is very different between murid and human orthologous genes. There was no statistically significant difference in overall intensity of purifying selection between humans and murids. Thus, purifying selection does not seem to be an important factor contributing to the observed differences in the rates of evolution between these two taxa.     

The effect of multifunctionality on the rate of evolution in yeast

Multifunctional genes are expected to evolve at lower rates because mutations in such genes that improve one function might often have deleterious effects on other functions. Here we tested for an association between multifunctionality and evolutionary rates in genes of Saccharomyces cerevisiae, and we find a highly significant negative correlation between the number of biological processes in which a gene is involved in and its rate of evolution. However, the magnitude of this effect is small, and the results do not support the notion that multifunctionality limits a gene's rate of evolution.     

Comparative Analysis of SSR Markers Developed in Exon,Intron, and Intergenic Regions and Distributed in Regions Controlling Fruit Quality Traits in <Emphasis Type="Italic">Prunus</Emphasis> Species: Genetic Diversity and Association Studies

Simple sequence repeats (SSRs) are genome domains located in both coding and non-coding regions in eukaryotic genomes. Although SSRs are often characterized by low polymorphism, their DNA-flanking sequences could be a useful source of DNA markers, which could help in genetic studies and breeding because they are associated with genes that control traits of interest. In this study, 56 genotypes from different Prunus species were used, including peach, apricot, plum, and almond (already phenotyped for several agronomical traits, including self-compatibility, flowering and ripening time, fruit type, skin and flesh color, and shell hardness). These Prunus genotypes were molecularly characterized using 28 SSR markers developed in exons, introns, and intergenic regions. All these genes were located in specific regions where quantitative trait loci (QTLs) for certain fruit quality traits were also located, including flowering and ripening times and fruit flesh and skin color. A sum of 309 SSR alleles were identified in the whole panel of analyzed cultivars, with expected heterozygosity values of 0.61 (upstream SSRs), 0.17 (exonic SSRs), 0.65 (intronic SSRs), and 0.58 (downstream SSRs). These values prove the low level of polymorphism of the exonic (gene-coding regions) markers. Cluster and structural analysis based on SSR data clearly differentiated the genotypes according to either specie (for the four species) and pedigree (apricot) or geographic origin (Japanese plum). In addition, some SSR markers mainly developed in intergenic regions could be associated with genes that control traits of interest in breeding and could therefore help in marker-assisted breeding. These findings highlight the importance of using molecular markers able to discriminate between the functional roles of the gene allelic variants.     

Comparative genomics reveals high rates of horizontal transfer and strong purifying selection on rhizobial symbiosis genes

Horizontal transfer (HT) alters the repertoire of symbiosis genes in rhizobial genomes and may play an important role in the on-going evolution of the rhizobia–legume symbiosis. To gain insight into the extent of HT of symbiosis genes with different functional roles (nodulation, N-fixation, host benefit and rhizobial fitness), we conducted comparative genomic and selection analyses of the full-genome sequences from 27 rhizobial genomes. We find that symbiosis genes experience high rates of HT among rhizobial lineages but also bear signatures of purifying selection (low Ka : Ks). HT and purifying selection appear to be particularly strong in genes involved in initiating the symbiosis (e.g. nodulation) and in genome-wide association candidates for mediating benefits provided to the host. These patterns are consistent with rhizobia adapting to the host environment through the loss and gain of symbiosis genes, but not with host-imposed positive selection driving divergence of symbiosis genes through recurring bouts of positive selection.     

Evolutionary conservation suggests a regulatory function of AUG triplets in 5'-UTRs of eukaryotic genes

By comparing sequences of human, mouse and rat orthologous genes, we show that in 5′-untranslated regions (5′-UTRs) of mammalian cDNAs but not in 3′-UTRs or coding sequences, AUG is conserved to a significantly greater extent than any of the other 63 nt triplets. This effect is likely to reflect, primarily, bona fide evolutionary conservation, rather than cDNA annotation artifacts, because the excess of conserved upstream AUGs (uAUGs) is seen in 5′-UTRs containing stop codons in-frame with the start AUG and many of the conserved AUGs are found in different frames, consistent with the location in authentic non-coding sequences. Altogether, conserved uAUGs are present in at least 20–30% of mammalian genes. Qualitatively similar results were obtained by comparison of orthologous genes from different species of the yeast genus Saccharomyces. Together with the observation that mammalian and yeast 5′-UTRs are significantly depleted in overall AUG content, these findings suggest that AUG triplets in 5′-UTRs are subject to the pressure of purifying selection in two opposite directions: the uAUGs that have no specific function tend to be deleterious and get eliminated during evolution, whereas those uAUGs that do serve a function are conserved. Most probably, the principal role of the conserved uAUGs is attenuation of translation at the initiation stage, which is often additionally regulated by alternative splicing in the mammalian 5′-UTRs. Consistent with this hypothesis, we found that open reading frames starting from conserved uAUGs are significantly shorter than those starting from non-conserved uAUGs, possibly, owing to selection for optimization of the level of attenuation.     

基于表达序列标签的微卫星标记(EST-SSRs)研究进展

作为一种新型分子标记,EST-SSR来自表达基因,因而除具备传统基因组来源的SSR标记所有优势外,可能与基因功能表达具有直接或间接关系,从而强化了SSR标记在遗传研究中的应用.本文综述了近几年来EST-SSR用于遗传图谱构建、基因定位、比较基因组学及重要基因筛选和发掘等方面的研究.     

MUSTANG is a novel family of domesticated transposase genes found in diverse angiosperms

While transposons have traditionally been viewed as genomic parasites or "junk DNA," the discovery of transposon-derived host genes has fueled an ongoing debate over the evolutionary role of transposons. In particular, while mobility-related open reading frames have been known to acquire host functions, the contribution of these types of events to the evolution of genes is not well understood. Here we report that genome-wide searches for Mutator transposase-derived host genes in Arabidopsis thaliana (Columbia-0) and Oryza sativa ssp. japonica (cv. Nipponbare) (domesticated rice) identified 121 sequences, including the taxonomically conserved MUSTANG1. Syntenic MUSTANG1 orthologs in such varied plant species as rice, poplar, Arabidopsis, and Medicago truncatula appear to be under purifying selection. However, despite the evidence of this pathway of gene evolution, MUSTANG1 belongs to one of only two Mutator-like gene families with members in both monocotyledonous and dicotyledonous plants, suggesting that Mutator-like elements seldom evolve into taxonomically widespread host genes.     

Intracellular symbionts and the evolution of uniparental cytoplasmic inheritance.

Uniparental inheritance of cytoplasmic elements is widespread among eukaryotic organisms and is achieved by a diverse range of mechanisms. This paper shows that the cytoplasmic genetic system would be expected to evolve towards uniparental inheritance, given the existence of deleterious symbionts capable of invading the host cytoplasm together with nuclear genes that lead to the elimination of cytoplasmic elements from one of the gamete types. The reason for this is that, under biparental inheritance, foreign symbionts with strong deleterious effects are able to spread through host populations. A nuclear modifier gene which leads to the loss of cytoplasmic elements from one gamete type gains a net advantage as a symbiont spreads, because the modifier sometimes gives rise to a symbiont-free zygote. Insofar as small gametes reduce the rate of symbiont transmission to the zygote, modifier genes causing small gamete size would tend to accumulate, so that cytoplasmic inheritance would become associated with maternal rather than paternal gametes. Once uniparental inheritance predominates in the host population, the population is protected from invasions by a large class of harmful symbionts, but at the same time those symbionts that benefit their hosts are still able to increase in frequency.     

Development and Application of Microsatellites in Candidate Genes Related to Wood Properties in the Chinese White Poplar (Populus tomentosa Carr.)

Gene-derived simple sequence repeats (genic SSRs), also known as functional markers, are often preferred over random genomic markers because they represent variation in gene coding and/or regulatory regions. We characterized 544 genic SSR loci derived from 138 candidate genes involved in wood formation, distributed throughout the genome of Populus tomentosa, a key ecological and cultivated wood production species. Of these SSRs, three-quarters were located in the promoter or intron regions, and dinucleotide (59.7%) and trinucleotide repeat motifs (26.5%) predominated. By screening 15 wild P. tomentosa ecotypes, we identified 188 polymorphic genic SSRs with 861 alleles, 2–7 alleles for each marker. Transferability analysis of 30 random genic SSRs, testing whether these SSRs work in 26 genotypes of five genus Populus sections (outgroup, Salix matsudana), showed that 72% of the SSRs could be amplified in Turanga and 100% could be amplified in Leuce. Based on genotyping of these 26 genotypes, a neighbour-joining analysis showed the expected six phylogenetic groupings. In silico analysis of SSR variation in 220 sequences that are homologous between P. tomentosa and Populus trichocarpa suggested that genic SSR variations between relatives were predominantly affected by repeat motif variations or flanking sequence mutations. Inheritance tests and single-marker associations demonstrated the power of genic SSRs in family-based linkage mapping and candidate gene-based association studies, as well as marker-assisted selection and comparative genomic studies of P. tomentosa and related species.     

Positive and negative selection on mammalian Y chromosomes

Y chromosomes are genetically degenerate in most organisms studied. The loss of genes from Y chromosomes is thought to be due to the inefficiency of purifying selection in nonrecombining regions, which leads to the accumulation of deleterious mutations via the processes of hitchhiking, background selection, and Muller's ratchet. As the severity of these processes depends on the number of functional genes linked together on the nonrecombining Y, it is not clear whether these processes are still at work on the old, gene-poor mammalian Y chromosomes. If purifying selection is indeed less efficient in the Y-linked, compared to the X-linked genes, deleterious nonsynonymous substitutions are expected to accumulate faster on the Y chromosome. However, positive selection on Y-linked genes could also increase the rate of amino acid-changing substitutions. Thus, the previous reports of an elevated nonsynonymous substitution rate in Y-linked genes are still open to interpretation. Here, we report evidence for positive selection in two out of three studied mammalian Y-linked genes, suggesting that adaptive Darwinian evolution may be common on mammalian Y chromosomes. Taking positive selection into account, we demonstrate that purifying selection is less efficient in mammalian Y-linked genes compared to their X-linked homologues, suggesting that these genes continue to degenerate.     

Neutral evolution of robustness in Drosophila microRNA precursors

Mutational robustness describes the extent to which a phenotype remains unchanged in the face of mutations. Theory predicts that the strength of direct selection for mutational robustness is at most the magnitude of the rate of deleterious mutation. As far as nucleic acid sequences are concerned, only long sequences in organisms with high deleterious mutation rates and large population sizes are expected to evolve mutational robustness. Surprisingly, recent studies have concluded that molecules that meet none of these conditions--the microRNA precursors (pre-miRNAs) of multicellular eukaryotes--show signs of selection for mutational and/or environmental robustness. To resolve the apparent disagreement between theory and these studies, we have reconstructed the evolutionary history of Drosophila pre-miRNAs and compared the robustness of each sequence to that of its reconstructed ancestor. In addition, we "replayed the tape" of pre-miRNA evolution via simulation under different evolutionary assumptions and compared these alternative histories with the actual one. We found that Drosophila pre-miRNAs have evolved under strong purifying selection against changes in secondary structure. Contrary to earlier claims, there is no evidence that these RNAs have been shaped by either direct or congruent selection for any kind of robustness. Instead, the high robustness of Drosophila pre-miRNAs appears to be mostly intrinsic and likely a consequence of selection for functional structures.     

Low microsatellite frequencies in neuron and brain expressed microRNAs

The locations of microsatellites in mammalian genomes are restricted by purifying selection in a number of ways. For example, with the exception of some trinucleotide repeats they are excluded from protein coding regions of genomes because of their tendency to cause frameshift mutations. Here we investigate whether purifying selection might affect the types and frequencies of microsatellites in microRNA (miRNA). We concentrate on miRNAs expressed in neurons and the brain (NB-miRNAs) as microsatellites in these genes might give rise to similar effects as disease-causing repeats in protein coding genes. We show that in human miRNAs in general AG and AT microsatellites are reduced in frequency compared to AC repeats and that NB-miRNA genes contain significantly fewer microsatellites than expected from frequencies of microsatellites in other miRNA genes. NB-miRNAs show lower levels of sequence divergence in comparisons of human-macaque orthologues and more often have detectable orthologues in non-human mammals than non-NB-miRNAs. This suggests that microsatellites in miRNAs may indeed be constrained by purifying selection and that the strength of this selection may differ between NB-miRNAs and non-NB-miRNAs. We identify a number of ways in which the potential disruption of pre-miRNA secondary structure might result in purifying selection. However other, non-selective forces could also play a role in generating the biases observed in miRNA microsatellites.     

The control of copy number of IS6110 in Mycobacterium tuberculosis

Insertion sequence (IS) elements are bacterial genes that are able to transpose to different locations in the genome. These elements are often used in molecular epidemiology as genetic markers that track the spread of pathogens. Transposable elements have frequently been described as "selfish DNA" because they facilitate their own transposition, causing damage when they insert into coding regions, while contributing little if anything to the bacterial host. According to this hypothesis, the expansion of copy number of insertion sequences is opposed by negative selection against high copy numbers. From an alternative point of view, we might expect IS elements to intrinsically regulate transposition within cells, thereby limiting damage to their bacterial host. Here, we report evidence that the copy number of IS6110 in Mycobacterium tuberculosis is controlled by selection against the element. We first construct 12 different models of marker change resulting from a combination of possible transposition functions and selective regimes. We then compute the Akaike Information Criterion for each model to identify the models that best explain data consisting of serial isolates of M. tuberculosis genotyped with IS6110. We find that the best performing models all include selection against the accumulation of copies. Specifically, our analysis points to the interaction of separate copies of the element causing lethal effects. We discuss the implications of these findings for genome evolution and molecular epidemiology.     

Ontogeny and phylogeny: molecular signatures of selection, constraint, and temporal pleiotropy in the development ofDrosophila

Background  

Karl Ernst Von Baer noted that species tend to show greater morphological divergence in later stages of development when compared to earlier stages. Darwin originally interpreted these observations via a selectionist framework, suggesting that divergence should be greatest during ontogenic stages in which organisms experienced varying 'conditions of existence' and opportunity for differential selection. Modern hypotheses have focused on the notion that genes and structures involved in early development will be under stronger purifying selection due to the deleterious pleiotropic effects of mutations propagating over the course of ontogeny, also known as the developmental constraint hypothesis.