Search

Long non-coding RNAs (lncRNAs) have begun to receive overdue attention for their regulatory roles in gene expression and other cellular processes. Although most lncRNAs are lowly expressed and tissue-specific, notable exceptions include MALAT1 and its genomic neighbor NEAT1, two highly and ubiquitously expressed oncogenes with roles in transcriptional regulation and RNA splicing. Previous studies have suggested that NEAT1 is found only in mammals, while MALAT1 is present in all gnathostomes (jawed vertebrates) except birds. Here we show that these assertions are incomplete, likely due to the challenges associated with properly identifying these two lncRNAs. Using phylogenetic analysis and structure-aware annotation of publicly available genomic and RNA-seq coverage data, we show that NEAT1 is a common feature of tetrapod genomes except birds and squamates. Conversely, we identify MALAT1 in representative species of all major gnathostome clades, including birds. Our in-depth examination of MALAT1, NEAT1, and their genomic context in a wide range of vertebrate species allows us to reconstruct the series of events that led to the formation of the locus containing these genes in taxa from cartilaginous fish to mammals. This evolutionary history includes the independent loss of NEAT1 in birds and squamates, since NEAT1 is found in the closest living relatives of both clades (crocodilians and tuataras, respectively). These data clarify the origins and relationships of MALAT1 and NEAT1 and highlight an opportunity to study the change and continuity in lncRNA structure and function over deep evolutionary time.

The contribution of alloresponses to mismatched HLA-DP in solid organ transplantation and hematopoietic stem cell transplantation (HCT) has been well documented. Exploring the regulatory mechanisms of DPB1 alleles has become an important question to be answered. In this study, our initial investigation focused on examining the correlation between the rs9277534G/A SNP and DPB1 mRNA expression. The result showed that there was a significant increase in DPB1 mRNA expression in B lymphoblastoid cell lines (BLCLs) with the rs9277534GG genotype compared to rs9277534AA genotype. In addition, B cells with the rs9277534GG exhibited significantly higher DP protein expression than those carrying the rs9277534AA genotype in primary B cells. Furthermore, we observed a significant upregulation of DP expression in B cells following treatment with Interleukin 13 (IL-13) compared to untreated B cells carrying rs9277534GG-linked DPB1 alleles. Fluorescence in situ hybridization (FISH) analysis of DPB1 in BLCL demonstrated significant differences in both the cytoplasmic (p=0.0003) and nuclear (p=0.0001) localization of DP mRNA expression comparing DPB1*04:01 (rs9277534AA) and DPB1*05:01 (rs9277534GG) homozygous cells. The study of the correlation between differential DPB1 expression and long non-coding RNAs (lncRNAs) showed that lnc-HLA-DPB1-13:1 is strongly associated with DP expression (r=0.85), suggesting the potential involvement of lncRNA in regulating DP expression. The correlation of DP donor specific antibody (DSA) with B cell flow crossmatch (B-FCXM) results showed a better linear correlation of DP DSA against GG and AG donor cells (R2 = 0.4243, p=0.0025 and R2 = 0.6172, p=0.0003, respectively), compared to DSA against AA donor cells (R2 = 0.0649, p=0.4244). This explained why strong DP DSA with a low expression DP leads to negative B-FCXM. In conclusion, this study provides evidence supporting the involvement of lncRNA in modulating HLA-DP expression, shedding lights on the intricate regulatory mechanisms of DP, particularly under inflammatory conditions in transplantation.

Proper activation of cytotoxic T cells via the T cell receptor and the costimulatory receptor CD28 is essential for adaptive immunity against viruses, intracellular bacteria, and cancers. Through biochemical analysis of RNA:protein interactions, we uncovered a non-coding RNA circuit regulating activation and differentiation of cytotoxic T cells composed of the long non-coding RNA Malat1 (Metastasis Associated Lung Adenocarcinoma Transcript 1) and the microRNA family miR-15/16. miR-15/16 is a widely and highly expressed tumor suppressor miRNA family important for cell proliferation and survival. miR-15/16 play important roles in T cell responses to viral infection, including the regulation of antigen-specific T cell expansion and memory. Comparative Argonaute-2 high-throughput sequencing of crosslinking immunoprecipitation (AHC) combined with gene expression profiling in normal and miR-15/16-deficient mouse T cells revealed a large network of hundreds of direct miR-15/16 target mRNAs, many with functional relevance for T cell activation, survival and memory formation. Among these targets, Malat1 contained the largest absolute magnitude miR-15/16-dependent AHC peak. This binding site was among the strongest lncRNA:miRNA interactions detected in the T cell transcriptome. We used CRISPR targeting with homology directed repair to generate mice with a 5-nucleotide mutation in the miR-15/16-binding site in Malat1. This mutation interrupted Malat1:miR-15/16 interaction, and enhanced the repression of other miR-15/16 target genes, including CD28. Interrupting Malat1 interaction with miR-15/16 decreased cytotoxic T cell activation, including the expression of interleukin 2 (IL-2) and a broader CD28-responsive gene program. Accordingly, Malat1 mutation diminished memory cell persistence in mice following LCMV Armstrong and Listeria monocytogenes infection. This study marks a significant advance in the study of long non-coding RNAs in the immune system by ascribing cell-intrinsic, sequence-specific in vivo function to Malat1. These findings have implications for T cell-mediated autoimmune diseases, antiviral and anti-tumor immunity, as well as lung adenocarcinoma and other malignancies where Malat1 is overexpressed.

The existence of latent cellular reservoirs is recognized as the major barrier to an HIV cure. Reactivating and eliminating "shock and kill" or permanently silencing "block and lock" the latent HIV reservoir, as well as gene editing, remain promising approaches, but so far have proven to be only partially successful. Moreover, using latency reversing agents or "block and lock" drugs pose additional considerations, including the ability to cause cellular toxicity, a potential lack of specificity for HIV, or low potency when each agent is used alone. RNA molecules, such as microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) are becoming increasingly recognized as important regulators of gene expression. RNA-based approaches for combatting HIV latency represent a promising strategy since both miRNAs and lncRNAs are more cell-type and tissue specific than protein coding genes. Thus, a higher specificity of targeting the latent HIV reservoir with less overall cellular toxicity can likely be achieved. In this review, we summarize current knowledge about HIV gene expression regulation by miRNAs and lncRNAs encoded in the human genome, as well as regulatory molecules encoded in the HIV genome. We discuss both the transcriptional and post-transcriptional regulation of HIV gene expression to align with the current definition of latency, and describe RNA molecules that either promote HIV latency or have anti-latency properties. Finally, we provide perspectives on using each class of RNAs as potential targets for combatting HIV latency, and describe the complexity of the interactions between different RNA molecules, their protein targets, and HIV.

