Noncoding RNA Loci Research Articles

Tumor suppressor miR-317 and lncRNA Peony are expressed from a polycistronic non-coding RNA locus that regulates germline differentiation and testis morphology.

This research focuses on investigating the impact of non-coding RNAs on stem cell biology and differentiation processes. We found that miR-317 plays a role in germline stem cell progeny differentiation. miR-317 and its neighbor, the lncRNA Peony, originate and are co-expressed from a singular polycistronic non-coding RNA locus. Alternative polyadenylation is implicated in regulation of their differential expression. While the increased expression of the lncRNA Peony results in the disruption of the muscle sheath covering the testis, the absence of miR-317 leads to the emergence of germline tumors in young flies. The deficiency of miR-317 increases Notch signaling activity in the somatic cyst cells, which drives germline tumorigenesis. Germline tumors also arise from upregulation of several predicted targets of miR-317 , among which are regulators of the Notch pathway. This implicates miR-317 as a novel tumor suppressor that modulates Notch signaling strength.

Precision run-on sequencing (PRO-seq) for microbiome transcriptomics.

Bacteria respond to environmental stimuli through precise regulation of transcription initiation and elongation. Bulk RNA sequencing primarily characterizes mature transcripts, so to identify actively transcribed loci we need to capture RNA polymerase (RNAP) complexed with nascent RNA. However, such capture methods have only previously been applied to culturable, genetically tractable organisms such as E. coli and B. subtilis. Here we apply precision run-on sequencing (PRO-seq) to profile nascent transcription in cultured E. coli and diverse uncultured bacteria. We demonstrate that PRO-seq can characterize the transcription of small, structured, or post-transcriptionally modified RNAs, which are often absent from bulk RNA-seq libraries. Applying PRO-seq to the human microbiome highlights taxon-specific RNAP pause motifs and pause-site distributions across non-coding RNA loci that reflect structure-coincident pausing. We also uncover concurrent transcription and cleavage of CRISPR guide RNAs and transfer RNAs. We demonstrate the utility of PRO-seq for exploring transcriptional dynamics in diverse microbial communities.

SAT432 Interactions Of Endocrine-disrupting Chemical (EDC) Exposure And Fear Recall On Small Noncoding RNA (sncRNA) Contents Of Rat Caput Extracellular Vesicles

Abstract Disclosure: R. Gillette: None. M. Aguirre: None. K. Bell: None. Y. Nadeem: None. L. Thompson: None. A.C. Gore: None. Extracellular vesicles (EVs) transfer a complex repertoire of small non-coding RNAs (sncRNA) to sperm as they transit the caput of the epididymis. These molecules are essential for the functional competence of sperm and are directly influenced by environmental factors such as stress, diet, and drug use. We hypothesized that prenatal exposure to EDCs would alter the sncRNA contents of caput extracellular vesicles, and that an exposure to a previously established contextual fear cue would exacerbate this effect. Pregnant dams were fed Nilla Wafers with either 3% DMSO (vehicle control) or Aroclor 1221 (1 mg/kg; a weakly estrogenic EDC) from embryonic day (E)8 until E18 and post-natal day (PND)1 to P21. In adulthood male rat offspring were subjected to a contextual fear conditioning paradigm in which each rat was exposed to an odorant paired with a shock (four 500 ms shocks at 0.7 mA) for three consecutive days. Two weeks later, all rats were returned to the conditioning chambers and exposed to the odorant for 5 minutes (no shocks). The next day, male rats were euthanized and the epididymis was dissected. EVs and total RNA isolated from the caput portion were sequenced on a NovaSeq 6000 SP in a single-end 75 bp format. Reads were checked for quality, trimmed for adapter sequences, aligned to the rat genome, and assigned an annotation designation across 10 different categories of sncRNA. We found that the majority of caput EV sncRNA were either tRNA fragments or piRNA and that unlike mice, miRNA accounted for a small portion of the total reads. The third largest category of aligned reads was in intergenic space and not associated with canonical small non-coding (snc) RNA loci. In-depth investigation determined these reads (∼19 nt) aligned strictly within the boundaries of known CpG islands (CpGi), which have not previously been reported to express any form of sncRNA. Re-exposure to a fear-associated scent 24 hours prior to sample collection substantially reduced the proportion tRNAs in EVs, which were supplanted by an increase in piRNA, CpGi sRNA, and lncRNA. This effect was exacerbated by prenatal EDC exposure and was particularly prevalent in the lncRNA and CpGi sRNA categories. Taken together our data show that the small RNA contents of caput epididymosomes of male rats is substantially different from previous observations reported in mice. We demonstrate that exposure to a contextual fear cue is a powerful influence on sncRNA composition of caput extracellular vesicles and that it is uniquely altered by prenatal EDC exposure. Finally, we provide evidence for what we believe may be a novel category of sncRNA. Presentation: Saturday, June 17, 2023

