Quantifying full-length circular RNAs in cancer.

Circular RNAs (circRNAs) are abundantly expressed in cancer. Their resistance to exonucleases enables them to have potentially stable interactions with different types of biomolecules. Alternative splicing can create different circRNA isoforms that have different sequences and unequal interaction potentials. The study of circRNA function thus requires knowledge of complete circRNA sequences. Here we describe psirc, a method that can identify full-length circRNA isoforms and quantify their expression levels from RNA sequencing data. We confirm the effectiveness and computational efficiency of psirc using both simulated and actual experimental data. Applying psirc on transcriptome profiles from nasopharyngeal carcinoma and normal nasopharynx samples, we discover and validate circRNA isoforms differentially expressed between the two groups. Compared with the assumed circular isoforms derived from linear transcript annotations, some of the alternatively spliced circular isoforms have 100 times higher expression and contain substantially fewer microRNA response elements, showing the importance of quantifying full-length circRNA isoforms.

CircRNA-miRNA interactions in atherogenesis.

Atherosclerosis is the major cause of coronary artery disease (CAD) which includes unstable angina, myocardial infarction, and heart failure. The onset of atherogenesis, a process of atherosclerotic lesion formation in the intima of arteries, is driven by lipid accumulation, a vicious cycle of reactive oxygen species (ROS)-induced oxidative stress and inflammatory reactions leading to endothelial cell (EC) dysfunction, vascular smooth muscle cell (VSMC) activation, and foam cell formation which further fuel plaque formation and destabilization. In recent years, there is a surge in the number of publications reporting the involvement of circular RNAs (circRNAs) in the pathogenesis of cardiovascular diseases, cancers, and metabolic syndromes. These studies have advanced our understanding on the biological functions of circRNAs. One of the most common mechanism of action of circRNAs reported is the sponging of microRNAs (miRNAs) by binding to the miRNAs response element (MRE), thereby indirectly increases the transcription of their target messenger RNAs (mRNAs). Individual networks of circRNA-miRNA-mRNA associated with atherogenesis have been extensively reported, however, there is a need to connect these findings for a complete overview. This review aims to provide an update on atherogenesis-related circRNAs and analyze the circRNA-miRNA-mRNA interactions in atherogenesis. The atherogenic mechanisms and clinical relevance of each atherogenesis-related circRNA were systematically discussed for better understanding of the knowledge gap in this area.

Inducible RasGEF1B circular RNA is a positive regulator of ICAM-1 in the TLR4/LPS pathway.

Circular RNAs (circRNAs) constitute a large class of RNA species formed by the back-splicing of co-linear exons, often within protein-coding transcripts. Despite much progress in the field, it remains elusive whether the majority of circRNAs are merely aberrant splicing by-products with unknown functions, or their production is spatially and temporally regulated to carry out specific biological functions. To date, the majority of circRNAs have been cataloged in resting cells. Here, we identify an LPS-inducible circRNA: mcircRasGEF1B, which is predominantly localized in cytoplasm, shows cell-type specific expression, and has a human homolog with similar properties, hcircRasGEF1B. We show that knockdown of the expression of mcircRasGEF1B reduces LPS-induced ICAM-1 expression. Additionally, we demonstrate that mcircRasGEF1B regulates the stability of mature ICAM-1 mRNAs. These findings expand the inventory of functionally characterized circRNAs with a novel RNA species that may play a critical role in fine-tuning immune responses and protecting cells against microbial infection.

Epigenetic modulation of the miR-200 family is associated with transition to a breast cancer stem-cell-like state.

The miR-200 family is a key regulator of the epithelial-mesenchymal transition, however, its role in controlling the transition between cancer stem-cell-like and non-stem-cell-like phenotypes is not well understood. We utilized immortalized human mammary epithelial (HMLE) cells to investigate the regulation of the miR-200 family during their conversion to a stem-like phenotype. HMLE cells were found to be capable of spontaneous conversion from a non-stem to a stem-like phenotype and this conversion was accompanied by the loss of miR-200 expression. Stem-like cell fractions isolated from metastatic breast cancers also displayed loss of miR-200 indicating similar molecular changes may occur during breast cancer progression. The phenotypic change observed in HMLE cells was directly controlled by miR-200 because restoration of its expression decreased stem-like properties while promoting a transition to an epithelial phenotype. Investigation of the mechanisms controlling miR-200 expression revealed both DNA methylation and histone modifications were significantly altered in the stem-like and non-stem phenotypes. In particular, in the stem-like phenotype, the miR-200b-200a-429 cluster was silenced primarily through polycomb group-mediated histone modifications whereas the miR-200c-141 cluster was repressed by DNA methylation. These results indicate that the miR-200 family plays a crucial role in the transition between stem-like and non-stem phenotypes and that distinct epigenetic-based mechanisms regulate each miR-200 gene in this process. Therapy targeted against miR-200 family members and epigenetic modifications might therefore be applicable to breast cancer.

