In silico identification of functional divergence between the multiple groEL gene paralogs in Chlamydiae

<p>Abstract</p> <p>Background</p> <p>Heat-shock proteins are specialized molecules performing different and essential roles in the cell including protein degradation, folding and trafficking. GroEL is a 60 Kda heat-shock protein ubiquitous in bacteria and has been regarded as an important molecule implicated in chronic inflammatory processes caused by <it>Chlamydiae </it>infections. GroEL in <it>Chlamydiae </it>became duplicated at the origin of the <it>Chlamydiae </it>lineage presenting three distinct molecular chaperones, namely the original protein GroEL1 (Ct110), and its paralogous proteins GroEL2 (Ct604) and GroEL3 (Ct755). These chaperones present differential and independent expressions during the different stages of <it>Chlamydiae </it>infections and have been suggested to present differential physiological and regulatory roles.</p> <p>Results</p> <p>In this comprehensive <it>in silico </it>study we show that GroEL protein paralogs have diverged functionally after the different gene duplication events and that this divergence has occurred mainly between GroEL3 and GroEL1. GroEL2 presents an intermediate functional divergence pattern from GroEL1. Our results point to the different protein-protein interaction patterns between GroEL paralogs and known GroEL protein clients supporting their functional divergence after <it>groEL </it>gene duplication. Analysis of selective constraints identifies periods of adaptive evolution after gene duplication that led to the fixation of amino acid replacements in GroEL protein domains involved in the interaction with GroEL protein clients.</p> <p>Conclusion</p> <p>We demonstrate that GroEL protein copies in <it>Chlamydiae </it>species have diverged functionally after the gene duplication events. We also show that functional divergence has occurred in important functional regions of these GroEL proteins and that very probably have affected the ancestral GroEL regulatory role and protein-protein interaction patterns with GroEL client proteins. Most of the amino acid replacements that have affected interaction with protein clients and that were responsible for the functional divergence between GroEL paralogs were fixed by adaptive evolution after the <it>groEL </it>gene duplication events.</p

Describing the structural robustness landscape of bacterial small RNAs

<p>Abstract</p> <p>Background</p> <p>The potential role of RNA molecules as gene expression regulators has led to a new perspective on the intracellular control and genome organization. Because secondary structures are crucial for their regulatory role, we sought to investigate their robustness to mutations and environmental changes.</p> <p>Results</p> <p>Here, we dissected the structural robustness landscape of the small non-coding RNAs (sncRNAs) encoded in the genome of the bacterium <it>Escherichia coli</it>. We found that bacterial sncRNAs are not significantly robust to both mutational and environmental perturbations when compared against artificial, unbiased sequences. However, we found that, on average, bacterial sncRNAs tend to be significantly plastic, and that mutational and environmental robustness strongly correlate. We further found that, on average, epistasis in bacterial sncRNAs is significantly antagonistic, and positively correlates with plasticity. Moreover, the evolution of robustness is likely dependent upon the environmental stability of the cell, with more fluctuating environments leading to the emergence and fixation of more robust molecules. Mutational robustness also appears to be correlated with structural functionality and complexity.</p> <p>Conclusion</p> <p>Our study provides a deep characterization of the structural robustness landscape of bacterial sncRNAs, suggesting that evolvability could be evolved as a consequence of selection for more plastic molecules. It also supports that environmental fluctuations could promote mutational robustness. As a result, plasticity emerges to link robustness, functionality and evolvability.</p

Unravelling Selection Shifts Among Foot-and-Mouth Disease Virus (FMDV) Serotypes

FMDV virus has been increasingly recognised as the most economically severe animal virus with a remarkable degree of antigenic diversity. Using an integrative evolutionary and computational approach we have compelling evidence for heterogeneity in the selection forces shaping the evolution of the seven different FMDV serotypes. Our results show that positive Darwinian selection has governed the evolution of the major antigenic regions of serotypes A, Asia1, O, SAT1 and SAT2, but not C or SAT3. Co-evolution between sites from antigenic regions under positive selection pinpoints their functional communication to generate immune-escape mutants while maintaining their ability to recognise the host-cell receptors. Neural network and functional divergence analyses strongly point to selection shifts between the different serotypes. Our results suggest that, unlike African FMDV serotypes, serotypes with wide geographical distribution have accumulated compensatory mutations as a strategy to ameliorate the effect of slightly deleterious mutations fixed by genetic drift. This strategy may have provided the virus by a flexibility to generate immune-escape mutants and yet recognise host-cell receptors. African serotypes presented no evidence for compensatory mutations. Our results support heterogeneous selective constraints affecting the different serotypes. This points to the possible accelerated rates of evolution diverging serotypes sharing geographical locations as to ameliorate the competition for the host

Why Should We Care About Molecular Coevolution?

Non-independent evolution of amino acid sites has become a noticeable limitation of most methods aimed at identifying selective constraints at functionally important amino acid sites or protein regions. The need for a generalised framework to account for non-independence of amino acid sites has fuelled the design and development of new mathematical models and computational tools centred on resolving this problem. Molecular coevolution is one of the most active areas of research, with an increasing rate of new models and methods being developed everyday. Both parametric and non-parametric methods have been developed to account for correlated variability of amino acid sites. These methods have been utilised for detecting phylogenetic, functional and structural coevolution as well as to identify surfaces of amino acid sites involved in protein-protein interactions. Here we discuss and briefly describe these methods, and identify their advantages and limitations