Myeloid leukemia vulnerabilities embedded in long noncoding RNA locus MYNRL15

The noncoding genome presents a largely untapped source of new biological insights, including thousands of long noncoding RNA (lncRNA) loci. While lncRNA dysregulation has been reported in myeloid malignancies, their functional relevance remains to be systematically interrogated. We performed CRISPRi screens of lncRNA signatures from normal and malignant hematopoietic cells and identified MYNRL15 as a myeloid leukemia dependency. Functional dissection suggests an RNA-independent mechanism mediated by two regulatory elements embedded in the locus. Genetic perturbation of these elements triggered a long-range chromatin interaction and downregulation of leukemia dependency genes near the gained interaction sites, as well as overall suppression of cancer dependency pathways. Thus, this study describes a new noncoding myeloid leukemia vulnerability and mechanistic concept for myeloid leukemia. Importantly, MYNRL15 perturbation caused strong and selective impairment of leukemia cells of various genetic backgrounds over normal hematopoietic stem and progenitor cells invitro, and depletion of patient-derived xenografts invivo.

The role of microRNAs in responses to drought and heat stress in peanut (Arachis hypogaea).

MicroRNAs (miRNAs) are 21-24nt small RNAs (sRNAs) that negatively regulate protein-coding genes and/or trigger phased small-interfering RNA (phasiRNA) production. Two thousand nine hundred miRNA families, of which ∼40 are deeply conserved, have been identified in ∼80 different plant species genomes. miRNA functions in response to abiotic stresses is less understood than their roles in development. Only seven peanut MIRNA families are documented in miRBase, yet a reference genome assembly is now published and over 480 plant-like MIRNA loci were predicted in the diploid peanut progenitor Arachis duranensis genome. We explored by computational analysis of a leaf sRNA library and publicly available sRNA, degradome, and transcriptome datasets the miRNA and phasiRNA space associated with drought and heat stresses in peanut. We characterized 33 novel candidate and 33 ancient conserved families of MIRNAs and present degradome evidence for their cleavage activities on mRNA targets, including several noncanonical targets and novel phasiRNA-producing noncoding and mRNA loci with validated novel targets such as miR1509 targeting serine/threonine-protein phosphatase7 and miRc20 and ahy-miR3514 targeting penta-tricopeptide repeats (PPRs), in contradistinction to other claims of miR1509/173/7122 superfamily miRNAs indirectly targeting PPRs via TAS-like noncoding RNA loci. We characterized the inverse correlations of significantly differentially expressed drought- and heat-regulated miRNAs, assayed by sRNA blots or transcriptome datasets, with target mRNA expressions in the same datasets. Meta-analysis of an expression atlas and over representation of miRNA target genes in co-expression networks suggest that miRNAs have functions in unique aspects of peanut gynophore development. Genome-wide MIRNA annotation of the published allopolyploid peanut genome can facilitate molecular breeding of value-added traits.

