Strong Selective Constraints Research Articles

Abundance and ribotypes of phosphate-solubilizing bacteria in Argentinean agricultural soils under no-till management

Phosphate-solubilizing bacteria (PSB) abundance and ribotypes were examined in the top soil (0–10 cm) of agricultural fields under no-till management sampled in two seasons, summer (February) and late winter (September), in 2010. No-till plots under sustainable agricultural practices (intense crop rotation) or under non-sustainable practices (soybean monocropping) were sampled at four different locations as replicates in a 400-km west–east transect in the most productive agricultural region of Argentina. Natural grasslands were selected close to the cultivated fields for comparative purposes. Culturable heterotrophic bacteria (CHB) were enumerated on nutrient agar plates, and PSB were counted on NBRIP agar plates containing Ca3(PO4)2. The PSB community structure was explored by ribotyping (16S rDNA PCR-restriction fragment length polymorphism). Quantitatively, data showed that the number of CHB and PSB per gram of dry soil was not statistically different among sampled sites, soil management programs or seasons. Qualitatively, ribotyping showed that the most abundant PSB species differed in their fingerprinting patterns among geographical sites, which suggests that local soil conditions impose strong selective constraints. The comparative analysis of PSB ribotypes revealed seasonal differences among February and September isolates for all sampling sites.

Divergent evolutionary rates in vertebrate and mammalian specific conserved non-coding elements (CNEs) in echolocating mammals.

BackgroundThe majority of DNA contained within vertebrate genomes is non-coding, with a certain proportion of this thought to play regulatory roles during development. Conserved Non-coding Elements (CNEs) are an abundant group of putative regulatory sequences that are highly conserved across divergent groups and thus assumed to be under strong selective constraint. Many CNEs may contain regulatory factor binding sites, and their frequent spatial association with key developmental genes – such as those regulating sensory system development – suggests crucial roles in regulating gene expression and cellular patterning. Yet surprisingly little is known about the molecular evolution of CNEs across diverse mammalian taxa or their role in specific phenotypic adaptations. We examined 3,110 vertebrate-specific and ~82,000 mammalian-specific CNEs across 19 and 9 mammalian orders respectively, and tested for changes in the rate of evolution of CNEs located in the proximity of genes underlying the development or functioning of auditory systems. As we focused on CNEs putatively associated with genes underlying the development/functioning of auditory systems, we incorporated echolocating taxa in our dataset because of their highly specialised and derived auditory systems.ResultsPhylogenetic reconstructions of concatenated CNEs broadly recovered accepted mammal relationships despite high levels of sequence conservation. We found that CNE substitution rates were highest in rodents and lowest in primates, consistent with previous findings. Comparisons of CNE substitution rates from several genomic regions containing genes linked to auditory system development and hearing revealed differences between echolocating and non-echolocating taxa. Wider taxonomic sampling of four CNEs associated with the homeobox genes Hmx2 and Hmx3 – which are required for inner ear development – revealed family-wise variation across diverse bat species. Specifically within one family of echolocating bats that utilise frequency-modulated echolocation calls varying widely in frequency and intensity high levels of sequence divergence were found.ConclusionsLevels of selective constraint acting on CNEs differed both across genomic locations and taxa, with observed variation in substitution rates of CNEs among bat species. More work is needed to determine whether this variation can be linked to echolocation, and wider taxonomic sampling is necessary to fully document levels of conservation in CNEs across diverse taxa.Electronic supplementary materialThe online version of this article (doi:10.1186/s12862-014-0261-5) contains supplementary material, which is available to authorized users.

Nucleotide sequence conservation in paramyxoviruses; the concept of codon constellation

The stability and conservation of the sequences of RNA viruses in the field and the high error rates measured in vitro are paradoxical. The field stability indicates that there are very strong selective constraints on sequence diversity. The nature of these constraints is discussed. Apart from constraints on variation in cis-acting RNA and the amino acid sequences of viral proteins, there are other ones relating to the presence of specific dinucleotides such CpG and UpA as well as the importance of RNA secondary structures and RNA degradation rates. Recent other constraints identified in other RNA viruses, such as effects of secondary RNA structure on protein folding or modification of cellular tRNA complements, are also discussed. Using the family Paramyxoviridae, I show that the codon usage pattern (CUP) is (i) specific for each virus species and (ii) that it is markedly different from the host - it does not vary even in vaccine viruses that have been derived by passage in a number of inappropriate host cells. The CUP might thus be an additional constraint on variation, and I propose the concept of codon constellation to indicate the informational content of the sequences of RNA molecules relating not only to stability and structure but also to the efficiency of translation of a viral mRNA resulting from the CUP and the numbers and position of rare codons.

