Cells release RNA-carrying vesicles and membrane-free RNA/protein complexes into the extracellular milieu. Horizontal vesicle-mediated transfer of such shuttle RNA between cells allows dissemination ...of genetically encoded messages, which may modify the function of target cells. Other studies used array analysis to establish the presence of microRNAs and mRNA in cell-derived vesicles from many sources. Here, we used an unbiased approach by deep sequencing of small RNA released by immune cells. We found a large variety of small non-coding RNA species representing pervasive transcripts or RNA cleavage products overlapping with protein coding regions, repeat sequences or structural RNAs. Many of these RNAs were enriched relative to cellular RNA, indicating that cells destine specific RNAs for extracellular release. Among the most abundant small RNAs in shuttle RNA were sequences derived from vault RNA, Y-RNA and specific tRNAs. Many of the highly abundant small non-coding transcripts in shuttle RNA are evolutionary well-conserved and have previously been associated to gene regulatory functions. These findings allude to a wider range of biological effects that could be mediated by shuttle RNA than previously expected. Moreover, the data present leads for unraveling how cells modify the function of other cells via transfer of specific non-coding RNA species.
MicroRNAs (miRNAs), non-coding RNAs regulating gene expression, are frequently aberrantly expressed in human cancers. Next-generation deep sequencing technology enables genome-wide expression ...profiling of known miRNAs and discovery of novel miRNAs at unprecedented quantitative and qualitative accuracy. Deep sequencing was performed on 11 fresh frozen clear cell renal cell carcinoma (ccRCC) and adjacent non-tumoral renal cortex (NRC) pairs, 11 additional frozen ccRCC tissues, and 2 ccRCC cell lines (n = 35). The 22 ccRCCs patients belonged to 3 prognostic sub-groups, i.e. those without disease recurrence, with recurrence and with metastatic disease at diagnosis. Thirty-two consecutive samples (16 ccRCC/NRC pairs) were used for stem-loop PCR validation. Novel miRNAs were predicted using 2 distinct bioinformatic pipelines. In total, 463 known miRNAs (expression frequency 1-150,000/million) were identified. We found that 100 miRNA were significantly differentially expressed between ccRCC and NRC. Differential expression of 5 miRNAs was confirmed by stem-loop PCR in the 32 ccRCC/NRC samples. With respect to RCC subgroups, 5 miRNAs discriminated between non-recurrent versus recurrent and metastatic disease, whereas 12 uniquely distinguished non-recurrent versus metastatic disease. Blocking overexpressed miR-210 or miR-27a in cell line SKCR-7 by transfecting specific antagomirs did not result in significant changes in proliferation or apoptosis. Twenty-three previously unknown miRNAs were predicted in silico. Quantitative genome-wide miRNA profiling accurately separated ccRCC from (benign) NRC. Individual differentially expressed miRNAs may potentially serve as diagnostic or prognostic markers or future therapeutic targets in ccRCC. The biological relevance of candidate novel miRNAs is unknown at present.
Genome sequencing projects are discovering millions of genetic variants in humans, and interpretation of their functional effects is essential for understanding the genetic basis of variation in ...human traits. Here we report sequencing and deep analysis of messenger RNA and microRNA from lymphoblastoid cell lines of 462 individuals from the 1000 Genomes Project--the first uniformly processed high-throughput RNA-sequencing data from multiple human populations with high-quality genome sequences. We discover extremely widespread genetic variation affecting the regulation of most genes, with transcript structure and expression level variation being equally common but genetically largely independent. Our characterization of causal regulatory variation sheds light on the cellular mechanisms of regulatory and loss-of-function variation, and allows us to infer putative causal variants for dozens of disease-associated loci. Altogether, this study provides a deep understanding of the cellular mechanisms of transcriptome variation and of the landscape of functional variants in the human genome.
MicroRNAs are small non-coding RNA transcripts that regulate post-transcriptional gene expression. The millions of short sequence reads generated by next generation sequencing technologies make this ...technique explicitly suitable for profiling of known and novel microRNAs. A modification to the small-RNA expression kit (SREK, Ambion) library preparation method for the SOLiD sequencing platform is described to generate microRNA sequencing libraries that are compatible with the Illumina Genome Analyzer.
High quality sequencing libraries can successfully be prepared from as little as 100 ng small RNA enriched RNA. An easy to use perl-based analysis pipeline called E-miR was developed to handle the sequencing data in several automated steps including data format conversion, 3' adapter removal, genome alignment and annotation to non-coding RNA transcripts. The sample preparation and E-miR pipeline were used to identify 37 cardiac enriched microRNAs in stage 16 chicken embryos. Isomir expression profiles between the heart and embryo were highly correlated for all miRNAs suggesting that tissue or cell specific miRNA modifications do not occur.
In conclusion, our alternative sample preparation method can successfully be applied to generate high quality miRNA sequencing libraries for the Illumina genome analyzer.
ABSTRACT
Cytochrome P450 2D6 (CYP2D6) is among the most important genes involved in drug metabolism. Specific variants are associated with changes in the enzyme's amount and activity. Multiple ...technologies exist to determine these variants, like the AmpliChip CYP450 test, Taqman qPCR, or Second‐Generation Sequencing, however, sequence homology between cytochrome P450 genes and pseudogene CYP2D7 impairs reliable CYP2D6 genotyping, and variant phasing cannot accurately be determined using these assays. To circumvent this, we sequenced CYP2D6 using the Pacific Biosciences RSII and obtained high‐quality, full‐length, phased CYP2D6 sequences, enabling accurate variant calling and haplotyping of the entire gene‐locus including exonic, intronic, and upstream and downstream regions. Unphased diplotypes (Roche AmpliChip CYP450 test) were confirmed for 24 of the 25 samples, including gene duplications. Cases with gene deletions required additional specific assays to resolve. In total, 61 unique variants were detected, including variants that had not previously been associated with specific haplotypes. To further aid genomic analysis using standard reference sequences, we have established an LOVD‐powered CYP2D6 gene‐variant database, and added all reference haplotypes and data reported here. We conclude that our CYP2D6 genotyping approach produces reliable CYP2D6 diplotypes and reveals information about additional variants, including phasing and copy‐number variation.
Specific variants in the CYP2D6 gene locus determine the enzyme's amount and activity. We sequenced CYP2D6 using the PacBios RSII and obtained high‐quality, full‐length, phased CYP2D6 sequences, enabling accurate variant calling and resolution of the separate haplotype sequences of the entire gene‐locus including exonic, intronic and up and downstream regions. Diplotypes used for CYP2D6 protein activity prediction were confirmed for 24 of 25 samples, including gene duplications, and, variants that had not previously been associated with specific haplotypes were detected.
RNA-seq is increasingly used for quantitative profiling of small RNAs (for example, microRNAs, piRNAs and snoRNAs) in diverse sample types, including isolated cells, tissues and cell-free biofluids. ...The accuracy and reproducibility of the currently used small RNA-seq library preparation methods have not been systematically tested. Here we report results obtained by a consortium of nine labs that independently sequenced reference, 'ground truth' samples of synthetic small RNAs and human plasma-derived RNA. We assessed three commercially available library preparation methods that use adapters of defined sequence and six methods using adapters with degenerate bases. Both protocol- and sequence-specific biases were identified, including biases that reduced the ability of small RNA-seq to accurately measure adenosine-to-inosine editing in microRNAs. We found that these biases were mitigated by library preparation methods that incorporate adapters with degenerate bases. MicroRNA relative quantification between samples using small RNA-seq was accurate and reproducible across laboratories and methods.
Differentiated derivatives of human pluripotent stem cells in culture are generally phenotypically immature compared to their adult counterparts. Their identity is often difficult to determine with ...certainty because little is known about their human fetal equivalents in vivo. Cellular identity and signaling pathways directing differentiation are usually determined by extrapolating information from either human adult tissue or model organisms, assuming conservation with humans. To resolve this, we generated a collection of human fetal transcriptional profiles at different developmental stages. Moreover, we developed an algorithm, KeyGenes, which uses this dataset to quantify the extent to which next-generation sequencing or microarray data resemble specific cell or tissue types in the human fetus. Using KeyGenes combined with the human fetal atlas, we identified multiple cell and tissue samples unambiguously on a limited set of features. We thus provide a flexible and expandable platform to monitor and evaluate the efficiency of differentiation in vitro.
Display omitted
•NGS-derived transcriptional profiles of human fetal tissues/organs are generated•Algorithm called KeyGenes uses a training set to predict the identity of a test set•KeyGenes using the fetal atlas identifies NGS- and microarray-derived data•KeyGenes is a flexible and expandable platform to monitor stem cell differentiations
In this article, Lopes and colleagues introduce a platform called KeyGenes, which uses a collection of human fetal transcriptional profiles to assign identity and developmental stages to NGS- and microarray-derived data. Furthermore, they show that KeyGenes is flexible and expandable to get the most meaningful output for a given experiment.
Human lymphoid tissues harbor, in addition to CD56
and CD56
natural killer (NK) cells, a third NK cell population: CD69
CXCR6
lymphoid tissue (lt)NK cells. The function and development of ltNK cells ...remain poorly understood. In this study, we performed RNA sequencing on the three NK cell populations derived from bone marrow (BM) and blood. In ltNK cells, 1,353 genes were differentially expressed compared to circulating NK cells. Several molecules involved in migration were downregulated in ltNK cells:
and
. By flow cytometry we confirmed that the expression profile of adhesion molecules (CD49e
, CD29
, CD81
, CD62L
, CD11c
) and transcription factors (Eomes
, Tbet
) of ltNK cells differed from their circulating counterparts. LtNK cells were characterized by enhanced expression of inhibitory receptors TIGIT and CD96 and low expression of DNAM1 and cytolytic molecules (
). Their proliferative capacity was reduced compared to the circulating NK cells. By performing gene set enrichment analysis, we identified DUSP6 and EGR2 as potential regulators of the ltNK cell transcriptome. Remarkably, comparison of the ltNK cell transcriptome to the published human spleen-resident memory CD8
T (Trm) cell transcriptome revealed an overlapping gene signature. Moreover, the phenotypic profile of ltNK cells resembled that of CD8
Trm cells in BM. Together, we provide transcriptional and phenotypic data that clearly distinguish ltNK cells from both the CD56
and CD56
NK cells and substantiate the view that ltNK cells are tissue-resident cells, which are functionally restrained in killing and have low proliferative activity.
RNA sequencing is an increasingly popular technology for genome-wide analysis of transcript sequence and abundance. However, understanding of the sources of technical and interlaboratory variation is ...still limited. To address this, the GEUVADIS consortium sequenced mRNAs and small RNAs of lymphoblastoid cell lines of 465 individuals in seven sequencing centers, with a large number of replicates. The variation between laboratories appeared to be considerably smaller than the already limited biological variation. Laboratory effects were mainly seen in differences in insert size and GC content and could be adequately corrected for. In small-RNA sequencing, the microRNA (miRNA) content differed widely between samples owing to competitive sequencing of rRNA fragments. This did not affect relative quantification of miRNAs. We conclude that distributing RNA sequencing among different laboratories is feasible, given proper standardization and randomization procedures. We provide a set of quality measures and guidelines for assessing technical biases in RNA-seq data.
The release and uptake of nano-sized extracellular vesicles (EV) is a highly conserved means of intercellular communication. The molecular composition of EV, and thereby their signaling function to ...target cells, is regulated by cellular activation and differentiation stimuli. EV are regarded as snapshots of cells and are, therefore, in the limelight as biomarkers for disease. Although research on EV-associated RNA has predominantly focused on microRNAs, the transcriptome of EV consists of multiple classes of small non-coding RNAs with potential gene-regulatory functions. It is not known whether environmental cues imposed on cells induce specific changes in a broad range of EV-associated RNA classes. Here, we investigated whether immune-activating or -suppressing stimuli imposed on primary dendritic cells affected the release of various small non-coding RNAs via EV. The small RNA transcriptomes of highly pure EV populations free from ribonucleoprotein particles were analyzed by RNA sequencing and RT-qPCR. Immune stimulus-specific changes were found in the miRNA, snoRNA, and Y-RNA content of EV from dendritic cells, whereas tRNA and snRNA levels were much less affected. Only part of the changes in EV-RNA content reflected changes in cellular RNA, which urges caution in interpreting EV as snapshots of cells. By comprehensive analysis of RNA obtained from highly purified EV, we demonstrate that multiple RNA classes contribute to genetic messages conveyed via EV. The identification of multiple RNA classes that display cell stimulation-dependent association with EV is the prelude to unraveling the function and biomarker potential of these EV-RNAs.