The long non-coding RNA (lncRNA) Metastasis-associated lung adenocarcinoma transcript 1 (MALAT1) maintains the integrity of the intestinal epithelial barrier and regulates local inflammation. However, its influences on intestinal microbial communities and tissue susceptibility to cancer development remain unexplored. Here, we report that MALAT1 regulates host anti-microbial response gene expression and the composition of mucosal-associated microbial communities in a region-specific manner. In the APC mutant mouse model of intestine tumorigenesis, knocking out MALAT1 results in higher polyp counts in the small intestine and colon. Interestingly, intestine polyps that developed in the absence of MALAT1 were smaller in size. These findings highlight the unexpected bivalent role of MALAT1 in restricting and promoting cancer progression at different disease stages. Among the 30 MALAT1-targets shared by both the small intestine and colon, ZNF638 and SENP8 levels are predictive of colon adenoma patient overall survival and disease-free survival. Genomic assays further revealed that MALAT1 modulates intestinal target expression and splicing through both direct and indirect mechanisms. This study expands the role of lncRNAs in regulating intestine homeostasis, microbial communities, and cancer pathogenesis.

Glioblastoma multiforme is the most common, highly aggressive malignant brain tumor which is marked by highest inter- and intra-tumoral heterogeneity. Despite, immunotherapy, and combination therapies developed; the clinical trials often result into large number of failures. Often cancer cells are known to communicate with surrounding cells in tumor microenvironment (TME). Extracellular vesicles (EVs) consisting of diverse cargo mediates this intercellular communication and is believed to modulate the immune function against GBM. Tumor-associated microglia (TAM), though being the resident innate immune cell of CNS, is known to attain pro-tumorigenic M2 phenotype, and this immunomodulation is aided by extracellular vesicle-mediated transfer of oncogenic, immunomodulatory molecules. Besides, oncogenic proteins, long non-coding RNAs (lncRNAs), are believed to carry oncogenic potential, and therefore, understanding the mechanism leading to microglial dysregulation mediated by GBM-derived extracellular vesicle (GDEV) lncRNAs becomes crucial. This review focuses on current understanding of role of GDEV and lncRNA in microglial dysfunction and its potential as a therapeutic target.

Domestic dogs (Canis lupus familiaris) are the most variable-sized mammalian species on Earth, displaying a 40-fold size difference between breeds.1 Although dogs of variable size are found in the archeological record,2-4 the most dramatic shifts in body size are the result of selection over the last two centuries, as dog breeders selected and propagated phenotypic extremes within closed breeding populations.5 Analyses of over 200 domestic breeds have identified approximately 20 body size genes regulating insulin processing, fatty acid metabolism, TGFβ signaling, and skeletal formation.6-10 Of these, insulin-like growth factor 1 (IGF1) predominates, controlling approximately 15% of body size variation between breeds.8 The identification of a functional mutation associated with IGF1 has thus far proven elusive.6,10,11 Here, to identify and elucidate the role of an ancestral IGF1 allele in the propagation of modern canids, we analyzed 1,431 genome sequences from 13 species, including both ancient and modern canids, thus allowing us to define the evolutionary history of both ancestral and derived alleles at this locus. We identified a single variant in an antisense long non-coding RNA (IGF1-AS) that interacts with the IGF1 gene, creating a duplex. While the derived mutation predominates in both modern gray wolves and large domestic breeds, the ancestral allele, which predisposes to small size, was common in small-sized breeds and smaller wild canids. Our analyses demonstrate that this major regulator of canid body size nearly vanished in Pleistocene wolves, before its recent resurgence resulting from human-imposed selection for small-sized breed dogs.

HOXA Transcript Antisense RNA, Myeloid-Specific 1 (HOTAIRM1) is a conserved long non-coding RNA (lncRNA) involved in myeloid and neural differentiation that is deregulated in acute myeloid leukemia and other cancers. Previous studies focused on the nuclear unspliced HOTAIRM1 transcript, however cytoplasmic splice variants exist whose roles have remained unknown. Here, we report novel functions of HOTAIRM1 in the kidney. HOTAIRM1 transcripts are induced during renal lineage differentiation of embryonic stem cells and required for expression of specific renal differentiation genes. We show that the major HOTAIRM1 transcript in differentiated cells is the spliced cytoplasmic HM1-3 isoform and that HM1-3 is downregulated in >90% of clear cell renal cell carcinomas (ccRCCs). Knockdown of HM1-3 in renal cells deregulates hypoxia-responsive and angiogenic genes, including ANGPTL4. Furthermore, HOTAIRM1 transcripts are downregulated by hypoxia-mimetic stress and knockdown of the cytoplasmic HM1-3 isoform in normoxic cells post-transcriptionally induces Hypoxia-Inducible Factor 1α (HIF1α) protein, a key activator of ANGPTL4. Our results demonstrate the pervasive downregulation of the specific HOTAIRM1 cytoplasmic isoform HM1-3 in ccRCC and suggest possible roles of HOTAIRM1 in kidney differentiation and suppression of HIF1-dependent angiogenic pathways.

Irritable bowel syndrome (IBS) is a brain-gut axis disorder characterized by abdominal pain and altered bowel habits. IBS is a multifactorial, stress-sensitive disorder with evidence for familial clustering attributed to genetic or shared environmental factors. However, there are weak genetic associations reported with IBS and a lack of evidence to suggest that major genetic factor(s) contribute to IBS pathophysiology. Studies on animal models of stress, including early life stress, suggest a role for environmental factors, specifically, stress associated with dysregulation of corticotropin releasing factor and hypothalamus-pituitary-adrenal (HPA) axis pathways in the pathophysiology of IBS. Recent evidence suggests that epigenetic mechanisms, which constitute molecular changes not driven by a change in gene sequence, can mediate environmental effects on central and peripheral function. Epigenetic alterations including DNA methylation changes, histone modifications, and differential expression of non-coding RNAs (microRNA [miRNA] and long non-coding RNA) have been associated with several diseases. The objective of this review is to elucidate the molecular factors in the pathophysiology of IBS with an emphasis on epigenetic mechanisms. Emerging evidence for epigenetic changes in IBS includes changes in DNA methylation in animal models of IBS and patients with IBS, and various miRNAs that have been associated with IBS and endophenotypes, such as increased visceral sensitivity and intestinal permeability. DNA methylation, in particular, is an emerging field in the realm of complex diseases and a promising mechanism which can provide important insights into IBS pathogenesis and identify potential targets for treatment.