Elucidating the genotype-phenotype relationships and network perturbations of human shared and specific disease genes from an evolutionary perspective.

To date, numerous studies have been attempted to determine the extent of variation in evolutionary rates between human disease and nondisease (ND) genes. In our present study, we have considered human autosomal monogenic (Mendelian) disease genes, which were classified into two groups according to the number of phenotypic defects, that is, specific disease (SPD) gene (one gene: one defect) and shared disease (SHD) gene (one gene: multiple defects). Here, we have compared the evolutionary rates of these two groups of genes, that is, SPD genes and SHD genes with respect to ND genes. We observed that the average evolutionary rates are slow in SHD group, intermediate in SPD group, and fast in ND group. Group-to-group evolutionary rate differences remain statistically significant regardless of their gene expression levels and number of defects. We demonstrated that disease genes are under strong selective constraint if they emerge through edgetic perturbation or drug-induced perturbation of the interactome network, show tissue-restricted expression, and are involved in transmembrane transport. Among all the factors, our regression analyses interestingly suggest the independent effects of 1) drug-induced perturbation and 2) the interaction term of expression breadth and transmembrane transport on protein evolutionary rates. We reasoned that the drug-induced network disruption is a combination of several edgetic perturbations and, thus, has more severe effect on gene phenotypes.

Introgression of mitochondrial DNA promoted by natural selection in the Japanese pipistrelle bat (Pipistrellus abramus).

Introgression of mitochondrial DNA (mtDNA) between closely related taxa can be promoted by either neutral processes or natural selection. Since mitochondrial gene-encoded proteins play critical roles in oxidative metabolism, mtDNA genes are commonly considered to experience strong selective constraint. However, metabolic requirements vary across climatic and ecological gradients, thus modifying potential selective pressures acting on mtDNA genes. Here we conducted tests to detect adaptive evolution occurring in two mtDNA genes (Cytb and ND5) in individuals of Japanese pipistrelle bat (Pipistrellus abramus) across the mainland of China and Hainan Island. Nuclear DNA markers identified two clades in both the mainland and Hainan Island populations, whereas each of these regions had a specific mtDNA clade. This cyto-nuclear discordance is most likely caused by introgression of the mtDNA by ruling out two other alternative scenarios (incomplete lineage sorting and sex-biased gene flow). Although population-based analyses revealed purifying selection acting on Cytb and neutrality in ND5, multiple nonsynonymous substitutions in both Cytb and ND5 were suggested to have been caused by positive selection by a divergence-based analysis. Our study supports the view that molecular adaptation can occur at genes under strong purifying selection if nonsynonymous substitutions cause radical changes in the physicochemical properties of amino acids.

Engineering microorganisms based on molecular evolutionary analysis: a succinate production case study.

Evolution has resulted in thousands of species possessing similar metabolic enzymes with identical functions that are, however, regulated by different mechanisms. It is thus difficult to select optimal gene to engineer novel or manipulated metabolic pathways. Here, we tested the ability of molecular evolutionary analysis to identify appropriate genes from various species. We calculated the fraction of synonymous substitution and the effective number of codons (ENC) for nine genes stemming from glycolysis. Our research indicated that an enzyme gene with a stronger selective constraint in synonymous sites would mainly regulate corresponding reaction flux through altering the concentration of the protein, whereas those with a more relaxed selective constraint would primarily affect corresponding reaction flux by changing kinetic properties of the enzyme. Further, molecular evolutionary analysis was investigated for three types of genes involved in succinate precursor supply by catalysis of pyruvate. In this model, overexpression of Corynebacterium glutamicum pyc should result in greater conversion of pyruvate. Succinate yields in two Escherichia coli strains that overexpressed each of the three types of genes supported the molecular evolutionary analysis. This approach may thus provide an alternative strategy for selecting genes from different species for metabolic engineering and synthetic biology.