Relationships of gag-pol diversity between Ty3/Gypsy and Retroviridae LTR retroelements and the three kings hypothesis

<p>Abstract</p> <p>Background</p> <p>The origin of vertebrate retroviruses (<it>Retroviridae</it>) is yet to be thoroughly investigated, but due to their similarity and identical gag-pol (and env) genome structure, it is accepted that they evolve from <it>Ty3/Gypsy </it>LTR retroelements the retrotransposons and retroviruses of plants, fungi and animals. These 2 groups of LTR retroelements code for 3 proteins rarely studied due to the high variability – gag polyprotein, protease and GPY/F module. In relation to 3 previously proposed <it>Retroviridae </it>classes I, II and II, investigation of the above proteins conclusively uncovers important insights regarding the ancient history of <it>Ty3/Gypsy </it>and <it>Retroviridae </it>LTR retroelements.</p> <p>Results</p> <p>We performed a comprehensive study of 120 non-redundant <it>Ty3/Gypsy </it>and <it>Retroviridae </it>LTR retroelements. Phylogenetic reconstruction inferred based on the concatenated analysis of the gag and pol polyproteins shows a robust phylogenetic signal regarding the clustering of OTUs. Evaluation of gag and pol polyproteins separately yields discordant information. While pol signal supports the traditional perspective (2 monophyletic groups), gag polyprotein describes an alternative scenario where each <it>Retroviridae </it>class can be distantly related with one or more <it>Ty3/Gypsy </it>lineages. We investigated more in depth this evidence through comparative analyses performed based on the gag polyprotein, the protease and the GPY/F module. Our results indicate that contrary to the traditional monophyletic view of the origin of vertebrate retroviruses, the <it>Retroviridae </it>class I is a molecular fossil, preserving features that were probably predominant among <it>Ty3/Gypsy </it>ancestors predating the split of plants, fungi and animals. In contrast, classes II and III maintain other phenotypes that emerged more recently during <it>Ty3/Gypsy </it>evolution.</p> <p>Conclusion</p> <p>The 3 <it>Retroviridae </it>classes I, II and III exhibit phenotypic differences that delineate a network never before reported between <it>Ty3/Gypsy </it>and <it>Retroviridae </it>LTR retroelements. This new scenario reveals how the diversity of vertebrate retroviruses is polyphyletically recurrent into the <it>Ty3/Gypsy </it>evolution, i.e. older than previously thought. The simplest hypothesis to explain this finding is that classes I, II and III trace back to at least 3 <it>Ty3/Gypsy </it>ancestors that emerged at different evolutionary times prior to protostomes-deuterostomes divergence. We have called this "the three kings hypothesis" concerning the origin of vertebrate retroviruses.</p

Reducing the false positive rate in the non-parametric analysis of molecular coevolution

<p>Abstract</p> <p>Background</p> <p>The strength of selective constraints operating on amino acid sites of proteins has a multifactorial nature. In fact, amino acid sites within proteins coevolve due to their functional and/or structural relationships. Different methods have been developed that attempt to account for the evolutionary dependencies between amino acid sites. Researchers have invested a significant effort to increase the sensitivity of such methods. However, the difficulty in disentangling functional co-dependencies from historical covariation has fuelled the scepticism over their power to detect biologically meaningful results. In addition, the biological parameters connecting linear sequence evolution to structure evolution remain elusive. For these reasons, most of the evolutionary studies aimed at identifying functional dependencies among protein domains have focused on the structural properties of proteins rather than on the information extracted from linear multiple sequence alignments (MSA). Non-parametric methods to detect coevolution have been reported to be especially susceptible to produce false positive results based on the properties of MSAs. However, no formal statistical analysis has been performed to definitively test the differential effects of these properties on the sensitivity of such methods.</p> <p>Results</p> <p>Here we test the effect that variations on the MSA properties have over the sensitivity of non-parametric methods to detect coevolution. We test the effect that the size of the MSA (number of sequences), mean pairwise amino acid distance per site and the strength of the coevolution signal have on the ability of non-parametric methods to detect coevolution. Our results indicate that all three factors have significant effects on the accuracy of non-parametric methods. Further, introducing statistical filters improves the sensitivity and increases the statistical power of the methods to detect functional coevolution. Statistical analysis of the physico-chemical properties of amino acid sites in the context of the protein structure reveals striking dependencies among amino acid sites. Results indicate a covariation trend in the hydrophobicities and molecular weight characteristics of amino acid sites when analysing a non-redundant set of 8000 protein structures. Using this biological information as filter in coevolutionary analyses minimises the false positive rate of these methods. Application of these filters to three different proteins with known functional domains supports the importance of using biological filters to detect coevolution.</p> <p>Conclusion</p> <p>Coevolutionary analyses using non-parametric methods have proved difficult and highly prone to provide spurious results depending on the properties of MSAs and on the strength of coevolution between amino acid sites. The application of statistical filters to the number of pairs detected as coevolving reduces significantly the number of artifactual results. Analysis of the physico-chemical properties of amino acid sites in the protein structure context reveals their structure-dependent covariation. The application of this known biological information to the analysis of covariation greatly enhances the functional coevolutionary signal and removes historical covariation. Simultaneous use of statistical and biological data is instrumental in the detection of functional amino acid sites dependencies and compensatory changes at the protein level.</p