Detection of Cytoplasmic and Nuclear Circular RNA via RT-qPCR.

Circular RNA (circRNA) is an intriguing class of non-coding RNA that exists as a continuous closed loop. With the improvements in high throughput sequencing, biochemical analysis, and bioinformatic algorithms, studies on circRNA expression became abundant in recent years. However, functional studies of circRNA are still limited. Subcellular localization of circRNA may provide some clues in elucidating its biological functions by performing subcellular fractionation assay. Notably, circRNAs that are predominantly found in the cytoplasm are more likely to be involved in post-transcriptional gene regulation, e.g., acting as micoRNA sponge, whereas nuclear-retained circRNAs are predicted to play a role in transcriptional regulation. Subcellular fractionation could help researchers to narrow down and prioritize downstream experiments. The majority of the currently available protocols describe the steps for subcellular fractionation followed by western blot analysis for protein molecules. Here, we present a protocol for the subcellular fractionation of cells to detect circRNA via RT-qPCR with divergent primers. Moreover, detailed steps for the generation of specific circRNAs-enriched cDNA included in this protocol will enhance the amplification and detection of low-abundance circRNAs. This will be useful for researchers studying low-abundance circRNAs. Key features This protocol builds upon the method developed by Gagnon et al. (2014) and extends its application to circRNA study. Protocol for amplification of low levels of circRNA expression. Analysis takes into consideration the ratio of cytoplasmic RNA concentration to nuclear RNA concentration.

In Silico Identification and Characterization of circRNAs During Host-Pathogen Interactions.

Circular RNAs (circRNAs) are a class of non-coding RNAs that are formed via back-splicing. These circRNAs are predominantly studied for their roles as regulators of various biological processes. Notably, emerging evidence demonstrates that host circRNAs can be differentially expressed (DE) upon infection with pathogens (e.g., influenza and coronaviruses), suggesting a role for circRNAs in regulating host innate immune responses. However, investigations on the role of circRNAs during pathogenic infections are limited by the knowledge and skills required to carry out the necessary bioinformatic analysis to identify DE circRNAs from RNA sequencing (RNA-seq) data. Bioinformatics prediction and identification of circRNAs is crucial before any verification, and functional studies using costly and time-consuming wet-lab techniques. To solve this issue, a step-by-step protocol of in silico prediction and characterization of circRNAs using RNA-seq data is provided in this manuscript. The protocol can be divided into four steps: 1) Prediction and quantification of DE circRNAs via the CIRIquant pipeline; 2) Annotation via circBase and characterization of DE circRNAs; 3) CircRNA-miRNA interaction prediction through Circr pipeline; 4) functional enrichment analysis of circRNA parental genes using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). This pipeline will be useful in driving future in vitro and in vivo research to further unravel the role of circRNAs in host-pathogen interactions.

Viruses join the circular RNA world.

Circular RNAs (circRNAs) are a recently discovered class of noncoding RNAs found in many species across the eukaryotic kingdom. These intriguing RNA species are formed through a unique mechanism that is known as back splicing in which the 5' and 3' termini are covalently joined. Recent research has revealed that viruses also encode a repertoire of circRNAs. Some of these viral circRNAs are abundantly expressed and are reported to play a role in disease pathogenesis. A growing number of studies also indicate that host circRNAs are involved in immune responses against virus infections with either an antiviral or proviral role. In this review, we briefly introduce circRNA, its biogenesis, and mechanism of action. We go on to summarize the latest research on the expression, regulation, and functions of viral and host-encoded circRNAs during the host-virus interaction, with the aim of highlighting the potential of viral and host circRNAs as a suitable target for diagnostic biomarker development and therapeutic treatment of viral-associated diseases. We conclude by discussing the current limitations in knowledge and significance of elucidating the roles of circRNAs in host-virus interactions, as well as future directions for this emerging field.

Identification and characterization of a novel Epstein-Barr Virus-encoded circular RNA from LMP-2 Gene.