AimsCausal transcripts at genomic loci associated with type 2 diabetes (T2D) are mostly unknown. The chr8p23.1 variant rs4841132, associated with an insulin-resistant diabetes risk phenotype, lies in the second exon of a long non-coding RNA (lncRNA) gene, LOC157273, located 175 kilobases from PPP1R3B, which encodes a key protein regulating insulin-mediated hepatic glycogen storage in humans. We hypothesized that LOC157273 regulates expression of PPP1R3B in human hepatocytes.MethodsWe tested our hypothesis using Stellaris fluorescent in situ hybridization to assess subcellular localization of LOC157273; small interfering RNA (siRNA) knockdown of LOC157273, followed by RT-PCR to quantify LOC157273 and PPP1R3B expression; RNA-seq to quantify the whole-transcriptome gene expression response to LOC157273 knockdown; and an insulin-stimulated assay to measure hepatocyte glycogen deposition before and after knockdown.ResultsWe found that siRNA knockdown decreased LOC157273 transcript levels by approximately 80%, increased PPP1R3B mRNA levels by 1.7-fold, and increased glycogen deposition by >50% in primary human hepatocytes. An A/G heterozygous carrier (vs. three G/G carriers) had reduced LOC157273 abundance due to reduced transcription of the A allele and increased PPP1R3B expression and glycogen deposition.ConclusionWe show that the lncRNA LOC157273 is a negative regulator of PPP1R3B expression and glycogen deposition in human hepatocytes and a causal transcript at an insulin-resistant T2D risk locus.

AimsCausal transcripts at genomic loci associated with type 2 diabetes (T2D) are mostly unknown. The chr8p23.1 variant rs4841132, associated with an insulin-resistant diabetes risk phenotype, lies in the second exon of a long non-coding RNA (lncRNA) gene, LOC157273, located 175 kilobases from PPP1R3B, which encodes a key protein regulating insulin-mediated hepatic glycogen storage in humans. We hypothesized that LOC157273 regulates expression of PPP1R3B in human hepatocytes.MethodsWe tested our hypothesis using Stellaris fluorescent in situ hybridization to assess subcellular localization of LOC157273; small interfering RNA (siRNA) knockdown of LOC157273, followed by RT-PCR to quantify LOC157273 and PPP1R3B expression; RNA-seq to quantify the whole-transcriptome gene expression response to LOC157273 knockdown; and an insulin-stimulated assay to measure hepatocyte glycogen deposition before and after knockdown.ResultsWe found that siRNA knockdown decreased LOC157273 transcript levels by approximately 80%, increased PPP1R3B mRNA levels by 1.7-fold, and increased glycogen deposition by >50% in primary human hepatocytes. An A/G heterozygous carrier (vs. three G/G carriers) had reduced LOC157273 abundance due to reduced transcription of the A allele and increased PPP1R3B expression and glycogen deposition.ConclusionWe show that the lncRNA LOC157273 is a negative regulator of PPP1R3B expression and glycogen deposition in human hepatocytes and a causal transcript at an insulin-resistant T2D risk locus.

BACKGROUND:Muscle-invasive bladder cancer (MIBC) is a heterogeneous disease, and gene expression profiling has identified several molecular subtypes with distinct biological and clinicopathological characteristics. While MIBC subtyping has primarily been based on messenger RNA (mRNA), long non-coding RNAs (lncRNAs) may provide additional resolution. METHODS:LncRNA expression was quantified from microarray data of a MIBC cohort treated with neoadjuvant chemotherapy (NAC) and radical cystectomy (RC) (n = 223). Unsupervised consensus clustering of highly variant lncRNAs identified a four-cluster solution, which was characterized using a panel of MIBC biomarkers, regulon activity profiles, gene signatures, and survival analysis. The four-cluster solution was confirmed in The Cancer Genome Atlas (TCGA) cohort (n = 405). A single-sample genomic classifier (GC) was trained using ridge-penalized logistic regression and validated in two independent cohorts (n = 255 and n = 94). RESULTS:NAC and TCGA cohorts both contained an lncRNA cluster (LC3) with favorable prognosis that was enriched with tumors of the luminal-papillary (LP) subtype. In both cohorts, patients with LP tumors in LC3 (LPL-C3) were younger and had organ-confined, node-negative disease. The LPL-C3 tumors had enhanced FGFR3, SHH, and wild-type p53 pathway activity. In the TCGA cohort, LPL-C3 tumors were enriched for FGFR3 mutations and depleted for TP53 and RB1 mutations. A GC trained to identify these LPL-C3 patients showed robust performance in two validation cohorts. CONCLUSIONS:Using lncRNA expression profiles, we identified a biologically distinct subgroup of luminal-papillary MIBC with a favorable prognosis. These data suggest that lncRNAs provide additional information for higher-resolution subtyping, potentially improving precision patient management.

BACKGROUND:Several studies have mined short-read RNA sequencing datasets to identify long non-coding RNAs (lncRNAs), and others have focused on the function of individual lncRNAs in abiotic stress response. However, our understanding of the complement, function and origin of lncRNAs - and especially transposon derived lncRNAs (TE-lncRNAs) - in response to abiotic stress is still in its infancy. RESULTS:We utilized a dataset of 127 RNA sequencing samples that included total RNA datasets and PacBio fl-cDNA data to discover lncRNAs in maize. Overall, we identified 23,309 candidate lncRNAs from polyA+ and total RNA samples, with a strong discovery bias within total RNA. The majority (65%) of the 23,309 lncRNAs had sequence similarity to transposable elements (TEs). Most had similarity to long-terminal-repeat retrotransposons from the Copia and Gypsy superfamilies, reflecting a high proportion of these elements in the genome. However, DNA transposons were enriched for lncRNAs relative to their genomic representation by ~ 2-fold. By assessing the fraction of lncRNAs that respond to abiotic stresses like heat, cold, salt and drought, we identified 1077 differentially expressed lncRNA transcripts, including 509 TE-lncRNAs. In general, the expression of these lncRNAs was significantly correlated with their nearest gene. By inferring co-expression networks across our large dataset, we found that 39 lncRNAs are as major hubs in co-expression networks that respond to abiotic stress, and 18 appear to be derived from TEs. CONCLUSIONS:Our results show that lncRNAs are enriched in total RNA samples, that most (65%) are derived from TEs, that at least 1077 are differentially expressed during abiotic stress, and that 39 are hubs in co-expression networks, including a small number that are evolutionary conserved. These results suggest that lncRNAs, including TE-lncRNAs, may play key regulatory roles in moderating abiotic responses.

BACKGROUND:Several studies have mined short-read RNA sequencing datasets to identify long non-coding RNAs (lncRNAs), and others have focused on the function of individual lncRNAs in abiotic stress response. However, our understanding of the complement, function and origin of lncRNAs - and especially transposon derived lncRNAs (TE-lncRNAs) - in response to abiotic stress is still in its infancy. RESULTS:We utilized a dataset of 127 RNA sequencing samples that included total RNA datasets and PacBio fl-cDNA data to discover lncRNAs in maize. Overall, we identified 23,309 candidate lncRNAs from polyA+ and total RNA samples, with a strong discovery bias within total RNA. The majority (65%) of the 23,309 lncRNAs had sequence similarity to transposable elements (TEs). Most had similarity to long-terminal-repeat retrotransposons from the Copia and Gypsy superfamilies, reflecting a high proportion of these elements in the genome. However, DNA transposons were enriched for lncRNAs relative to their genomic representation by ~ 2-fold. By assessing the fraction of lncRNAs that respond to abiotic stresses like heat, cold, salt and drought, we identified 1077 differentially expressed lncRNA transcripts, including 509 TE-lncRNAs. In general, the expression of these lncRNAs was significantly correlated with their nearest gene. By inferring co-expression networks across our large dataset, we found that 39 lncRNAs are as major hubs in co-expression networks that respond to abiotic stress, and 18 appear to be derived from TEs. CONCLUSIONS:Our results show that lncRNAs are enriched in total RNA samples, that most (65%) are derived from TEs, that at least 1077 are differentially expressed during abiotic stress, and that 39 are hubs in co-expression networks, including a small number that are evolutionary conserved. These results suggest that lncRNAs, including TE-lncRNAs, may play key regulatory roles in moderating abiotic responses.

