Gene-pseudogene evolution: a probabilistic approach.
Ontology highlight
ABSTRACT:
SUBMITTER: Mahmudi O
PROVIDER: S-EPMC4602177 | biostudies-literature | 2015
REPOSITORIES: biostudies-literature
Publications
Gene-pseudogene evolution: a probabilistic approach.
Mahmudi Owais O Sennblad Bengt B Arvestad Lars L Nowick Katja K Lagergren Jens J
BMC genomics 20151002
Over the last decade, methods have been developed for the reconstruction of gene trees that take into account the species tree. Many of these methods have been based on the probabilistic duplication-loss model, which describes how a gene-tree evolves over a species-tree with respect to duplication and losses, as well as extension of this model, e.g., the DLRS (Duplication, Loss, Rate and Sequence evolution) model that also includes sequence evolution under relaxed molecular clock. A disjoint, al ...[more]
Similar Datasets
Project description:Pseudogenes are fossil relatives of genes. Pseudogenes have long been thought of as "junk DNAs", since they do not code proteins in normal tissues. Although most of the human pseudogenes do not have noticeable functions, ?20% of them exhibit transcriptional activity. There has been evidence showing that some pseudogenes adopted functions as lncRNAs and work as regulators of gene expression. Furthermore, pseudogenes can even be "reactivated" in some conditions, such as cancer initiation. Some pseudogenes are transcribed in specific cancer types, and some are even translated into proteins as observed in several cancer cell lines. All the above have shown that pseudogenes could have functional roles or potentials in the genome. Evaluating the relationships between pseudogenes and their gene counterparts could help us reveal the evolutionary path of pseudogenes and associate pseudogenes with functional potentials. It also provides an insight into the regulatory networks involving pseudogenes with transcriptional and even translational activities.In this study, we develop a novel approach integrating graph analysis, sequence alignment and functional analysis to evaluate pseudogene-gene relationships, and apply it to human gene homologs and pseudogenes. We generated a comprehensive set of 445 pseudogene-gene (PGG) families from the original 3,281 gene families (13.56%). Of these 438 (98.4% PGG, 13.3% total) were non-trivial (containing more than one pseudogene). Each PGG family contains multiple genes and pseudogenes with high sequence similarity. For each family, we generate a sequence alignment network and phylogenetic trees recapitulating the evolutionary paths. We find evidence supporting the evolution history of olfactory family (both genes and pseudogenes) in human, which also supports the validity of our analysis method. Next, we evaluate these networks in respect to the gene ontology from which we identify functions enriched in these pseudogene-gene families and infer functional impact of pseudogenes involved in the networks. This demonstrates the application of our PGG network database in the study of pseudogene function in disease context.
| S-EPMC5744670 | biostudies-literature
Project description:The gene encoding the cytosolic β-glucosidase GBA3 shows pseudogenization due to a truncated allele (rs358231) that is polymorphic in humans. Since this enzyme is involved in the transformation of many plant β-glycosides, this particular case of gene loss may have been influenced by dietary adaptations during evolution. In humans, apart from the inactivating allele, we found that GBA3 accumulated additional damaging mutations, implying an extensive GBA3 loss. The allelic distribution of loss-of-function alleles revealed significant differences between human populations which can be partially related with their staple diet. The analysis of mammalian orthologs disclosed that GBA3 underwent at least nine pseudogenization events. Most events of pseudogenization occurred in carnivorous lineages, suggesting a possible link to a β-glycoside poor diet. However, GBA3 was also lost in omnivorous and herbivorous species, hinting that the physiological role of GBA3 is not fully understood and other unknown causes may underlie GBA3 pseudogenization. Such possibility relies upon a putative role in sialic acid biology, where GBA3 participates in a cellular network involving NEU2 and CMAH. Overall, our data shows that the recurrent loss of GBA3 in mammals is likely to represent an evolutionary endpoint of the relaxation of selective constraints triggered by diet-related factors.
| S-EPMC7360587 | biostudies-literature
Project description:BackgroundMost models of genome evolution concern either genetic sequences, gene content or gene order. They sometimes integrate two of the three levels, but rarely the three of them. Probabilistic models of gene order evolution usually have to assume constant gene content or adopt a presence/absence coding of gene neighborhoods which is blind to complex events modifying gene content.ResultsWe propose a probabilistic evolutionary model for gene neighborhoods, allowing genes to be inserted, duplicated or lost. It uses reconciled phylogenies, which integrate sequence and gene content evolution. We are then able to optimize parameters such as phylogeny branch lengths, or probabilistic laws depicting the diversity of susceptibility of syntenic regions to rearrangements. We reconstruct a structure for ancestral genomes by optimizing a likelihood, keeping track of all evolutionary events at the level of gene content and gene synteny. Ancestral syntenies are associated with a probability of presence.
| S-EPMC4603630 | biostudies-literature
Project description:BackgroundPseudogenes are an integral component of the human genome. Little attention, however, has so far been paid to the phenomenon that some pseudogenes are transcriptionally active. Recently, we demonstrated that the human ortholog of the rodent testis-specific ATP-binding cassette (ABC) transporter Abca17 is a ubiquitously transcribed pseudogene (ABCA17P). The aim of the present study was to establish a complete inventory of all ABC transporter pseudogenes in the human genome and to identify transcriptionally active ABC transporter pseudogenes. Moreover, we tested the hypothesis that a regulatory interdependency exists between ABC transporter pseudogenes and their parental protein coding equivalents.ResultsSystematic bioinformatic analysis revealed the existence of 22 ABC transporter pseudogenes within the human genome. We identified two clusters on chromosomes 15 and 16, respectively, which harbor almost half of all pseudogenes (n = 10). Available information from EST and mRNA databases and RT-PCR expression profiling indicate that a large portion of the ABC transporter pseudogenes (45%, n = 10) are transcriptionally active and some of them are expressed as alternative splice variants. We demonstrate that both pseudogenes of the pseudoxanthoma elasticum gene ABCC6, ABCC6P1 and ABCC6P2, are transcribed. ABCC6P1 and ABCC6 possess near-identical promoter sequences and their tissue-specific expression profiles are strikingly similar raising the possibility that they form a gene-pseudogene dual transcription unit. Intriguingly, targeted knockdown of the transcribed pseudogene ABCC6P1 resulted in a significant reduction of ABCC6 mRNA expression levels.ConclusionThe human genome contains a surprisingly small number of ABC transporter pseudogenes relative to other known gene families. They are unevenly distributed across the chromosomes. Importantly, a significant portion of the ABC transporter pseudogenes is transcriptionally active. The downregulation of ABCC6 mRNA levels by targeted suppression of the expression of its pseudogene ABCC6P1 provides evidence, for the first time, for a regulatory interdependence of a transcribed pseudogene and its protein coding counterpart in the human genome.
| S-EPMC2329642 | biostudies-literature
Project description:We performed evolution of ∆menF∆entC strain of Escherichia coli K12 MG1655 to study how the system adapt to the loss of siderophore biosynthesis. RNA-Seq was performed to examine the underlying transcriptional rewiring.
2018-12-20 | GSE122779 | GEO
Project description:Regulation by microRNAs (miRNAs) and modulation of miRNA activity are critical components of diverse cellular processes. Recent research has shown that miRNA-based regulation of the tumor suppressor gene PTEN can be modulated by the expression of other miRNA targets acting as competing endogenous RNAs (ceRNAs). However, the key sequence-based features enabling a transcript to act as an effective ceRNA are not well understood and a quantitative model associating statistical significance to such features is currently lacking. To identify and assess features characterizing target recognition by PTEN-regulating miRNAs, we analyze multiple datasets from PAR-CLIP experiments in conjunction with RNA-Seq data. We consider a set of miRNAs known to regulate PTEN and identify high-confidence binding sites for these miRNAs on the 3' UTR of protein coding genes. Based on the number and spatial distribution of these binding sites, we calculate a set of probabilistic features that are used to make predictions for novel ceRNAs of PTEN. Using a series of experiments in human prostate cancer cell lines, we validate the highest ranking prediction (TNRC6B) as a ceRNA of PTEN. The approach developed can be applied to map ceRNA networks of critical cellular regulators and to develop novel insights into crosstalk between different pathways involved in cancer.
| S-EPMC5552881 | biostudies-literature
Project description:BackgroundPairwise stochastic context-free grammars (Pair SCFGs) are powerful tools for evolutionary analysis of RNA, including simultaneous RNA sequence alignment and secondary structure prediction, but the associated algorithms are intensive in both CPU and memory usage. The same problem is faced by other RNA alignment-and-folding algorithms based on Sankoff's 1985 algorithm. It is therefore desirable to constrain such algorithms, by pre-processing the sequences and using this first pass to limit the range of structures and/or alignments that can be considered.ResultsWe demonstrate how flexible classes of constraint can be imposed, greatly reducing the computational costs while maintaining a high quality of structural homology prediction. Any score-attributed context-free grammar (e.g. energy-based scoring schemes, or conditionally normalized Pair SCFGs) is amenable to this treatment. It is now possible to combine independent structural and alignment constraints of unprecedented general flexibility in Pair SCFG alignment algorithms. We outline several applications to the bioinformatics of RNA sequence and structure, including Waterman-Eggert N-best alignments and progressive multiple alignment. We evaluate the performance of the algorithm on test examples from the RFAM database.ConclusionA program, Stemloc, that implements these algorithms for efficient RNA sequence alignment and structure prediction is available under the GNU General Public License.
| S-EPMC1090553 | biostudies-literature
Project description:BackgroundThe CRISPR/Cas system is known to act as an adaptive and heritable immune system in Eubacteria and Archaea. Immunity is encoded in an array of spacer sequences. Each spacer can provide specific immunity to invasive elements that carry the same or a similar sequence. Even in closely related strains, spacer content is very dynamic and evolves quickly. Standard models of nucleotide evolution cannot be applied to quantify its rate of change since processes other than single nucleotide changes determine its evolution.MethodsWe present probabilistic models that are specific for spacer content evolution. They account for the different processes of insertion and deletion. Insertions can be constrained to occur on one end only or are allowed to occur throughout the array. One deletion event can affect one spacer or a whole fragment of adjacent spacers. Parameters of the underlying models are estimated for a pair of arrays by maximum likelihood using explicit ancestor enumeration.ResultsSimulations show that parameters are well estimated on average under the models presented here. There is a bias in the rate estimation when including fragment deletions. The models also estimate times between pairs of strains. But with increasing time, spacer overlap goes to zero, and thus there is an upper bound on the distance that can be estimated. Spacer content similarities are displayed in a distance based phylogeny using the estimated times.We use the presented models to analyze different Yersinia pestis data sets and find that the results among them are largely congruent. The models also capture the variation in diversity of spacers among the data sets. A comparison of spacer-based phylogenies and Cas gene phylogenies shows that they resolve very different time scales for this data set.ConclusionsThe simulations and data analyses show that the presented models are useful for quantifying spacer content evolution and for displaying spacer content similarities of closely related strains in a phylogeny. This allows for comparisons of different CRISPR arrays or for comparisons between CRISPR arrays and nucleotide substitution rates.
| S-EPMC3704272 | biostudies-literature
Project description:Comparative sequence analyses, including such fundamental bioinformatics techniques as similarity searching, sequence alignment and phylogenetic inference, have become a mainstay for researchers studying type 1 Human Immunodeficiency Virus (HIV-1) genome structure and evolution. Implicit in comparative analyses is an underlying model of evolution, and the chosen model can significantly affect the results. In general, evolutionary models describe the probabilities of replacing one amino acid character with another over a period of time. Most widely used evolutionary models for protein sequences have been derived from curated alignments of hundreds of proteins, usually based on mammalian genomes. It is unclear to what extent these empirical models are generalizable to a very different organism, such as HIV-1-the most extensively sequenced organism in existence. We developed a maximum likelihood model fitting procedure to a collection of HIV-1 alignments sampled from different viral genes, and inferred two empirical substitution models, suitable for describing between-and within-host evolution. Our procedure pools the information from multiple sequence alignments, and provided software implementation can be run efficiently in parallel on a computer cluster. We describe how the inferred substitution models can be used to generate scoring matrices suitable for alignment and similarity searches. Our models had a consistently superior fit relative to the best existing models and to parameter-rich data-driven models when benchmarked on independent HIV-1 alignments, demonstrating evolutionary biases in amino-acid substitution that are unique to HIV, and that are not captured by the existing models. The scoring matrices derived from the models showed a marked difference from common amino-acid scoring matrices. The use of an appropriate evolutionary model recovered a known viral transmission history, whereas a poorly chosen model introduced phylogenetic error. We argue that our model derivation procedure is immediately applicable to other organisms with extensive sequence data available, such as Hepatitis C and Influenza A viruses.
| S-EPMC1876811 | biostudies-literature
Project description:RNA-Sequencing has been the leading technology to quantify expression of thousands of genes simultaneously. The data analysis of an RNA-Seq experiment starts from aligning short reads to the reference genome/transcriptome or reconstructed transcriptome. However, current aligners lack the sensitivity to distinguish reads that come from homologous regions of an genome. One group of these homologies is the paralog pseudogenes. Pseudogenes arise from duplication of a set of protein coding genes, and have been considered as degraded paralogs in the genome due to their lost of functionality. Recent studies have provided evidence to support their novel regulatory roles in biological processes. With the growing interests in quantifying the expression level of pseudogenes at different tissues or cell lines, it is critical to have a sensitive method that can correctly align ambiguous reads and accurately estimate the expression level among homologous genes. Previously in PseudoLasso, we proposed a linear regression approach to learn read alignment behaviors, and to leverage this knowledge for abundance estimation and alignment correction. In this paper, we extend the work of PseudoLasso by grouping the homologous genomic regions into different communities using a community detection algorithm, followed by building a linear regression model separately for each community. The results show that this approach is able to retain the same accuracy as PseudoLasso. By breaking the genome into smaller homologous communities, the running time is improved from quadratic growth to linear with respect to the number of genes.
| S-EPMC5514313 | biostudies-literature