Complex regulation of Eomes levels mediated through distinct functional features of the Meteor long non-coding RNA locus

Long non-coding RNAs (lncRNAs) are implicated in a plethora of cellular processes, but an in-depth understanding of their functional features or their mechanisms of action is currently lacking. Here we study Meteor, a lncRNA transcribed near the gene encoding EOMES, a pleiotropic transcription factor implicated in various processes throughout development and in adult tissues. Using a wide array of perturbation techniques, we show that transcription elongation through the Meteor locus is required for Eomes activation in mouse embryonic stem cells, with Meteor repression linked to a change in the subpopulation primed to differentiate to the mesoderm lineage. We further demonstrate that a distinct functional feature of the locus-namely, the underlying DNA element-is required for suppressing Eomes expression following neuronal differentiation. Our results demonstrate the complex regulation that can be conferred by a single locus and emphasize the importance of careful selection of perturbation techniques when studying lncRNA loci.

A lncRNA from an inflammatory bowel disease risk locus maintains intestinal host-commensal homeostasis

Inflammatory bowel diseases (IBD) are known to have complex, genetically influenced etiologies, involving dysfunctional interactions between the intestinal immune system and the microbiome. Here, we characterized how the RNA transcript from an IBD-associated long non-coding RNA locus (“CARINH-Colitis Associated IRF1 antisense Regulator of Intestinal Homeostasis”) protects against IBD. We show that CARINH and its neighboring gene coding for the transcription factor IRF1 together form a feedforward loop in host myeloid cells. The loop activation is sustained by microbial factors, and functions to maintain the intestinal host-commensal homeostasis via the induction of the anti-inflammatory factor IL-18BP and anti-microbial factors called guanylate-binding proteins (GBPs). Extending these mechanistic insights back to humans, we demonstrate that the function of the CARINH/IRF1 loop is conserved between mice and humans. Genetically, the T allele of rs2188962, the most probable causal variant of IBD within the CARINH locus from the human genetics study, impairs the inducible expression of the CARINH/IRF1 loop and thus increases genetic predisposition to IBD. Our study thus illustrates how an IBD-associated lncRNA maintains intestinal homeostasis and protects the host against colitis.

De novo genes with an lncRNA origin encode unique human brain developmental functionality

Human de novo genes can originate from neutral long non-coding RNA (lncRNA) loci and are evolutionarily significant in general, yet how and why this all-or-nothing transition to functionality happens remains unclear. Here, in 74 human/hominoid-specific de novo genes, we identified distinctive U1 elements and RNA splice-related sequences accounting for RNA nuclear export, differentiating mRNAs from lncRNAs, and driving the origin of de novo genes from lncRNA loci. The polymorphic sites facilitating the lncRNA–mRNA conversion through regulating nuclear export are selectively constrained, maintaining a boundary that differentiates mRNAs from lncRNAs. The functional new genes actively passing through it thus showed a mode of pre-adaptive origin, in that they acquire functions along with the achievement of their coding potential. As a proof of concept, we verified the regulations of splicing and U1 recognition on the nuclear export efficiency of one of these genes, the ENSG00000205704, in human neural progenitor cells. Notably, knock-out or over-expression of this gene in human embryonic stem cells accelerates or delays the neuronal maturation of cortical organoids, respectively. The transgenic mice with ectopically expressed ENSG00000205704 showed enlarged brains with cortical expansion. We thus demonstrate the key roles of nuclear export in de novo gene origin. These newly originated genes should reflect the novel uniqueness of human brain development.

The role of the long non-coding RNA gene PVT1 polymorphism in the development of multiple sclerosis

