Annotation of nuclear lncRNAs based on chromatin interactions

The human genome is pervasively transcribed and produces a wide variety of long non-coding RNAs (lncRNAs), constituting the majority of transcripts across human cell types. Some specific nuclear lncRNAs have been shown to be important regulatory components acting locally. As RNA-chromatin interaction and Hi-C chromatin conformation data showed that chromatin interactions of nuclear lncRNAs are determined by the local chromatin 3D conformation, we used Hi-C data to identify potential target genes of lncRNAs. RNA-protein interaction data suggested that nuclear lncRNAs act as scaffolds to recruit regulatory proteins to target promoters and enhancers. Nuclear lncRNAs may therefore play a role in directing regulatory factors to locations spatially close to the lncRNA gene. We provide the analysis results through an interactive visualization web portal at https://fantom.gsc.riken.jp/zenbu/reports/#F6_3D_lncRNA.

Use of clinical variables for preoperative prediction of lymph node metastasis in endometrial cancer

OBJECTIVE

Endometrial cancer is the most common gynaecological cancer, and most patients are identified during early disease stages. Noninvasive evaluation of lymph node metastasis likely will improve the quality of clinical treatment, for example, by omitting unnecessary lymphadenectomy.

METHODS

The study population comprised 611 patients with endometrial cancer who underwent lymphadenectomy at four types of institutions, comprising seven hospitals in total. We systematically assessed the association of 18 preoperative clinical variables with postoperative lymph node metastasis. We then constructed statistical models for preoperative lymph node metastasis prediction and assessed their performance with a previously proposed system, in which the score was determined by counting the number of high-risk variables among the four predefined ones.

RESULTS

Of the preoperative 18 variables evaluated, 10 were significantly associated with postoperative lymph node metastasis. A logistic regression model achieved an area under the curve of 0.85 in predicting lymph node metastasis; this value is significantly higher than that from the previous system (area under the curve, 0.74). When we set the false-negative rate to ~1%, the new predictive model increased the rate of true negatives to 21%, compared with 6.8% from the previous one. We also provide a spreadsheet-based tool for further evaluation of its ability to predict lymph node metastasis in endometrial cancer.

CONCLUSIONS

Our new lymph node metastasis prediction method, which was based solely on preoperative clinical variables, performed significantly better than the previous method. Although additional evaluation is necessary for its clinical use, our noninvasive system may help improve the clinical treatment of endometrial cancer, complementing minimally invasive sentinel lymph node biopsy.

ROR2 expression predicts human induced pluripotent stem cell differentiation into neural stem/progenitor cells and GABAergic neurons

Despite the development of various in vitro differentiation protocols for the efficient derivation of specific cell types, human induced pluripotent stem cell (hiPSC) lines have varing ability to differentiate into specific lineages. Therefore, surrogate markers for accurately predicting the differentiation propensity of hiPSC lines may facilitate cell-based therapeutic product development and manufacture. We attempted to identify marker genes that could predict the differentiation propensity of hiPSCs into neural stem/progenitor cells (NS/PCs). Using Spearman's rank correlation coefficients, we investigated genes in the undifferentiated state, the expression levels of which were significantly correlated with the neuronal differentiation propensity of several hiPSC lines. Among genes significantly correlated with NS/PC differentiation (P < 0.01), we identified ROR2 as a novel predictive marker. ROR2 expression in hiPSCs was negatively correlated with NS/PC differentiation tendency, regardless of the differentiation method, whereas its knockdown enhanced differentiation. ROR2 regulates NS/PC differentiation, suggesting that ROR2 is functionally essential for NS/PC differentiation. Selecting cell lines with relatively low ROR2 expression facilitated identification of hiPSCs that can differentiate into NS/PCs. Cells with ROR2 knockdown showed increased efficiency of differentiation into forebrain GABAergic neurons compared to controls. These findings suggest that ROR2 is a surrogate marker for selecting hiPSC lines appropriate for NS/PC and GABAergic neuronal differentiations.

Single-Cell RNA-Seq Reveals LRRC75A-Expressing Cell Population Involved in VEGF Secretion of Multipotent Mesenchymal Stromal/Stem Cells Under Ischemia

Human multipotent mesenchymal stromal/stem cells (MSCs) have been utilized in cell therapy for various diseases and their clinical applications are expected to increase in the future. However, the variation in MSC-based product quality due to the MSC heterogeneity has resulted in significant constraints in the clinical utility of MSCs. Therefore, we hypothesized that it might be important to identify and ensure/enrich suitable cell subpopulations for therapies using MSC-based products. In this study, we aimed to identify functional cell subpopulations to predict the efficacy of angiogenic therapy using bone marrow-derived MSCs (BM-MSCs). To assess its angiogenic potency, we observed various levels of vascular endothelial growth factor (VEGF) secretion among 11 donor-derived BM-MSC lines under in vitro ischemic culture conditions. Next, by clarifying the heterogeneity of BM-MSCs using single-cell RNA-sequencing analysis, we identified a functional cell subpopulation that contributed to the overall VEGF production in BM-MSC lines under ischemic conditions. We also found that leucine-rich repeat-containing 75A (LRRC75A) was more highly expressed in this cell subpopulation than in the others. Importantly, knockdown of LRRC75A using small interfering RNA resulted in significant inhibition of VEGF secretion in ischemic BM-MSCs, indicating that LRRC75A regulates VEGF secretion under ischemic conditions. Therefore, LRRC75A may be a useful biomarker to identify cell subpopulations that contribute to the angiogenic effects of BM-MSCs. Our work provides evidence that a strategy based on single-cell transcriptome profiles is effective for identifying functional cell subpopulations in heterogeneous MSC-based products.

CALML5 is a novel diagnostic marker for differentiating thymic squamous cell carcinoma from type B3 thymoma

BACKGROUND

