Small genomic insertions form enhancers that misregulate oncogenes
Ontology highlight
ABSTRACT: ChIP-Seq data from 102 cell lines were processed to identify insertions relative to the human reference genome. To confirm the ability of a computational pipeline to identify insertions that are actually present in these genomes, we PCR-amplified and deeply sequenced regions where insertions were predicted to exist in MOLT4 T-cell acute lymphoblastic leukemia cells. The current study release makes available sequences of amplicons from a tumor cell genome (MOLT4) where a computational pipeline predicted small insertions from orthogonal data.
SECONDARY ACCESSION(S): PRJNA353149PRJNA353150
REPOSITORIES: dbGaP
Similar Datasets
Project description:The noncoding regions of tumor cell genomes harbor a considerable fraction of total variation, but the functional contribution of such variants to tumorigenesis is ill-defined. Among these noncoding variants, somatic insertions are among the least well-characterized due to challenges with interpreting short-read DNA sequences. To address this, we have identified enhancer-associated small insertion variants across a spectrum of tumor cells by using ChIP-seq to enrich and sequence enhancer DNA, and a computational approach with multiple alignment procedures to identify small insertions. Among the 82 tumor genomes studied here, small insertions were frequently observed in enhancer DNA sequences located near known oncogenes. Further study of one such insertion, found in primary leukemia tumor genomes, revealed that it nucleates formation of an active enhancer that drives expression of the LMO2 oncogene. The catalogue of enhancer-associated small insertion variants described here provides a foundation for further investigating the functions of this class of variants across a range of human cancers.
2016-07-08 | GSE76783 | GEO
Project description:Purpose: The goals of this study are to identify Mutator insertion sites and determine the affect of these insertions on nearby genes Methods: Wiseq-based Mu element profiling was performed as described previously using B73/Mo17 F1 hybrid progeny seedlings. The B73 parent was carried Mutator activity that had been introgressed into the B73 genetic background. The Mo17 parent lacked active Mu elements. Thus, all new insertions were into the B73 genome. Genomic DNA was extracted from 6-day-old seedlings of B73/Mo17 hybrid plants. Amplicon-based enrichment of Mu flanking DNA was then performed. The purified PCR products were subject to Miseq-based Wideseq pipeline at Purdue Genomics Core Facility (https://www.purdue.edu/hla/sites/genomics/wideseq-2/). Wideseq reads were mapped to the B73 reference genome as described previously. By identifying the Mu target site duplications (TSDs), a set of genes targeted by Mu insertions that segregated in hybrid progeny was obtained and those containing B73/Mo17 SNPs in their mRNA sequences were used for allele-specific expression analysis. Because new insertions were into the B73 genome, the effect of these insertions would be expected to be specific in all cases to the B73 allele. To quantify the allele frequency, we performed RT-PCR followed by Wideseq from the identical shoot tissues of the hybrid seedlings mentioned above. Total RNA was extracted using the RNA Extraction Kit (Zymo) and cDNAs were synthesized using Promega M-MLV Reverse Transcriptase. For for a subset of 16 genes that carried SNPs, RNA fragments containing B73/Mo17 SNPs were amplified by RT-PCR. Libraries were prepared using an mRNA preparation kit (Illumina) following a standard protocol (https://www.purdue.edu/hla/sites/genomics/wi-deseq-2/). The RT-PCR products were then sequenced by the “WideSeq” pipeline. Results: A knockdown index was deduced by normalizing the observed ratio to that observed in genes that lacked Mu insertions in both genetic backgrounds for each insertion. We found that none of the four promoter insertions changed the expression of nearby genes. A quarter of 5’ UTR insertions (5 out of 20) caused knockout or strong knockdown effects and one of seven intronic insertions resulted in a knockout effect. Collectively, of a total of 33 Mu insertions, all of which were within 200 bp of genes, only 11 significantly reduced gene expression, and only two knocked it out completely. Conclusions: These results indicate that that Mu element insertions near TSSs of host genes are often associated with quantitative and in many cases neglectable functional consequences on nearby gene expression.
2020-05-09 | GSE146647 | GEO
Project description:The wild grass Brachypodium distachyon has emerged as a model system for temperate grasses and biofuel plants. However, the global analysis of miRNAs, molecules known to be key for eukaryotic gene regulation, has been limited in B. distachyon to studies examining a few samples or that rely on computational predictions. Similarly an in-depth global analysis of miRNA-mediated target cleavage using Parallel Analysis of RNA Ends (PARE) data is lacking in B. distachyon. B. distachyon small RNAs were cloned and deeply sequenced from 17 libraries that represent different tissues and stresses. Using a computational pipeline, we identified 116 miRNAs including not only conserved miRNAs that have not been reported in B. distachyon, but also non-conserved miRNAs that were not found in other plants. To investigate miRNA-mediated cleavage function, four PARE libraries were constructed from key tissues and sequenced to a total depth of approximately 70 million sequences. The roughly 5 million distinct genome-matched sequences that resulted represent an extensive dataset to analyze small RNA-guided cleavage events. Analysis of the PARE and miRNA data provided experimental evidence for miRNA-mediated cleavage of 264 sites in predicted miRNA targets. In addition, PARE analysis revealed that differentially expressed miRNAs in the same family guide specific target RNA cleavage in a correspondingly tissue-preferential manner. B. distachyon miRNAs and target RNAs were experimentally identified and analyzed. Knowledge gained from this study should provide insights into the roles of miRNAs and the regulation of their targets in B. distachyon and related plants. Examination of various tissues and stresses in Brachypodium by high throughput sequencing for small RNA profiling and PARE (Parallel Analysis of RNA Ends)
2013-12-06 | E-GEOD-52441 | biostudies-arrayexpress
Project description:There is increasing evidence that transcripts or transcript regions annotated as non-coding can harbor functional short open reading frames (sORFs). Loss-of-function experiments have identified essential developmental or physiological roles for a few of the encoded peptides (micropeptides), but genome-wide experimental or computational identification of functional sORFs remains challenging. In this project we expand our previously developed method and present results of an integrated computational pipeline for the identification of conserved sORFs in human, mouse, zebrafish, fruit fly and the nematode C. elegans. This dataset represents the peaklist of identified novel and annotated mouse micropeptides by mass spectrometry.
2015-07-28 | PXD002583 | Pride
Project description:Recent advances in transcriptome sequencing have enabled the discovery of thousands of long non-coding RNAs (lncRNAs) across many species. Though several lncRNAs have been shown to play important roles in diverse biological processes, the functions and mechanisms of most lncRNAs remain unknown. Two significant obstacles lie between transcriptome sequencing and functional characterization of lncRNAs: identifying truly non-coding genes from de novo reconstructed transcriptomes, and prioritizing the hundreds of resulting putative lncRNAs for downstream experimental interrogation. We present slncky, a computational lncRNA discovery tool that produces a high-quality set of lncRNAs from RNA-sequencing data and further uses evolutionary constraint to prioritize lncRNAs that are likely to be functionally important. Our automated filtering pipeline is comparable to manual curation efforts and more sensitive than previously published computational approaches. Furthermore, we develop a sensitive alignment pipeline for aligning lncRNA loci and propose new evolutionary metrics relevant for analyzing sequence and transcript evolution. Our analysis reveals that evolutionary selection acts in several distinct patterns, and uncovers two notable classes of intergenic lncRNAs: one showing strong purifying selection on RNA sequence and another where constraint is restricted to the regulation but not the sequence of the transcript. To study a comprehensive and comparable set of lncRNAs expressed in the pluripotent state, we analyzed RNA-Seq data from pluripotent cells derived from several strains and species, and grown in two physiological conditions. First we derived ânaïveâ ES cells (ESCs) from each of three different mice strains: 129SvEv, NOD, and Mus musculus castaneus (cast) mouse, a wild mouse subspecies originally from Thailand, as well as naïve induced pluripotent stem (iPS) cells from rat and human. Next, to facilitate comparisons across pluripotent cells from different species, we also cultured mouse and rat cells in âprimedâ epiblast stem cell (EpiSC) culture conditions, since human iPS cells in culture are much more similar molecularly and functionally to mouse primed EpiSCs rather than to a ground state naïve ESCs. We polyA selected RNA from each cell type and deeply sequenced on HiSeq2500
2016-01-07 | E-GEOD-64818 | biostudies-arrayexpress
Project description:We analyze whether tRNA quantifications based on small RNA-seq data are informative enough to distinguish between different human cell lines covering multiple tissue types. We therefore apply both small RNA-seq and Hydro-tRNAseq to HEK293 (kidney), HCT116 (colon), HeLa (cervix), MDA-MB-231 (breast), and BJ fibroblasts. First, the correlations between the two methods of identical samples and computational mapping pipeline range between 0.93 and 0.96 for all cell lines. tRNA quantifications from both protocols are compared and significantly higher Spearman correlations are obtained within matching samples versus mismatching cell lines. In consequence, we demonstrate that small RNA-seq quantifications of sample-specific tRNA profiles show a good agreement with conventional tRNA-seq.
2020-02-28 | GSE137834 | GEO
Project description:MicroRNAs are a class of small non-coding RNAs that regulate mRNA expression at the post-transcriptional level and thereby many fundamental biological processes. A number of methods, such as multiplex polymerase chain reaction, microarrays have been developed for profiling levels of known miRNAs. These methods lack ability to identify novel miRNAs and accurately determine expression at a range of concentration. Deep or massively parallel sequencing methods are providing suitable platforms for genome wide transcriptome analysis and have the ability to identify novel transcripts. The results of analysis of small RNA sequences obtained by Solexa technology of normal peripheral blood mononuclear cells, tumor cell lines K562 (chronic myelogenous leukemia) and HL60 (acute promyelogenous leukemia) are presented. Custom computation pipelines were used to generate expression profiles of known and for identification of novel miRNAs. Some of the highly expressed miRNAs in the leukocytes include several members of the let 7 family, mir-21, 103, 185, 191 and 320a. Comparison of the miRNA profiles of normal versus K562 cells or HL60 revealed a specific set of differentially expressed molecules. Correlation of the miRNA with that of mRNA expression profiles, obtained by microarray, revealed a set of target genes showing inverse correlation with miRNA levels. Our computational pipeline also predicted a number of novel miRNAs. Some of the predictions were validated by realtime RT-PCR and or RNAase protection assay. The small RNA population from four samples - Two Normal Peripheral blood mononuclear cells (PBMC) samples, K562 cell line (This cell line is used as a model to study Chronic Myelogenous Leukemia), HL60 (This cell line is used to study acute promyelogenous leukemia) was sequenced using Solexa technology.
2010-05-13 | E-GEOD-19812 | biostudies-arrayexpress
Project description:CRC Organoid and corresponding PDX data are part of a larger project that includes proteomeXchange accessions PXD006616, PXD008032, and PXD009995. We present a computational pipeline called Integrative Inferred Kinase Activity (INKA) scoring that integrates the observed abundance of phosphopeptides derived from (i) kinases as a whole, (ii) their kinase activation loop segments, (iii) proteins deduced to be kinase substrates based on prior experimental knowledge of kinase-substrate relationships, and/or (iv) kinase substrates predicted by the sequence motif- and network-based NetworKIN algorithm.
2019-04-18 | PXD012565 | Pride
Project description:We describe a series of computational pipelines for the in silico analysis of small RNAs (sRNA) produced in response to viral infections in plants. Our workflow is primarily focused on the analysis of sRNA populations derived from known or previously undescribed viruses infecting host plants. Furthermore, we provide an additional pipeline to examine host-specific endogenous sRNAs activated or specifically expressed during viral infections in plants. We present some key points for a successful and cost-efficient processing of next generation sequencing sRNA libraries, from purification of high quality RNA to guidance for library preparation and sequencing strategies. We report a series of free available tools and programs as well as in-house Perl scripts to perform in-house sRNA-seq data mining. A multi-step analysis pipeline is extensively detailed so previous bioinformatic background is not required, but experience with basic Unix commands is desirable.
2017-10-12 | GSE104799 | GEO
Project description:The wild grass Brachypodium distachyon has emerged as a model system for temperate grasses and biofuel plants. However, the global analysis of miRNAs, molecules known to be key for eukaryotic gene regulation, has been limited in B. distachyon to studies examining a few samples or that rely on computational predictions. Similarly an in-depth global analysis of miRNA-mediated target cleavage using Parallel Analysis of RNA Ends (PARE) data is lacking in B. distachyon. B. distachyon small RNAs were cloned and deeply sequenced from 17 libraries that represent different tissues and stresses. Using a computational pipeline, we identified 116 miRNAs including not only conserved miRNAs that have not been reported in B. distachyon, but also non-conserved miRNAs that were not found in other plants. To investigate miRNA-mediated cleavage function, four PARE libraries were constructed from key tissues and sequenced to a total depth of approximately 70 million sequences. The roughly 5 million distinct genome-matched sequences that resulted represent an extensive dataset to analyze small RNA-guided cleavage events. Analysis of the PARE and miRNA data provided experimental evidence for miRNA-mediated cleavage of 264 sites in predicted miRNA targets. In addition, PARE analysis revealed that differentially expressed miRNAs in the same family guide specific target RNA cleavage in a correspondingly tissue-preferential manner. B. distachyon miRNAs and target RNAs were experimentally identified and analyzed. Knowledge gained from this study should provide insights into the roles of miRNAs and the regulation of their targets in B. distachyon and related plants.
2013-12-06 | GSE52441 | GEO