Long non-coding RNA locus PVT1 produces multiple regulatory RNA molecules via alternative splicing, and therefore is a promising object for the development of targeted therapies of different disorders. Performing associationanalysis of genetic variants rs4410871 andrs759648 at the PVT1 locus with multiple sclerosis (MS) in the ethnic group of Russians, Tatars, and Bashkirs from the Volga-Ural region of Russia, we established an association between rs759648*С and MS in the group of Tatars (OR=1.42, PBonf0.046), confirmed by the results of meta-analysisin the three ethnic groups (OR=1.28, P = 0.015). Haplotypic analysis has revealed an association of the rs759648*С/rs4410871*С haplotype with MS in the total study group (χ2=6.083, Рperm =0.012), whilers4410871*C/C + rs759648* C/C combination conferred an increased risk of MS (OR = 2.71, Pperm = 0.027) in this groupaccording to the data yielded by the multilocus analysis. The results of our analysis indicate the significance of the rs759648 polymorphism located near the DNA sequence producing miR-1208. Further studies will help to elucidate the molecular mechanisms of the involvement of PVT1 gene in the etiopathogenesis of MS.

From genotype to phenotype: genetics of mammalian long non-coding RNAs in vivo.

Genome-wide sequencing has led to the discovery of thousands of long non-coding RNA (lncRNA) loci in the human genome, but evidence of functional significance has remained controversial for many lncRNAs. Genetically engineered model organisms are considered the gold standard for linking genotype to phenotype. Recent advances in CRISPR-Cas genome editing have led to a rapid increase in the use of mouse models to more readily survey lncRNAs for functional significance. Here, we review strategies to investigate the physiological relevance of lncRNA loci by highlighting studies that have used genetic mouse models to reveal key in vivo roles for lncRNAs, from fertility to brain development. We illustrate how an investigative approach, starting with whole-gene deletion followed by transcription termination and/or transgene rescue strategies, can provide definitive evidence for the in vivo function of mammalian lncRNAs.

Myeloid Leukemia Dependencies at CTCF-Enriched Long Noncoding RNA Loci

The noncoding genome presents a largely untapped source of biological insights, including tens of thousands of long noncoding RNA (lncRNA) loci. While some produce bona fide lncRNAs, others exert transcript-independent cis-regulatory effects, and the lack of predictive features renders their mechanistic dissection highly challenging. Here, we describe CTCF-enriched lncRNA loci (C-LNC) as a putative new subclass of functional genetic elements exemplified by MYNRL15 - myeloid leukemia noncoding regulatory locus on chromosome 15.Initially identified by an expression-guided CRISPRi screen of hematopoietic stem and progenitor (HSPC) / acute myeloid leukemia (AML) lncRNA signatures (480 genes, 1545 sgRNAs), we found MYNRL15 dependency in myeloid leukemia cells of diverse genetic backgrounds. Interestingly, cis and trans perturbation approaches revealed both the MYNRL15 transcript and its flanking protein-coding genes to be dispensable. High density CRISPR tiling of a 15 kb area centered on MYNRL15 (1613 sgRNAs) instead uncovered two crucial, candidate cis-regulatory DNA elements in the locus, which drive the MYNRL15 perturbation phenotype. To determine the molecular basis of MYNRL15 dependence, we performed transcriptome, chromatin conformation, chromatin accessibility, and CTCF profiling. RNA-sequencing established MYNRL15's involvement in maintaining key cancer dependency pathways (e.g. cell cycle, ribosome, spliceosome). Further, MYNRL15 perturbation associated with the coordinated dysregulation of several chromosome 15 neighbourhoods, and formation of a long-range chromatin interaction between the locus and the base of a distal loop, as detected via next-generation Capture-C. The gained interaction was accompanied by diffuse gains in chromatin accessibility across the distal interaction sites (ATAC-seq) as well as reduced CTCF occupancy at the MYRNL15 locus (CTCF CUT&RUN), altogether indicating the 3D re-organization of chromosome 15 following MYNRL15 perturbation. Integrative analysis of the chromatin conformation and transcriptome data, combined with a small CRISPR-Cas9 knockout screen of protein-coding genes from the gained interaction region (29 genes, 149 sgRNAs), pinpointed two potent cancer dependency genes that are located in the region and downregulated following MYNRL15 perturbation: namely, WDR61 and IMP3. Individual knockout of both genes robustly depleted myeloid leukemia cells, recapitulating the MYNRL15 perturbation phenotype and positioning WDR61 and IMP3 as its regulatory targets. Importantly, in primary cells, MYNRL15 perturbation eradicated AML blasts while sparing 50-60% of CD34 + HSPCs in vitro, and reduced patient-derived AML xenografts up to 10-fold in vivo, indicating a potential therapeutic window.Having implicated MYNRL15 in 3D genome organization and demonstrated its role in myeloid leukemia cells, we explored whether MYNRL15 may belong to a sub-category of biologically relevant lncRNA loci that have thus far been overlooked due to their lack of transcript-specific functions. Remarkably, elevated CTCF density (e.g. number of CTCF binding sites per kb of gene length) distinguishes MYNRL15 and 531 other lncRNA loci in K562 cells, of which 43-54% associate with genetic subgroups and/or survival in AML patient cohorts, and 18.4% are functionally required for leukemia maintenance as determined by CRISPR-Cas9 screening. The latter hit identification rate represents a substantial improvement over typical lncRNA essentiality screens (which range from 2-6%) - illustrating the effectiveness of CTCF density metrics in refining functional lncRNA candidate lists, and underlining the relevance such loci hold for AML and cancer pathophysiology in general. Curated C-LNC catalogs in other cell types will facilitate the search for noncoding oncogenic vulnerabilities in AML and other malignancies. [Display omitted] DisclosuresReinhardt: Celgene Corporation: Consultancy; Novartis: Consultancy; Bluebird Bio: Consultancy; Janssen: Consultancy; CLS Behring: Research Funding; Roche: Research Funding. Klusmann: Bluebird Bio: Consultancy; Novartis: Consultancy; Roche: Consultancy; Jazz Pharmaceuticals: Consultancy.

