Enhanced Performance of DNA Methylation Markers by Simultaneous Measurement of Sense and Antisense DNA Strands after Cytosine Conversion.

BACKGROUND: Most existing DNA methylation-based methods for detection of circulating tumor DNA (ctDNA) are based on conversion of unmethylated cytosines to uracil. After conversion, the 2 DNA strands are no longer complementary; therefore, targeting only 1 DNA strand merely utilizes half of the available input DNA. We investigated whether the sensitivity of methylation-based ctDNA detection strategies could be increased by targeting both DNA strands after bisulfite conversion. METHODS: Dual-strand digital PCR assays were designed for the 3 colorectal cancer (CRC)-specific methylation markers KCNQ5, C9orf50, and CLIP4 and compared with previously reported single-strand assays. Performance was tested in tumor and leukocyte DNA, and the ability to detect ctDNA was investigated in plasma from 43 patients with CRC stages I to IV and 42 colonoscopy-confirmed healthy controls. RESULTS: Dual-strand assays quantified close to 100% of methylated control DNA input, whereas single-strand assays quantified approximately 50%. Furthermore, dual-strand assays showed a 2-fold increase in the number of methylated DNA copies detected when applied to DNA purified from tumor tissue and plasma from CRC patients. When the results of the 3 DNA methylation markers were combined into a ctDNA detection test and applied to plasma, the dual-strand assay format detected 86% of the cancers compared with 74% for the single-strand assay format. The specificity was 100% for both the dual- and single-strand test formats. CONCLUSION: Dual-strand assays enabled more sensitive detection of methylated ctDNA than single-strand assays.

Putative cis-regulatory drivers in colorectal cancer.

The cis-regulatory effects responsible for cancer development have not been as extensively studied as the perturbations of the protein coding genome in tumorigenesis. To better characterize colorectal cancer (CRC) development we conducted an RNA-sequencing experiment of 103 matched tumour and normal colon mucosa samples from Danish CRC patients, 90 of which were germline-genotyped. By investigating allele-specific expression (ASE) we show that the germline genotypes remain important determinants of allelic gene expression in tumours. Using the changes in ASE in matched pairs of samples we discover 71 genes with excess of somatic cis-regulatory effects in CRC, suggesting a cancer driver role. We correlate genotypes and gene expression to identify expression quantitative trait loci (eQTLs) and find 1,693 and 948 eQTLs in normal samples and tumours, respectively. We estimate that 36% of the tumour eQTLs are exclusive to CRC and show that this specificity is partially driven by increased expression of specific transcription factors and changes in methylation patterns. We show that tumour-specific eQTLs are more enriched for low CRC genome-wide association study (GWAS) P values than shared eQTLs, which suggests that some of the GWAS variants are tumour specific regulatory variants. Importantly, tumour-specific eQTL genes also accumulate more somatic mutations when compared to the shared eQTL genes, raising the possibility that they constitute germline-derived cancer regulatory drivers. Collectively the integration of genome and the transcriptome reveals a substantial number of putative somatic and germline cis-regulatory cancer changes that may have a role in tumorigenesis.

Natural RNA circles function as efficient microRNA sponges.

MicroRNAs (miRNAs) are important post-transcriptional regulators of gene expression that act by direct base pairing to target sites within untranslated regions of messenger RNAs. Recently, miRNA activity has been shown to be affected by the presence of miRNA sponge transcripts, the so-called competing endogenous RNA in humans and target mimicry in plants. We previously identified a highly expressed circular RNA (circRNA) in human and mouse brain. Here we show that this circRNA acts as a miR-7 sponge; we term this circular transcript ciRS-7 (circular RNA sponge for miR-7). ciRS-7 contains more than 70 selectively conserved miRNA target sites, and it is highly and widely associated with Argonaute (AGO) proteins in a miR-7-dependent manner. Although the circRNA is completely resistant to miRNA-mediated target destabilization, it strongly suppresses miR-7 activity, resulting in increased levels of miR-7 targets. In the mouse brain, we observe overlapping co-expression of ciRS-7 and miR-7, particularly in neocortical and hippocampal neurons, suggesting a high degree of endogenous interaction. We further show that the testis-specific circRNA, sex-determining region Y (Sry), serves as a miR-138 sponge, suggesting that miRNA sponge effects achieved by circRNA formation are a general phenomenon. This study serves as the first, to our knowledge, functional analysis of a naturally expressed circRNA.