Epstein-Barr virus (EBV) has been recently found to generate novel circular RNAs (circRNAs) through backsplicing. However, comprehensive catalogs of EBV circRNAs in other cell lines and their functional characterization are still lacking. In this study, we have identified a list of putative EBV circRNAs in GM12878, an EBV-transformed lymphoblastoid cell line, with a significant majority encoded from the EBV latent genes. A novel EBV circRNA derived from the exon 5 of LMP-2 gene which exhibited highest prevalence, was further validated using RNase R assay and Sanger sequencing. This circRNA, which we term circLMP-2_e5, can be universally detected in a panel of EBV-positive cell lines modelling different latency programs. It ranges from lower expression in nasopharyngeal carcinoma (NPC) cells to higher expression in B cells, and is localized to both the cytoplasm and the nucleus. We provide evidence that circLMP-2_e5 is expressed concomitantly with its cognate linear LMP-2 RNA upon EBV lytic reactivation, and may be produced as a result of exon skipping, with its circularization possibly occurring without the involvement of cis elements in the short flanking introns. Furthermore, we show that circLMP-2_e5 is not involved in regulating cell proliferation, host innate immune response, its linear parental transcripts, or EBV lytic reactivation. Taken together, our study expands the current repertoire of putative EBV circRNAs, broadens our understanding of the biology of EBV circRNAs, and lays the foundation for further investigation of their function in the EBV life cycle and disease development.

JMJD8 is a novel endoplasmic reticulum protein with a JmjC domain.

Jumonji C (JmjC) domain-containing proteins have been shown to regulate cellular processes by hydroxylating or demethylating histone and non-histone targets. JMJD8 belongs to the JmjC domain-only family that was recently shown to be involved in angiogenesis and TNF-induced NF-κB signaling. Here, we employed bioinformatic analysis and immunofluorescence microscopy to examine the physiological properties of JMJD8. We demonstrated that JMJD8 localizes to the lumen of endoplasmic reticulum and that JMJD8 forms dimers or oligomers in vivo. Furthermore, we identified potential JMJD8-interacting proteins that are known to regulate protein complex assembly and protein folding. Taken together, this work demonstrates that JMJD8 is the first JmjC domain-containing protein found in the lumen of the endoplasmic reticulum that may function in protein complex assembly and protein folding.

Transcriptomic analysis of the role of RasGEF1B circular RNA in the TLR4/LPS pathway.

Circular RNAs (circRNAs) have recently emerged as a large class of novel non-coding RNA species. However, the detailed functional significance of the vast majority of them remains to be elucidated. Most functional characterization studies targeting circRNAs have been limited to resting cells, leaving their role in dynamic cellular responses to stimuli largely unexplored. In this study, we focus on the LPS-induced cytoplasmic circRNA, mcircRasGEF1B, and combine targeted mcircRasGEF1B depletion with high-throughput transcriptomic analysis to gain insight into its function during the cellular response to LPS stimulation. We show that knockdown of mcircRasGEF1B results in altered expression of a wide array of genes. Pathway analysis revealed an overall enrichment of genes involved in cell cycle progression, mitotic division, active metabolism, and of particular interest, NF-κB, LPS signaling pathways, and macrophage activation. These findings expand the set of functionally characterized circRNAs and support the regulatory role of mcircRasGEF1B in immune response during macrophage activation and protection against microbial infections.

Identification of 5-Methoxy-2-(Diformylmethylidene)-3,3-Dimethylindole as an Anti-Influenza A Virus Agent.

Influenza virus is estimated to cause 3-5 million severe complications and about 250-500 thousand deaths per year. Different kinds of anti-influenza virus drugs have been developed. However, the emergence of drug resistant strains has presented a big challenge for efficient antiviral therapy. Indole derivatives have been shown to exhibit both antiviral and anti-inflammatory activities. In this study, we adopted a cell-based system to screen for potential anti-IAV agents. Four indole derivatives (named 525A, 526A, 527A and 528A) were subjected to the antiviral screening, of which 526A was selected for further investigation. We reported that pre-treating cells with 526A protects cells from IAV infection. Furthermore, 526A inhibits IAV replication by inhibiting the expression of IAV genes. Interestingly, 526A suppresses the activation of IRF3 and STAT1 in host cells and thus represses the production of type I interferon response and cytokines in IAV-infected cells. Importantly, 526A also partially blocks the activation of RIG-I pathway. Taken together, these results suggest that 526A may be a potential anti-influenza A virus agent.

SNP variants associated with non-Hodgkin lymphoma (NHL) correlate with human leukocyte antigen (HLA) class II expression.