Genome-wide CRISPR interference screen identifies long non-coding RNA loci required for differentiation and pluripotency.

Although many long non-coding RNAs (lncRNAs) exhibit lineage-specific expression, the vast majority remain functionally uncharacterized in the context of development. Here, we report the first described human embryonic stem cell (hESC) lines to repress (CRISPRi) or activate (CRISPRa) transcription during differentiation into all three germ layers, facilitating the modulation of lncRNA expression during early development. We performed an unbiased, genome-wide CRISPRi screen targeting thousands of lncRNA loci expressed during endoderm differentiation. While dozens of lncRNA loci were required for proper differentiation, most differentially expressed lncRNAs were not, supporting the necessity for functional screening instead of relying solely on gene expression analyses. In parallel, we developed a clustering approach to infer mechanisms of action of lncRNA hits based on a variety of genomic features. We subsequently identified and validated FOXD3-AS1 as a functional lncRNA essential for pluripotency and differentiation. Taken together, the cell lines and methodology described herein can be adapted to discover and characterize novel regulators of differentiation into any lineage.

Genome-Wide Analysis of the Regulatory Landscape in Different Pre-B-ALL Subtypes

BackgroundAcute pre-B leukemia results from a failure in cell differentiation. In normal development, cells commit to their lineage and differentiated phenotype state by the activation of transcription factors (TFs), or TF modules, in a series of decision points at which a choice is made between alternative lineage fates. Previously, we showed that the ETV6-RUNX1 gene translocation that occurs in 25% of pediatric leukemia patients silences regulatory regions containing RUNX binding sites (enhancers) at target gene loci implicated in lymphoid lineage differentiation (Teppo et al. 2016). However, the genome-wide changes in enhancer activities at binding sites of other TFs functioning in B-cell differentiation remain poorly understood. Moreover, the genome-wide identification of enhancers altered during leukemogenesis is incomplete.Material & MethodsTo measure transcripts generated both at enhancer and gene regions, global run-on sequencing (GRO-seq) technique was used to measure nascent RNAs (Heinäniemi et al. 2016) from human pre-B-ALL patient samples (N=35) and several cell lines representing different subtypes. For validation, additional samples were generated from cell line models with overexpression or silencing of leukemic TFs. To this collection of 151 nascent transcriptomes, additional 92 GRO-seq samples representing different cell lineages (Liivrand et al 2017) were used to compare the leukemic cells to primary cells and other cancers. Enhancer RNA (eRNA) levels were compared at TF binding motif centered enhancer regions. TF ChIP-seq data from leukemic cells and HSC were acquired from the NCBI Gene Expression Omnibus for validation of this approach.ResultsTo gain insight into the derailed regulatory landscape in leukemic cells, we identified eRNA regions both from normal and leukemic samples and performed statistical analysis distinguishing regulatory regions with differential activity between pre-B-ALL-subtypes or compared to normal CD34+ bone marrow precursors. Next, association of the significantly altered enhancers with gene transcripts was utilized to build regulatory maps for differentially expressed gene loci. For this purpose, we developed a transcript identification and quantification pipeline suitable for detecting altered expression at different coding and non-coding RNA loci from nascent transcriptomes (Liivrand et al 2017; unpublished). To globally evaluate altered TF activities across pre-B-ALL subtypes underlying the differences observed, we present here an unbiased approach based on eRNA quantification at TF motif-centered enhancers. The analysis of eRNA levels at binding sites of TFs with established role in leukemia was used as a benchmark: we could confirm low eRNA transcription at RUNX1 binding sites in the ETV6-RUNX1 subtype, elevated eRNA levels at HOXA9 and MEIS1 binding sites in the MLL rearranged subtype, while TCF3, BCL6 and PBX1 bound sites were most active in TCF3-PBX1 samples.ConclusionsOur results provide the first genome-wide transcriptional regulatory landscape of different pre-B-ALL subtypes. We identified novel subtype-specific active enhancers and annotated regulatory non-coding regions that can be utilized to evaluate the functional impact of genetic alterations outside gene regions. DisclosuresNo relevant conflicts of interest to declare.