Thymic squamous cell carcinoma and type B3 thymoma are primary neoplasms of the anterior mediastinum that are sometimes difficult to differentiate from one another histologically. However, only a few immunohistochemical markers are available for the differential diagnosis. The purpose of this study was to discover a novel marker for differentiating between thymic squamous cell carcinoma and type B3 thymoma.

METHODS

We used histological samples of thymic carcinomas (n = 26) and type B3 thymomas (n = 38) which were resected between 1986 and 2017. To search for candidates of differential markers, gene expression levels were evaluated in samples using promoter analysis by cap analysis of gene expression (CAGE) sequencing.

RESULTS

Promoter level expression of CALML5 genes was significantly higher in thymic carcinomas than in type B3 thymomas. We further validated the results of the CAGE analysis in all 26 thymic carcinomas and 38 type B3 thymomas by immunohistochemistry (IHC). CALML5 was strongly expressed in the cytoplasm in 19 of 26 cases with thymic carcinoma, whereas positivity at the protein level was shown in two of 38 type B3 thymomas. Thus, the sensitivity (73.1%) and specificity (94.7%) of CALML5 as markers for immunohistochemical diagnosis of thymic carcinoma were extremely high.

CONCLUSION

We identified CALML5 as a potential marker for differentiating thymic squamous cell carcinoma from type B3 thymoma. It is assumed that future clinical use of CALML5 may improve the diagnostic accuracy of differentiating between these two diseases.

Antisense-oligonucleotide-mediated perturbation of long non-coding RNA reveals functional features in stem cells and across cell types

Within the scope of the FANTOM6 consortium, we perform a large-scale knockdown of 200 long non-coding RNAs (lncRNAs) in human induced pluripotent stem cells (iPSCs) and systematically characterize their roles in self-renewal and pluripotency. We find 36 lncRNAs (18%) exhibiting cell growth inhibition. From the knockdown of 123 lncRNAs with transcriptome profiling, 36 lncRNAs (29.3%) show molecular phenotypes. Integrating the molecular phenotypes with chromatin-interaction assays further reveals cis- and trans-interacting partners as potential primary targets. Additionally, cell-type enrichment analysis identifies lncRNAs associated with pluripotency, while the knockdown of LINC02595, CATG00000090305.1, and RP11-148B6.2 modulates colony formation of iPSCs. We compare our results with previously published fibroblasts phenotyping data and find that 2.9% of the lncRNAs exhibit a consistent cell growth phenotype, whereas we observe 58.3% agreement in molecular phenotypes. This highlights that molecular phenotyping is more comprehensive in revealing affected pathways.

Deep sequencing of short capped RNAs reveals novel families of noncoding RNAs

In eukaryotes, capped RNAs include long transcripts such as messenger RNAs and long noncoding RNAs, as well as shorter transcripts such as spliceosomal RNAs, small nucleolar RNAs, and enhancer RNAs. Long capped transcripts can be profiled using cap analysis gene expression (CAGE) sequencing and other methods. Here, we describe a sequencing library preparation protocol for short capped RNAs, apply it to a differentiation time course of the human cell line THP-1, and systematically compare the landscape of short capped RNAs to that of long capped RNAs. Transcription initiation peaks associated with genes in the sense direction have a strong preference to produce either long or short capped RNAs, with one out of six peaks detected in the short capped RNA libraries only. Gene-associated short capped RNAs have highly specific 3' ends, typically overlapping splice sites. Enhancers also preferentially generate either short or long capped RNAs, with 10% of enhancers observed in the short capped RNA libraries only. Enhancers producing either short or long capped RNAs show enrichment for GWAS-associated disease SNPs. We conclude that deep sequencing of short capped RNAs reveals new families of noncoding RNAs and elucidates the diversity of transcripts generated at known and novel promoters and enhancers.

High quality genome assembly of the anhydrobiotic midge provides insights on a single chromosome-based emergence of extreme desiccation tolerance

Non-biting midges (Chironomidae) are known to inhabit a wide range of environments, and certain species can tolerate extreme conditions, where the rest of insects cannot survive. In particular, the sleeping chironomid Polypedilum vanderplanki is known for the remarkable ability of its larvae to withstand almost complete desiccation by entering a state called anhydrobiosis. Chromosome numbers in chironomids are higher than in other dipterans and this extra genomic resource might facilitate rapid adaptation to novel environments. We used improved sequencing strategies to assemble a chromosome-level genome sequence for P. vanderplanki for deep comparative analysis of genomic location of genes associated with desiccation tolerance. Using whole genome-based cross-species and intra-species analysis, we provide evidence for the unique functional specialization of Chromosome 4 through extensive acquisition of novel genes. In contrast to other insect genomes, in the sleeping chironomid a uniquely high degree of subfunctionalization in paralogous anhydrobiosis genes occurs in this chromosome, as well as pseudogenization in a highly duplicated gene family. Our findings suggest that the Chromosome 4 in Polypedilum is a site of high genetic turnover, allowing it to act as a ‘sandbox’ for evolutionary experiments, thus facilitating the rapid adaptation of midges to harsh environments.

The choice of negative control antisense oligonucleotides dramatically impacts downstream analysis depending on the cellular background

Background

The lymphatic and the blood vasculature are closely related systems that collaborate to ensure the organism’s physiological function. Despite their common developmental origin, they present distinct functional fates in adulthood that rely on robust lineage-specific regulatory programs. The recent technological boost in sequencing approaches unveiled long noncoding RNAs (lncRNAs) as prominent regulatory players of various gene expression levels in a cell-type-specific manner.

Results

To investigate the potential roles of lncRNAs in vascular biology, we performed antisense oligonucleotide (ASO) knockdowns of lncRNA candidates specifically expressed either in human lymphatic or blood vascular endothelial cells (LECs or BECs) followed by Cap Analysis of Gene Expression (CAGE-Seq). Here, we describe the quality control steps adopted in our analysis pipeline before determining the knockdown effects of three ASOs per lncRNA target on the LEC or BEC transcriptomes. In this regard, we especially observed that the choice of negative control ASOs can dramatically impact the conclusions drawn from the analysis depending on the cellular background.