Negative purifying selection drives prion and doppel protein evolution.

The prion protein (PrP) when misfolded into the pathogenic conformer PrP(Sc) is the major causative agent of several lethal transmissible spongiform encephalopathies in mammals. Studies of evolutionary pressure on the corresponding gene using different datasets have yielded conflicting results. In addition, putative PrP or PrP interacting partners with strong similarity to PrP such as the doppel protein have not been examined to determine if the same evolutionary mechanisms apply to prion paralogs or if there are coselected sites that might indicate how and where the proteins interact. We examined several taxonomic groups that contain model organisms of prion diseases focusing on primates, bovids, and an expanded dataset of rodents for selection pressure on the prion gene (PRNP) and doppel gene (PRND) individually and for coevolving sites within. Overall, the results clearly indicate that both proteins are under strong selective constraints with relaxed selection on amino acid residues connecting α-helices 1 and 2.

Remarkable selective constraints on exonic dinucleotide repeats.

Long dinucleotide repeats found in exons present a substantial mutational hazard: mutations at these loci occur often and generate frameshifts. Here, we provide clear and compelling evidence that exonic dinucleotides experience strong selective constraint. In humans, only 18 exonic dinucleotides have repeat lengths greater than six, which contrasts sharply with the genome-wide distribution of dinucleotides. We genotyped each of these dinucleotides in 200 humans from eight 1000 Genomes Project populations and found a near-absence of polymorphism. More remarkably, divergence data demonstrate that repeat lengths have been conserved across the primate phylogeny in spite of what is likely considerable mutational pressure. Coalescent simulations show that even a very low mutation rate at these loci fails to explain the anomalous patterns of polymorphism and divergence. Our data support two related selective constraints on the evolution of exonic dinucleotides: a short-term intolerance for any change to repeat length and a long-term prevention of increases to repeat length. In general, our results implicate purifying selection as the force that eliminates new, deleterious mutants at exonic dinucleotides. We briefly discuss the evolution of the longest exonic dinucleotide in the human genome--a 10 x CA repeat in fibroblast growth factor receptor-like 1 (FGFRL1)--that should possess a considerably greater mutation rate than any other exonic dinucleotide and therefore generate a large number of deleterious variants.

Rapid birth-and-death evolution of imprinted snoRNAs in the Prader-Willi syndrome locus: implications for neural development in Euarchontoglires.

Imprinted small nucleolar RNAs (snoRNAs) are only found in eutherian genomes and closely related to brain functions. A complex human neurological disease, Prader-Willi syndrome (PWS), is primarily attributed to the deletion of imprinted snoRNAs in chromosome 15q11-q13. Here we investigated the snoRNA repertoires in the PWS locus of 12 mammalian genomes and their evolution processes. A total of 613 imprinted snoRNAs were identified in the PWS homologous loci and the gene number was highly variable across lineages, with a peak in Euarchontoglires. Lineage-specific gene gain and loss events account for most extant genes of the HBII-52 (SNORD115) and the HBII-85 (SNORD116) gene family, and remarkable high gene-birth rates were observed in the primates and the rodents. Meanwhile, rapid sequence substitution occurred only in imprinted snoRNA genes, rather than their flanking sequences or the protein-coding genes located in the same imprinted locus. Strong selective constraints on the functional elements of these imprinted snoRNAs further suggest that they are subjected to birth-and-death evolution. Our data suggest that the regulatory role of HBII-52 on 5-HT2CR pre-mRNA might originate in the Euarchontoglires through adaptive process. We propose that the rapid evolution of PWS-related imprinted snoRNAs has contributed to the neural development of Euarchontoglires.

Phenotypic Characterization of Salmonella RyhB-1 Mutations that Modulate Target Regulation

