Small RNA-sequencing (RNA-Seq) is being increasingly used for profiling of circulating microRNAs (miRNAs), a new group of promising biomarkers. Unfortunately, small RNA-Seq protocols are prone to biases limiting quantification accuracy, which motivated development of several novel methods. Here, we present comparison of all small RNA-Seq library preparation approaches that are commercially available for quantification of miRNAs in biofluids. Using synthetic and human plasma samples, we compared performance of traditional two-adaptor ligation protocols (Lexogen, Norgen), as well as methods using randomized adaptors (NEXTflex), polyadenylation (SMARTer), circularization (RealSeq), capture probes (EdgeSeq), or unique molecular identifiers (QIAseq). There was no single protocol outperforming others across all metrics. Limited overlap of measured miRNA profiles was documented between methods largely owing to protocol-specific biases. Methods designed to minimize bias largely differ in their performance, and contributing factors were identified. Usage of unique molecular identifiers has rather negligible effect and, if designed incorrectly, can even introduce spurious results. Together, these results identify strengths and weaknesses of all current methods and provide guidelines for applications of small RNA-Seq in biomarker research.
The role of alternative promoter usage in tissue-specific gene expression has been well established; however, its role in complex diseases is poorly understood. We performed cap analysis of gene expression (CAGE) sequencing from the left ventricle of a rat model of hypertension, the spontaneously hypertensive rat (SHR), and a normotensive strain, Brown Norway to understand the role of alternative promoter usage in complex disease. We identified 26,560 CAGE-defined transcription start sites in the rat left ventricle, including 1,970 novel cardiac transcription start sites. We identified 28 genes with alternative promoter usage between SHR and Brown Norway, which could lead to protein isoforms differing at the amino terminus between two strains and 475 promoter switching events altering the length of the 5' UTR. We found that the shift in Insr promoter usage was significantly associated with insulin levels and blood pressure within a panel of HXB/BXH recombinant inbred rat strains, suggesting that hyperinsulinemia due to insulin resistance might lead to hypertension in SHR. Our study provides a preliminary evidence of alternative promoter usage in complex diseases.
- MeSH
- genetická transkripce genetika MeSH
- hypertenze * genetika metabolismus MeSH
- krysa rodu rattus MeSH
- potkani inbrední SHR MeSH
- promotorové oblasti (genetika) genetika MeSH
- sekvenční analýza RNA metody MeSH
- stanovení celkové genové exprese metody MeSH
- zvířata MeSH
- Check Tag
- krysa rodu rattus MeSH
- mužské pohlaví MeSH
- ženské pohlaví MeSH
- zvířata MeSH
- Publikační typ
- časopisecké články MeSH
- MeSH
- alternativní sestřih * MeSH
- DNA vazebné proteiny genetika MeSH
- fúzní onkogenní proteiny genetika MeSH
- lidé MeSH
- lipoblastom genetika patologie MeSH
- sekvenční analýza RNA metody MeSH
- transkripční faktor Zeb2 genetika MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- dopisy MeSH
- kazuistiky MeSH
- práce podpořená grantem MeSH
Extracellular RNAs present in biofluids have emerged as potential biomarkers for disease. Where most studies focus on blood-derived fluids, other biofluids may be more informative. We present an atlas of messenger, circular, and small RNA transcriptomes of a comprehensive collection of 20 human biofluids. By means of synthetic spike-in controls, we compare RNA content across biofluids, revealing a 10,000-fold difference in concentration. The circular RNA fraction is increased in most biofluids compared to tissues. Each biofluid transcriptome is enriched for RNA molecules derived from specific tissues and cell types. Our atlas enables an informed selection of the most relevant biofluid to monitor particular diseases. To verify the biomarker potential in these biofluids, four validation cohorts representing a broad spectrum of diseases were profiled, revealing numerous differential RNAs between case and control subjects. Spike-normalized data are publicly available in the R2 web portal for further exploration.
- MeSH
- biologické markery * MeSH
- kohortové studie MeSH
- lidé MeSH
- RNA genetika metabolismus MeSH
- sekvenční analýza RNA metody MeSH
- stanovení celkové genové exprese metody MeSH
- tělesné tekutiny metabolismus MeSH
- transkriptom * MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
- Research Support, N.I.H., Extramural MeSH
BACKGROUND/AIM: Rectal cancer accounts for approximately one-third of all colorectal cancers. Currently, the standard treatment for locally advanced rectal cancer (LARC) is neoadjuvant chemoradiotherapy (CRT) with capecitabine or 5-fluorouracil followed by curative surgery. Unfortunately, only 20% of patients with LARC present complete pathological response after CRT, whereas in 20-40% cases the response is poor or absent. The aim of our study was to evaluate whether microRNAs (miRNAs) in tumor biopsy specimen have the potential to predict therapeutic response in LARC patients. PATIENTS AND METHODS: In total 87 LARC patients treated by CRT were enrolled in our prospective study. To identify predictive miRNAs, we used small RNA sequencing in 40 tumor biopsy samples of LARC patients (20 responders, 20 non-responders) and qPCR validation of selected miRNA candidates. RESULTS: In the discovery phase of the study, we identified 69 miRNAs to have significantly different expression between the group of responders (TRG 1,2) and a group of non-responders (TRG 4,5) to neoadjuvant CRT. Among these miRNAs, 48 showed a lower expression and 21 showed higher expression in tumor tissues from poorly responding LARC patients. Five miRNAs were selected for validation, but only miR-487a-3p was confirmed to have a significantly higher expression in the tumor biopsy specimens of non-responders to neoadjuvant CRT (p<0.0006, AUC=0.766). Gene Ontology (GO) clustering and pathway enrichment analysis of the miR-487a-3p mRNA targets, revealed potential mechanisms behind miR-487a-3p roles in chemoradioresistance (e.g. TGF-beta signaling pathway, protein kinase activity, double-stranded DNA binding, or microRNAs in cancer). CONCLUSION: By combination of miRNA expression profiling and integrative computational biology we identified miR-487a-3p as a potential predictive biomarker of CRT response in LARC patients.
- MeSH
- adenokarcinom genetika patologie terapie MeSH
- chemoradioterapie MeSH
- dospělí MeSH
- lidé středního věku MeSH
- lidé MeSH
- malá nekódující RNA genetika MeSH
- mikro RNA genetika MeSH
- nádorové biomarkery genetika MeSH
- nádory rekta genetika patologie terapie MeSH
- neoadjuvantní terapie MeSH
- prospektivní studie MeSH
- ROC křivka MeSH
- sekvenční analýza RNA metody MeSH
- senioři nad 80 let MeSH
- senioři MeSH
- výsledek terapie MeSH
- Check Tag
- dospělí MeSH
- lidé středního věku MeSH
- lidé MeSH
- mužské pohlaví MeSH
- senioři nad 80 let MeSH
- senioři MeSH
- ženské pohlaví MeSH
- Publikační typ
- časopisecké články MeSH
Recent technological advances have made next-generation sequencing (NGS) a popular and financially accessible technique allowing a broad range of analyses to be done simultaneously. A huge amount of newly generated NGS data, however, require advanced software support to help both in analyzing the data and biologically interpreting the results. In this article, we describe SATrans (Software for Annotation of Transcriptome), a software package providing fast and robust functional annotation of novel sequences obtained from transcriptome sequencing. Moreover, it performs advanced gene ontology analysis of differentially expressed genes, thereby helping to interpret biologically-and in a user-friendly form-the quantitative changes in gene expression. The software is freely available and provides the possibility to work with thousands of sequences using a standard personal computer or notebook running on the Linux operating system.
Uterine leiomyosarcomas (uLMS) are rare, aggressive malignancies for which limited treatment options are available. To gain novel molecular insights into uLMS and identify potential novel therapeutic targets, we characterized 84 uLMS samples for genome-wide somatic copy number alterations, mutations, gene fusions and gene expression and performed a data integration analysis. We found that alterations affecting TP53, RB1, PTEN, MED12, YWHAE and VIPR2 were present in the majority of uLMS. Pathway analyses additionally revealed that the PI3K/AKT/mTOR, estrogen-mediated S-phase entry and DNA damage response signaling pathways, for which inhibitors have already been developed and approved, frequently harbored genetic changes. Furthermore, a significant proportion of uLMS was characterized by amplifications and overexpression of known oncogenes (CCNE1, TDO2), as well as deletions and reduced expression of tumor suppressor genes (PTEN, PRDM16). Overall, it emerged that the most frequently affected gene in our uLMS samples was VIPR2 (96%). Interestingly, VIPR2 deletion also correlated with unfavorable survival in uLMS patients (multivariate analysis; HR = 4.5, CI = 1.4-14.3, p = 1.2E-02), while VIPR2 protein expression was reduced in uLMS vs. normal myometrium. Moreover, stimulation of VIPR2 with its natural agonist VIP decreased SK-UT-1 uLMS cell proliferation in a dose-dependent manner. These data suggest that VIPR2, which is a negative regulator of smooth muscle cell proliferation, might be a novel tumor suppressor gene in uLMS. Our work further highlights the importance of integrative molecular analyses, through which we were able to uncover the genes and pathways most frequently affected by somatic alterations in uLMS.
- MeSH
- dospělí MeSH
- genomika metody MeSH
- Kaplanův-Meierův odhad MeSH
- karcinogeneze genetika MeSH
- leiomyosarkom genetika mortalita patologie terapie MeSH
- lidé středního věku MeSH
- lidé MeSH
- myometrium patologie MeSH
- nádory dělohy genetika mortalita patologie terapie MeSH
- onkogeny genetika MeSH
- proliferace buněk genetika MeSH
- receptory vazoaktivního intestinálního peptidu typu II genetika MeSH
- regulace genové exprese u nádorů genetika MeSH
- sekvenční analýza RNA metody MeSH
- sekvenování celého genomu metody MeSH
- senioři nad 80 let MeSH
- senioři MeSH
- signální transdukce genetika MeSH
- tumor supresorové geny MeSH
- vysoce účinné nukleotidové sekvenování metody MeSH
- Check Tag
- dospělí MeSH
- lidé středního věku MeSH
- lidé MeSH
- senioři nad 80 let MeSH
- senioři MeSH
- ženské pohlaví MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
Motivation: Sanger sequencing is still being employed for sequence variant detection by many laboratories, especially in a clinical setting. However, chromatogram interpretation often requires manual inspection and in some cases, considerable expertise. Results: We present GLASS, a web-based Sanger sequence trace viewer, editor, aligner and variant caller, built to assist with the assessment of variations in 'curated' or user-provided genes. Critically, it produces a standardized variant output as recommended by the Human Genome Variation Society. Availability and implementation: GLASS is freely available at http://bat.infspire.org/genomepd/glass/ with source code at https://github.com/infspiredBAT/GLASS. Contact: nikos.darzentas@gmail.com or malcikova.jitka@fnbrno.cz. Supplementary information: Supplementary data are available at Bioinformatics online.