Conclusion

In conclusion, the comparison of negative control ASO effects on the targeted cell type transcriptomes highlights the essential need to select a proper control set of multiple negative control ASO based on the investigated cell types.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-021-00992-1.

Nanopore sequencing reveals TACC2 locus complexity and diversity of isoforms transcribed from an intronic promoter

Gene expression is controlled at the transcriptional and post-transcriptional levels. The TACC2 gene was known to be associated with tumors but the control of its expression is unclear. We have reported that activity of the intronic promoter p10 of TACC2 in primary lesion of endometrial cancer is indicative of lymph node metastasis among a low-risk patient group. Here, we analyze the intronic promoter derived isoforms in JHUEM-1 endometrial cancer cells, and primary tissues of endometrial cancers and normal endometrium. Full-length cDNA amplicons are produced by long-range PCR and subjected to nanopore sequencing followed by computational error correction. We identify 16 stable, 4 variable, and 9 rare exons including 3 novel exons validated independently. All variable and rare exons reside N-terminally of the TACC domain and contribute to isoform variety. We found 240 isoforms as high-confidence, supported by more than 20 reads. The large number of isoforms produced from one minor promoter indicates the post-transcriptional complexity coupled with transcription at the TACC2 locus in cancer and normal cells.

Cancer-associated fibroblast migration in non-small cell lung cancers is modulated by increased integrin α11 expression

Cancer‐associated fibroblasts (CAFs) regulate cancer progression through the modulation of extracellular matrix (ECM) and cancer cell adhesion. While undergoing a series of phenotypic changes, CAFs control cancer–stroma interactions through integrin receptor signaling. Here, we isolated CAFs from patients with non‐small‐cell lung cancer (NSCLC) and examined their gene expression profiles. We identified collagen type XI α1 (COL11A1), integrin α11 (ITGA11), and the ITGA11 major ligand collagen type I α1 (COL1A1) among the 390 genes that were significantly enriched in NSCLC‐associated CAFs. Increased ITGA11 expression in cancer stroma was correlated with a poor clinical outcome in patients with NSCLC. Increased expression of fibronectin and collagen type I induced ITGA11 expression in CAFs. The cellular migration of CAFs toward collagen type I and fibronectin was promoted via ERK1/2 signaling, independently of the fibronectin receptor integrin α5β1. Additionally, ERK1/2 signaling induced ITGA11 and COL11A1 expression in cancer stroma. We, therefore, propose that targeting ITGA11 and COL11A1 expressing CAFs to block cancer–stroma interactions may serve as a novel, promising anti‐tumor strategy.

Ex vivo conditioning of peripheral blood mononuclear cells of diabetic patients promotes vasculogenic wound healing

The quality and quantity of endothelial progenitor cells (EPCs) are impaired in patients with diabetes mellitus patients, leading to reduced tissue repair during autologous EPC therapy. This study aimed to address the limitations of the previously described serum-free Quantity and Quality Control Culture System (QQc) using CD34+ cells by investigating the therapeutic potential of a novel mononuclear cell (MNC)-QQ. MNCs were isolated from 50 mL of peripheral blood of patients with diabetes mellitus and healthy volunteers (n = 13 each) and subjected to QQc for 7 days in serum-free expansion media with VEGF, Flt-3 ligand, TPO, IL-6, and SCF. The vascular regeneration capability of MNC-QQ cells pre- or post-QQc was evaluated with an EPC colony-forming assay, FACS, EPC culture, tube formation assay, and quantitative real time PCR. For in vivo assessment, 1 × 10⁴ pre- and post-MNC-QQc cells from diabetic donors were injected into a murine wound-healing model using Balb/c nude mice. The percentage of wound closure and angio-vasculogenesis was then assessed. This study revealed vasculogenic, anti-inflammatory, and wound-healing effects of MNC-QQ therapy in both in vitro and in vivo models. This system addresses the low efficiency and efficacy of the current naïve MNC therapy for wound-healing in diabetic patients. As this technique requires a simple blood draw, isolation, and peripheral blood MNC suspension culture for only a week, it can be used as a simple and effective outpatient-based vascular and regenerative therapy for patients with diabetes mellitus.

CREB3L1 overexpression as a potential diagnostic marker of Philadelphia chromosome-negative myeloproliferative neoplasms

Discrimination of Philadelphia-negative myeloproliferative neoplasms (Ph-MPNs) from reactive hypercytosis and myelofibrosis requires a constellation of testing including driver mutation analysis and bone marrow biopsies. We searched for a biomarker that can more easily distinguish Ph-MPNs from reactive hypercytosis and myelofibrosis by using RNA-seq analysis utilizing platelet-rich plasma (PRP)-derived RNAs from patients with essential thrombocythemia (ET) and reactive thrombocytosis, and CREB3L1 was found to have an extremely high impact in discriminating the two disorders. To validate and further explore the result, expression levels of CREB3L1 in PRP were quantified by reverse-transcription quantitative PCR and compared among patients with ET, other Ph-MPNs, chronic myeloid leukemia (CML), and reactive hypercytosis and myelofibrosis. A CREB3L1 expression cutoff value determined based on PRP of 18 healthy volunteers accurately discriminated 150 driver mutation-positive Ph-MPNs from other entities (71 reactive hypercytosis and myelofibrosis, 6 CML, and 18 healthy volunteers) and showed both sensitivity and specificity of 1.0000. Importantly, CREB3L1 expression levels were significantly higher in ET compared with reactive thrombocytosis (P < .0001), and polycythemia vera compared with reactive erythrocytosis (P < .0001). Pathology-affirmed triple-negative ET (TN-ET) patients were divided into a high- and low-CREB3L1-expression group, and some patients in the low-expression group achieved a spontaneous remission during the clinical course. In conclusion, CREB3L1 analysis has the potential to single-handedly discriminate driver mutation-positive Ph-MPNs from reactive hypercytosis and myelofibrosis, and also may identify a subgroup within TN-ET showing distinct clinical features including spontaneous remission.