Large consortia efforts and genome-wide association studies (GWASs) have linked a number of genetic variants within the 6p21 chromosomal region to non-Hodgkin lymphoma (NHL). Complementing these efforts, we genotyped previously reported SNPs in the human leukocyte antigen (HLA) class I (rs6457327) and class II (rs9271100, rs2647012 and rs10484561) regions in a total of 1,145 subjects (567 NHL cases and 578 healthy controls) from two major ethnic groups in Malaysia, the Malays and the Chinese. We identified a NHL-associated (PNHL_add = 0.0008; ORNHL_add = 0.54; 95% CI = 0.37-0.77) and B-cell associated (PBcell_add = 0.0007; ORBcell_add = 0.51; 95% CI = 0.35-0.76) SNP rs2647012 in the Malaysian Malays. In silico cis-eQTL analysis of rs2647012 suggests potential regulatory function of nearby HLA class II molecules. Minor allele rs2647012-T is linked to higher expression of HLA-DQB1, rendering a protective effect to NHL risk. Our findings suggest that the HLA class II region plays an important role in NHL etiology.

JMJD8 is a positive regulator of TNF-induced NF-κB signaling.

TNF-induced signaling mediates pleiotropic biological consequences including inflammation, immunity, cell proliferation and apoptosis. Misregulation of TNF signaling has been attributed as a major cause of chronic inflammatory diseases and cancer. Jumonji domain-containing protein 8 (JMJD8) belongs to the JmjC family. However, only part of the family members has been described as hydroxylase enzymes that function as histone demethylases. Here, we report that JMJD8 positively regulates TNF-induced NF-κB signaling. Silencing the expression of JMJD8 using RNA interference (RNAi) greatly suppresses TNF-induced expression of several NF-κB-dependent genes. Furthermore, knockdown of JMJD8 expression reduces RIP ubiquitination, IKK kinase activity, delays IκBα degradation and subsequently blocks nuclear translocation of p65. In addition, JMJD8 deficiency enhances TNF-induced apoptosis. Taken together, these findings indicate that JMJD8 functions as a positive regulator of TNF-induced NF-κB signaling.

Inclusion bodies of recombinant Epstein-Barr virus capsid antigen p18 as potential immobilized antigens in enzyme immunoassays for detection of nasopharyngeal carcinoma.

BACKGROUND: Development of indirect enzyme-linked immunosorbent assays (ELISAs) often utilizes synthetic peptides or recombinant proteins from Escherichia coli as immobilized antigens. Because inclusion bodies (IBs) formed during recombinant protein expression in E. coli are commonly thought as misfolded aggregates, only refolded proteins from IBs are used to develop new or in-house diagnostic assays. However, the promising utilities of IBs as nanomaterials and immobilized enzymes as shown in recent studies have led us to explore the potential use of IBs of recombinant Epstein-Barr virus viral capsid antigen p18 (VCA p18) as immobilized antigens in ELISAs for serologic detection of nasopharyngeal carcinoma (NPC). METHODS: Thioredoxin fusion VCA p18 (VCA-Trx) and IBs of VCA p18 without fusion tags (VCA-IBs) were purified from E. coli. The diagnostic performances of IgG/VCA-IBs, IgG/VCA-Denat-IBs (using VCA-IBs coated in 8mol/l urea), IgG/VCA-Trx, and IgG/VCA-Peptide assays were compared by screening 100 NPC case-control pairs. RESULTS: The IgG/VCA-Denat-IBs assay showed the best area under the receiver operating characteristic curve (AUC: 0.802; p<0.05), while the AUCs for the IgG/VCA-IBs, IgG/VCA-Trx, and IgG/VCA-Peptide assays were comparable (AUC: 0.740, 0.727, and 0.741, respectively). CONCLUSION: We improved the diagnostic performance of the ELISA significantly using IBs of recombinant VCA p18.

Analysis of interleukin-10 promoter single nucleotide polymorphisms and risk of non-Hodgkin lymphoma in a Malaysian population.

We evaluated the association of two IL10 single nucleotide polymorphisms (SNPs) (rs1800896 and rs1800871) with non-Hodgkin lymphoma (NHL) risk in the three major races of the Malaysian population (Malay, Chinese and Indian; 317 cases and 330 controls). Our initial screening demonstrated that rs1800871 but not rs1800896 was significantly associated with increased NHL risk in Malays (pMalay-Rec = 0.007) and Chinese only (pChinese-Rec = 0.039). Subsequent combined analysis of the Malay and Chinese revealed significant association of rs1800871 with all (ALL) NHL subtypes (pMeta-ALL-NHL-Rec = 0.001), ALL B-cell subtypes (pMeta-ALL-B-cell-Rec = 0.003), diffuse large B-cell lymphoma (DLBCL) subtype (pMeta-DLBCL-Rec = 0.002) and ALL T-cell subtypes (pMeta-ALL-T-cell-Rec = 0.031). SNP rs1800896 showed increased risk only in follicular lymphoma (FL) (pMeta-FL-Dom = 0.0004). We also detected a male-specific association of rs1800871 with increased NHL risk (pMeta-Male-ALL-NHL-Rec = 0.006) in the combined analysis. To our knowledge, this is the first report on the association of IL10 promoter SNPs with NHL susceptibility in the three major races of Malaysia.