Intergenic RNA mainly derives from nascent transcripts of known genes

BackgroundEukaryotic genomes undergo pervasive transcription, leading to the production of many types of stable and unstable RNAs. Transcription is not restricted to regions with annotated gene features but includes almost any genomic context. Currently, the source and function of most RNAs originating from intergenic regions in the human genome remain unclear.ResultsWe hypothesize that many intergenic RNAs can be ascribed to the presence of as-yet unannotated genes or the “fuzzy” transcription of known genes that extends beyond the annotated boundaries. To elucidate the contributions of these two sources, we assemble a dataset of more than 2.5 billion publicly available RNA-seq reads across 5 human cell lines and multiple cellular compartments to annotate transcriptional units in the human genome. About 80% of transcripts from unannotated intergenic regions can be attributed to the fuzzy transcription of existing genes; the remaining transcripts originate mainly from putative long non-coding RNA loci that are rarely spliced. We validate the transcriptional activity of these intergenic RNAs using independent measurements, including transcriptional start sites, chromatin signatures, and genomic occupancies of RNA polymerase II in various phosphorylation states. We also analyze the nuclear localization and sensitivities of intergenic transcripts to nucleases to illustrate that they tend to be rapidly degraded either on-chromatin by XRN2 or off-chromatin by the exosome.ConclusionsWe provide a curated atlas of intergenic RNAs that distinguishes between alternative processing of well-annotated genes from independent transcriptional units based on the combined analysis of chromatin signatures, nuclear RNA localization, and degradation pathways.

Chromosome-level genome assembly of the hard-shelled mussel Mytilus coruscus, a widely distributed species from the temperate areas of East Asia.