Characteristics of Mandibular Canal Branches Related to Nociceptive Marker

It is difficult to correlate the direction of mandibular canal branches (MCBs) with altered sensation in dental treatments. In contrast, calcitonin gene-related peptide (CGRP) is related to vasodilation, bone formation, and the interaction with the peripheral nervous system. Therefore, we investigated the detailed morphological characteristics of MCBs using cone-beam computed tomography (CBCT) and observation of the CGRP distribution around the MCB. The MCB measurements were evaluated using principal component analysis (PCA) to identify morphological correlations. A total of 168 sides of mandibles from 84 cadavers were analyzed in this study. Most of the MCBs were primarily in the direction of the clock model from X to XI in sagittal sections and XII to I in coronal sections of the mandible. The structure of the MCB was divided into the fine canal branch (60.4%, 223/369), partial branch (24.4%, 90/369), and no canal branch (15.2%, 56/369). PCA indicated that the measurement element with the MCB and its structures were correlated in contrast to tooth factors. Positive CGRP reactions were clearly observed in the no-canal branch group compared to other groups. These data provide useful suggestions for MCB dynamics and information for clinical dental treatment.

Low Quantity Single Strand CAGE (LQ-ssCAGE) Maps Regulatory Enhancers and Promoters

The Cap Analysis of Gene Expression (CAGE) is a powerful method to identify Transcription Start Sites (TSSs) of capped RNAs while simultaneously measuring transcripts expression level. CAGE allows mapping at single nucleotide resolution at all active promoters and enhancers. Large CAGE datasets have been produced over the years from individual laboratories and consortia, including the Encyclopedia of DNA Elements (ENCODE) and Functional Annotation of the Mammalian Genome (FANTOM) consortia. These datasets constitute open resource for TSS annotations and gene expression analysis. Here, we provide an experimental protocol for the most recent CAGE method called Low Quantity (LQ) single strand (ss) CAGE "LQ-ssCAGE", which enables cost-effective profiling of low quantity RNA samples. LQ-ssCAGE is especially useful for samples derived from cells cultured in small volumes, cellular compartments such as nuclear RNAs or for samples from developmental stages. We demonstrate the reproducibility and effectiveness of the method by constructing 240 LQ-ssCAGE libraries from 50 ng of THP-1 cell extracted RNAs and discover lowly expressed novel enhancer and promoter-derived lncRNAs.

Scanning single-molecule counting system for Eprobe with highly simple and effective approach

Here, we report a rapid and ultra-sensitive detection technique for fluorescent molecules called scanning single molecular counting (SSMC). The method uses a fluorescence-based digital measurement system to count single molecules in a solution. In this technique, noise is reduced by conforming the signal shape to the intensity distribution of the excitation light via a circular scan of the confocal region. This simple technique allows the fluorescent molecules to freely diffuse into the solution through the confocal region and be counted one by one and does not require statistical analysis. Using this technique, 28 to 62 aM fluorescent dye was detected through measurement for 600 s. Furthermore, we achieved a good signal-to-noise ratio (S/N = 2326) under the condition of 100 pM target nucleic acid by only mixing a hybridization-sensitive fluorescent probe, called Eprobe, into the target oligonucleotide solution. Combination of SSMC and Eprobe provides a simple, rapid, amplification-free, and high-sensitive target nucleic acid detection system. This method is promising for future applications to detect particularly difficult to design primers for amplification as miRNAs and other short oligo nucleotide biomarkers by only hybridization with high sensitivity.

Development of a Tet-On Inducible Expression System for the Anhydrobiotic Cell Line, Pv11

The Pv11 cell line established from an African chironomid, Polypedilum vanderplanki, is the only cell line tolerant to complete desiccation. In Pv11 cells, a constitutive expression system for Pv11 cells was previously exploited and several reporter genes were successfully expressed. Here we report the identification of an effective minimal promoter for Pv11 cells and its application to the Tet-On inducible expression system. First, using a luciferase reporter assay, we showed that a 202 bp deletion fragment derived from the constitutively active 121-promoter functions in Pv11 cells as an appropriate minimal promoter with the Tet-On inducible expression system. The AcGFP1 (Aequorea coerulescens green fluorescent protein) was also successfully expressed in Pv11 cells using the inducible system. In addition to these reporter genes, the avian myeloblastosis virus reverse transcriptase α subunit (AMV RTα), which is one of the most widely commercially available RNA-dependent DNA polymerases, was successfully expressed through the inducible expression system and its catalytic activity was verified. These results demonstrate the establishment of an inducible expression system in cells that can be preserved in the dry state and highlight a possible application to the production of large and complex proteins.

Comparative transcriptomics of primary cells in vertebrates

Gene expression profiles in homologous tissues have been observed to be different between species, which may be due to differences between species in the gene expression program in each cell type, but may also reflect differences in cell type composition of each tissue in different species. Here, we compare expression profiles in matching primary cells in human, mouse, rat, dog, and chicken using Cap Analysis Gene Expression (CAGE) and short RNA (sRNA) sequencing data from FANTOM5. While we find that expression profiles of orthologous genes in different species are highly correlated across cell types, in each cell type many genes were differentially expressed between species. Expression of genes with products involved in transcription, RNA processing, and transcriptional regulation was more likely to be conserved, while expression of genes encoding proteins involved in intercellular communication was more likely to have diverged during evolution. Conservation of expression correlated positively with the evolutionary age of genes, suggesting that divergence in expression levels of genes critical for cell function was restricted during evolution. Motif activity analysis showed that both promoters and enhancers are activated by the same transcription factors in different species. An analysis of expression levels of mature miRNAs and of primary miRNAs identified by CAGE revealed that evolutionary old miRNAs are more likely to have conserved expression patterns than young miRNAs. We conclude that key aspects of the regulatory network are conserved, while differential expression of genes involved in cell-to-cell communication may contribute greatly to phenotypic differences between species.

Functional annotation of human long noncoding RNAs via molecular phenotyping