Inhibition of euchromatic histone methyltransferase 1 and 2 sensitizes chronic myeloid leukemia cells to interferon treatment.

BACKGROUND: H3K9 methylation is one of the essential histone post-translational modifications for heterochromatin formation and transcriptional repression. Recently, several studies have demonstrated that H3K9 methylation negatively regulates the type I interferon response. RESULTS: We report the application of EHMT1 and EHMT2 specific chemical inhibitors to sensitize CML cell lines to interferon and imatinib treatments. Inhibition of EHMT1 and EHMT2 with BIX01294 enhances the cytotoxicity of IFNα2a in four CML cell lines, K562, KCL22, BV173 and KT1 cells. Chromatin immunoprecipitation assay shows that BIX01294 treatment enhances type I interferon response by reducing H3K9me2 at the promoters of interferon-stimulated genes. Additionally, BIX01294 treatment augments IFNα2a- and imatinib-mediated apoptosis in CML cell lines. Moreover, our data suggest that the expression level of EHMT1 and EHMT2 inversely correlates with the type I interferon responsiveness in CML cell lines. CONCLUSIONS: Our study sheds light on the role of EHMT1 and EHMT2 as potential targets in improving the efficacy of standard treatments of CML.

Identification of an enhancer that increases miR-200b~200a~429 gene expression in breast cancer cells.

The miR-200b~200a~429 gene cluster is a key regulator of EMT and cancer metastasis, however the transcription-based mechanisms controlling its expression during this process are not well understood. We have analyzed the miR-200b~200a~429 locus for epigenetic modifications in breast epithelial and mesenchymal cell lines using chromatin immunoprecipitation assays and DNA methylation analysis. We discovered a novel enhancer located approximately 5.1kb upstream of the miR-200b~200a~429 transcriptional start site. This region was associated with the active enhancer chromatin signature comprising H3K4me1, H3K27ac, RNA polymerase II and CpG dinucleotide hypomethylation. Luciferase reporter assays revealed the upstream enhancer stimulated the transcription of the miR-200b~200a~429 minimal promoter region approximately 27-fold in breast epithelial cells. Furthermore, we found that a region of the enhancer was transcribed, producing a short, GC-rich, mainly nuclear, non-polyadenylated RNA transcript designated miR-200b eRNA. Over-expression of miR-200b eRNA had little effect on miR-200b~200a~429 promoter activity and its production did not correlate with miR-200b~200a~429 gene expression. While additional investigations of miR-200b eRNA function will be necessary, it is possible that miR-200b eRNA may be involved in the regulation of miR-200b~200a~429 gene expression and silencing. Taken together, these findings reveal the presence of a novel enhancer, which contributes to miR-200b~200a~429 transcriptional regulation in epithelial cells.

An autocrine TGF-beta/ZEB/miR-200 signaling network regulates establishment and maintenance of epithelial-mesenchymal transition.

Epithelial-mesenchymal transition (EMT) is a form of cellular plasticity that is critical for embryonic development and tumor metastasis. A double-negative feedback loop involving the miR-200 family and ZEB (zinc finger E-box-binding homeobox) transcription factors has been postulated to control the balance between epithelial and mesenchymal states. Here we demonstrate using the epithelial Madin Darby canine kidney cell line model that, although manipulation of the ZEB/miR-200 balance is able to repeatedly switch cells between epithelial and mesenchymal states, the induction and maintenance of a stable mesenchymal phenotype requires the establishment of autocrine transforming growth factor-ß (TGF-ß) signaling to drive sustained ZEB expression. Furthermore, we show that prolonged autocrine TGF-ß signaling induced reversible DNA methylation of the miR-200 loci with corresponding changes in miR-200 levels. Collectively, these findings demonstrate the existence of an autocrine TGF-ß/ZEB/miR-200 signaling network that regulates plasticity between epithelial and mesenchymal states. We find a strong correlation between ZEBs and TGF-ß and negative correlations between miR-200 and TGF-ß and between miR-200 and ZEBs, in invasive ductal carcinomas, consistent with an autocrine TGF-ß/ZEB/miR-200 signaling network being active in breast cancers.