Independent evolution of transcript abundance and gene regulatory dynamics.
Ontology highlight
ABSTRACT:
SUBMITTER: Krieger G
PROVIDER: S-EPMC7397873 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:A high-quality, barley gene reference transcript dataset (BaRTv1.0), was used to quantify gene and transcript abundances from 22 RNA-seq experiments, covering 843 separate samples. Using the abundance data we developed a Barley Expression Database (EORNA*) to underpin a visualisation tool that displays comparative gene and transcript abundance data on demand as transcripts per million (TPM) across all samples and all the genes. EORNA provides gene and transcript models for all of the transcripts contained in BaRTV1.0, and these can be conveniently identified through either BaRT or HORVU gene names, or by direct BLAST of query sequences. Browsing the quantification data reveals cultivar, tissue and condition specific gene expression and shows changes in the proportions of individual transcripts that have arisen via alternative splicing. TPM values can be easily extracted to allow users to determine the statistical significance of observed transcript abundance variation among samples or perform meta analyses on multiple RNA-seq experiments. * Eòrna is the Scottish Gaelic word for Barley.
| S-EPMC7994555 | biostudies-literature
Project description:BACKGROUND:Acute Kawasaki disease (KD) is difficult to distinguish from other illnesses that involve acute rash or fever, in part because the etiologic agent(s) and pathophysiology remain poorly characterized. As a result, diagnosis and critical therapies may be delayed. METHODS:We used DNA microarrays to identify possible diagnostic features of KD. We compared gene expression patterns in the blood of 23 children with acute KD and 18 age-matched febrile children with 3 illnesses that resemble KD. RESULTS:Genes associated with platelet and neutrophil activation were expressed at higher levels in patients with KD than in patients with acute adenovirus infections or systemic adverse drug reactions, but levels in patients with KD were not higher than those in patients with scarlet fever. Genes associated with B cell activation were also expressed at higher levels in patients with KD than in control subjects. A striking absence of interferon-stimulated gene expression in patients with KD was confirmed in an independent cohort of patients with KD. Using a set of 38 gene transcripts, we successfully predicted the diagnosis for 21 of 23 patients with KD and 7 of 8 patients with adenovirus infection. CONCLUSIONS:These findings provide insight into the molecular features that distinguish KD from other febrile illnesses and support the feasibility of developing novel diagnostic reagents for KD based on the host response.
| S-EPMC2878183 | biostudies-literature
Project description:Flavonoids are secondary metabolites present in all terrestrial plants. The flavonoid pathway has been extensively studied, and many of the involved genes and metabolites have been described in the literature. Despite this extensive knowledge, the functioning of the pathway in vivo is still poorly understood. Here, we study the flavonoid pathway using both experiments and mathematical models. We measured flavonoid metabolite dynamics in two tissues, hypocotyls and cotyledons, during tomato seedling development. Interestingly, the same backbone of interactions leads to very different accumulation patterns in the different tissues. Initially, we developed a mathematical model with constant enzyme concentrations that described the metabolic networks separately in both tissues. This model was unable to fit the measured flavonoid dynamics in the hypocotyls, even if we allowed unrealistic parameter values. This suggested us to investigate the effect of transcript abundance on flavonoid accumulation. We found that the expression of candidate flavonoid genes varies considerably with time. Variation in transcript abundance results in enzymatic variation, which could have a large effect on metabolite accumulation. Candidate transcript abundance was included in the mathematical model as representative for enzyme concentration. We fitted the resulting model to the flavonoid dynamics in the cotyledons, and tested it by applying it to the data from hypocotyls. When transcript abundance is included, we are indeed able to explain flavonoid dynamics in both tissues. Importantly, this is possible under the biologically relevant restriction that the enzymatic properties estimated by the model are conserved between the tissues.
| S-EPMC3724892 | biostudies-literature
Project description:Developmental dynamics in Boolean models of gene networks self-organize, either into point attractors (stable repeating patterns of gene expression) or limit cycles (stable repeating sequences of patterns), depending on the network interactions specified by a genome of evolvable bits. Genome specifications for dynamics that can map specific gene expression patterns in early development onto specific point attractor patterns in later development are essentially impossible to discover by chance mutation alone, even for small networks. We show that selection for approximate mappings, dynamically maintained in the states comprising limit cycles, can accelerate evolution by at least an order of magnitude. These results suggest that self-organizing dynamics that occur within lifetimes can, in principle, guide natural selection across lifetimes.
| S-EPMC6856163 | biostudies-literature
Project description:BACKGROUND: Isatis indigotica Fort. is one of the most commonly used traditional Chinese medicines. Its antiviral compound is a kind of lignan, which is formed with the action of dirigent proteins (DIR). DIR proteins are members of a large family of proteins which impart stereoselectivity on the phenoxy radical-coupling reaction, yielding optically active lignans from two molecules of E-coniferyl alcohol. They exist in almost every vascular plant. However, the DIR and DIR-like protein gene family in I. indigotica has not been analyzed in detail yet. This study focuses on discovery and analysis of this protein gene family in I. indigotica for the first time. RESULTS: Analysis of transcription profiling database from I. indigotica revealed a family of 19 full-length unique DIR and DIR-like proteins. Sequence analysis found that I. indigotica DIR and DIR-like proteins (IiDIR) were all-beta strand proteins, with a signal peptide at the N-terminus. Phylogenetic analysis of the 19 proteins indicated that the IiDIR genes cluster into three distinct subfamilies, DIR-a, DIR-b/d, and DIR-e, of a larger plant DIR and DIR-like gene family. Gene-specific primers were designed for 19 unique IiDIRs and were used to evaluate patterns of constitutive expression in different organs. It showed that most IiDIR genes were expressed comparatively higher in roots and flowers than stems and leaves. CONCLUSIONS: New DIR and DIR-like proteins were discovered from the transcription profiling database of I. indigotica through bioinformatics methods for the first time. Sequence characteristics and transcript abundance of these new genes were analyzed. This study will provide basic data necessary for further studies.
| S-EPMC4052678 | biostudies-literature
Project description:Suspended animation (e.g. hibernation, diapause) allows organisms to survive extreme environments. But the mechanisms underlying the evolution of suspended animation states are unknown. The African turquoise killifish has evolved diapause as a form of suspended development to survive the complete drought that occurs every summer. Here, we show that gene duplicates – paralogs – exhibit specialized expression in diapause compared to normal development in the African turquoise killifish. Surprisingly, paralogs with specialized expression in diapause are evolutionarily very ancient and are present even in vertebrates that do not exhibit diapause. To determine if evolution of diapause is due to the regulatory landscape rewiring at ancient paralogs, we assessed chromatin accessibility genome-wide in fish species with or without diapause. This analysis revealed an evolutionary recent increase in chromatin accessibility at very ancient paralogs in African turquoise killifish. The increase in chromatin accessibility is linked to the presence of new binding sites for transcription factors, likely due to de novo mutations and transposable element (TE) insertion. Interestingly, accessible chromatin regions in diapause are enriched for lipid metabolism genes, and our lipidomics studies uncover a striking difference in lipid species in African turquoise killifish diapause, which could be critical for long-term survival. Together, our results show that diapause likely originated by repurposing pre-existing gene programs via recent changes in the regulatory landscape. This work raises the possibility that suspended animation programs could be reactivated in other species for long-term preservation via transcription factor remodeling and suggests a mechanism for how complex adaptations evolve in nature.
2021-08-05 | ST001898 | MetabolomicsWorkbench
Project description:Codon usage bias (CUB) refers to the phenomena that synonymous codons are used in different frequencies in most genes and organisms. The general assumption is that codon biases reflect a balance between mutational biases and natural selection. Today we understand that the codon content is related and can affect all gene expression steps. Starting from the 1980s, codon-based indices have been used for answering different questions in all biomedical fields, including systems biology, agriculture, medicine, and biotechnology. In general, codon usage bias indices weigh each codon or a small set of codons to estimate the fitting of a certain coding sequence to a certain phenomenon (e.g., bias in codons, adaptation to the tRNA pool, frequencies of certain codons, transcription elongation speed, etc.) and are usually easy to implement. Today there are dozens of such indices; thus, this paper aims to review and compare the different codon usage bias indices, their applications, and advantages. In addition, we perform analysis that demonstrates that most indices tend to correlate even though they aim to capture different aspects. Due to the centrality of codon usage bias on different gene expression steps, it is important to keep developing new indices that can capture additional aspects that are not modeled with the current indices.
| S-EPMC8122159 | biostudies-literature
Project description:BackgroundConvergent and parallel evolution provide unique insights into the mechanisms of natural selection. Some of the most striking convergent and parallel (collectively recurrent) amino acid substitutions in proteins are adaptive, but there are also many that are selectively neutral. Accordingly, genome-wide assessment has shown that recurrent sequence evolution in orthologs is chiefly explained by nearly neutral evolution. For paralogs, more frequent functional change is expected because additional copies are generally not retained if they do not acquire their own niche. Yet, it is unknown to what extent recurrent sequence differentiation is discernible after independent gene duplications in different eukaryotic taxa.ResultsWe develop a framework that detects patterns of recurrent sequence evolution in duplicated genes. This is used to analyze the genomes of 90 diverse eukaryotes. We find a remarkable number of families with a potentially predictable functional differentiation following gene duplication. In some protein families, more than ten independent duplications show a similar sequence-level differentiation between paralogs. Based on further analysis, the sequence divergence is found to be generally asymmetric. Moreover, about 6% of the recurrent sequence evolution between paralog pairs can be attributed to recurrent differentiation of subcellular localization. Finally, we reveal the specific recurrent patterns for the gene families Hint1/Hint2, Sco1/Sco2 and vma11/vma3.ConclusionsThe presented methodology provides a means to study the biochemical underpinning of functional differentiation between paralogs. For instance, two abundantly repeated substitutions are identified between independently derived Sco1 and Sco2 paralogs. Such identified substitutions allow direct experimental testing of the biological role of these residues for the repeated functional differentiation. We also uncover a diverse set of families with recurrent sequence evolution and reveal trends in the functional and evolutionary trajectories of this hitherto understudied phenomenon.
| S-EPMC7414715 | biostudies-literature
Project description:Understanding the mechanisms leading to new traits or additional features in organisms is a fundamental goal of evolutionary biology. We show that HOXDB regulatory changes have been used repeatedly in different fish genera to alter the length and number of the prominent dorsal spines used to classify stickleback species. In Gasterosteus aculeatus (typically 'three-spine sticklebacks'), a variant HOXDB allele is genetically linked to shortening an existing spine and adding an additional spine. In Apeltes quadracus (typically 'four-spine sticklebacks'), a variant HOXDB allele is associated with lengthening a spine and adding an additional spine in natural populations. The variant alleles alter the same non-coding enhancer region in the HOXDB locus but do so by diverse mechanisms, including single-nucleotide polymorphisms, deletions and transposable element insertions. The independent regulatory changes are linked to anterior expansion or contraction of HOXDB expression. We propose that associated changes in spine lengths and numbers are partial identity transformations in a repeating skeletal series that forms major defensive structures in fish. Our findings support the long-standing hypothesis that natural Hox gene variation underlies key patterning changes in wild populations and illustrate how different mutational mechanisms affecting the same region may produce opposite gene expression changes with similar phenotypic outcomes.
| S-EPMC9525239 | biostudies-literature
Project description:BackgroundThe accuracy of transcript quantification using RNA-seq data depends on many factors, such as the choice of alignment or mapping method and the quantification model being adopted. While the choice of quantification model has been shown to be important, considerably less attention has been given to comparing the effect of various read alignment approaches on quantification accuracy.ResultsWe investigate the influence of mapping and alignment on the accuracy of transcript quantification in both simulated and experimental data, as well as the effect on subsequent differential expression analysis. We observe that, even when the quantification model itself is held fixed, the effect of choosing a different alignment methodology, or aligning reads using different parameters, on quantification estimates can sometimes be large and can affect downstream differential expression analyses as well. These effects can go unnoticed when assessment is focused too heavily on simulated data, where the alignment task is often simpler than in experimentally acquired samples. We also introduce a new alignment methodology, called selective alignment, to overcome the shortcomings of lightweight approaches without incurring the computational cost of traditional alignment.ConclusionWe observe that, on experimental datasets, the performance of lightweight mapping and alignment-based approaches varies significantly, and highlight some of the underlying factors. We show this variation both in terms of quantification and downstream differential expression analysis. In all comparisons, we also show the improved performance of our proposed selective alignment method and suggest best practices for performing RNA-seq quantification.
| S-EPMC7487471 | biostudies-literature