Long noncoding RNAs (lncRNAs) constitute the majority of transcripts in the mammalian genomes, and yet, their functions remain largely unknown. As part of the FANTOM6 project, we systematically knocked down the expression of 285 lncRNAs in human dermal fibroblasts and quantified cellular growth, morphological changes, and transcriptomic responses using Capped Analysis of Gene Expression (CAGE). Antisense oligonucleotides targeting the same lncRNAs exhibited global concordance, and the molecular phenotype, measured by CAGE, recapitulated the observed cellular phenotypes while providing additional insights on the affected genes and pathways. Here, we disseminate the largest-to-date lncRNA knockdown data set with molecular phenotyping (over 1000 CAGE deep-sequencing libraries) for further exploration and highlight functional roles for ZNF213-AS1 and lnc-KHDC3L-2.

Recounting the FANTOM CAGE-Associated Transcriptome

Long noncoding RNAs (lncRNAs) have emerged as key coordinators of biological and cellular processes. Characterizing lncRNA expression across cells and tissues is key to understanding their role in determining phenotypes, including human diseases. We present here FC-R2, a comprehensive expression atlas across a broadly defined human transcriptome, inclusive of over 109,000 coding and noncoding genes, as described in the FANTOM CAGE-Associated Transcriptome (FANTOM-CAT) study. This atlas greatly extends the gene annotation used in the original recount2 resource. We demonstrate the utility of the FC-R2 atlas by reproducing key findings from published large studies and by generating new results across normal and diseased human samples. In particular, we (a) identify tissue-specific transcription profiles for distinct classes of coding and noncoding genes, (b) perform differential expression analysis across thirteen cancer types, identifying novel noncoding genes potentially involved in tumor pathogenesis and progression, and (c) confirm the prognostic value for several enhancer lncRNAs expression in cancer. Our resource is instrumental for the systematic molecular characterization of lncRNA by the FANTOM6 Consortium. In conclusion, comprised of over 70,000 samples, the FC-R2 atlas will empower other researchers to investigate functions and biological roles of both known coding genes and novel lncRNAs.

Cap Analysis of Gene Expression (CAGE): A Quantitative and Genome-Wide Assay of Transcription Start Sites

Cap analysis of gene expression (CAGE) is an approach to identify and monitor the activity (transcription initiation frequency) of transcription start sites (TSSs) at single base-pair resolution across the genome. It has been effectively used to identify active promoter and enhancer regions in cancer cells, with potential utility to identify key factors to immunotherapy. Here, we overview a series of CAGE protocols and describe detailed experimental steps of the latest protocol based on the Illumina sequencing platform; both experimental steps (see Subheadings 3.1-3.11) and computational processing steps (see Subheadings 3.12-3.20) are described.

Analysis of the Beaufort Gyre Freshwater Content in 2003-2018

Hydrographic data collected from research cruises, bottom-anchored moorings, drifting Ice-Tethered Profilers, and satellite altimetry in the Beaufort Gyre region of the Arctic Ocean document an increase of more than 6,400 km3 of liquid freshwater content from 2003 to 2018: a 40% growth relative to the climatology of the 1970s. This fresh water accumulation is shown to result from persistent anticyclonic atmospheric wind forcing (1997-2018) accompanied by sea ice melt, a wind-forced redirection of Mackenzie River discharge from predominantly eastward to westward flow, and a contribution of low salinity waters of Pacific Ocean origin via Bering Strait. Despite significant uncertainties in the different observations, this study has demonstrated the synergistic value of having multiple diverse datasets to obtain a more comprehensive understanding of Beaufort Gyre freshwater content variability. For example, Beaufort Gyre Observational System (BGOS) surveys clearly show the interannual increase in freshwater content, but without satellite or Ice-Tethered Profiler measurements, it is not possible to resolve the seasonal cycle of freshwater content, which in fact is larger than the year-to-year variability, or the more subtle interannual variations.

Identification of novel cerebellar developmental transcriptional regulators with motif activity analysis

Background

The work of the FANTOM5 Consortium has brought forth a new level of understanding of the regulation of gene transcription and the cellular processes involved in creating diversity of cell types. In this study, we extended the analysis of the FANTOM5 Cap Analysis of Gene Expression (CAGE) transcriptome data to focus on understanding the genetic regulators involved in mouse cerebellar development.

Results

We used the HeliScopeCAGE library sequencing on cerebellar samples over 8 embryonic and 4 early postnatal times. This study showcases temporal expression pattern changes during cerebellar development. Through a bioinformatics analysis that focused on transcription factors, their promoters and binding sites, we identified genes that appear as strong candidates for involvement in cerebellar development. We selected several candidate transcriptional regulators for validation experiments including qRT-PCR and shRNA transcript knockdown. We observed marked and reproducible developmental defects in Atf4, Rfx3, and Scrt2 knockdown embryos, which support the role of these genes in cerebellar development.

Conclusions

The successful identification of these novel gene regulators in cerebellar development demonstrates that the FANTOM5 cerebellum time series is a high-quality transcriptome database for functional investigation of gene regulatory networks in cerebellar development.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-6063-9) contains supplementary material, which is available to authorized users.

RNA sequencing analysis revealed the induction of CCL3 expression in human intracranial aneurysms

Intracranial aneurysm (IA) is a socially important disease as a major cause of subarachnoid hemorrhage. Recent experimental studies mainly using animal models have revealed a crucial role of macrophage-mediated chronic inflammatory responses in its pathogenesis. However, as findings from comprehensive analysis of unruptured human IAs are limited, factors regulating progression and rupture of IAs in humans remain unclear. Using surgically dissected human unruptured IA lesions and control arterial walls, gene expression profiles were obtained by RNA sequence analysis. RNA sequencing analysis was done with read count about 60~100 million which yielded 6~10 billion bases per sample. 79 over-expressed and 329 under-expressed genes in IA lesions were identified. Through Gene Ontology analysis, ‘chemokine activity’, ‘defense response’ and ‘extracellular region’ were picked up as over-represented terms which included CCL3 and CCL4 in common. Among these genes, quantitative RT-PCR analysis using another set of samples reproduced the above result. Finally, increase of CCL3 protein compared with that in control arterial walls was clarified in IA lesions. Findings of the present study again highlight importance of macrophage recruitment via CCL3 in the pathogenesis of IA progression.