BACKGROUND:Muscle-invasive bladder cancer (MIBC) is a heterogeneous disease, and gene expression profiling has identified several molecular subtypes with distinct biological and clinicopathological characteristics. While MIBC subtyping has primarily been based on messenger RNA (mRNA), long non-coding RNAs (lncRNAs) may provide additional resolution. METHODS:LncRNA expression was quantified from microarray data of a MIBC cohort treated with neoadjuvant chemotherapy (NAC) and radical cystectomy (RC) (n = 223). Unsupervised consensus clustering of highly variant lncRNAs identified a four-cluster solution, which was characterized using a panel of MIBC biomarkers, regulon activity profiles, gene signatures, and survival analysis. The four-cluster solution was confirmed in The Cancer Genome Atlas (TCGA) cohort (n = 405). A single-sample genomic classifier (GC) was trained using ridge-penalized logistic regression and validated in two independent cohorts (n = 255 and n = 94). RESULTS:NAC and TCGA cohorts both contained an lncRNA cluster (LC3) with favorable prognosis that was enriched with tumors of the luminal-papillary (LP) subtype. In both cohorts, patients with LP tumors in LC3 (LPL-C3) were younger and had organ-confined, node-negative disease. The LPL-C3 tumors had enhanced FGFR3, SHH, and wild-type p53 pathway activity. In the TCGA cohort, LPL-C3 tumors were enriched for FGFR3 mutations and depleted for TP53 and RB1 mutations. A GC trained to identify these LPL-C3 patients showed robust performance in two validation cohorts. CONCLUSIONS:Using lncRNA expression profiles, we identified a biologically distinct subgroup of luminal-papillary MIBC with a favorable prognosis. These data suggest that lncRNAs provide additional information for higher-resolution subtyping, potentially improving precision patient management.

Growing evidence suggests that long non-coding RNAs NORAD and miR-205 play a significant role in regulating cancer progression and metastasis. In this study, high expression of NORAD was firstly observed in melanoma tissues and human malignant melanoma cell lines, our aim was to study the interaction of them in the process of invasion and migration of malignant melanoma cells. NORAD, miR-205, and EGLN2 mRNA level in MM cells was detected by qRT-PCR. In situ hybridization (ISH) was performed to detect NORAD expression in MM tissues specimens. Effects of NORAD and miR-205 on Prolyl hydroxylase 2 (EGLN2) expression was explored by western blot in MM cells line. Dual-luciferase reporter assay was performed to verify the interaction relationship between NORAD and miR-205, as well as, miR-205 and EGLN2. Transwell assay was conducted to explore the effects of NORAD and miR-205 in vitro. Xenografts in nude mice experiment were used to confirm the role of NORAD and miR-205 in vivo. In vitro, NORAD knockdown significantly inhibited migration and invasion of malignant melanoma cells and elevated the expression of miR-205, there was an interaction between miR-205 and NORAD in the RNA-induced silencing complex. Upregulation of miR-205 induced significant inhibition of migratory and invasive ability compared with the scrambled control. However, downregulating NORAD largely reversed this effect. Furthermore, the regulatory effects of miR-205 on EGLN2 levels and the induction of endoplasmic reticulum stress were reversed by NORAD. In vivo, deletion of miR-205 induced tumor growth in nude mice. NORAD may play critical roles in tumorigenesis and progression of malignant melanoma by regulating of the miR-205-EGLN2 pathway, and may serve as a new therapeutic target.

Growing evidence suggests that long non-coding RNAs NORAD and miR-205 play a significant role in regulating cancer progression and metastasis. In this study, high expression of NORAD was firstly observed in melanoma tissues and human malignant melanoma cell lines, our aim was to study the interaction of them in the process of invasion and migration of malignant melanoma cells. NORAD, miR-205, and EGLN2 mRNA level in MM cells was detected by qRT-PCR. In situ hybridization (ISH) was performed to detect NORAD expression in MM tissues specimens. Effects of NORAD and miR-205 on Prolyl hydroxylase 2 (EGLN2) expression was explored by western blot in MM cells line. Dual-luciferase reporter assay was performed to verify the interaction relationship between NORAD and miR-205, as well as, miR-205 and EGLN2. Transwell assay was conducted to explore the effects of NORAD and miR-205 in vitro. Xenografts in nude mice experiment were used to confirm the role of NORAD and miR-205 in vivo. In vitro, NORAD knockdown significantly inhibited migration and invasion of malignant melanoma cells and elevated the expression of miR-205, there was an interaction between miR-205 and NORAD in the RNA-induced silencing complex. Upregulation of miR-205 induced significant inhibition of migratory and invasive ability compared with the scrambled control. However, downregulating NORAD largely reversed this effect. Furthermore, the regulatory effects of miR-205 on EGLN2 levels and the induction of endoplasmic reticulum stress were reversed by NORAD. In vivo, deletion of miR-205 induced tumor growth in nude mice. NORAD may play critical roles in tumorigenesis and progression of malignant melanoma by regulating of the miR-205-EGLN2 pathway, and may serve as a new therapeutic target.

Non-coding RNA (ncRNA) has been shown to regulate diverse cellular processes and functions through controlling gene expression. Long non-coding RNAs (lncRNAs) act as a competing endogenous RNAs (ceRNAs) where microRNAs (miRNAs) and lncRNAs regulate each other through their biding sites. Interactions of miRNAs and lncRNAs have been reported to trigger decay of the targeted lncRNAs and have important roles in target gene regulation. These interactions form complicated and intertwined networks. Certain lncRNAs encode miRNAs and small nucleolar RNAs (snoRNAs), and may regulate expression of these small RNAs as precursors. SnoRNAs have also been reported to be precursors for PIWI-interacting RNAs (piRNAs) and thus may regulate the piRNAs as a precursor. These miRNAs and piRNAs target messenger RNAs (mRNAs) and regulate gene expression. In this review, we will present and discuss these interactions, cross-talk, and co-regulation of ncRNAs and gene regulation due to these interactions.

