Jian Gu

MicroRNAs and their isomiRs function cooperatively to target common biological pathways

BackgroundVariants of microRNAs (miRNAs), called isomiRs, are commonly reported in deep-sequencing studies; however, the functional significance of these variants remains controversial. Observational studies show that isomiR patterns are non-random, hinting that these molecules could be regulated and therefore functional, although no conclusive biological role has been demonstrated for these molecules.ResultsTo assess the biological relevance of isomiRs, we have performed ultra-deep miRNA-seq on ten adult human tissues, and created an analysis pipeline called miRNA-MATE to align, annotate, and analyze miRNAs and their isomiRs. We find that isomiRs share sequence and expression characteristics with canonical miRNAs, and are generally strongly correlated with canonical miRNA expression. A large proportion of isomiRs potentially derive from AGO2 cleavage independent of Dicer. We isolated polyribosome-associated mRNA, captured the mRNA-bound miRNAs, and found that isomiRs and canonical miRNAs are equally associated with translational machinery. Finally, we transfected cells with biotinylated RNA duplexes encoding isomiRs or their canonical counterparts and directly assayed their mRNA targets. These studies allow us to experimentally determine genome-wide mRNA targets, and these experiments showed substantial overlap in functional mRNA networks suppressed by both canonical miRNAs and their isomiRs.ConclusionsTogether, these results find isomiRs to be biologically relevant and functionally cooperative partners of canonical miRNAs that act coordinately to target pathways of functionally related genes. This work exposes the complexity of the miRNA-transcriptome, and helps explain a major miRNA paradox: how specific regulation of biological processes can occur when the specificity of miRNA targeting is mediated by only 6 to 11 nucleotides.

Tumor Transcriptome Sequencing Reveals Allelic Expression Imbalances Associated with Copy Number Alterations

Due to growing throughput and shrinking cost, massively parallel sequencing is rapidly becoming an attractive alternative to microarrays for the genome-wide study of gene expression and copy number alterations in primary tumors. The sequencing of transcripts (RNA-Seq) should offer several advantages over microarray-based methods, including the ability to detect somatic mutations and accurately measure allele-specific expression. To investigate these advantages we have applied a novel, strand-specific RNA-Seq method to tumors and matched normal tissue from three patients with oral squamous cell carcinomas. Additionally, to better understand the genomic determinants of the gene expression changes observed, we have sequenced the tumor and normal genomes of one of these patients. We demonstrate here that our RNA-Seq method accurately measures allelic imbalance and that measurement on the genome-wide scale yields novel insights into cancer etiology. As expected, the set of genes differentially expressed in the tumors is enriched for cell adhesion and differentiation functions, but, unexpectedly, the set of allelically imbalanced genes is also enriched for these same cancer-related functions. By comparing the transcriptomic perturbations observed in one patient to his underlying normal and tumor genomes, we find that allelic imbalance in the tumor is associated with copy number mutations and that copy number mutations are, in turn, strongly associated with changes in transcript abundance. These results support a model in which allele-specific deletions and duplications drive allele-specific changes in gene expression in the developing tumor.

The Complete Exosome Workflow Solution: From Isolation to Characterization of RNA Cargo

Exosomes are small (30–150 nm) vesicles containing unique RNA and protein cargo, secreted by all cell types in culture. They are also found in abundance in body fluids including blood, saliva, and urine. At the moment, the mechanism of exosome formation, the makeup of the cargo, biological pathways, and resulting functions are incompletely understood. One of their most intriguing roles is intercellular communication—exosomes function as the messengers, delivering various effector or signaling macromolecules between specific cells. There is an exponentially growing need to dissect structure and the function of exosomes and utilize them for development of minimally invasive diagnostics and therapeutics. Critical to further our understanding of exosomes is the development of reagents, tools, and protocols for their isolation, characterization, and analysis of their RNA and protein contents. Here we describe a complete exosome workflow solution, starting from fast and efficient extraction of exosomes from cell culture media and serum to isolation of RNA followed by characterization of exosomal RNA content using qRT-PCR and next-generation sequencing techniques. Effectiveness of this workflow is exemplified by analysis of the RNA content of exosomes derived from HeLa cell culture media and human serum, using Ion Torrent PGM as a sequencing platform.

The Human Salivary RNA Transcriptome Revealed by Massively Parallel Sequencing

BACKGROUND Evaluation of the salivary transcriptome is an emerging diagnostic technology with discriminatory power for disease detection. This study explored massively parallel sequencing for providing nucleotide-level sequence information for each RNA in saliva. METHODS Transcriptome profiling with the SOLiD™ system was applied to RNA isolated from unstimulated cell-free saliva (CFS) and whole saliva (WS) from healthy human volunteers. Sequenced reads were aligned to human genome build 18 and the Human Oral Microbiome Database (HOMD). RESULTS Massively parallel sequencing enabled the acquisition of complete sequence information for each nucleotide position of the human salivary transcriptome through alignment to multiple sequence databases. Approximately 20%-25% of the sequenced reads from CFS aligned to the human genome, and approximately 30% of the sequenced reads aligned to the HOMD. We detected the expression of >4000 coding and noncoding genes in CFS and WS. Gene expression at different genome loci showed that the structural integrity of the transcripts for the annotated genes was preserved in saliva. CONCLUSIONS A single measurement provided RNA sequence information of gene transcript abundance for both coding and noncoding RNAs and identified sequences from >400 different microbial species within a single sample. Contrary to previous data suggesting that salivary RNA is highly fragmented, in our study the structural integrity of RNA was preserved. The high degree of sequence alignment to annotated exons and introns for each of the respective reference genomes, with sequence coverage spanning the full length of the genes, provides strong evidence that the salivary transcriptome is a complex RNA network.

RNA-Seq Mapping and Detection of Gene Fusions with a Suffix Array Algorithm

High-throughput RNA sequencing enables quantification of transcripts (both known and novel), exon/exon junctions and fusions of exons from different genes. Discovery of gene fusions–particularly those expressed with low abundance– is a challenge with short- and medium-length sequencing reads. To address this challenge, we implemented an RNA-Seq mapping pipeline within the LifeScope software. We introduced new features including filter and junction mapping, annotation-aided pairing rescue and accurate mapping quality values. We combined this pipeline with a Suffix Array Spliced Read (SASR) aligner to detect chimeric transcripts. Performing paired-end RNA-Seq of the breast cancer cell line MCF-7 using the SOLiD system, we called 40 gene fusions among over 120,000 splicing junctions. We validated 36 of these 40 fusions with TaqMan assays, of which 25 were expressed in MCF-7 but not the Human Brain Reference. An intra-chromosomal gene fusion involving the estrogen receptor alpha gene ESR1, and another involving the RPS6KB1 (Ribosomal protein S6 kinase beta-1) were recurrently expressed in a number of breast tumor cell lines and a clinical tumor sample.