Regulation of the cohesin-loading factor NIPBL: Role of the lncRNA NIPBL-AS1 and identification of a distal enhancer element

Identification of Human Pathways Acting on Nuclear Non-Coding RNAs Using the Mirror Forward Genetic Approach

Despite critical roles in diseases, human pathways acting on strictly nuclear non-coding RNAs have been refractory to forward genetics. To enable their forward genetic discovery, we developed a single-cell approach that "Mirrors" activities of nuclear pathways with cytoplasmic fluorescence. Application of Mirror to two nuclear pathways targeting MALAT1's 3' end, the pathway of its maturation and the other, the degradation pathway blocked by the triple-helical Element for Nuclear Expression (ENE), identified nearly all components of three complexes: Ribonuclease P and the RNA Exosome, including nuclear DIS3, EXOSC10, and C1D, as well as the Nuclear Exosome Targeting (NEXT) complex. Additionally, Mirror identified DEAD-box helicase DDX59 associated with the genetic disorder Oral-Facial-Digital syndrome (OFD), yet lacking known substrates or roles in nuclear RNA degradation. Knockout of DDX59 exhibits stabilization of the full-length MALAT1 with a stability-compromised ENE and increases levels of 3'-extended forms of small nuclear RNAs. It also exhibits extensive retention of minor introns, including in OFD-associated genes, suggesting a mechanism for DDX59 association with OFD. Mirror efficiently identifies pathways acting on strictly nuclear non-coding RNAs, including essential and indirectly-acting components, and, as a result, uncovers unexpected links to human disease.

Molecular models of bidirectional promoter regulation

Approximately 11% of human genes are transcribed by a bidirectional promoter (BDP), defined as two genes with <1 kb between their transcription start sites. Despite their evolutionary conservation and enrichment for housekeeping genes and oncogenes, the regulatory role of BDPs remains unclear. BDPs have been suggested to facilitate gene coregulation and/or decrease expression noise. This review discusses these potential regulatory functions through the context of six prospective underlying mechanistic models: a single nucleosome free region, shared transcription factor/regulator binding, cooperative negative supercoiling, bimodal histone marks, joint activation by enhancer(s), and RNA-mediated recruitment of regulators. These molecular mechanisms may act independently and/or cooperatively to facilitate the coregulation and/or decreased expression noise predicted of BDPs.

Profiling and Bioinformatics Analyses of Hypoxia-Induced Differential Expression of Long Non-coding RNA in Glioblastoma Multiforme Cells

Hypoxic microenvironments are intricately linked to malignant characteristics of glioblastoma multiforme (GBM). Long non-coding ribonucleic acids (lncRNAs) have been reported to be involved in the progression of GBM and closely associated with hypoxia. Nevertheless, the differential expression profiles as well as functional roles of lncRNAs in GBM cells under hypoxic conditions remain largely obscure. We explored the expression profiles of lncRNAs in hypoxic U87 cells as well as T98G cells using sequencing analysis. The effect of differentially expressed lncRNAs (DElncRNAs) was assessed through bioinformatic analysis. Furthermore, the expression of lncRNAs significantly dysregulated in both U87 and T98G cells was further validated using quantitative reverse-transcription polymerase chain reaction (qRT-PCR). Relevant cell functional experiments were also conducted. We used predicted RNA-binding proteins (RBPs) to construct an interaction network via the interaction prediction module. U87 and T98G cells showed dysregulation of 1115 and 597 lncRNAs, respectively. Gene Ontology (GO) analysis indicated that altered lncRNA expression was associated with nucleotide-excision repair and cell metabolism in GBM cells. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed the association between dysregulated lncRNAs and the Hippo signaling pathway under hypoxia. The dysregulation of six selected lncRNAs (ENST00000371192, uc003tnq.3, ENST00000262952, ENST00000609350, ENST00000610036, and NR_046262) was validated by qRT-PCR. Investigation of lncRNA-microRNA (miRNA)-mRNA networks centered on HIF-1α demonstrated cross-talk between the six validated lncRNAs and 16 related miRNAs. Functional experiments showed the significant inhibition of GBM cell proliferation, invasion, and migration by the knockdown of uc003tnq.3 in vitro. Additionally, uc003tnq.3 was used to construct a comprehensive RBP-transcription factor (TF)-miRNA interaction network. The expression of LncRNAs was dysregulated in GBM cells under hypoxic conditions. The identified six lncRNAs might exert important effect on the development of GBM under hypoxic microenvironment.

Live-cell imaging of chromatin contacts opens a new window into chromatin dynamics