BackgroundThe hard-shelled mussel (Mytilus coruscus) is widely distributed in the temperate seas of East Asia and is an important commercial bivalve in China. Chromosome-level genome information of this species will contribute not only to the development of hard-shelled mussel genetic breeding but also to studies on larval ecology, climate change biology, marine biology, aquaculture, biofouling, and antifouling.FindingsWe applied a combination of Illumina sequencing, Oxford Nanopore Technologies sequencing, and high-throughput chromosome conformation capture technologies to construct a chromosome-level genome of the hard-shelled mussel, with a total length of 1.57 Gb and a median contig length of 1.49 Mb. Approximately 90.9% of the assemblies were anchored to 14 linkage groups. We assayed the genome completeness using BUSCO. In the metazoan dataset, the present assemblies have 89.4% complete, 1.9% incomplete, and 8.7% missing BUSCOs. Gene modeling enabled the annotation of 37,478 protein-coding genes and 26,917 non-coding RNA loci. Phylogenetic analysis showed that M. coruscus is the sister taxon to the clade including Modiolus philippinarum and Bathymodiolus platifrons. Conserved chromosome synteny was observed between hard-shelled mussel and king scallop, suggesting that this is shared ancestrally. Transcriptomic profiling indicated that the pathways of catecholamine biosynthesis and adrenergic signaling in cardiomyocytes might be involved in metamorphosis.ConclusionsThe chromosome-level assembly of the hard-shelled mussel genome will provide novel insights into mussel genome evolution and serve as a fundamental platform for studies regarding the planktonic-sessile transition, genetic diversity, and genomic breeding of this bivalve.

Genomic Space of MGMT in Human Glioma Revisited: Novel Motifs, Regulatory RNAs, NRF1, 2, and CTCF Involvement in Gene Expression.

Background: The molecular regulation of increased MGMT expression in human brain tumors, the associated regulatory elements, and linkages of these to its epigenetic silencing are not understood. Because the heightened expression or non-expression of MGMT plays a pivotal role in glioma therapeutics, we applied bioinformatics and experimental tools to identify the regulatory elements in the MGMT and neighboring EBF3 gene loci. Results: Extensive genome database analyses showed that the MGMT genomic space was rich in and harbored many undescribed RNA regulatory sequences and recognition motifs. We extended the MGMT’s exon-1 promoter to 2019 bp to include five overlapping alternate promoters. Consensus sequences in the revised promoter for (a) the transcriptional factors CTCF, NRF1/NRF2, GAF, (b) the genetic switch MYC/MAX/MAD, and (c) two well-defined p53 response elements in MGMT intron-1, were identified. A putative protein-coding or non-coding RNA sequence was located in the extended 3′ UTR of the MGMT transcript. Eleven non-coding RNA loci coding for miRNAs, antisense RNA, and lncRNAs were identified in the MGMT-EBF3 region and six of these showed validated potential for curtailing the expression of both MGMT and EBF3 genes. ChIP analysis verified the binding site in MGMT promoter for CTCF which regulates the genomic methylation and chromatin looping. CTCF depletion by a pool of specific siRNA and shRNAs led to a significant attenuation of MGMT expression in human GBM cell lines. Computational analysis of the ChIP sequence data in ENCODE showed the presence of NRF1 in the MGMT promoter and this occurred only in MGMT-proficient cell lines. Further, an enforced NRF2 expression markedly augmented the MGMT mRNA and protein levels in glioma cells. Conclusions: We provide the first evidence for several new regulatory components in the MGMT gene locus which predict complex transcriptional and posttranscriptional controls with potential for new therapeutic avenues.

The Missing Link Between Cancer-Associated Variants and LncRNAs

Despite studies indicating that long noncoding RNAs, or lncRNAs, can act as proto-oncogenes, the implications of large numbers of cancer-associated variants found within noncoding RNA loci remain largely unknown. Here, we draw upon emerging studies to speculate on how variants of lncRNAs might play a role in cancer development.