Recently, we reported the global regulatory roles of the two small RNAs in Salmonella typhimurium, ryhB-1 and ryhB-2. However, the genetic basis of the sRNA-target interactions remains unknown. To identify the nucleotides of RyhB-1 that are functionally important for its regulatory actions, we introduced random single-point mutations into ryhB-1 gene on the chromosome of Salmonella typhimurium carrying a sodB-lacZ translational fusion by an error-prone PCR method. We reasoned that mutants expressing variant RyhB-1 with weakened interaction with sodB transcript would produce a higher level of SodB when compared to wild type, leading to darker blue colonies on X-gal agar plates. Five mutants displaying a significant increase in β-galactosidase activity under the condition inducing RyhB-1 expression were isolated. Quantitative real-time PCR analysis showed that the expression levels of eight target mRNAs in these five mutants were significantly changed when compared to the parent strain. Interestingly, two mutations affected growth and cell survival under H2O2-stressed conditions. The results suggest that there are strong selective constraints against mutational changes in ryhB-1 gene sequence, leading to high levels of nucleotide conservation in ryhB-1 gene sequences among the genus of Salmonella.

Sequence Diversity of Pan troglodytes Subspecies and the Impact of WFDC6 Selective Constraints in Reproductive Immunity

Recent efforts have attempted to describe the population structure of common chimpanzee, focusing on four subspecies: Pan troglodytes verus, P. t. ellioti, P. t. troglodytes, and P. t. schweinfurthii. However, few studies have pursued the effects of natural selection in shaping their response to pathogens and reproduction. Whey acidic protein (WAP) four-disulfide core domain (WFDC) genes and neighboring semenogelin (SEMG) genes encode proteins with combined roles in immunity and fertility. They display a strikingly high rate of amino acid replacement (dN/dS), indicative of adaptive pressures during primate evolution. In human populations, three signals of selection at the WFDC locus were described, possibly influencing the proteolytic profile and antimicrobial activities of the male reproductive tract. To evaluate the patterns of genomic variation and selection at the WFDC locus in chimpanzees, we sequenced 17 WFDC genes and 47 autosomal pseudogenes in 68 chimpanzees (15 P. t. troglodytes, 22 P. t. verus, and 31 P. t. ellioti). We found a clear differentiation of P. t. verus and estimated the divergence of P. t. troglodytes and P. t. ellioti subspecies in 0.173 Myr; further, at the WFDC locus we identified a signature of strong selective constraints common to the three subspecies in WFDC6—a recent paralog of the epididymal protease inhibitor EPPIN. Overall, chimpanzees and humans do not display similar footprints of selection across the WFDC locus, possibly due to different selective pressures between the two species related to immune response and reproductive biology.

Conserved Proteins Are Fragile

Levels of selective constraint vary among proteins. Although strong constraint on a protein is often attributed to its functional importance, evolutionary rate may also be limited if a protein is fragile, such that a large proportion of amino acid replacements reduce its fitness. To determine the relative contributions of essentiality and fragility to selective constraint, we compared relationships of selection against nonsense mutations (snon) and selection against missense mutations (smis) to protein sequence conservation (Ka). As expected, snon is greater than smis; however, the correlation between smis and Ka is nearly three times stronger than the correlation between snon and Ka. Moreover, examination of relationships to gene expression level, tissue specificity, and number of protein-protein interactions shows that smis is more strongly correlated than snon to all three measures of biological function. Thus, our analysis reveals that slowly evolving proteins are under strong selective constraint primarily because they are fragile, and that this association likely exists because allowing a protein to function improperly, rather than removing it from a biological network, can negatively affect the functions of other molecules it interacts with and their downstream products.

Characterization of pseudogenes in members of the order Frankineae

Pseudogenes are defined as non-functional relatives of genes whose protein-coding abilities are lost and are no longer expressed within cells. They are an outcome of accumulation of mutations within a gene whose end product is not essential for survival. Proper investigation of the procedure of pseudogenization is relevant for estimating occurrence of duplications in genomes. Frankineae houses an interesting group of microorganisms, carving a niche in the microbial world. This study was undertaken with the objective of determining the abundance of pseudogenes, understanding strength of purifying selection, investigating evidence of pseudogene expression, and analysing their molecular nature, their origin, evolution and deterioration patterns amongst domain families. Investigation revealed the occurrence of 956 core pFAM families sharing common characteristics indicating co-evolution. WD40, Rve_3, DDE_Tnp_IS240 and phage integrase core domains are larger families, having more pseudogenes, signifying a probability of harmful foreign genes being disabled within transposable elements. High selective pressure depicted that gene families rapidly duplicating and evolving undoubtedly facilitated creation of a number of pseudogenes in Frankineae. Codon usage analysis between protein-coding genes and pseudogenes indicated a wide degree of variation with respect to different factors. Moreover, the majority of pseudogenes were under the effect of purifying selection. Frankineae pseudogenes were under stronger selective constraints, indicating that they were functional for a very long time and became pseudogenes abruptly. The origin and deterioration of pseudogenes has been attributed to selection and mutational pressure acting upon sequences for adapting to stressed soil environments.