Our understanding of the organization of the chromatin fiber within the cell nucleus has made great progress in the last few years. High-resolution techniques based on next-generation sequencing as well as optical imaging that can investigate chromatin conformations down to the single cell level have revealed that chromatin structure is highly heterogeneous at the level of the individual allele. While TAD boundaries and enhancer-promoter pairs emerge as hotspots of 3D proximity, the spatiotemporal dynamics of these different types of chromatin contacts remain largely unexplored. Investigation of chromatin contacts in live single cells is necessary to close this knowledge gap and further enhance the current models of 3D genome organization and enhancer-promoter communication. In this review, we first discuss the potential of single locus labeling to study architectural and enhancer-promoter contacts and provide an overview of the available single locus labeling techniques such as FROS, TALE, CRISPR-dCas9 and ANCHOR, and discuss the latest developments and applications of these systems.

Deep intronic NIPBL de novo mutations and differential diagnoses revealed by whole genome and RNA sequencing in Cornelia de Lange syndrome patients

Cornelia de Lange syndrome (CdLS; MIM# 122470) is a rare developmental disorder. Pathogenic variants in 5 genes explain approximately 50% cases, leaving the other 50% unsolved. We performed whole genome sequencing (WGS) ± RNA sequencing (RNA-seq) in 5 unsolved trios fulfilling the following criteria: (i) clinical diagnosis of classic CdLS, (ii) negative gene panel sequencing from blood and saliva-isolated DNA, (iii) unaffected parents' DNA samples available and (iv) proband's blood-isolated RNA available. A pathogenic de novo mutation (DNM) was observed in a CdLS differential diagnosis gene in 3/5 patients, namely POU3F3, SPEN, and TAF1. In the other two, we identified two distinct deep intronic DNM in NIPBL predicted to create a novel splice site. RT-PCRs and RNA-Seq showed aberrant transcripts leading to the creation of a novel frameshift exon. Our findings suggest the relevance of WGS in unsolved suspected CdLS cases and that deep intronic variants may account for a proportion of them.

Case report: A novel heterozygous synonymous variant in deep exon region of NIPBL gene generating a non-canonical splice donor in a patient with cornelia de lange syndrome

Cornelia de Lange syndrome (CdLS) is an autosomal dominant or X-linked genetic disease with significant genetic heterogeneity. Variants of the NIPBL gene are responsible for CdLS in 60% of patients. Herein, we report the case of a patient with CdLS showing distinctive facial features, microcephaly, developmental delay, and growth retardation. Whole exome sequencing was performed for the patient, and a novel de novo heterozygous synonymous variant was identified in the deep region of exon 40 in the NIPBL gene (NM_133433.4: c. 6819G > T, p. Gly2273 = ). The clinical significance of the variant was uncertain according to the ACMG/AMP guidelines; however, based on in silico analysis, it was predicted to alter mRNA splicing. To validate the prediction, a reverse transcriptase-polymerase chain reaction was conducted. The variant activated a cryptic splice donor, generating a short transcript of NIPBL. A loss of 137 bp at the 3' end of NIPBL exon 40 was detected, which potentially altered the open reading frame by inserting multiple premature termination codons. Quantitative real-time PCR analysis showed that the ratio of the transcription level of the full-length transcript to that of the altered short transcript in the patient was 5:1, instead of 1:1. These findings may explain the relatively mild phenotype of the patient, regardless of the loss of function of the truncated protein due to a frameshift in the mRNA. To the best of our knowledge, this study is the first to report a synonymous variant in the deep exon regions of the NIPBL gene responsible for CdLS. The identified variant expands the mutational spectrum of the NIPBL gene. Furthermore, synonymous variations may be pathogenic, which should not be ignored in the clinical and genetic diagnosis of the disease.

The Cohesin Complex and Its Interplay with Non-Coding RNAs

The cohesin complex is a multi-subunit protein complex initially discovered for its role in sister chromatid cohesion. However, cohesin also has several other functions and plays important roles in transcriptional regulation, DNA double strand break repair, and chromosome architecture thereby influencing gene expression and development in organisms from yeast to man. While most of these functions rely on protein–protein interactions, post-translational protein, as well as DNA modifications, non-coding RNAs are emerging as additional players that facilitate and modulate the function or expression of cohesin and its individual components. This review provides a condensed overview about the architecture as well as the function of the cohesin complex and highlights its multifaceted interplay with both short and long non-coding RNAs.

Identification of co‑expression modules and hub genes of retinoblastoma via co‑expression analysis and protein‑protein interaction networks

Retinoblastoma is a common intraocular malignant tumor in children. However, the molecular and genetic mechanisms of retinoblastoma remain unclear. The gene expression dataset GSE110811 was retrieved from Gene Expression Omnibus. After preprocessing, coexpression modules were constructed by weighted gene coexpression network analysis (WGCNA), and modules associated with clinical traits were identified. In addition, functional enrichment analysis was performed for genes in the indicated modules, and protein-protein interaction (PPI) networks and subnetworks were constructed based on these genes. Eight coexpression modules were constructed through WGCNA. Of these, the yellow module had the highest association with severity and age (r=0.82 and P=3e-07; r=0.72 and P=3e-05). The turquoise module had the highest association with months (r=−0.63 and P=5e-04). The genes in the two modules participate in multiple pathways of retinoblastoma, and by combining the PPI network and subnetworks; 10 hub genes were identified in the two modules. The present study identified coexpression modules and hub genes associated with clinical traits of retinoblastoma, providing novel insight into retinoblastoma progression.

MAU2 and NIPBL Variants Impair the Heterodimerization of the Cohesin Loader Subunits and Cause Cornelia de Lange Syndrome

The NIPBL/MAU2 heterodimer loads cohesin onto chromatin. Mutations in NIPBL account for most cases of the rare developmental disorder Cornelia de Lange syndrome (CdLS). Here we report a MAU2 variant causing CdLS, a deletion of seven amino acids that impairs the interaction between MAU2 and the NIPBL N terminus. Investigating this interaction, we discovered that MAU2 and the NIPBL N terminus are largely dispensable for normal cohesin and NIPBL function in cells with a NIPBL early truncating mutation. Despite a predicted fatal outcome of an out-of-frame single nucleotide duplication in NIPBL, engineered in two different cell lines, alternative translation initiation yields a form of NIPBL missing N-terminal residues. This form cannot interact with MAU2, but binds DNA and mediates cohesin loading. Altogether, our work reveals that cohesin loading can occur independently of functional NIPBL/MAU2 complexes and highlights a novel mechanism protective against out-of-frame mutations that is potentially relevant for other genetic conditions.

Pan-tissue transcriptome analysis of long noncoding RNAs in the American beaver Castor canadensis

BACKGROUND

Long noncoding RNAs (lncRNAs) have roles in gene regulation, epigenetics, and molecular scaffolding and it is hypothesized that they underlie some mammalian evolutionary adaptations. However, for many mammalian species, the absence of a genome assembly precludes the comprehensive identification of lncRNAs. The genome of the American beaver (Castor canadensis) has recently been sequenced, setting the stage for the systematic identification of beaver lncRNAs and the characterization of their expression in various tissues. The objective of this study was to discover and profile polyadenylated lncRNAs in the beaver using high-throughput short-read sequencing of RNA from sixteen beaver tissues and to annotate the resulting lncRNAs based on their potential for orthology with known lncRNAs in other species.

RESULTS

Using de novo transcriptome assembly, we found 9528 potential lncRNA contigs and 187 high-confidence lncRNA contigs. Of the high-confidence lncRNA contigs, 147 have no known orthologs (and thus are putative novel lncRNAs) and 40 have mammalian orthologs. The novel lncRNAs mapped to the Oregon State University (OSU) reference beaver genome with greater than 90% sequence identity. While the novel lncRNAs were on average shorter than their annotated counterparts, they were similar to the annotated lncRNAs in terms of the relationships between contig length and minimum free energy (MFE) and between coverage and contig length. We identified beaver orthologs of known lncRNAs such as XIST, MEG3, TINCR, and NIPBL-DT. We profiled the expression of the 187 high-confidence lncRNAs across 16 beaver tissues (whole blood, brain, lung, liver, heart, stomach, intestine, skeletal muscle, kidney, spleen, ovary, placenta, castor gland, tail, toe-webbing, and tongue) and identified both tissue-specific and ubiquitous lncRNAs.

CONCLUSIONS

To our knowledge this is the first report of systematic identification of lncRNAs and their expression atlas in beaver. LncRNAs-both novel and those with known orthologs-are expressed in each of the beaver tissues that we analyzed. For some beaver lncRNAs with known orthologs, the tissue-specific expression patterns were phylogenetically conserved. The lncRNA sequence data files and raw sequence files are available via the web supplement and the NCBI Sequence Read Archive, respectively.

Analysis of Survival-Related lncRNA Landscape Identifies A Role for LINC01537 in Energy Metabolism and Lung Cancer Progression

Many long non-coding RNAs (lncRNAs) have emerged as good biomarkers and potential therapeutic targets for various cancers. We aimed to get a detailed understanding of the lncRNA landscape that is associated with lung cancer survival. A comparative analysis between our RNA sequencing (RNA-seq) data and TCGA datasets was conducted to reveal lncRNAs with significant correlations with lung cancer survival and then the association of the most promising lncRNA was validated in a cohort of 243 lung cancer patients. Comparing RNA-seq data with TCGA ones, 84 dysregulated lncRNAs were identified in lung cancer tissues, among which 10 lncRNAs were significantly associated with lung cancer survival. LINC01537 was the most significant one (p = 2.95 × 10⁻⁶). Validation analysis confirmed the downregulation of LINC01537 in lung cancer. LINC01537 was observed to inhibit tumor growth and metastasis. It also increased cellular sensitivity to nilotinib. PDE2A (phosphodiesterase 2A) was further identified to be a target of LINC01537 and it was seen that LINC01537 promoted PDE2A expression via RNA–RNA interaction to stabilize PDE2A mRNA and thus echoed effects of PDE2A on energy metabolism including both Warburg effect and mitochondrial respiration. Other regulators of tumor energy metabolism were also affected by LINC01537. These results elucidate a suppressed role of LINC01537 in lung cancer development involving tumor metabolic reprogramming, and we believe that it might be a biomarker for cancer survival prediction and therapy.

Emerging strategies to target the dysfunctional cohesin complex in cancer

INTRODUCTION

Mutations in cohesin genes have been described in numerous solid cancers and hematologic malignancies; subsequent experimental evidence has linked these mutations with carcinogenesis. Areas covered: In this review, we present current information about the physiological role of the cohesin complex in normal and malignant cells and describe current therapeutic strategies that are being explored in cohesin-mutated cancers. We discuss a range of targets and strategies that should be explored to develop targeted therapies for patients with aberrant cohesin. Expert opinion: Targeting of the cohesin complex is an underexplored area of drug development. There is a high frequency of cohesin mutations in multiple cancers, hence specific targeting strategies should be explored. Cohesins play a crucial role in cellular organization; therefore, we expect a narrow therapeutic window of direct inhibitors of cohesin components. Exploiting experimental approaches that correct dysfunctional cohesins and coupling them with current therapeutic strategies can provide novel, innovative and more effective treatment regimens.

Next generation sequencing identified two novel mutations in NIPBL and a frame shift mutation in CREBBP in three Chinese children

BACKGROUND

Cornelia de Lange syndrome (CdLS) and Rubinstein-Taybi syndrome (RSTS) are both rare congenital multiple malformation disorders caused by genes associated with transcription. They share a number of similar features clinically. In addition, it is difficult to make a molecular diagnosis rapidly and detect the mosaic mutation when only sanger sequencing is taken. This study aims to report three novel mutations in three Chinese children identified by next generation sequencing.

RESULTS

We describe patient 1 and patient 2 presenting with characteristics of CdLS with mutations in NIPBL and patient 3 with a frame shift mutation in CREBBP who can be diagnosed as RSTS clinically and also have similar symptoms with CdLS to some extent. The splicing site c.4321-1G > A transversion in NIPBL is a mosaic mutation and produces an abnormal transcript bearing the loss of exon 20. The nonsense mutation c.218C > A in NIPBL and the frame shift c.1715delC mutation in CREBBP generate stop codon and yield the premature termination of proteins.

CONCLUSIONS

In general, we detect three novel heterozygous mutations including a splicing mutation and a nonsense mutation in NIPBL and a frame shift in CREBBP. And several similar features observed in patients indicate the clinical complexity and clinically overlapping of CdLS and RSTS termed "transcriptomopathies", suggest the underlying molecular mechanism and emphasize the utilization of next generation sequencing technologies.

Roles of NIPBL in maintenance of genome stability

A cohesin-loading factor (NIPBL) is one of important regulatory factors in the maintenance of 3D genome organization and function, by interacting with a large number of factors, e.g. cohesion, CCCTC-binding factor (CTCF) or cohesin complex component. The present article overviews the critical and regulatory roles of NIBPL in cohesion loading on chromotin and in gene expression and transcriptional signaling. We explore molecular mechanisms by which NIPBL recruits endogenous histone deacetylase (HDAC) to induce histone deacetylation and influence multi-dimensions of genome, through which NIPBL "hop" movement in chromatin regulates gene expression and alters genome folding. NIPBL regulates the process of CTCF and cohesion into chromatin loops and topologically associated domains, binding of cohesion and H3K4mes3 through interaction among promoters and enhancers. HP1 recruits NIPBL to DNA damage site through RNF8/RNF168 ubiquitylation pathway. NIPBL contributes to regulation of genome-controlled gene expression through the influence of cohesin in chromosome structure. NIPBL interacts with cohesin and then increases transcriptional activities of REC8 promoter, leading to up-regulation of gene expression. NIPBL movement among chromosomal loops regulates gene expression through dynamic alterations of genome organization. Thus, we expect a new and deep insight to understand dynamics of chromosome and explore potential strategies of therapiesc on basis of NIPBL.