Mutational dynamics of murine angiogenin duplicates

<p>Abstract</p> <p>Background</p> <p>Angiogenin (Ang) is a protein involved in angiogenesis by inducing the formation of blood vessels. The biomedical importance of this protein has come from findings linking mutations in Ang to cancer progression and neurodegenerative diseases. These findings highlight the evolutionary constrain on Ang amino acid sequence. However, previous studies comparing human Angiogenin with homologs from other phylogenetically related organisms have led to the conclusion that Ang presents a striking variability. Whether this variability has an adaptive value <it>per se </it>remains elusive. Understanding why many functional Ang paralogs have been preserved in mouse and rat and identifying functional divergence mutations at these copies may explain the relationship between mutations and function. In spite of the importance of testing this hypothesis from the evolutionarily and biomedical perspectives, this remains yet unaccomplished. Here we test the main mutational dynamics driving the evolution and function of Ang paralogs in mammals.</p> <p>Results</p> <p>We analysed the phylogenetic asymmetries between the different Ang gene copies in mouse and rat in the context of vertebrate Ang phylogeny. This analysis shows strong evidence in support of accelerated evolution in some Ang murine copies (mAng). This acceleration is not due to non-functionalisation because constraints on amino acid replacements remain strong. We identify many of the amino acid sites involved in signal localization and nucleotide binding by Ang to have evolved under diversifying selection. Compensatory effects of many of the mutations at these paralogs and their key structural location in or nearby important functional regions support a possible functional shift (functional divergence) in many Ang copies. Similarities between 3D-structural models for mAng copies suggest that their divergence is mainly functional.</p> <p>Conclusions</p> <p>We identify the main evolutionary dynamics shaping the variability of Angiogenin in vertebrates and highlight the plasticity of this protein after gene duplication. Our results suggest functional divergence among mAng paralogs. This puts forward mAng as a good system candidate for testing functional plasticity of such an important protein while stresses caution when using mouse as a model to infer the consequences of mutations in the single Ang copy of humans.</p

Evidence of Positively Selected Sites in Mammalian a-Defensins

Defensins are a family of mammalian antimicrobial peptides that exhibit variable activity against a panel of microbes, including bacteria, fungi, and enveloped viruses. We have employed a maximum-likelihood approach to detect evidence of positive selection (adaptive evolution) in the evolution of these important molecules of the innate immune response. We have identified 14 amino acid sites that are predicted to be subject to positive selection. Furthermore, we show that all these sites are located in the mature antimicrobial peptide and not in the prepropeptide region of the molecule, implying that they are of functional importance. These results suggest that mammalian a-defensins have been under selective pressure to evolve in response to potentially infectious challenges by fast-evolving microbes

Evidence of Positively Selected Sites in Mammalian a-Defensins

Defensins are a family of mammalian antimicrobial peptides that exhibit variable activity against a panel of microbes, including bacteria, fungi, and enveloped viruses. We have employed a maximum-likelihood approach to detect evidence of positive selection (adaptive evolution) in the evolution of these important molecules of the innate immune response. We have identified 14 amino acid sites that are predicted to be subject to positive selection. Furthermore, we show that all these sites are located in the mature antimicrobial peptide and not in the prepropeptide region of the molecule, implying that they are of functional importance. These results suggest that mammalian a-defensins have been under selective pressure to evolve in response to potentially infectious challenges by fast-evolving microbes

Evidence from comparative genomics for a complete sexual cycle in the 'asexual' pathogenic yeast Candida glabrata

BACKGROUND: Candida glabrata is a pathogenic yeast of increasing medical concern. It has been regarded as asexual since it was first described in 1917, yet phylogenetic analyses have revealed that it is more closely related to sexual yeasts than other Candida species. We show here that the C. glabrata genome contains many genes apparently involved in sexual reproduction. RESULTS: By genome survey sequencing, we find that genes involved in mating and meiosis are as numerous in C. glabrata as in the sexual species Kluyveromyces delphensis, which is its closest known relative. C. glabrata has a putative mating-type (MAT) locus and a pheromone gene (MFALPHA2), as well as orthologs of at least 31 other Saccharomyces cerevisiae genes that have no known roles apart from mating or meiosis, including FUS3, IME1 and SMK1. CONCLUSIONS: We infer that C. glabrata is likely to have an undiscovered sexual stage in its life cycle, similar to that recently proposed for C. albicans. The two Candida species represent two distantly related yeast lineages that have independently become both pathogenic and 'asexual'. Parallel evolution in the two lineages as they adopted mammalian hosts resulted in separate but analogous switches from overtly sexual to cryptically sexual life cycles, possibly in response to defense by the host immune system