Integration of genetics and miRNA-target gene network identified disease biology implicated in tissue specificity

MicroRNAs (miRNAs) modulate the post-transcriptional regulation of target genes and are related to biology of complex human traits, but genetic landscape of miRNAs remains largely unknown. Given the strikingly tissue-specific miRNA expression profiles, we here expand a previous method to quantitatively evaluate enrichment of genome-wide association study (GWAS) signals on miRNA-target gene networks (MIGWAS) to further estimate tissue-specific enrichment. Our approach integrates tissue-specific expression profiles of miRNAs (∼1800 miRNAs in 179 cells) with GWAS to test whether polygenic signals enrich in miRNA-target gene networks and whether they fall within specific tissues. We applied MIGWAS to 49 GWASs (nTotal = 3 520 246), and successfully identified biologically relevant tissues. Further, MIGWAS could point miRNAs as candidate biomarkers of the trait. As an illustrative example, we performed differentially expressed miRNA analysis between rheumatoid arthritis (RA) patients and healthy controls (n = 63). We identified novel biomarker miRNAs (e.g. hsa-miR-762) by integrating differentially expressed miRNAs with MIGWAS results for RA, as well as novel associated loci with significant genetic risk (rs56656810 at MIR762 at 16q11; n = 91 482, P = 3.6 × 10-8). Our result highlighted that miRNA-target gene network contributes to human disease genetics in a cell type-specific manner, which could yield an efficient screening of miRNAs as promising biomarkers.

CXCL4/PF4 is a predictive biomarker of cardiac differentiation potential of human induced pluripotent stem cells

Selection of human induced pluripotent stem cell (hiPSC) lines with high cardiac differentiation potential is important for regenerative therapy and drug screening. We aimed to identify biomarkers for predicting cardiac differentiation potential of hiPSC lines by comparing the gene expression profiles of six undifferentiated hiPSC lines with different cardiac differentiation capabilities. We used three platforms of gene expression analysis, namely, cap analysis of gene expression (CAGE), mRNA array, and microRNA array to efficiently screen biomarkers related to cardiac differentiation of hiPSCs. Statistical analysis revealed candidate biomarker genes with significant correlation between the gene expression levels in the undifferentiated hiPSCs and their cardiac differentiation potential. Of the candidate genes, PF4 was validated as a biomarker expressed in undifferentiated hiPSCs with high potential for cardiac differentiation in 13 additional hiPSC lines. Our observations suggest that PF4 may be a useful biomarker for selecting hiPSC lines appropriate for the generation of cardiomyocytes.

Antisense Transcription in Loci Associated to Hereditary Neurodegenerative Diseases

Natural antisense transcripts are common features of mammalian genes providing additional regulatory layers of gene expression. A comprehensive description of antisense transcription in loci associated to familial neurodegenerative diseases may identify key players in gene regulation and provide tools for manipulating gene expression. We take advantage of the FANTOM5 sequencing datasets that represent the largest collection to date of genome-wide promoter usage in almost 2000 human samples. Transcription start sites (TSSs) are mapped at high resolution by the use of a modified protocol of cap analysis of gene expression (CAGE) for high-throughput single molecule next-generation sequencing with Helicos (hCAGE). Here we present the analysis of antisense transcription at 17 loci associated to hereditary Alzheimer’s disease, Frontotemporal Dementia, Parkinson’s disease, Amyotrophic Lateral Sclerosis, and Huntington’s disease. We focused our analysis on libraries derived from brain tissues and primary cells. We also screened libraries from total blood and blood cell populations in the quest for peripheral biomarkers of neurodegenerative diseases. We identified 63 robust promoters in antisense orientation to genes associated to familial neurodegeneration. When applying a less stringent cutoff, this number increases to over 400. A subset of these promoters represents alternative TSSs for 24 FANTOM5 annotated long noncoding RNA (lncRNA) genes, in antisense orientation to 13 of the loci analyzed here, while the remaining contribute to the expression of additional transcript variants. Intersection with GWAS studies, sample ontology, and dynamic expression reveals association to specific genetic traits as well as cell and tissue types, not limited to neurodegenerative diseases. Antisense transcription was validated for a subset of genes, including those encoding for Microtubule-Associated Protein Tau, α-synuclein, Parkinsonism-associated deglycase DJ-1, and Leucin-Rich Repeat Kinase 2. This work provides evidence for the existence of additional regulatory mechanisms of the expression of neurodegenerative disease-causing genes by previously not-annotated and/or not-validated antisense long noncoding RNAs.

Effective multiple sideband generation using an electro-optic modulator for a multiple isotope magneto-optical trap

Herein, we report an effective method for the generation of radio-frequency (RF) sidebands in an electro-optic modulator for the simultaneous magneto-optical trapping of two isotopes. This is achieved by switching the RF signals alternately, which suppresses the generation of unwanted frequency signals and improves the laser power per sideband. The generated sidebands are successfully applied to a dual-rubidium-isotope magneto-optical trap (MOT), which results in an increased number of trapped atoms. This simple, flexible, and robust technique can be implemented in experiments that require a large number of atoms in multiple-isotope MOTs and for various applications.

Promoter Usage and Dynamics in Vascular Smooth Muscle Cells Exposed to Fibroblast Growth Factor-2 or Interleukin-1β