Non-coding RNA (ncRNA) has been shown to regulate diverse cellular processes and functions through controlling gene expression. Long non-coding RNAs (lncRNAs) act as a competing endogenous RNAs (ceRNAs) where microRNAs (miRNAs) and lncRNAs regulate each other through their biding sites. Interactions of miRNAs and lncRNAs have been reported to trigger decay of the targeted lncRNAs and have important roles in target gene regulation. These interactions form complicated and intertwined networks. Certain lncRNAs encode miRNAs and small nucleolar RNAs (snoRNAs), and may regulate expression of these small RNAs as precursors. SnoRNAs have also been reported to be precursors for PIWI-interacting RNAs (piRNAs) and thus may regulate the piRNAs as a precursor. These miRNAs and piRNAs target messenger RNAs (mRNAs) and regulate gene expression. In this review, we will present and discuss these interactions, cross-talk, and co-regulation of ncRNAs and gene regulation due to these interactions.

Esophageal squamous cell carcinoma (ESCC) is one of the most lethal malignant cancers with high incidence and mortality. Current reliable effective diagnostic and prognostic biomarkers are very limited in clinic. Emerging evidence indicates that dysregulated expression of the long non-coding RNAs (lncRNAs) was examined in various types of cancer including ESCC. ESCC associated lncRNA _1 (ESCCAL_1) was first time identified to be increased expression in ESCC, and therefore named by our research team. However, its potential function in the progression of ESCC remains unclear. In this study, we investigated the effect of ESCCAL_1 knockdown on ESCC tumorigenicity using a xenograft mouse model and explored the underlying molecular mechanism. Here we showed that ESCCAL_1 knockdown significantly inhibited EC9706 cell growth in nude mice. Interestingly, we also found that reduced expression of ESCCAL_1 resulted in distinct alterations of relative phosphorylation level of kinases (p-p38α, p-JNK, p-FAK and p-Src), and significant changes of the expression level of apoptosis-related proteins (p53, BAX, Bcl-2 and Caspase-3). In summary, our results suggest that lncRNA ESCCAL_1 is a potential diagnostic and prognostic target of ESCC.

BackgroundPolycomb Repressive Complex 2 (PRC2) catalyzes histone methylation at H3 Lys27, and plays crucial roles during development and diseases in numerous systems. Its catalytic subunit EZH2 represents a key nuclear target for long non-coding RNAs (lncRNAs) that emerging to be a novel class of epigenetic regulator and participate in diverse cellular processes. LncRNAs are characterized by high tissue-specificity; however, little is known about the tissue profile of the EZH2- interacting lncRNAs.ObjectiveHere we performed a global screening for EZH2-binding lncRNAs in tissues including brain, lung, heart, liver, kidney, intestine, spleen, testis, muscle and blood by combining RNA immuno- precipitation and RNA sequencing. We identified 1328 EZH2-binding lncRNAs, among which 470 were shared in at least two tissues while 858 were only detected in single tissue. An RNA motif with specific secondary structure was identified in a number of lncRNAs, albeit not in all EZH2-binding lncRNAs. The EZH2-binding lncRNAs fell into four categories including intergenic lncRNA, antisense lncRNA, intron-related lncRNA and promoter-related lncRNA, suggesting diverse regulations of both cis and trans-mechanisms. A promoter-related lncRNA Hnf1aos1 bound to EZH2 specifically in the liver, a feature same as its paired coding gene Hnf1a, further confirming the validity of our study. In addition to the well known EZH2-binding lncRNAs like Kcnq1ot1, Gas5, Meg3, Hotair and Malat1, majority of the lncRNAs were firstly reported to be associated with EZH2.ConclusionOur findings provide a profiling view of the EZH2-interacting lncRNAs across different tissues, and suggest critical roles of lncRNAs during cell differentiation and maturation.

BackgroundPolycomb Repressive Complex 2 (PRC2) catalyzes histone methylation at H3 Lys27, and plays crucial roles during development and diseases in numerous systems. Its catalytic subunit EZH2 represents a key nuclear target for long non-coding RNAs (lncRNAs) that emerging to be a novel class of epigenetic regulator and participate in diverse cellular processes. LncRNAs are characterized by high tissue-specificity; however, little is known about the tissue profile of the EZH2- interacting lncRNAs.ObjectiveHere we performed a global screening for EZH2-binding lncRNAs in tissues including brain, lung, heart, liver, kidney, intestine, spleen, testis, muscle and blood by combining RNA immuno- precipitation and RNA sequencing. We identified 1328 EZH2-binding lncRNAs, among which 470 were shared in at least two tissues while 858 were only detected in single tissue. An RNA motif with specific secondary structure was identified in a number of lncRNAs, albeit not in all EZH2-binding lncRNAs. The EZH2-binding lncRNAs fell into four categories including intergenic lncRNA, antisense lncRNA, intron-related lncRNA and promoter-related lncRNA, suggesting diverse regulations of both cis and trans-mechanisms. A promoter-related lncRNA Hnf1aos1 bound to EZH2 specifically in the liver, a feature same as its paired coding gene Hnf1a, further confirming the validity of our study. In addition to the well known EZH2-binding lncRNAs like Kcnq1ot1, Gas5, Meg3, Hotair and Malat1, majority of the lncRNAs were firstly reported to be associated with EZH2.ConclusionOur findings provide a profiling view of the EZH2-interacting lncRNAs across different tissues, and suggest critical roles of lncRNAs during cell differentiation and maturation.

The innate immune system protects against infections by initiating an inducible inflammatory response. NF-κB is one of the critical transcription factors controlling this complex response, but some aspects of its regulation remain unclear. For example, although long non-coding RNAs (lncRNAs) have been shown to critically regulate gene expression, only a fraction of these have been functionally characterized, and the extent to which lncRNAs control NF-κB expression is unknown. Here, we describe the generation of a GFP-based NF-κB reporter system in immortalized murine bone marrow-derived macrophages (iBMDM). Activation of this reporter, using Toll-like receptor ligands, resulted in GFP expression, which could be monitored by flow cytometry. We also established a CRISPR/Cas9 gene deletion system in this NF-κB reporter line, enabling us to screen for genes that regulate NF-κB signaling. Our deletion-based approach identified two long intergenic non-coding(linc)RNAs, lincRNA-Cox2 and lincRNA-AK170409, that control NF-κB signaling. We demonstrate a potential novel role for lincRNA-Cox2 in promoting IκBα degradation in the cytoplasm. For lincRNA-AK170409, we provide evidence that this nuclearly-localized lincRNA regulates a number of inflammation-related genes. In conclusion, we have established an NF-κB-GFP iBMDM reporter cell line and a line that stably expresses Cas9. Our approach enabled the identification of lincRNA-Cox2 and lincRNA-AK170409 as NF-κB regulators, and this tool will be useful for identifying additional genes involved in regulating this transcription factor critical for immune function.

