Gene set enrichment tests (a.k.a. functional enrichment analysis) are among the most frequently used methods in computational biology. Despite this popularity, there are concerns that these methods ...are being applied incorrectly and the results of some peer-reviewed publications are unreliable. These problems include the use of inappropriate background gene lists, lack of false discovery rate correction and lack of methodological detail. To ascertain the frequency of these issues in the literature, we performed a screen of 186 open-access research articles describing functional enrichment results. We find that 95% of analyses using over-representation tests did not implement an appropriate background gene list or did not describe this in the methods. Failure to perform p-value correction for multiple tests was identified in 43% of analyses. Many studies lacked detail in the methods section about the tools and gene sets used. An extension of this survey showed that these problems are not associated with journal or article level bibliometrics. Using seven independent RNA-seq datasets, we show misuse of enrichment tools alters results substantially. In conclusion, most published functional enrichment studies suffered from one or more major flaws, highlighting the need for stronger standards for enrichment analysis.
Inference of biological pathway activity via gene set enrichment analysis is frequently used in the interpretation of clinical and other omics data. With the proliferation of new omics profiling ...approaches and ever-growing size of data sets generated, there is a lack of tools available to perform and visualise gene set enrichments in analyses involving multiple contrasts.
To address this, we developed mitch, an R package for multi-contrast gene set enrichment analysis. It uses a rank-MANOVA statistical approach to identify sets of genes that exhibit joint enrichment across multiple contrasts. Its unique visualisation features enable the exploration of enrichments in up to 20 contrasts. We demonstrate the utility of mitch with case studies spanning multi-contrast RNA expression profiling, integrative multi-omics, tool benchmarking and single-cell RNA sequencing. Using simulated data we show that mitch has similar accuracy to state of the art tools for single-contrast enrichment analysis, and superior accuracy in identifying multi-contrast enrichments.
mitch is a versatile tool for rapidly and accurately identifying and visualising gene set enrichments in multi-contrast omics data. Mitch is available from Bioconductor ( https://bioconductor.org/packages/mitch ).
BACKGROUND:Dietary intake of fruit and vegetables is associated with lower incidence of hypertension, but the mechanisms involved have not been elucidated. Here, we evaluated the effect of a ...high-fiber diet and supplementation with the short-chain fatty acid acetate on the gut microbiota and the prevention of cardiovascular disease.
METHODS:Gut microbiome, cardiorenal structure/function, and blood pressure were examined in sham and mineralocorticoid excess–treated mice with a control diet, high-fiber diet, or acetate supplementation. We also determined the renal and cardiac transcriptome of mice treated with the different diets.
RESULTS:We found that high consumption of fiber modified the gut microbiota populations and increased the abundance of acetate-producing bacteria independently of mineralocorticoid excess. Both fiber and acetate decreased gut dysbiosis, measured by the ratio of Firmicutes to Bacteroidetes, and increased the prevalence of Bacteroides acidifaciens. Compared with mineralocorticoid-excess mice fed a control diet, both high-fiber diet and acetate supplementation significantly reduced systolic and diastolic blood pressures, cardiac fibrosis, and left ventricular hypertrophy. Acetate had similar effects and markedly reduced renal fibrosis. Transcriptome analyses showed that the protective effects of high fiber and acetate were accompanied by the downregulation of cardiac and renal Egr1, a master cardiovascular regulator involved in cardiac hypertrophy, cardiorenal fibrosis, and inflammation. We also observed the upregulation of a network of genes involved in circadian rhythm in both tissues and downregulation of the renin-angiotensin system in the kidney and mitogen-activated protein kinase signaling in the heart.
CONCLUSIONS:A diet high in fiber led to changes in the gut microbiota that played a protective role in the development of cardiovascular disease. The favorable effects of fiber may be explained by the generation and distribution of one of the main metabolites of the gut microbiota, the short-chain fatty acid acetate. Acetate effected several molecular changes associated with improved cardiovascular health and function.
Gene name errors: Lessons not learned Abeysooriya, Mandhri; Soria, Megan; Kasu, Mary Sravya ...
PLOS computational biology/PLoS computational biology,
07/2021, Letnik:
17, Številka:
7
Journal Article
Recenzirano
Odprti dostop
Erroneous conversion of gene names into other dates and other data types has been a frustration for computational biologists for years. We hypothesized that such errors in supplementary files might ...diminish after a report in 2016 highlighting the extent of the problem. To assess this, we performed a scan of supplementary files published in PubMed Central from 2014 to 2020. Overall, gene name errors continued to accumulate unabated in the period after 2016. An improved scanning software we developed identified gene name errors in 30.9% (3,436/11,117) of articles with supplementary Excel gene lists; a figure significantly higher than previously estimated. This is due to gene names being converted not just to dates and floating-point numbers, but also to internal date format (five-digit numbers). These findings further reinforce that spreadsheets are ill-suited to use with large genomic data.
BACKGROUND:High blood pressure (BP) continues to be a major, poorly controlled but modifiable risk factor for cardiovascular death. Among key Western lifestyle factors, a diet poor in fiber is ...associated with prevalence of high BP. The impact of lack of prebiotic fiber and the associated mechanisms that lead to higher BP are unknown. Here we show that lack of prebiotic dietary fiber leads to the development of a hypertensinogenic gut microbiota, hypertension and its complications, and demonstrate a role for G-protein coupled-receptors (GPCRs) that sense gut metabolites.
METHODS:One hundred seventy-nine mice including C57BL/6J, gnotobiotic C57BL/6J, and knockout strains for GPR41, GPR43, GPR109A, and GPR43/109A were included. C57BL/6J mice were implanted with minipumps containing saline or a slow-pressor dose of angiotensin II (0.25 mg·kg·d). Mice were fed diets lacking prebiotic fiber with or without addition of gut metabolites called short-chain fatty acids (SCFA) produced during fermentation of prebiotic fiber in the large intestine), or high prebiotic fiber diets. Cardiac histology and function, BP, sodium and potassium excretion, gut microbiome, flow cytometry, catecholamines and methylation-wide changes were determined.
RESULTS:Lack of prebiotic fiber predisposed mice to hypertension in the presence of a mild hypertensive stimulus, with resultant pathological cardiac remodeling. Transfer of a hypertensinogenic microbiota to gnotobiotic mice recapitulated the prebiotic-deprived hypertensive phenotype, including cardiac manifestations. Reintroduction of SCFAs to fiber-depleted mice had protective effects on the development of hypertension, cardiac hypertrophy, and fibrosis. The cardioprotective effect of SCFAs were mediated via the cognate SCFA receptors GPR43/GPR109A, and modulated L-3,4-dihydroxyphenylalanine levels and the abundance of T regulatory cells regulated by DNA methylation.
CONCLUSIONS:The detrimental effects of low fiber Westernized diets may underlie hypertension, through deficient SCFA production and GPR43/109A signaling. Maintaining a healthy, SCFA-producing microbiota is important for cardiovascular health.
Abstract
Background
RNA sequencing (RNA-seq) is an indispensable tool in the study of gene regulation. While the technology has brought with it better transcript coverage and quantification, there ...remain considerable barriers to entry for the computational biologist to analyse large data sets. There is a real need for a repository of uniformly processed RNA-seq data that is easy to use.
Findings
To address these obstacles, we developed Digital Expression Explorer 2 (DEE2), a web-based repository of RNA-seq data in the form of gene-level and transcript-level expression counts. DEE2 contains >5.3 trillion assigned reads from 580,000 RNA-seq data sets including species Escherichia coli, yeast, Arabidopsis, worm, fruit fly, zebrafish, rat, mouse, and human. Base-space sequence data downloaded from the National Center for Biotechnology Information Sequence Read Archive underwent quality control prior to transcriptome and genome mapping using open-source tools. Uniform data processing methods ensure consistency across experiments, facilitating fast and reproducible meta-analyses.
Conclusions
The web interface allows users to quickly identify data sets of interest using accession number and keyword searches. The data can also be accessed programmatically using a specifically designed R package. We demonstrate that DEE2 data are compatible with statistical packages such as edgeR or DESeq. Bulk data are also available for download. DEE2 can be found at http://dee2.io.
Aquaporins (AQPs) are major intrinsic proteins (MIPs) that mediate bidirectional flux of water and other substrates across cell membranes, and play critical roles in plant-water relations, ...dehydration stress responses and crop productivity. However, limited data are available as yet on the contributions of these proteins to the physiology of the major crop barley (Hordeum vulgare). The present work reports the identification and expression analysis of the barley MIP family. A comprehensive search of publicly available leaf mRNA-seq data, draft barley genome data, GenBank transcripts and sixteen new annotations together revealed that the barley MIP family is comprised of at least forty AQPs. Alternative splicing events were likely in two plasma membrane intrinsic protein (PIP) AQPs. Analyses of the AQP signature sequences and specificity determining positions indicated a potential of several putative AQP isoforms to transport non-aqua substrates including physiological important substrates, and respond to abiotic stresses. Analysis of our publicly available leaf mRNA-seq data identified notable differential expression of HvPIP1;2 and HvTIP4;1 under salt stress. Analyses of other gene expression resources also confirmed isoform-specific responses in different tissues and/or in response to salinity, as well as some potentially inter-cultivar differences. The work reports systematic and comprehensive analysis of most, if not all, barley AQP genes, their sequences, expression patterns in different tissues, potential transport and stress response functions, and a strong framework for selection and/or development of stress tolerant barley varieties. In addition, the barley data would be highly valuable for genetic studies of the evolutionarily closely related wheat (Triticum aestivum L.).
Objectives
Focal cortical dysplasia (FCD) is a major cause of drug‐resistant focal epilepsy in children, and the clinicopathological classification remains a challenging issue in daily practice. With ...the recent progress in DNA methylation–based classification of human brain tumors we examined whether genomic DNA methylation and gene expression analysis can be used to also distinguish human FCD subtypes.
Methods
DNA methylomes and transcriptomes were generated from massive parallel sequencing in 15 surgical FCD specimens, matched with 5 epilepsy and 6 nonepilepsy controls.
Results
Differential hierarchical cluster analysis of DNA methylation distinguished major FCD subtypes (ie, Ia, IIa, and IIb) from patients with temporal lobe epilepsy patients and nonepileptic controls. Targeted panel sequencing identified a novel likely pathogenic variant in DEPDC5 in a patient with FCD type IIa. However, no enrichment of differential DNA methylation or gene expression was observed in mechanistic target of rapamycin (mTOR) pathway–related genes.
Significance
Our studies extend the evidence for disease‐specific methylation signatures toward focal epilepsies in favor of an integrated clinicopathologic and molecular classification system of FCD subtypes incorporating genomic methylation.
Seaweed extracts are important sources of plant biostimulants that boost agricultural productivity to meet current world demand. The ability of seaweed extracts based on either of the Phaeophyceaean ...species
or
to enhance plant growth or suppress plant disease have recently been shown. However, very limited information is available on the mechanisms of suppression of plant disease by such extracts. In addition, there is no information on the ability of a combination of extracts from
and
to suppress a plant pathogen or to induce plant defense. The present study has explored the transcriptome, using RNA-seq, of
following treatment with extracts from the two species, or a mixture of both, prior to inoculation with the root pathogen
. Following inoculation, five time points (0-24 h post-inoculation) that represented early stages in the interaction of the pathogen with its host were assessed for each treatment and compared with their respective water controls. Wide scale transcriptome reprogramming occurred predominantly related to phytohormone biosynthesis and signaling, changes in metabolic processes and cell wall biosynthesis, there was a broad induction of proteolysis pathways, a respiratory burst and numerous defense-related responses were induced. The induction by each seaweed extract of defense-related genes coincident with the time of inoculation showed that the plants were primed for defense prior to infection. Each seaweed extract acted differently in inducing plant defense-related genes. However, major systemic acquired resistance (SAR)-related genes as well as salicylic acid-regulated marker genes (
,
, and
) and auxin associated genes were found to be commonly up-regulated compared with the controls following treatment with each seaweed extract. Moreover, each seaweed extract suppressed
growth within the roots of inoculated
by the early induction of defense pathways and likely through ROS-based signaling pathways that were linked to production of ROS. Collectively, the RNA-seq transcriptome analysis revealed the induction by seaweed extracts of suites of genes that are associated with direct or indirect plant defense in addition to responses that require cellular energy to maintain plant growth during biotic stress.
Despite increasing knowledge about the factors involved in the progression of diabetic complications, diabetic kidney disease (DKD) continues to be a major health burden. Current therapies only slow ...but do not prevent the progression of DKD. Thus, there is an urgent need to develop novel therapy to halt the progression of DKD and improve disease prognosis. In our preclinical study where we administered a histone deacetylase (HDAC) inhibitor, valproic acid, to streptozotocin-induced diabetic mice, albuminuria and glomerulosclerosis were attenuated. Furthermore, we discovered that valproic acid attenuated diabetes-induced upregulation of complement C5a receptors, with a concomitant reduction in markers of cellular senescence and senescence-associated secretory phenotype. Interestingly, further examination of mice lacking the C5a receptor 1 (C5aR1) gene revealed that cellular senescence was attenuated in diabetes. Similar results were observed in diabetic mice treated with a C5aR1 inhibitor, PMX53. RNA-sequencing analyses showed that PMX53 significantly regulated genes associated with cell cycle pathways leading to cellular senescence. Collectively, these results for the first time demonstrated that complement C5a mediates cellular senescence in diabetic kidney disease. Cellular senescence has been implicated in the pathogenesis of diabetic kidney disease, thus therapies to inhibit cellular senescence such as complement inhibitors present as a novel therapeutic option to treat diabetic kidney disease.