miRNA-dependent gene silencing involving Ago2-mediated cleavage of a circular antisense RNA.

MicroRNAs (miRNAs) are ∼22 nt non-coding RNAs that typically bind to the 3' UTR of target mRNAs in the cytoplasm, resulting in mRNA destabilization and translational repression. Here, we report that miRNAs can also regulate gene expression by targeting non-coding antisense transcripts in human cells. Specifically, we show that miR-671 directs cleavage of a circular antisense transcript of the Cerebellar Degeneration-Related protein 1 (CDR1) locus in an Ago2-slicer-dependent manner. The resulting downregulation of circular antisense has a concomitant decrease in CDR1 mRNA levels, independently of heterochromatin formation. This study provides the first evidence for non-coding antisense transcripts as functional miRNA targets, and a novel regulatory mechanism involving a positive correlation between mRNA and antisense circular RNA levels.

Exploring the biology of ctDNA release in colorectal cancer.

BACKGROUND: Circulating tumor DNA (ctDNA) has emerged as a promising tool for early cancer detection and minimal residual disease monitoring. However, the biology underlying ctDNA release and its variation across cancer types and histologies remains poorly understood. This study investigated the biology behind ctDNA shedding in colorectal cancer. METHODS: The study included a local cohort of 747 stage I-III colorectal cancer patients. All patients had ctDNA measurement prior to treatment and extensive clinical data. Primary tumor RNA sequencing and whole exome sequencing was performed in 95 and 652 patients respectively. Additionally, the study evaluated 89 non-small cell lung cancer patients from the TRACERx cohort, comprising primary tumor RNA sequencing and ctDNA measurement. RESULTS: We found tumor size and proliferative capacity to be key factors associated with ctDNA shedding in colorectal cancer. Furthermore, we found that the secretory and CMS3 colorectal cancer subtypes exhibited lower ctDNA shedding, while microsatellite instability (MSI) tumors had higher levels of ctDNA. Mutational analysis did not reveal any genes or pathways associated with ctDNA shedding in colorectal cancer. A comparison of transcriptomic profiles across multiple cancer types demonstrated that colorectal cancer and lung squamous cell carcinoma tumors shared a high-proliferative ctDNA shedding phenotype, while lung adenocarcinoma tumors displayed a distinct low-proliferative subgroup. Additionally, proliferation levels correlated with ctDNA detection sensitivity across multiple cancer types. CONCLUSION: These findings suggest that tumor size and proliferative capacity are drivers of ctDNA release in colorectal cancer and provide insights into the biology of ctDNA shedding on a pan-cancer level.

In vivo knockdown of Brachyury results in skeletal defects and urorectal malformations resembling caudal regression syndrome.

The T-box transcription factor BRACHYURY (T) is a key regulator of mesoderm formation during early development. Complete loss of T has been shown to lead to embryonic lethality around E10.0. Here we characterize an inducible miRNA-based in vivo knockdown mouse model of T, termed KD3-T, which exhibits a hypomorphic phenotype. KD3-T embryos display axial skeletal defects caused by apoptosis of paraxial mesoderm, which is accompanied by urorectal malformations resembling the murine uro-recto-caudal syndrome and human caudal regression syndrome phenotypes. We show that there is a reduction of T in the notochord of KD3-T embryos which results in impaired notochord differentiation and its subsequent loss, whereas levels of T in the tailbud are sufficient for axis extension and patterning. Furthermore, the notochord in KD3-T embryos adopts a neural character and loses its ability to act as a signaling center. Since KD3-T animals survive until birth, they are useful for examining later roles for T in the development of urorectal tissues.

Spatio-temporal regulation of ADAR editing during development in porcine neural tissues.