The innate immune system protects against infections by initiating an inducible inflammatory response. NF-κB is one of the critical transcription factors controlling this complex response, but some aspects of its regulation remain unclear. For example, although long non-coding RNAs (lncRNAs) have been shown to critically regulate gene expression, only a fraction of these have been functionally characterized, and the extent to which lncRNAs control NF-κB expression is unknown. Here, we describe the generation of a GFP-based NF-κB reporter system in immortalized murine bone marrow-derived macrophages (iBMDM). Activation of this reporter, using Toll-like receptor ligands, resulted in GFP expression, which could be monitored by flow cytometry. We also established a CRISPR/Cas9 gene deletion system in this NF-κB reporter line, enabling us to screen for genes that regulate NF-κB signaling. Our deletion-based approach identified two long intergenic non-coding(linc)RNAs, lincRNA-Cox2 and lincRNA-AK170409, that control NF-κB signaling. We demonstrate a potential novel role for lincRNA-Cox2 in promoting IκBα degradation in the cytoplasm. For lincRNA-AK170409, we provide evidence that this nuclearly-localized lincRNA regulates a number of inflammation-related genes. In conclusion, we have established an NF-κB-GFP iBMDM reporter cell line and a line that stably expresses Cas9. Our approach enabled the identification of lincRNA-Cox2 and lincRNA-AK170409 as NF-κB regulators, and this tool will be useful for identifying additional genes involved in regulating this transcription factor critical for immune function.

BackgroundEfforts to resolve the transcribed sequences in the equine genome have focused on protein-coding RNA. The transcription of the intergenic regions, although detected via total RNA sequencing (RNA-seq), has yet to be characterized in the horse. The most recent equine transcriptome based on RNA-seq from several tissues was a prime opportunity to obtain a concurrent long non-coding RNA (lncRNA) database.ResultsThis lncRNA database has a breadth of eight tissues and a depth of over 20 million reads for select tissues, providing the deepest and most expansive equine lncRNA database. Utilizing the intergenic reads and three categories of novel genes from a previously published equine transcriptome pipeline, we better describe these groups by annotating the lncRNA candidates. These lncRNA candidates were filtered using an approach adapted from human lncRNA annotation, which removes transcripts based on size, expression, protein-coding capability and distance to the start or stop of annotated protein-coding transcripts.ConclusionOur equine lncRNA database has 20,800 transcripts that demonstrate characteristics unique to lncRNA including low expression, low exon diversity and low levels of sequence conservation. These candidate lncRNA will serve as a baseline lncRNA annotation and begin to describe the RNA-seq reads assigned to the intergenic space in the horse.

Silibinin is the major active constituent of silymarin, an extract of milk thistle seeds. Silibinin has been shown to have significant anti-cancer effects in a variety of malignancies. However, the molecular mechanisms of silibinin action in bladder cancer have not been studied extensively. In the present study, we found that silibinin (10 μM) significantly suppressed proliferation, migration, invasion and induced apoptosis of T24 and UM-UC-3 human bladder cancer cells. Silibinin down-regulated the actin cytoskeleton and phosphatidylinositide 3-kinase (PI3K)/Akt signaling pathways in these cancer cell lines. These pathways were found to crosstalk through RAS cascades. We found that silibinin suppressed levels of trimethylated histone H3 lysine 4 and acetylated H3 at the KRAS promoter. Furthermore, silibinin targets long non-coding RNA: HOTAIR and ZFAS1, which are known to play roles as oncogenic factors in various cancers. This study shows that silibinin exerts anti-cancer effects through down-regulation of actin cytoskeleton and PI3K/Akt pathways and thus suppresses bladder cancer growth and progression.

Silibinin is the major active constituent of silymarin, an extract of milk thistle seeds. Silibinin has been shown to have significant anti-cancer effects in a variety of malignancies. However, the molecular mechanisms of silibinin action in bladder cancer have not been studied extensively. In the present study, we found that silibinin (10 μM) significantly suppressed proliferation, migration, invasion and induced apoptosis of T24 and UM-UC-3 human bladder cancer cells. Silibinin down-regulated the actin cytoskeleton and phosphatidylinositide 3-kinase (PI3K)/Akt signaling pathways in these cancer cell lines. These pathways were found to crosstalk through RAS cascades. We found that silibinin suppressed levels of trimethylated histone H3 lysine 4 and acetylated H3 at the KRAS promoter. Furthermore, silibinin targets long non-coding RNA: HOTAIR and ZFAS1, which are known to play roles as oncogenic factors in various cancers. This study shows that silibinin exerts anti-cancer effects through down-regulation of actin cytoskeleton and PI3K/Akt pathways and thus suppresses bladder cancer growth and progression.

BackgroundEfforts to resolve the transcribed sequences in the equine genome have focused on protein-coding RNA. The transcription of the intergenic regions, although detected via total RNA sequencing (RNA-seq), has yet to be characterized in the horse. The most recent equine transcriptome based on RNA-seq from several tissues was a prime opportunity to obtain a concurrent long non-coding RNA (lncRNA) database.ResultsThis lncRNA database has a breadth of eight tissues and a depth of over 20 million reads for select tissues, providing the deepest and most expansive equine lncRNA database. Utilizing the intergenic reads and three categories of novel genes from a previously published equine transcriptome pipeline, we better describe these groups by annotating the lncRNA candidates. These lncRNA candidates were filtered using an approach adapted from human lncRNA annotation, which removes transcripts based on size, expression, protein-coding capability and distance to the start or stop of annotated protein-coding transcripts.ConclusionOur equine lncRNA database has 20,800 transcripts that demonstrate characteristics unique to lncRNA including low expression, low exon diversity and low levels of sequence conservation. These candidate lncRNA will serve as a baseline lncRNA annotation and begin to describe the RNA-seq reads assigned to the intergenic space in the horse.

Execution of pluripotency requires progression from the naïve status represented by mouse embryonic stem cells (ESCs) to a state capacitated for lineage specification. This transition is coordinated at multiple levels. Non-coding RNAs may contribute to this regulatory orchestra. We identified a rodent-specific long non-coding RNA (lncRNA) linc1281, hereafter Ephemeron (Eprn), that modulates the dynamics of exit from naïve pluripotency. Eprn deletion delays the extinction of ESC identity, an effect associated with perduring Nanog expression. In the absence of Eprn, Lin28a expression is reduced which results in persistence of let-7 microRNAs, and the up-regulation of de novo methyltransferases Dnmt3a/b is delayed. Dnmt3a/b deletion retards ES cell transition, correlating with delayed Nanog promoter methylation and phenocopying loss of Eprn or Lin28a. The connection from lncRNA to miRNA and DNA methylation facilitates the acute extinction of naïve pluripotency, a pre-requisite for rapid progression from preimplantation epiblast to gastrulation in rodents. Eprn illustrates how lncRNAs may introduce species-specific network modulations.