Extensive nuclear gyration and pervasive non-genic transcription during primordial germ cell development in zebrafish.

ABSTRACTPrimordial germ cells (PGCs) are the precursors of germ cells, which migrate to the genital ridge during early development. Relatively little is known about PGCs after their migration. We studied this post-migratory stage using microscopy and sequencing techniques, and found that many PGC-specific genes, including genes known to induce PGC fate in the mouse, are only activated several days after migration. At this same time point, PGC nuclei become extremely gyrated, displaying general broad opening of chromatin and high levels of intergenic transcription. This is accompanied by changes in nuage morphology, expression of large loci (PGC-expressed non-coding RNA loci, PERLs) that are enriched for retro-transposons and piRNAs, and a rise in piRNA biogenesis signatures. Interestingly, no nuclear Piwi protein could be detected at any time point, indicating that the zebrafish piRNA pathway is fully cytoplasmic. Our data show that the post-migratory stage of zebrafish PGCs holds many cues to both germ cell fate establishment and piRNA pathway activation.

Aberrant expression of a novel circular RNA in pancreatic cancer.

Circular RNAs (circRNAs) are single-stranded, covalently closed RNA molecules that are produced from pre-mRNAs through a process known as back-splicing. Although circRNAs are expressed under specific conditions, current understanding of their comprehensive expression status is still limited. Here, we performed a large-scale circRNA profiling analysis in human pancreatic ductal adenocarcinoma (PDAC) tissues, using circular RNA-specific RNA sequencing. We identified more than 40,000 previously unknown circRNAs, some of which were upregulated in PDAC tissues, compared with normal pancreatic tissues. We determined the full-length sequence of a circRNA upregulated in PDAC, which was derived from two noncoding RNA loci on chromosome 12. The novel circRNA, named circPDAC RNA, was not expressed in normal human cells, but was expressed in PDAC and other carcinoma cells. While postulated biological functions, such as peptide production from the circPDAC RNA, were not detected, its aberrant expression was confirmed in other PDAC tissues and in serum from a PDAC patient. These results demonstrate that comprehensive studies are necessary to reveal the expression status of circRNAs and that the circPDAC RNA identified here might serve as a novel biomarker for cancers, including PDAC.

CARMAL Is a Long Non-coding RNA Locus That Regulates MFGE8 Expression.

Genome-wide association studies have identified several genetic loci linked to coronary artery disease (CAD) most of them located in non-protein coding regions of the genome. One such locus is the CAD Associated Region between MFGE8 and ABHD2 (CARMA), a ∼18 kb haplotype that was recently shown to regulate vicinal protein coding genes. Here, we further investigate the region by examining a long non-coding RNA gene locus (CARMAL/RP11-326A19.4/AC013565) abutting the CARMA region. Expression-genotype correlation analyses of public databases indicate that CARMAL levels are influenced by CAD associated variants suggesting that it might have cardioprotective functions. We found CARMAL to be stably expressed at relatively low levels and enriched in the cytosol. CARMAL function was investigated by several gene targeting approaches in HEK293T: inactive CRISPR fusion proteins, antisense, overexpression and inactivation by CRISPR-mediated knock-out. Modest increases in CARMAL (3–4×) obtained via CRISPRa using distinct single-guided RNAs did not result in consistent transcriptome effects. By contrast, CARMAL deletion or reduced CARMAL expression via CRISPRi increased MFGE8 levels, suggesting that CARMAL is contributing to reduce MFGE8 expression under basal conditions. While future investigations are required to clarify the mechanism(s) by which CARMAL acts on MFGE8, integrative bioinformatic analyses of the transcriptome of CARMAL deleted cells suggest that this locus may also be involved in leucine metabolism, splicing, transcriptional regulation and Shwachman-Bodian-Diamond syndrome protein function.