Editing by ADAR enzymes is essential for mammalian life. Still, knowledge of the spatio-temporal editing patterns in mammals is limited. By use of 454 amplicon sequencing we examined the editing status of 12 regionally extracted mRNAs from porcine developing brain encompassing a total of 64 putative ADAR editing sites. In total 24 brain tissues, dissected from up to five regions from embryonic gestation day 23, 42, 60, 80, 100 and 115, were examined for editing.   Generally, editing increased during embryonic development concomitantly with an increase in ADAR2 mRNA level. Notably, the Gria2 (GluR-B) Q/R site, reported to be ~100% edited in previous studies, is only 54% edited at embryonic day 23. Transcripts with multiple editing sites in close proximity to each other exhibit coupled editing and an extraordinary incidence of long-range coupling of editing events more than 32 kb apart is observed for the kainate glutamate receptor 2 transcript, Grik2. Our study reveals complex spatio-temporal ADAR editing patterns of coordinated editing events that may play important roles in the development of the mammalian brain.

A screen of chemical modifications identifies position-specific modification by UNA to most potently reduce siRNA off-target effects.

Small interfering RNAs (siRNAs) are now established as the preferred tool to inhibit gene function in mammalian cells yet trigger unintended gene silencing due to their inherent miRNA-like behavior. Such off-target effects are primarily mediated by the sequence-specific interaction between the siRNA seed regions (position 2-8 of either siRNA strand counting from the 5'-end) and complementary sequences in the 3'UTR of (off-) targets. It was previously shown that chemical modification of siRNAs can reduce off-targeting but only very few modifications have been tested leaving more to be identified. Here we developed a luciferase reporter-based assay suitable to monitor siRNA off-targeting in a high throughput manner using stable cell lines. We investigated the impact of chemically modifying single nucleotide positions within the siRNA seed on siRNA function and off-targeting using 10 different types of chemical modifications, three different target sequences and three siRNA concentrations. We found several differently modified siRNAs to exercise reduced off-targeting yet incorporation of the strongly destabilizing unlocked nucleic acid (UNA) modification into position 7 of the siRNA most potently reduced off-targeting for all tested sequences. Notably, such position-specific destabilization of siRNA-target interactions did not significantly reduce siRNA potency and is therefore well suited for future siRNA designs especially for applications in vivo where siRNA concentrations, expectedly, will be low.

Coordinated epigenetic repression of the miR-200 family and miR-205 in invasive bladder cancer.

MicroRNAs (miRNA) are small noncoding RNAs commonly deregulated in cancer. The miR-200 family (miR-200a, -200b, -200c, -141 and -429) and miR-205 are frequently silenced in advanced cancer and have been implicated in epithelial to mesenchymal transition (EMT) and tumor invasion by targeting the transcriptional repressors of E-cadherin, ZEB1 and ZEB2. ZEB1 is also known to repress miR-200c-141 transcription in a negative feedback loop, but otherwise little is known about the transcriptional regulation of the miR-200 family and miR-205. Recently, miR-200 silencing was also reported in cancer stem cells, implying that miR-200 deregulation is a key event in multiple levels of tumor biology. However, what prevents miR-200 expression remains largely unanswered. Here we report concerted transcriptional regulation of the miR-200 and miR-205 loci in bladder tumors and bladder cell lines. Using a combination of miRNA expression arrays, qPCR assays and mass spectrometry DNA methylation analyses, we show that the miR-200 and miR-205 loci are specifically silenced and gain promoter hypermethylation and repressive chromatin marks in muscle invasive bladder tumors and undifferentiated bladder cell lines. Moreover, we report that miR-200c expression is significantly correlated with early stage T1 bladder tumor progression, and propose miR-200 and miR-205 silencing and DNA hypermethylation as possible prognostic markers in bladder cancer. In addition, we observe that the mesoderm transcription factor TWIST1 and miR-200 expression are inversely correlated in bladder tumor samples and cell lines. TWIST1 associates directly with the miR-200 and miR-205 promoters, and may act as a repressor of miR-200 and miR-205 expression.

Enhancing miRNA annotation confidence in miRBase by continuous cross dataset analysis.

The immaculate annotation of all microRNAs (miRNAs) is a prerequisite to study their biological function on a genome-wide scale. However, the original criteria for proper miRNA annotation seem unsuited for the automated analysis of the immense number of small RNA reads available in next generation sequencing (NGS) datasets. Here we analyze the confidence of past miRNA annotation in miRBase by cross-analyzing publicly available NGS datasets using strengthened annotation requirements. Our analysis highlights that a large number of annotated human miRNAs in miRBase seems to require more experimental validation to be confidently annotated. Notably, our dataset analysis also identified almost 300 currently non-annotated miRNA*s and 28 novel miRNAs. These observations hereby greatly increase the confidence of past miRNA annotation in miRBase but also illustrate the usefulness of continuous re-evaluating NGS datasets in the identification of novel miRNAs.

A large-scale chemical modification screen identifies design rules to generate siRNAs with high activity, high stability and low toxicity.

The use of chemically synthesized short interfering RNAs (siRNAs) is currently the method of choice to manipulate gene expression in mammalian cell culture, yet improvements of siRNA design is expectably required for successful application in vivo. Several studies have aimed at improving siRNA performance through the introduction of chemical modifications but a direct comparison of these results is difficult. We have directly compared the effect of 21 types of chemical modifications on siRNA activity and toxicity in a total of 2160 siRNA duplexes. We demonstrate that siRNA activity is primarily enhanced by favouring the incorporation of the intended antisense strand during RNA-induced silencing complex (RISC) loading by modulation of siRNA thermodynamic asymmetry and engineering of siRNA 3'-overhangs. Collectively, our results provide unique insights into the tolerance for chemical modifications and provide a simple guide to successful chemical modification of siRNAs with improved activity, stability and low toxicity.

MethCORR infers gene expression from DNA methylation and allows molecular analysis of ten common cancer types using fresh-frozen and formalin-fixed paraffin-embedded tumor samples.

BACKGROUND: Transcriptional analysis is widely used to study the molecular biology of cancer and hold great biomarker potential for clinical patient stratification. Yet, accurate transcriptional profiling requires RNA of a high quality, which often cannot be retrieved from formalin-fixed, paraffin-embedded (FFPE) tumor tissue that is routinely collected and archived in clinical departments. To overcome this roadblock to clinical testing, we previously developed MethCORR, a method that infers gene expression from DNA methylation data, which is robustly retrieved from FFPE tissue. MethCORR was originally developed for colorectal cancer and with this study, we aim to: (1) extend the MethCORR method to 10 additional cancer types and (2) to illustrate that the inferred gene expression is accurate and clinically informative. RESULTS: Regression models to infer gene expression information from DNA methylation were developed for ten common cancer types using matched RNA sequencing and DNA methylation profiles (HumanMethylation450 BeadChip) from The Cancer Genome Atlas Project. Robust and accurate gene expression profiles were inferred for all cancer types: on average, the expression of 11,000 genes was modeled with good accuracy and an intra-sample correlation of R2 = 0.90 between inferred and measured gene expression was observed. Molecular pathway analysis and transcriptional subtyping were performed for breast, prostate, and lung cancer samples to illustrate the general usability of the inferred gene expression profiles: overall, a high correlation of r = 0.96 (Pearson) in pathway enrichment scores and a 76% correspondence in molecular subtype calls were observed when using measured and inferred gene expression as input. Finally, inferred expression from FFPE tissue correlated better with RNA sequencing data from matched fresh-frozen tissue than did RNA sequencing data from FFPE tissue (P < 0.0001; Wilcoxon rank-sum test). CONCLUSIONS: In all cancers investigated, MethCORR enabled DNA methylation-based transcriptional analysis, thus enabling future analysis of cancer in situations where high-quality DNA, but not RNA, is available. Here, we provide the framework and resources for MethCORR modeling of ten common cancer types, thereby widely expanding the possibilities for transcriptional studies of archival FFPE material.

MethCORR modelling of methylomes from formalin-fixed paraffin-embedded tissue enables characterization and prognostication of colorectal cancer.

Transcriptional characterization and classification has potential to resolve the inter-tumor heterogeneity of colorectal cancer and improve patient management. Yet, robust transcriptional profiling is difficult using formalin-fixed, paraffin-embedded (FFPE) samples, which complicates testing in clinical and archival material. We present MethCORR, an approach that allows uniform molecular characterization and classification of fresh-frozen and FFPE samples. MethCORR identifies genome-wide correlations between RNA expression and DNA methylation in fresh-frozen samples. This information is used to infer gene expression information in FFPE samples from their methylation profiles. MethCORR is here applied to methylation profiles from 877 fresh-frozen/FFPE samples and comparative analysis identifies the same two subtypes in four independent cohorts. Furthermore, subtype-specific prognostic biomarkers that better predicts relapse-free survival (HR = 2.66, 95%CI [1.67-4.22], P value < 0.001 (log-rank test)) than UICC tumor, node, metastasis (TNM) staging and microsatellite instability status are identified and validated using DNA methylation-specific PCR. The MethCORR approach is general, and may be similarly successful for other cancer types.

Publisher Correction: MethCORR modelling of methylomes from formalin-fixed paraffin-embedded tissue enables characterization and prognostication of colorectal cancer.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Improved silencing properties using small internally segmented interfering RNAs.

RNA interference is mediated by small interfering RNAs (siRNAs) that upon incorporation into the RNA-induced silencing complex (RISC) can target complementary mRNA for degradation. Standard siRNA design usually feature a 19-27 base pair contiguous double-stranded region that is believed to be important for RISC incorporation. Here, we describe a novel siRNA design composed of an intact antisense strand complemented with two shorter 10-12 nt sense strands. This three-stranded construct, termed small internally segmented interfering RNA (sisiRNA), is highly functional demonstrating that an intact sense strand is not a prerequisite for RNA interference. Moreover, when using the sisiRNA design only the antisense strand is functional in activated RISC thereby completely eliminating unintended mRNA targeting by the sense strand. Interestingly, the sisiRNA design supports the function of chemically modified antisense strands, which are non-functional within the context of standard siRNA designs. This suggests that the sisiRNA design has a clear potential of improving the pharmacokinetic properties of siRNA in vivo.

Synthesis of 2'-O-modified adenosine building blocks and application for RNA interference.

RNA-interference has been recognized as a powerful tool to control gene function and has been used for gene silencing by knocking down mRNA. Chemically modified RNAs, especially 2'-O-modification, successfully improved their physicochemical and pharmaceutical properties such as stability, nuclease resistance and delivery. Here, we report the synthesis of adenosine building blocks with different 2'-tethered modifications like aminoethyl and guanidinoethyl and show that they are compatible with RNAi function. They enhance the half life of the siRNA in serum suggesting that these modifications can enhance the pharmacokinetic properties and knock down activity of siRNAs in vivo.

Characterization of genetic intratumor heterogeneity in colorectal cancer and matching patient-derived spheroid cultures.

Patient-derived in vitro cultures of colorectal cancer (CRC) may help guide treatment strategies prior to patient treatment. However, most previous studies have been performed on a single biopsy per tumor. The purpose of this study was to analyze multiple spatially distinct biopsies from CRCs and see how well intratumor heterogeneity (ITH) was recapitulated in matching patient-derived spheroids. Three to five biopsies were collected from six CRC tumors. Each biopsy was split in two; one half was used for spheroid culturing, while the other half was used for DNA and RNA purification. For two patients, lymph node metastases were analyzed. Somatic mutations were called from whole exome sequencing data. Each tumor contained mutations shared across all biopsies and spheroids, including major CRC drivers such as APC, KRAS, and TP53. At the same time, all tumors exhibited ITH on both mutation and copy number level. The concordance between biopsies and spheroids ranged between 40 and 70% for coding mutations. For three patients, the biopsy and spheroid from matching areas clustered together, meaning that the spheroid resembled the area of origin more than the other areas. However, all biopsies and spheroids contained private mutations. Therefore, multiple cultures from spatially distinct sites of the tumor increase the insight into the genetic profile of the entire tumor. Molecular subtypes were called from RNA sequencing data. When based on transcripts from both cancer and noncancerous cells, the subtypes were largely independent of sampling site. In contrast, subtyping based on cancer cell transcripts alone was dependent on sample site and genetic ITH. In conclusion, all examined CRC tumors showed genetic ITH. Spheroid cultures partly reflected this ITH, and having multiple cultures from distinct tumor sites improved the representation of the genetic tumor subclones. This should be taken into account when establishing patient-derived models for drug screening.

An Optimized Shotgun Strategy for the Rapid Generation of Comprehensive Human Proteomes.

This study investigates the challenge of comprehensively cataloging the complete human proteome from a single-cell type using mass spectrometry (MS)-based shotgun proteomics. We modify a classical two-dimensional high-resolution reversed-phase peptide fractionation scheme and optimize a protocol that provides sufficient peak capacity to saturate the sequencing speed of modern MS instruments. This strategy enables the deepest proteome of a human single-cell type to date, with the HeLa proteome sequenced to a depth of ∼584,000 unique peptide sequences and ∼14,200 protein isoforms (∼12,200 protein-coding genes). This depth is comparable with next-generation RNA sequencing and enables the identification of post-translational modifications, including ∼7,000 N-acetylation sites and ∼10,000 phosphorylation sites, without the need for enrichment. We further demonstrate the general applicability and clinical potential of this proteomics strategy by comprehensively quantifying global proteome expression in several different human cancer cell lines and patient tissue samples.

The non-coding variant rs1800734 enhances DCLK3 expression through long-range interaction and promotes colorectal cancer progression.

Genome-wide association studies have identified a great number of non-coding risk variants for colorectal cancer (CRC). To date, the majority of these variants have not been functionally studied. Identification of allele-specific transcription factor (TF) binding is of great importance to understand regulatory consequences of such variants. A recently developed proteome-wide analysis of disease-associated SNPs (PWAS) enables identification of TF-DNA interactions in an unbiased manner. Here we perform a large-scale PWAS study to comprehensively characterize TF-binding landscape that is associated with CRC, which identifies 731 allele-specific TF binding at 116 CRC risk loci. This screen identifies the A-allele of rs1800734 within the promoter region of MLH1 as perturbing the binding of TFAP4 and consequently increasing DCLK3 expression through a long-range interaction, which promotes cancer malignancy through enhancing expression of the genes related to epithelial-to-mesenchymal transition.

SNHG16 is regulated by the Wnt pathway in colorectal cancer and affects genes involved in lipid metabolism.

It is well established that lncRNAs are aberrantly expressed in cancer where they have been shown to act as oncogenes or tumor suppressors. RNA profiling of 314 colorectal adenomas/adenocarcinomas and 292 adjacent normal colon mucosa samples using RNA-sequencing demonstrated that the snoRNA host gene 16 (SNHG16) is significantly up-regulated in adenomas and all stages of CRC. SNHG16 expression was positively correlated to the expression of Wnt-regulated transcription factors, including ASCL2, ETS2, and c-Myc. In vitro abrogation of Wnt signaling in CRC cells reduced the expression of SNHG16 indicating that SNHG16 is regulated by the Wnt pathway. Silencing of SNHG16 resulted in reduced viability, increased apoptotic cell death and impaired cell migration. The SNHG16 silencing particularly affected expression of genes involved in lipid metabolism. A connection between SNHG16 and genes involved in lipid metabolism was also observed in clinical tumors. Argonaute CrossLinking and ImmunoPrecipitation (AGO-CLIP) demonstrated that SNHG16 heavily binds AGO and has 27 AGO/miRNA target sites along its length, indicating that SNHG16 may act as a competing endogenous RNA (ceRNA) "sponging" miRNAs off their cognate targets. Most interestingly, half of the miRNA families with high confidence targets on SNHG16 also target the 3'UTR of Stearoyl-CoA Desaturase (SCD). SCD is involved in lipid metabolism and is down-regulated upon SNHG16 silencing. In conclusion, up-regulation of SNHG16 is a frequent event in CRC, likely caused by deregulated Wnt signaling. In vitro analyses demonstrate that SNHG16 may play an oncogenic role in CRC and that it affects genes involved in lipid metabolism, possible through ceRNA related mechanisms.