Long non-coding RNAs (lncRNAs) play an important role in gene regulation and are involving in diverse cellular processes. However, their roles in reprogramming of gene expression profiles during lineage commitment and maturation of mesenchymal stem cells (MSCs) remain poorly understood. In the current study, we characterize the expression of a lncRNA, HoxA-AS3, during the differentiation of MSCs. We showed that HoxA-AS3 is increased upon adipogenic induction of MSCs, while HoxA-AS3 remains unaltered during osteogenic induction. Silencing of HoxA-AS3 in MSCs resulted in decreased adipogenesis and expression of adipogenic markers, PPARG, CEBPA, FABP4 and ADIPOQ. Conversely, knockdown of HoxA-AS3 expression in MSCs exhibited an enhanced osteogenesis and osteogenic markers expression, including RUNX2, SP7, COL1A1, IBSP, BGLAP and SPP1. Mechanistically, HoxA-AS3 interacts with Enhancer Of Zeste 2 (EZH2) and is required for H3 lysine-27 trimethylation (H3K27me3) of key osteogenic transcription factor Runx2. Our data reveal that HoxA-AS3 acts as an epigenetic switch that determines the lineage specification of MSC.

Long non-coding RNAs (lncRNAs) play an important role in gene regulation and are involving in diverse cellular processes. However, their roles in reprogramming of gene expression profiles during lineage commitment and maturation of mesenchymal stem cells (MSCs) remain poorly understood. In the current study, we characterize the expression of a lncRNA, HoxA-AS3, during the differentiation of MSCs. We showed that HoxA-AS3 is increased upon adipogenic induction of MSCs, while HoxA-AS3 remains unaltered during osteogenic induction. Silencing of HoxA-AS3 in MSCs resulted in decreased adipogenesis and expression of adipogenic markers, PPARG, CEBPA, FABP4 and ADIPOQ. Conversely, knockdown of HoxA-AS3 expression in MSCs exhibited an enhanced osteogenesis and osteogenic markers expression, including RUNX2, SP7, COL1A1, IBSP, BGLAP and SPP1. Mechanistically, HoxA-AS3 interacts with Enhancer Of Zeste 2 (EZH2) and is required for H3 lysine-27 trimethylation (H3K27me3) of key osteogenic transcription factor Runx2. Our data reveal that HoxA-AS3 acts as an epigenetic switch that determines the lineage specification of MSC.

ObjectivesThe role of long non-coding RNA (lncRNA) expression in human head and neck squamous cell carcinoma (HNSCC) is still poorly understood. In this study, we aimed at establishing the onco-lncRNAome profiling of HNSCC and to identify lncRNAs correlating with prognosis and patient survival.Materials and methodsThe Atlas of Noncoding RNAs in Cancer (TANRIC) database was employed to retrieve the lncRNA expression information generated from The Cancer Genome Atlas (TCGA) HNSCC RNA-sequencing data. RNA-sequencing data from HNSCC cell lines were also considered for this study. Bioinformatics approaches, such as differential gene expression analysis, survival analysis, principal component analysis, and Co-LncRNA enrichment analysis were performed.ResultsUsing TCGA HNSCC RNA-sequencing data from 426 HNSCC and 42 adjacent normal tissues, we found 728 lncRNA transcripts significantly and differentially expressed in HNSCC. Among the 728 lncRNAs, 55 lncRNAs were significantly associated with poor prognosis, such as overall survival and/or disease-free survival. Next, we found 140 lncRNA transcripts significantly and differentially expressed between Human Papilloma Virus (HPV) positive tumors and HPV negative tumors. Thirty lncRNA transcripts were differentially expressed between TP53 mutated and TP53 wild type tumors. Co-LncRNA analysis suggested that protein-coding genes that are co-expressed with these deregulated lncRNAs might be involved in cancer associated molecular events. With consideration of differential expression of lncRNAs in a HNSCC cell lines panel (n=22), we found several lncRNAs that may represent potential targets for diagnosis, therapy and prevention of HNSCC.ConclusionLncRNAs profiling could provide novel insights into the potential mechanisms of HNSCC oncogenesis.

FBXW7, an E3-ubiquitin protein ligase in SCFs (SKP1-cullin-F-box) complex, is a major human tumor suppressor gene, and understanding mechanisms by which FBXW7 contributes to tumorigenesis is critical for the treatment of human cancers with FBXW7 deficiency. Long non-coding RNAs (lncRNAs) have emerged as key regulators of various biological processes. Here we have identified a set of lncRNAs that are associated with Fbxw7 deficiency. The correlation network and functional annotation analysis revealed that Fbxw7-associated lncRNAs regulate genes involved in cell cycle, DNA repair, metabolic process, and cell communication and adhesion. The number of coding genes that correlated with individual lncRNAs varied largely. A lncRNA on chromosome 15 (A_30_P01032978), which was upregulated in tumors from Fbxw7 deficient mice was positively correlated with 15 coding genes. High expression of this 15-gene signature was associated with poor prognosis in two independent human breast cancer studies. Our results open possible new avenues to understand mechanisms by which Fbxw7 deficiency increases tumor susceptibility via the alteration of lncRNAs.

FBXW7, an E3-ubiquitin protein ligase in SCFs (SKP1-cullin-F-box) complex, is a major human tumor suppressor gene, and understanding mechanisms by which FBXW7 contributes to tumorigenesis is critical for the treatment of human cancers with FBXW7 deficiency. Long non-coding RNAs (lncRNAs) have emerged as key regulators of various biological processes. Here we have identified a set of lncRNAs that are associated with Fbxw7 deficiency. The correlation network and functional annotation analysis revealed that Fbxw7-associated lncRNAs regulate genes involved in cell cycle, DNA repair, metabolic process, and cell communication and adhesion. The number of coding genes that correlated with individual lncRNAs varied largely. A lncRNA on chromosome 15 (A_30_P01032978), which was upregulated in tumors from Fbxw7 deficient mice was positively correlated with 15 coding genes. High expression of this 15-gene signature was associated with poor prognosis in two independent human breast cancer studies. Our results open possible new avenues to understand mechanisms by which Fbxw7 deficiency increases tumor susceptibility via the alteration of lncRNAs.

Adult-onset autosomal dominant leukodystrophy (ADLD) is a progressive and fatal hereditary demyelination disorder characterized initially by autonomic dysfunction and loss of myelin in the central nervous system (CNS). Majority of ADLD is caused by a genomic duplication of the nuclear lamin B1 gene (LMNB1) encoding lamin B1 protein, resulting in increased gene dosage in brain tissue. In vitro, excessive lamin B1 at the cellular level reduces transcription of myelin genes, leading to premature arrest of oligodendrocyte differentiation. Murine models of ADLD overexpressing LMNB1 exhibited age-dependent motor deficits and myelin defects, which are associated with reduced occupancy of the Yin Yang 1 transcription factor at the promoter region of the proteolipid protein gene. Lamin B1 overexpression mediates oligodendrocyte cell-autonomous neuropathology in ADLD and suggests lamin B1 as an important regulator of myelin formation and maintenance during aging. Identification of microRNA-23 (miR-23) as a negative regulator of lamin B1 can ameliorate the consequences of excessive lamin B1 at the cellular level. miR-23a-overexpressing mice display enhanced oligodendrocyte differentiation and myelin synthesis. miR-23a targets include a protein coding transcript PTEN (phosphatase and tensin homolog on chromosome 10), and a long noncoding RNA (2700046G09Rik), indicating a unique role for miR-23a in the coordination of proteins and noncoding RNAs in generating and maintaining healthy myelin. Here, we provide a concise review of the current literature on clinical presentations of ADLD and how lamin B1 affects myelination and other developmental processes. Moreover, we address the emerging role of non-coding RNAs (ncRNAs) in modulating gene networks, specifically investigating miR-23 as a potential target for the treatment of ADLD and other demyelinating disorders.