Smooth muscle cells (SMC) in blood vessels are normally growth quiescent and transcriptionally inactive. Our objective was to understand promoter usage and dynamics in SMC acutely exposed to a prototypic growth factor or pro-inflammatory cytokine. Using cap analysis gene expression (FANTOM5 project) we report differences in promoter dynamics for immediate-early genes (IEG) and other genes when SMC are exposed to fibroblast growth factor-2 or interleukin-1β. Of the 1871 promoters responding to FGF2 or IL-1β considerably more responded to FGF2 (68.4%) than IL-1β (18.5%) and 13.2% responded to both. Expression clustering reveals sets of genes induced, repressed or unchanged. Among IEG responding rapidly to FGF2 or IL-1β were FOS, FOSB and EGR-1, which mediates human SMC migration. Motif activity response analysis (MARA) indicates most transcription factor binding motifs in response to FGF2 were associated with a sharp induction at 1 h, whereas in response to IL-1β, most motifs were associated with a biphasic change peaking generally later. MARA revealed motifs for FOS_FOS{B,L1}_JUN{B,D} and EGR-1..3 in the cluster peaking 1 h after FGF2 exposure whereas these motifs were in clusters peaking 1 h or later in response to IL-1β. Our findings interrogating CAGE data demonstrate important differences in promoter usage and dynamics in SMC exposed to FGF2 or IL-1β.

Conserved temporal ordering of promoter activation implicates common mechanisms governing the immediate early response across cell types and stimuli

The promoters of immediate early genes (IEGs) are rapidly activated in response to an external stimulus. These genes, also known as primary response genes, have been identified in a range of cell types, under diverse extracellular signals and using varying experimental protocols. Whereas genomic dissection on a case-by-case basis has not resulted in a comprehensive catalogue of IEGs, a rigorous meta-analysis of eight genome-wide FANTOM5 CAGE (cap analysis of gene expression) time course datasets reveals successive waves of promoter activation in IEGs, recapitulating known relationships between cell types and stimuli: we obtain a set of 57 (42 protein-coding) candidate IEGs possessing promoters that consistently drive a rapid but transient increase in expression over time. These genes show significant enrichment for known IEGs reported previously, pathways associated with the immediate early response, and include a number of non-coding RNAs with roles in proliferation and differentiation. Surprisingly, we also find strong conservation of the ordering of activation for these genes, such that 77 pairwise promoter activation orderings are conserved. Using the leverage of comprehensive CAGE time series data across cell types, we also document the extensive alternative promoter usage by such genes, which is likely to have been a barrier to their discovery until now. The common activation ordering of the core set of early-responding genes we identify may indicate conserved underlying regulatory mechanisms. By contrast, the considerably larger number of transiently activated genes that are specific to each cell type and stimulus illustrates the breadth of the primary response.

Shared activity patterns arising at genetic susceptibility loci reveal underlying genomic and cellular architecture of human disease

Genetic variants underlying complex traits, including disease susceptibility, are enriched within the transcriptional regulatory elements, promoters and enhancers. There is emerging evidence that regulatory elements associated with particular traits or diseases share similar patterns of transcriptional activity. Accordingly, shared transcriptional activity (coexpression) may help prioritise loci associated with a given trait, and help to identify underlying biological processes. Using cap analysis of gene expression (CAGE) profiles of promoter- and enhancer-derived RNAs across 1824 human samples, we have analysed coexpression of RNAs originating from trait-associated regulatory regions using a novel quantitative method (network density analysis; NDA). For most traits studied, phenotype-associated variants in regulatory regions were linked to tightly-coexpressed networks that are likely to share important functional characteristics. Coexpression provides a new signal, independent of phenotype association, to enable fine mapping of causative variants. The NDA coexpression approach identifies new genetic variants associated with specific traits, including an association between the regulation of the OCT1 cation transporter and genetic variants underlying circulating cholesterol levels. NDA strongly implicates particular cell types and tissues in disease pathogenesis. For example, distinct groupings of disease-associated regulatory regions implicate two distinct biological processes in the pathogenesis of ulcerative colitis; a further two separate processes are implicated in Crohn's disease. Thus, our functional analysis of genetic predisposition to disease defines new distinct disease endotypes. We predict that patients with a preponderance of susceptibility variants in each group are likely to respond differently to pharmacological therapy. Together, these findings enable a deeper biological understanding of the causal basis of complex traits.

Restricted Presence of POU6F2 in Human Corneal Endothelial Cells Uncovered by Extension of the Promoter-level Expression Atlas

Corneal endothelial cells (CECs) are essential for maintaining the clarity of the cornea. Because CECs have limited proliferative ability, interest is growing in their potentially therapeutic regeneration from pluripotent stem cells. However, the molecular mechanisms of human CEC differentiation remain largely unknown. To determine the key regulators of CEC characteristics, here we generated a comprehensive promoter-level expression profile of human CECs, using cap analysis of gene expression (CAGE) with a single molecule sequencer. Integration with the FANTOM5 promoter-level expression atlas, which includes transcriptome profiles of various human tissues and cells, enabled us to identify 45 promoters at 28 gene loci that are specifically expressed in CECs. We further discovered that the expression of transcription factor POU class 6 homeobox 2 (POU6F2) is restricted to CECs, and upregulated during human CEC differentiation, suggesting that POU6F2 is pivotal to terminal differentiation of CECs. These CEC-specific promoters would be useful for the assessment of fully differentiated CECs derived from pluripotent stem cells. These findings promote the development of corneal regenerative medicine.

•

We comprehensively profiled promoter-level expression of human corneal endothelial cells.

•

Integrative transcriptome analysis identified 28 corneal endothelial cell-specific marker genes.

•

POU6F2 expression is restricted to corneal endothelial cells, and upregulated during differentiation.

Corneal endothelial cells (CECs) are essential for maintaining corneal transparency. Owing to the high prevalence of corneal endothelial disorders, CECs are important targets in regenerative medicine. However, it has been difficult to evaluate the final CEC products owing to the lack of appropriate CEC-specific markers. In this study, we identified 28 CEC-specific marker genes by integrative transcriptome analysis. One gene of particular interest, POU6F2, is expressed almost exclusively in CECs, and upregulated during differentiation. These markers would be useful for the assessment of CECs derived from pluripotent stem cells, and this study will facilitate the translation of corneal regenerative medicine.

