MTR3D: identifying regions within protein tertiary structures under purifying selection.
Ontology highlight
ABSTRACT:
SUBMITTER: Silk M
PROVIDER: S-EPMC8265191 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Advances in genomic sequencing have enormous potential to revolutionize personalized medicine, however distinguishing disease-causing from benign variants remains a challenge. The increasing number of human genome and exome sequences available has revealed areas where unfavourable variation is removed through purifying selection. Here, we present the MTR-Viewer, a web-server enabling easy visualization at the gene or variant level of the Missense Tolerance Ratio (MTR), a measure of regional intolerance to missense variation calculated using variation from 240 000 exome and genome sequences. The MTR-Viewer enables exploration of MTR calculations, using different sliding windows, for over 18 000 human protein-coding genes and 85 000 alternative transcripts. Users can also view MTR scores calculated for specific ethnicities, to enable easy exploration of regions that may be under different selective pressure. The spatial distribution of population and known disease variants is also displayed on the protein's domain structure. Intolerant regions were found to be highly enriched for ClinVar pathogenic and COSMIC somatic missense variants (Mann-Whitney U test P < 2.2 × 10-16). As the MTR is not biased by known domains and protein features, it can highlight functionally important regions within genes overlooked or inaccessible by traditional methods. MTR-Viewer is freely available via a user friendly web-server at http://biosig.unimelb.edu.au/mtr-viewer/.
| S-EPMC6602522 | biostudies-literature
Project description:The phenotypic effect of an allele at one genetic site may depend on alleles at other sites, a phenomenon known as epistasis. Epistasis can profoundly influence the process of evolution in populations and shape the patterns of protein divergence across species. Whereas epistasis between adaptive substitutions has been studied extensively, relatively little is known about epistasis under purifying selection. Here we use computational models of thermodynamic stability in a ligand-binding protein to explore the structure of epistasis in simulations of protein sequence evolution. Even though the predicted effects on stability of random mutations are almost completely additive, the mutations that fix under purifying selection are enriched for epistasis. In particular, the mutations that fix are contingent on previous substitutions: Although nearly neutral at their time of fixation, these mutations would be deleterious in the absence of preceding substitutions. Conversely, substitutions under purifying selection are subsequently entrenched by epistasis with later substitutions: They become increasingly deleterious to revert over time. Our results imply that, even under purifying selection, protein sequence evolution is often contingent on history and so it cannot be predicted by the phenotypic effects of mutations assayed in the ancestral background.
| S-EPMC4485141 | biostudies-literature
Project description:Placental mammals (Placentalia) are a very successful group that, today, comprise 94% of all mammalian species. Recent phylogenetic analyses, coupled with new, quite complete fossils, suggest that the crown orders were all established rapidly from a common ancestor just after the Cretaceous/Tertiary (K/T) boundary 65 million years ago. Extensive molecular and morphologic evidence has led to a description of the common ancestor of all Placentalia in which a two-horned uterus and a hemochorial placenta are present. Thus, the process of placentation in which the placenta invades and anchors to the uterine epithelium was already established. One factor that has been suggested as a crucial component of this process is placenta-specific protein 1 (PLAC1). A phylogenetic analysis of the PLAC1 protein in 25 placental mammal species, representing nine of the sixteen crown orders of the Placentalia, suggests that this protein was present in the placental common ancestor in the form we see it today, that it evolved in the Placentalia and has been subject to the effects of purifying selection since its appearance.
| S-EPMC4142310 | biostudies-literature
Project description:An excess of nonsynonymous over synonymous substitution at individual amino acid sites is an important indicator that positive selection has affected the evolution of a protein between the extant sequences under study and their most recent common ancestor. Several methods exist to detect the presence, and sometimes location, of positively selected sites in alignments of protein-coding sequences. This article describes the "sitewise likelihood-ratio" (SLR) method for detecting nonneutral evolution, a statistical test that can identify sites that are unusually conserved as well as those that are unusually variable. We show that the SLR method can be more powerful than currently published methods for detecting the location of positive selection, especially in difficult cases where the strength of selection is low. The increase in power is achieved while relaxing assumptions about how the strength of selection varies over sites and without elevated rates of false-positive results that have been reported with some other methods. We also show that the SLR method performs well even under circumstances where the results from some previous methods can be misleading.
| S-EPMC1449526 | biostudies-literature
Project description:Diverse life forms are driven by the evolution of gene regulatory programs including changes in regulator proteins and cis-regulatory elements. Alterations of cis-regulatory elements are likely to dominate the evolution of the gene regulatory networks, as they are subjected to smaller selective constraints compared with proteins and hence may evolve quickly to adapt the environment. Prior studies on cis-regulatory element evolution focus primarily on sequence substitutions of known transcription factor-binding motifs. However, evolutionary models for the dynamics of motif occurrence are relatively rare, and comprehensive characterization of the evolution of all possible motif sequences has not been pursued. In the present study, we propose an algorithm to estimate the strength of purifying selection of a motif sequence based on an evolutionary model capturing the birth and death of motif occurrences on promoters. We term this measure as the 'evolutionary retention coefficient', as it is related yet distinct from the canonical definition of selection coefficient in population genetics. Using this algorithm, we estimate and report the evolutionary retention coefficients of all possible 10-nucleotide sequences from the aligned promoter sequences of 27 748. orthologous gene families in 34 mammalian species. Intriguingly, the evolutionary retention coefficients of motifs are intimately associated with their functional relevance. Top-ranking motifs (sorted by evolutionary retention coefficients) are significantly enriched with transcription factor-binding sequences according to the curated knowledge from the TRANSFAC database and the ChIP-seq data generated from the ENCODE Consortium. Moreover, genes harbouring high-scoring motifs on their promoters retain significantly coherent expression profiles, and those genes are over-represented in the functional classes involved in gene regulation. The validation results reveal the dependencies between natural selection and functions of cis-regulatory elements and shed light on the evolution of gene regulatory networks.
| S-EPMC3575792 | biostudies-literature
Project description:Like many other large double-stranded DNA (dsDNA) viruses, herpesviruses are known to capture host genes to evade host defenses. Little is known about the detailed natural history of such genes, nor do we fully understand their evolutionary dynamics. A major obstacle is that they are often highly divergent, maintaining very low sequence similarity to host homologs. Here we use the herpesvirus genus Rhadinovirus as a model system to develop an analytical approach that combines complementary evolutionary and bioinformatic techniques, offering results that are both detailed and robust for a range of genes. Using a systematic phylogenetic strategy, we identify the original host lineage of viral genes with high confidence. We show that although host immunomodulatory genes evolve rapidly compared to other host genes, they undergo a clear increase in purifying selection once captured by a virus. To characterize this shift in detail, we developed a novel technique to identify changes in selection pressure that can be attributable to particular domains. These findings will inform us on how viruses develop strategies to evade the immune system, and our synthesis of techniques can be reapplied to other viruses or biological systems with similar analytical challenges.IMPORTANCE Viruses and hosts have been shown to capture genes from one another as part of the evolutionary arms race. Such genes offer a natural experiment on the effects of evolutionary pressure, since the same gene exists in vastly different selective environments. However, sequences of viral homologs often bear little similarity to the original sequence, complicating the reconstruction of their shared evolutionary history with host counterparts. In this study, we use a genus of herpesviruses as a model system to comprehensively investigate the evolution of host-derived viral genes, using a synthesis of genomics, phylogenetics, selection analysis, and nucleotide and amino acid modeling.
| S-EPMC6146799 | biostudies-other
Project description:With only 2 % of the human genome consisting of protein coding genes, functionality across the rest of the genome has been the subject of much debate. This has gained further impetus in recent years due to a rapidly growing catalogue of genomic elements, based primarily on biochemical signatures (e.g. the ENCODE project). While the assessment of functionality is a complex task, the presence of selection acting on a genomic region is a strong indicator of importance. In this study, we apply population genetic methods to investigate signals overlaying several classes of regulatory elements.We disentangle signals of purifying selection acting directly on regulatory elements from the confounding factors of demography and purifying selection linked to e.g. nearby protein coding regions. We confirm the importance of regulatory regions proximal to coding sequence, while also finding differential levels of selection at distal regions. We note differences in purifying selection among transcription factor families. Signals of constraint at some genomic classes were also strongly dependent on their physical location relative to coding sequence. In addition, levels of selection efficacy across genomic classes differed between African and non-African populations.In order to assign a valid signal of selection to a particular class of genomic sequence, we show that it is crucial to isolate the signal by accounting for the effects of demography and linked-purifying selection. Our study highlights the intricate interplay of factors affecting signals of selection on functional elements.
| S-EPMC5787233 | biostudies-literature
Project description:Bindin is a sperm protein that mediates attachment and membrane fusion of gametes. The mode of bindin evolution varies across sea urchin genera studied to date. In three genera it evolves under positive selection, in four under mostly purifying selection, and in one, results have been mixed. We studied bindin evolution in the pantropical sea urchin Diadema, which split from other studied genera 250 million years ago. We found that Diadema bindin is structurally similar to that of other genera, but much longer (418 amino acids). In seven species of Diadema, bindin evolves under purifying selection, more slowly than in any other sea urchin genus. Only bindin of the recently rediscovered D. clarki shows evidence of positive selection. As D. clarki is sympatric with D. setosum and D. savignyi, positive selection could arise from avoidance of maladaptive hybridization. However, D. setosum and D. savignyi overlap in the Indo-West Pacific, yet their bindins show no evidence of positive selection, possibly because the two species spawn at different times. Bindin in the East Pacific D. mexicanum, the West Atlantic D. antillarum, the East Atlantic D. africanum, and the Indo-Pacific D. paucispinum also evolves slowly under purifying selection.
| S-EPMC7299941 | biostudies-literature
Project description:Identification of cancer driver genes using somatic mutation patterns indicative of positive selection has become a major goal in cancer genomics. However, cancer cells additionally depend on a large number of genes involved in basic cellular processes. While such genes should in theory be subject to strong purifying (negative) selection against damaging somatic mutations, these patterns have been elusive and purifying selection remains inadequately explored in cancer. Here, we hypothesized that purifying selection should be evident in hemizygous genomic regions, where damaging mutations cannot be compensated for by healthy alleles. Using a 7,781-sample pan-cancer dataset, we first confirmed this in POLR2A, an essential gene where hemizygous deletions are known to confer elevated sensitivity to pharmacological suppression. We next used this principle to identify several genes and pathways that show patterns indicative of purifying selection to avoid deleterious mutations. These include the POLR2A interacting protein INTS10 as well as genes involved in mRNA splicing, nonsense-mediated mRNA decay and other RNA processing pathways. Many of these genes belong to large protein complexes, and strong overlaps were observed with recent functional screens for gene essentiality in human cells. Our analysis supports that purifying selection acts to preserve the remaining function of many hemizygously deleted essential genes in tumors, indicating vulnerabilities that might be exploited by future therapeutic strategies.
| S-EPMC5189942 | biostudies-literature
Project description:A critical question in biology is the identification of functionally important amino acid sites in proteins. Because functionally important sites are under stronger purifying selection, site-specific substitution rates tend to be lower than usual at these sites. A large number of phylogenetic models have been developed to estimate site-specific substitution rates in proteins and the extraordinarily low substitution rates have been used as evidence of function. Most of the existing tools, e.g. Rate4Site, assume that site-specific substitution rates are independent across sites. However, site-specific substitution rates may be strongly correlated in the protein tertiary structure, since functionally important sites tend to be clustered together to form functional patches. We have developed a new model, GP4Rate, which incorporates the Gaussian process model with the standard phylogenetic model to identify slowly evolved regions in protein tertiary structures. GP4Rate uses the Gaussian process to define a nonparametric prior distribution of site-specific substitution rates, which naturally captures the spatial correlation of substitution rates. Simulations suggest that GP4Rate can potentially estimate site-specific substitution rates with a much higher accuracy than Rate4Site and tends to report slowly evolved regions rather than individual sites. In addition, GP4Rate can estimate the strength of the spatial correlation of substitution rates from the data. By applying GP4Rate to a set of mammalian B7-1 genes, we found a highly conserved region which coincides with experimental evidence. GP4Rate may be a useful tool for the in silico prediction of functionally important regions in the proteins with known structures.
| S-EPMC3894161 | biostudies-literature