Evolution von Genomen am Beispiel pathogener Pilze

The evolutionary interplay of fungal virulence and host resistance imposes strong selective constraints that shape fungal pathogens genome structure. Genomes of pathogenic fungi are highly diverse in architecture and gene content. Moreover, analyses of these genomes unveiled pivotal insights into structural diversity and evolutionary mechanisms of genomes as well as into functional repertoires and speciation processes of pathogens. Here, we highlight key features and evolutionary mechanisms of fungal pathogen genomes.

Different selective pressures shape the evolution of Toll-like receptors in human and African great ape populations

The study of the genetic and selective landscape of immunity genes across primates can provide insight into the existing differences in susceptibility to infection observed between human and non-human primates. Here, we explored how selection has driven the evolution of a key family of innate immunity receptors, the Toll-like receptors (TLRs), in African great ape species. We sequenced the 10 TLRs in various populations of chimpanzees and gorillas, and analysed these data jointly with a human data set. We found that purifying selection has been more pervasive in great apes than in humans. Furthermore, in chimpanzees and gorillas, purifying selection has targeted TLRs irrespectively of whether they are endosomal or cell surface, in contrast to humans where strong selective constraints are restricted to endosomal TLRs. These observations suggest important differences in the relative importance of TLR-mediated pathogen sensing, such as that of recognition of flagellated bacteria by TLR5, between humans and great apes. Lastly, we used a population genetics-phylogenetics method that jointly analyses polymorphism and divergence data to detect fine-scale variation in selection pressures at specific codons within TLR genes. We identified different codons at different TLRs as being under positive selection in each species, highlighting that functional variation at these genes has conferred a selective advantage in immunity to infection to specific primate species. Overall, this study showed that the degree of selection driving the evolution of TLRs has largely differed between human and non-human primates, increasing our knowledge on their respective biological contribution to host defence in the natural setting.

Strong Purifying Selection at Synonymous Sites in D. melanogaster

Synonymous sites are generally assumed to be subject to weak selective constraint. For this reason, they are often neglected as a possible source of important functional variation. We use site frequency spectra from deep population sequencing data to show that, contrary to this expectation, 22% of four-fold synonymous (4D) sites in Drosophila melanogaster evolve under very strong selective constraint while few, if any, appear to be under weak constraint. Linking polymorphism with divergence data, we further find that the fraction of synonymous sites exposed to strong purifying selection is higher for those positions that show slower evolution on the Drosophila phylogeny. The function underlying the inferred strong constraint appears to be separate from splicing enhancers, nucleosome positioning, and the translational optimization generating canonical codon bias. The fraction of synonymous sites under strong constraint within a gene correlates well with gene expression, particularly in the mid-late embryo, pupae, and adult developmental stages. Genes enriched in strongly constrained synonymous sites tend to be particularly functionally important and are often involved in key developmental pathways. Given that the observed widespread constraint acting on synonymous sites is likely not limited to Drosophila, the role of synonymous sites in genetic disease and adaptation should be reevaluated.

Characterization of cys-loop receptor genes involved in inhibitory amine neurotransmission in parasitic and free living nematodes

