Jiannis Ragoussis

The small-nucleolar RNAs commonly used for microRNA normalisation correlate with tumour pathology and prognosis.

Background:To investigate small-nucleolar RNAs (snoRNAs) as reference genes when measuring miRNA expression in tumour samples, given emerging evidence for their role in cancer.Methods:Four snoRNAs, commonly used for normalisation, RNU44, RNU48, RNU43 and RNU6B, and miRNA known to be associated with pathological factors, were measured by real-time polymerase chain reaction in two patient series: 219 breast cancer and 46 head and neck squamous cell carcinoma (HNSCC). SnoRNA and miRNA were then correlated with clinicopathological features and prognosis.Results:Small-nucleolar RNA expression was as variable as miRNA expression (miR-21, miR-210, miR-10b). Normalising miRNA PCR expression data to these recommended snoRNAs introduced bias in associations between miRNA and pathology or outcome. Low snoRNA expression correlated with markers of aggressive pathology. Low levels of RNU44 were associated with a poor prognosis. RNU44 is an intronic gene in a cluster of highly conserved snoRNAs in the growth arrest specific 5 (GAS5) transcript, which is normally upregulated to arrest cell growth under stress. Low-tumour GAS5 expression was associated with a poor prognosis. RNU48 and RNU43 were also identified as intronic snoRNAs within genes that are dysregulated in cancer.Conclusion:Small-nucleolar RNAs are important in cancer prognosis, and their use as reference genes can introduce bias when determining miRNA expression.

Tumor hypoxia induces nuclear paraspeckle formation through HIF-2α dependent transcriptional activation of NEAT1 leading to cancer cell survival

Activation of cellular transcriptional responses, mediated by hypoxia-inducible factor (HIF), is common in many types of cancer, and generally confers a poor prognosis. Known to induce many hundreds of protein-coding genes, HIF has also recently been shown to be a key regulator of the non-coding transcriptional response. Here, we show that NEAT1 long non-coding RNA (lncRNA) is a direct transcriptional target of HIF in many breast cancer cell lines and in solid tumors. Unlike previously described lncRNAs, NEAT1 is regulated principally by HIF-2 rather than by HIF-1. NEAT1 is a nuclear lncRNA that is an essential structural component of paraspeckles and the hypoxic induction of NEAT1 induces paraspeckle formation in a manner that is dependent upon both NEAT1 and on HIF-2. Paraspeckles are multifunction nuclear structures that sequester transcriptionally active proteins as well as RNA transcripts that have been subjected to adenosine-to-inosine (A-to-I) editing. We show that the nuclear retention of one such transcript, F11R (also known as junctional adhesion molecule 1, JAM1), in hypoxia is dependent upon the hypoxic increase in NEAT1, thereby conferring a novel mechanism of HIF-dependent gene regulation. Induction of NEAT1 in hypoxia also leads to accelerated cellular proliferation, improved clonogenic survival and reduced apoptosis, all of which are hallmarks of increased tumorigenesis. Furthermore, in patients with breast cancer, high tumor NEAT1 expression correlates with poor survival. Taken together, these results indicate a new role for HIF transcriptional pathways in the regulation of nuclear structure and that this contributes to the pro-tumorigenic hypoxia-phenotype in breast cancer.

miR-210 is a target of hypoxia-inducible factors 1 and 2 in renal cancer, regulates ISCU and correlates with good prognosis.

Background:Clear cell renal cancer frequently harbours von Hippel–Lindau (VHL) gene mutations, leading to stabilisation of the hypoxia-inducible factors (HIFs) and expression of their target genes. We investigated HIF-1 and HIF-2 in the regulation of microRNA-210 (miR-210), and its clinical relevance in renal tumours.Methods:RCC4 and 786-O renal cancer cell lines transfected with either an empty vector or functional VHL and incubated in normoxia or hypoxia were examined for miR-210 expression. Hypoxia-inducible factor siRNAs were used to examine their regulation of miR-210. Seventy-one clear cell renal tumours were sequenced for VHL mutations. Expression of miR-210, VHL, CA9, ISCU and Ki-67 were determined by immunohistochemistry and qRT–PCR.Results:In addition to HIF-1 regulating miR-210 in renal cancer, HIF-2 can regulate this microRNA in the absence of HIF-1. MicroRNA-210 is upregulated in renal cancer compared with normal renal cortex tissue. MicroRNA-210 correlates negatively with its gene target ISCU at the protein and mRNA level. MicroRNA-210 correlated with positive outcome variables and negatively with Ki-67.Conclusion:We provide further evidence of miR-210 activity in vivo, and show that high miR-210 expression is associated with better clinico-pathological prognostic factors.

Extensive regulation of the non‐coding transcriptome by hypoxia: role of HIF in releasing paused RNApol2