Kaposi's sarcoma-associated herpesvirus (KSHV/human herpesvirus 8) is a γ-herpesvirus linked to Kaposi's sarcoma (KS) and two lymphoproliferative disorders, primary effusion lymphoma (PEL or body-cavity B-lymphoma [BCBL]) and a subset of Multicentric Castleman's Disease. During lytic growth, pervasive viral transcription generating a variety of transcripts with uncertain protein-coding potential has been described on a genome-wide scale in β- and γ-herpesviruses. One class of such RNAs is called long non-coding RNAs (lncRNAs). KSHV encodes a viral lncRNA known as polyadenylated nuclear RNA (PAN RNA), a copious early gene product. PAN RNA has been implicated in KSHV gene expression, replication, and immune modulation. PAN RNA expression is required for optimal expression of the entire KSHV lytic gene expression program. Latent KSHV episomes are coated with viral latency-associated nuclear antigen (LANA). LANA rapidly dissociates from episomes during reactivation. Here we review recent studies suggesting that PAN RNA may function as a viral lncRNA, including a role in the facilitation of LANA-episomal dissociation during lytic replication.

Kaposi's sarcoma-associated herpesvirus (KSHV/human herpesvirus 8) is a γ-herpesvirus linked to Kaposi's sarcoma (KS) and two lymphoproliferative disorders, primary effusion lymphoma (PEL or body-cavity B-lymphoma [BCBL]) and a subset of Multicentric Castleman's Disease. During lytic growth, pervasive viral transcription generating a variety of transcripts with uncertain protein-coding potential has been described on a genome-wide scale in β- and γ-herpesviruses. One class of such RNAs is called long non-coding RNAs (lncRNAs). KSHV encodes a viral lncRNA known as polyadenylated nuclear RNA (PAN RNA), a copious early gene product. PAN RNA has been implicated in KSHV gene expression, replication, and immune modulation. PAN RNA expression is required for optimal expression of the entire KSHV lytic gene expression program. Latent KSHV episomes are coated with viral latency-associated nuclear antigen (LANA). LANA rapidly dissociates from episomes during reactivation. Here we review recent studies suggesting that PAN RNA may function as a viral lncRNA, including a role in the facilitation of LANA-episomal dissociation during lytic replication.

The human transcriptome is complex, comprising multiple transcript types, mostly in the form of non-coding RNA (ncRNA). The majority of ncRNA is of the long form (lncRNA, ≥ 200 bp), which plays an important role in gene regulation through multiple mechanisms including epigenetics, chromatin modification, control of transcription factor binding, and regulation of alternative splicing. Both mRNA and ncRNA exhibit additional variability in the form of alternative splicing and RNA editing. All aspects of the human transcriptome can potentially be dysregulated by environmental exposures. Next-generation RNA sequencing (RNA-Seq) is the best available methodology to measure this although it has limitations, including experimental bias. The third phase of the MicroArray Quality Control Consortium project (MAQC-III), also called Sequencing Quality Control (SeQC), aims to address these limitations through standardization of experimental and bioinformatic methodologies. A limited number of toxicogenomic studies have been conducted to date using RNA-Seq. This review describes the complexity of the human transcriptome, the application of transcriptomics by RNA-Seq or microarray in molecular epidemiology studies, and limitations of these approaches including the type of cell or tissue analyzed, experimental variation, and confounding. By using good study designs with precise, individual exposure measurements, sufficient power and incorporation of phenotypic anchors, studies in human populations can identify biomarkers of exposure and/or early effect and elucidate mechanisms of action underlying associated diseases, even at low doses. Analysis of datasets at the pathway level can compensate for some of the limitations of RNA-Seq and, as more datasets become available, will increasingly elucidate the exposure-disease continuum.

Kaposi sarcoma-associated herpesvirus (KSHV) ORF57, also known as Mta (mRNA transcript accumulation), enhances viral intron-less transcript accumulation and promotes splicing of intron-containing viral RNA transcripts. In this study, we identified KSHV PAN, a long non-coding polyadenylated nuclear RNA as a main target of ORF57 by a genome-wide CLIP (cross-linking and immunoprecipitation) approach. KSHV genome lacking ORF57 expresses only a minimal amount of PAN. In cotransfection experiments, ORF57 alone increased PAN expression by 20-30-fold when compared to vector control. This accumulation function of ORF57 was dependent on a structured RNA element in the 5' PAN, named MRE (Mta responsive element), but not much so on an ENE (expression and nuclear retention element) in the 3' PAN previously reported by other studies. We showed that the major function of the 5' PAN MRE is increasing the RNA half-life of PAN in the presence of ORF57. Further mutational analyses revealed a core motif consisting of 9 nucleotides in the MRE-II , which is responsible for ORF57 interaction and function. The 9-nt core in the MRE-II also binds cellular PABPC1, but not the E1B-AP5 which binds another region of the MRE-II. In addition, we found that PAN RNA is partially exportable in the presence of ORF57. Together, our data provide compelling evidence as to how ORF57 functions to accumulate a non-coding viral RNA in the course of virus lytic infection.

Kaposi sarcoma-associated herpesvirus (KSHV) ORF57, also known as Mta (mRNA transcript accumulation), enhances viral intron-less transcript accumulation and promotes splicing of intron-containing viral RNA transcripts. In this study, we identified KSHV PAN, a long non-coding polyadenylated nuclear RNA as a main target of ORF57 by a genome-wide CLIP (cross-linking and immunoprecipitation) approach. KSHV genome lacking ORF57 expresses only a minimal amount of PAN. In cotransfection experiments, ORF57 alone increased PAN expression by 20-30-fold when compared to vector control. This accumulation function of ORF57 was dependent on a structured RNA element in the 5' PAN, named MRE (Mta responsive element), but not much so on an ENE (expression and nuclear retention element) in the 3' PAN previously reported by other studies. We showed that the major function of the 5' PAN MRE is increasing the RNA half-life of PAN in the presence of ORF57. Further mutational analyses revealed a core motif consisting of 9 nucleotides in the MRE-II , which is responsible for ORF57 interaction and function. The 9-nt core in the MRE-II also binds cellular PABPC1, but not the E1B-AP5 which binds another region of the MRE-II. In addition, we found that PAN RNA is partially exportable in the presence of ORF57. Together, our data provide compelling evidence as to how ORF57 functions to accumulate a non-coding viral RNA in the course of virus lytic infection.