Transcription start site profiling of 15 anatomical regions of the Macaca mulatta central nervous system

Rhesus macaque was the second non-human primate whose genome has been fully sequenced and is one of the most used model organisms to study human biology and disease, thanks to the close evolutionary relationship between the two species. But compared to human, where several previously unknown RNAs have been uncovered, the macaque transcriptome is less studied. Publicly available RNA expression resources for macaque are limited, even for brain, which is highly relevant to study human cognitive abilities. In an effort to complement those resources, FANTOM5 profiled 15 distinct anatomical regions of the aged macaque central nervous system using Cap Analysis of Gene Expression, a high-resolution, annotation-independent technology that allows monitoring of transcription initiation events with high accuracy. We identified 25,869 CAGE peaks, representing bona fide promoters. For each peak we provide detailed annotation, expanding the landscape of ‘known’ macaque genes, and we show concrete examples on how to use the resulting data. We believe this data represents a useful resource to understand the central nervous system in macaque.

Promoter-level transcriptome in primary lesions of endometrial cancer identified biomarkers associated with lymph node metastasis

For endometrial cancer patients, lymphadenectomy is recommended to exclude rarely metastasized cancer cells. This procedure is performed even in patients with low risk of recurrence despite the risk of complications such as lymphedema. A method to accurately identify cases with no lymph node metastases (LN-) before lymphadenectomy is therefore highly required. We approached this clinical problem by examining primary lesions of endometrial cancers with CAGE (Cap Analysis Gene Expression), which quantifies promoter-level expression across the genome. Fourteen profiles delineated distinct transcriptional networks between LN + and LN- cases, within those classified as having the low or intermediate risk of recurrence. Subsequent quantitative reverse transcription polymerase chain reaction (qRT-PCR) analyses of 115 primary tumors showed SEMA3D mRNA and TACC2 isoforms expressed through a novel promoter as promising biomarkers with high accuracy (area under the receiver operating characteristic curve, 0.929) when used in combination. Our high-resolution transcriptome provided evidence of distinct molecular profiles underlying LN + /LN- status in endometrial cancers, raising the possibility of preoperative diagnosis to reduce unnecessary operations in patients with minimum recurrence risk.

Systematic analysis of transcription start sites in avian development

Cap Analysis of Gene Expression (CAGE) in combination with single-molecule sequencing technology allows precision mapping of transcription start sites (TSSs) and genome-wide capture of promoter activities in differentiated and steady state cell populations. Much less is known about whether TSS profiling can characterize diverse and non-steady state cell populations, such as the approximately 400 transitory and heterogeneous cell types that arise during ontogeny of vertebrate animals. To gain such insight, we used the chick model and performed CAGE-based TSS analysis on embryonic samples covering the full 3-week developmental period. In total, 31,863 robust TSS peaks (>1 tag per million [TPM]) were mapped to the latest chicken genome assembly, of which 34% to 46% were active in any given developmental stage. ZENBU, a web-based, open-source platform, was used for interactive data exploration. TSSs of genes critical for lineage differentiation could be precisely mapped and their activities tracked throughout development, suggesting that non-steady state and heterogeneous cell populations are amenable to CAGE-based transcriptional analysis. Our study also uncovered a large set of extremely stable housekeeping TSSs and many novel stage-specific ones. We furthermore demonstrated that TSS mapping could expedite motif-based promoter analysis for regulatory modules associated with stage-specific and housekeeping genes. Finally, using Brachyury as an example, we provide evidence that precise TSS mapping in combination with Clustered Regularly Interspaced Short Palindromic Repeat (CRISPR)-on technology enables us, for the first time, to efficiently target endogenous avian genes for transcriptional activation. Taken together, our results represent the first report of genome-wide TSS mapping in birds and the first systematic developmental TSS analysis in any amniote species (birds and mammals). By facilitating promoter-based molecular analysis and genetic manipulation, our work also underscores the value of avian models in unravelling the complex regulatory mechanism of cell lineage specification during amniote development.

Integrative CAGE and DNA Methylation Profiling Identify Epigenetically Regulated Genes in NSCLC

Lung cancer is the leading cause of cancer-related deaths worldwide. The majority of cancer driver mutations have been identified; however, relevant epigenetic regulation involved in tumorigenesis has only been fragmentarily analyzed. Epigenetically regulated genes have a great theranostic potential, especially in tumors with no apparent driver mutations. Here, epigenetically regulated genes were identified in lung cancer by an integrative analysis of promoter-level expression profiles from Cap Analysis of Gene Expression (CAGE) of 16 non-small cell lung cancer (NSCLC) cell lines and 16 normal lung primary cell specimens with DNA methylation data of 69 NSCLC cell lines and 6 normal lung epithelial cells. A core set of 49 coding genes and 10 long noncoding RNAs (lncRNA), which are upregulated in NSCLC cell lines due to promoter hypomethylation, was uncovered. Twenty-two epigenetically regulated genes were validated (upregulated genes with hypomethylated promoters) in the adenocarcinoma and squamous cell cancer subtypes of lung cancer using The Cancer Genome Atlas data. Furthermore, it was demonstrated that multiple copies of the REP522 DNA repeat family are prominently upregulated due to hypomethylation in NSCLC cell lines, which leads to cancer-specific expression of lncRNAs, such as RP1-90G24.10, AL022344.4, and PCAT7. Finally, Myeloma Overexpressed (MYEOV) was identified as the most promising candidate. Functional studies demonstrated that MYEOV promotes cell proliferation, survival, and invasion. Moreover, high MYEOV expression levels were associated with poor prognosis.Implications: This report identifies a robust list of 22 candidate driver genes that are epigenetically regulated in lung cancer; such genes may complement the known mutational drivers.Visual Overview: http://mcr.aacrjournals.org/content/molcanres/15/10/1354/F1.large.jpg Mol Cancer Res; 15(10); 1354-65. ©2017 AACR.