Hypoxia is central to both ischaemic and neoplastic diseases. However, the non‐coding transcriptional response to hypoxia is largely uncharacterized. We undertook integrated genomic analyses of both non‐coding and coding transcripts using massively parallel sequencing and interfaced this data with pan‐genomic analyses of hypoxia‐inducible factor (HIF) and RNApol2 binding in hypoxic cells. These analyses revealed that all classes of RNA are profoundly regulated by hypoxia and implicated HIF as a major direct regulator of both the non‐coding and coding transcriptome, acting predominantly through release of pre‐bound promoter‐paused RNApol2. These findings indicate that the transcriptional response to hypoxia is substantially more extensive than previously considered.

Biomarker discovery: quantification of microRNAs and other small non-coding RNAs using next generation sequencing

BackgroundSmall ncRNAs (sncRNAs) offer great hope as biomarkers of disease and response to treatment. This has been highlighted in the context of several medical conditions such as cancer, liver disease, cardiovascular disease, and central nervous system disorders, among many others. Here we assessed several steps involved in the development of an ncRNA biomarker discovery pipeline, ranging from sample preparation to bioinformatic processing of small RNA sequencing data.MethodsA total of 45 biological samples were included in the present study. All libraries were prepared using the Illumina TruSeq Small RNA protocol and sequenced using the HiSeq2500 or MiSeq Illumina sequencers. Small RNA sequencing data was validated using qRT-PCR. At each stage, we evaluated the pros and cons of different techniques that may be suitable for different experimental designs. Evaluation methods included quality of data output in relation to hands-on laboratory time, cost, and efficiency of processing.ResultsOur results show that good quality sequencing libraries can be prepared from small amounts of total RNA and that varying degradation levels in the samples do not have a significant effect on the overall quantification of sncRNAs via NGS. In addition, we describe the strengths and limitations of three commercially available library preparation methods: (1) Novex TBE PAGE gel; (2) Pippin Prep automated gel system; and (3) AMPure XP beads. We describe our bioinformatics pipeline, provide recommendations for sequencing coverage, and describe in detail the expression and distribution of all sncRNAs in four human tissues: whole-blood, brain, heart and liver.ConclusionsUltimately this study provides tools and outcome metrics that will aid researchers and clinicians in choosing an appropriate and effective high-throughput sequencing quantification method for various study designs, and overall generating valuable information that can contribute to our understanding of small ncRNAs as potential biomarkers and mediators of biological functions and disease.

An oncogenic role of eIF3e/INT6 in human breast cancer

Altered expression of the eukaryotic translation initiation factor 3 (eIF3) subunit eIF3e/INT6 has been described in various types of human cancer, but the nature of its involvement in tumorigenesis is not yet clear. Using immunohistochemical analysis of 81 primary breast cancers, we found that high tumor grade correlated significantly with elevated cytoplasmic eIF3e level in epithelial tumor cells. Analysis of protein synthesis after siRNA-mediated knockdown in breast cancer cell lines indicated that eIF3e is not required for bulk translation. Microarray analysis of total and polysomal RNAs nonetheless identified distinct sets of mRNAs regulated either positively or negatively by eIF3e; functional classification of these revealed a marked enrichment of genes involved in cell proliferation, invasion and apoptosis. Validated mRNA targets regulated positively at the translational level by eIF3e included urokinase-type plasminogen activator and apoptotic regulator BCL-XL, whereas synthesis of proteins including the mitotic checkpoint component MAD2L1 was negatively regulated. Finally, eIF3e-depleted breast carcinoma cells showed reduced in vitro invasion and proliferation. Taken together, our study data suggest that eIF3e has a positive role in breast cancer progression. It regulates the translation, and in some cases abundance, of mRNAs involved in key aspects of cancer cell biology.

The role of hypoxia regulated microRNAs in cancer.

The molecular response of cancer cells to hypoxia is the focus of intense research. In the last decade, research into microRNAs (miRNAs), small RNAs which have a role in regulation of mRNA and translation, has grown exponentially. miR-210 has emerged as the predominant miRNA regulated by hypoxia. Elucidation of its targets points to a variety of roles for this, and other hypoxia-regulated miRNAs (HRMs), in tumour growth and survival. miR-210 expression correlates with poor survival in cancer patients, and shows promise for future use as a tumour marker or therapeutic agent. The role of miR-210 and other HRMs in cancer biology is the subject of this review.

Benchmarking of the Oxford Nanopore MinION sequencing for quantitative and qualitative assessment of cDNA populations