We have isolated two genes, Hco-lgc-53 and Hco-mod-1, from the parasitic nematode Haemonchus contortus, which are orthologs of previously characterized genes that encode dopamine and serotonin-gated chloride channels, respectively, in Caenorhabditis elegans. A search of transcriptome data for the filarial nematode parasites Loa loa, Brugia malayi, and Wucheria bancrofti revealed predicted coding sequences for orthologs of acetylcholine, serotonin and dopamine-gated chloride channels, which correspond to the C. elegans clades acc-1, mod-1 and ggr-3, respectively. Genome data for the more distantly related nematode parasite, Trichinella spiralis, contain genes predicted to encode members of the acc-1 clade only, but all three clades were absent from the trematode Schistosoma mansoni. Analysis of the ratio of non-synonymous to synonymous substitutions (ω) for receptor subunit sequences revealed strong selective constraint over the entire protein, consistent with the known highly conserved 3D structure of cys-loop receptors. This constraint was significantly greater for binding loop residues that are predicted to contact bound ligand and residues of the transmembrane domains. The substitution rate for ligand binding residues was significantly higher for branches leading to the acc-1 and mod-1 clades, where the convergent evolution for binding acetylcholine and serotonin, respectively, is thought to have occurred. Homology models of both Hco-MOD-1 and Hco-LGC-53 channels revealed the presence of binding structures typical of the cys-loop receptor family, including the presence of an aromatic box that is important for the formation of the binding pocket. Both receptors contain a tryptophan in loop C that appears to be a key residue important for the binding of amines to ligand-gated chloride channels. As additional ligand-gated chloride-channel sequences become available for a wider range of species the combination of molecular modeling and analysis of sequence evolution should provide an effective tool to understand the wide diversity of neurotransmitters that bind to this unique group of receptors.

The Molecular Evolutionary Patterns of the Insulin/FOXO Signaling Pathway

The insulin/insulin growth factor-1(IGF1)/FOXO (IIF) signal transduction pathway plays a core role in the endocrine system. Although the components of this pathway have been well characterized, the evolutionary pattern remains poorly understood. Here, we perform a comprehensive analysis to study whether the differences of signaling transduction elements exist as well as to determine whether the genes are subject to equivalent evolutionary forces and how natural selection shapes the evolution pattern of proteins in an interacting system. Our results demonstrate that most IIF pathway components are present throughout all animal phyla investigated here, and they are under strong selective constraint. Remarkably, we detect that the components in the middle of the pathway undergo stronger purifying selection, which is different from previous similar reports. We also find that the dN/dS may be influenced by quite complicated factors including codon bias, protein length among others.

Selective Constraint on Copy Number Variation in Human Piwi-Interacting RNA Loci

Piwi-interacting RNAs (piRNAs) are a recently discovered class of small non-coding RNA found in animals. PiRNAs are primarily expressed in the germline where their best understood function is to repress transposable elements. Unlike previous studies that investigated the evolution of piRNA-generating loci at the level of nucleotide substitutions, here we studied the evolution of piRNA-generating loci at the level of copy number variation (i.e. duplications and deletions) using genome-wide copy number variation data from three human populations. Our analysis shows that at the level of copy number variation there is strong selective constraint and a very high mutation rate in human piRNA-generating loci. Our results differ from a model of positive selection on copy number variation in piRNA-generating loci previously proposed in rodents. We discuss possible reasons for this difference based on the transposable element insertion histories in the rodent and primate lineages.

DEVELOPMENTAL BASIS OF TOOTHLESSNESS IN TURTLES: INSIGHT INTO CONVERGENT EVOLUTION OF VERTEBRATE MORPHOLOGY

The tooth is a major component of the vertebrate feeding apparatus and plays a crucial role in species survival, thus subjecting tooth developmental programs to strong selective constraints. However, irrespective of their functional importance, teeth have been lost in multiple lineages of tetrapod vertebrates independently. To understand both the generality and the diversity of developmental mechanisms that cause tooth agenesis in tetrapods, we investigated expression patterns of a series of tooth developmental genes in the lower jaw of toothless turtles and compared them to that of toothed crocodiles and the chicken as a representative of toothless modern birds. In turtle embryos, we found impairment of Shh signaling in the oral epithelium and early-stage arrest of odontoblast development caused by termination of Msx2 expression in the dental mesenchyme. Our data indicate that such changes underlie tooth agenesis in turtles and suggest that the mechanism that leads to early-stage odontogenic arrest differs between birds and turtles. Our results demonstrate that the cellular and molecular mechanisms that regulate early-stage arrest of tooth development are diverse in tetrapod lineages, and odontogenic developmental programs may respond to changes in upstream molecules similarly thereby evolving convergently with feeding morphology.