To assess the performance of the Oxford Nanopore Technologies MinION sequencing platform, cDNAs from the External RNA Controls Consortium (ERCC) RNA Spike-In mix were sequenced. This mix mimics mammalian mRNA species and consists of 92 polyadenylated transcripts with known concentration. cDNA libraries were generated using a template switching protocol to facilitate the direct comparison between different sequencing platforms. The MinION performance was assessed for its ability to sequence the cDNAs directly with good accuracy in terms of abundance and full length. The abundance of the ERCC cDNA molecules sequenced by MinION agreed with their expected concentration. No length or GC content bias was observed. The majority of cDNAs were sequenced as full length. Additionally, a complex cDNA population derived from a human HEK-293 cell line was sequenced on an Illumina HiSeq 2500, PacBio RS II and ONT MinION platforms. We observed that there was a good agreement in the measured cDNA abundance between PacBio RS II and ONT MinION (rpearsonu2009=u20090.82, isoforms with length more than 700bp) and between Illumina HiSeq 2500 and ONT MinION (rpearsonu2009=u20090.75). This indicates that the ONT MinION can sequence quantitatively both long and short full length cDNA molecules.

High-sensitivity sequencing reveals multi-organ somatic mosaicism causing DICER1 syndrome

Background Somatic mosaicism is being increasingly recognised as an important cause of non-Mendelian presentations of hereditary syndromes. A previous whole-exome sequencing study using DNA derived from peripheral blood identified mosaic mutations in DICER1 in two children with overgrowth and developmental delay as well as more typical phenotypes of germline DICER1 mutation. However, very-low-frequency mosaicism is difficult to detect, and thus, causal mutations can go unnoticed. Highly sensitive, cost-effective approaches are needed to molecularly diagnose these persons. We studied four children with multiple primary tumours known to be associated with the DICER1 syndrome, but in whom germline DICER1 mutations were not detected by conventional mutation detection techniques. Methods and results We observed the same missense mutation within the DICER1 RNase IIIb domain in multiple tumours from different sites in each patient, raising suspicion of somatic mosaicism. We implemented three different targeted-capture technologies, including the novel HaloPlexHS (Agilent Technologies), followed by deep sequencing, and confirmed that the identified mutations are mosaic in origin in three patients, detectable in 0.24–31% of sequencing reads in constitutional DNA. The mosaic origin of patient 4s mutation remains to be unequivocally established. We also discovered likely pathogenic second somatic mutations or loss of heterozygosity (LOH) in tumours from all four patients. Conclusions Mosaic DICER1 mutations are an important cause of the DICER1 syndrome in patients with severe phenotypes and often appear to be accompanied by second somatic truncating mutations or LOH in the associated tumours. Furthermore, the molecular barcode-containing HaloPlexHS provides the sensitivity required for detection of such low-level mosaic mutations and could have general applicability.

Olive fly transcriptomics analysis implicates energy metabolism genes in spinosad resistance

BackgroundThe olive fly, Bactrocera oleae, is the most devastating pest of cultivated olives. Its control has been traditionally based on insecticides, mainly organophosphates and pyrethroids. In recent years, the naturalyte spinosad is used against the olive fly. As with other insecticides, spinosad is subject to selection pressures that have led to resistance development. Mutations in the α6 subunit of the nicotinic acetylcholine receptor (nAChR) have been implicated in spinosad resistance in several species (e.g., Drosophila melanogaster) but excluded in others (e.g., Musca domestica). Yet, additional mechanisms involving enhanced metabolism of detoxification enzymes (such as P450 monooxygenases or mixed function oxidases) have also been reported. In order to clarify the spinosad resistance mechanisms in the olive fly, we searched for mutations in the α6-subunit of the nAChR and for up-regulated genes in the entire transcriptome of spinosad resistant olive flies.ResultsThe olive fly α6-subunit of the nAChR was cloned from the laboratory sensitive strain and a spinosad selected resistant line. The differences reflected silent nucleotide substitutions or conserved amino acid changes. Additionally, whole transcriptome analysis was performed in the two strains in order to reveal any underlying resistance mechanisms. Comparison of over 13,000 genes showed that in spinosad resistant flies nine genes were significantly over-expressed, whereas ~40 were under-expressed. Further functional analyses of the nine over-expressed and eleven under-expressed loci were performed. Four of these loci (Yolk protein 2, ATP Synthase FO subunit 6, Low affinity cationic amino acid transporter 2 and Serine protease 6) showed consistently higher expression both in the spinosad resistant strain and in wild flies from a resistant California population. On the other side, two storage protein genes (HexL1 and Lsp1) and two heat-shock protein genes (Hsp70 and Hsp23) were unfailingly under-expressed in resistant flies.ConclusionThe observed nucleotide differences in the nAChR-α6 subunit between the sensitive and spinosad resistant olive fly strains did not advocate for the involvement of receptor mutations in spinosad resistance. Instead, the transcriptome comparison between the two strains indicated that several immune system loci as well as elevated energy requirements of the resistant flies might be necessary to lever the detoxification process.