Cap Analysis Of Gene Expression Research Articles

Enhancers regulate genes linked to severe and mild childhood asthma

BackgroundChildren with severe asthma suffer from recurrent symptoms and impaired quality of life despite advanced treatment. Underlying causes of severe asthma are not completely understood, although genetic mechanisms are known to be important. ObjectiveThe aim of this study was to identify gene regulatory enhancers in leukocytes, to describe the role of these enhancers in regulating genes related to severe and mild asthma in children, and to identify known asthma-related SNPs situated in proximity to enhancers. MethodsGene enhancers were identified and expression of enhancers and genes were measured by Cap Analysis Gene Expression (CAGE) data from peripheral blood leukocytes from children with severe asthma (n = 13), mild asthma (n = 15), and age-matched controls (n = 9). ResultsFrom a comprehensive set of 8,289 identified enhancers, we further defined a robust sub-set of the high-confidence and most highly expressed 4,738 enhancers. Known single nucleotide polymorphisms, SNPs, related to asthma coincided with enhancers in general as well as with specific enhancer-gene interactions. Blocks of enhancer clusters were associated with genes including TGF-beta, PPAR and IL-11 signaling as well as genes related to vitamin A and D metabolism. A signature of 91 enhancers distinguished between children with severe and mild asthma as well as controls. ConclusionsGene regulatory enhancers were identified in leukocytes with potential roles related to severe and mild asthma in children. Enhancers hosting known SNPs give the opportunity to formulate mechanistic hypotheses about the functions of these SNPs.

Transcriptome analysis of long non-coding RNAs in Mycobacterium avium complex-infected macrophages.

Mycobacterium avium complex (MAC) is a non-tuberculous mycobacterium widely distributed in the environment. Even though MAC infection is increasing in older women and immunocompromised patients, to our knowledge there has been no comprehensive analysis of the MAC-infected host-cell transcriptome-and particularly of long non-coding RNAs (lncRNAs). By using in vitro-cultured primary mouse bone-marrow-derived macrophages (BMDMs) and Cap analysis of gene expression, we analyzed the transcriptional and kinetic landscape of macrophage genes, with a focus on lncRNAs, during MAC infection. MAC infection of macrophages induced the expression of immune/inflammatory response genes and other genes similar to those involved in M1 macrophage activation, consistent with previous reports, although Nos2 (M1 activation) and Arg1 (M2 activation) had distinct expression profiles. We identified 31 upregulated and 30 downregulated lncRNA promoters corresponding respectively to 18 and 26 lncRNAs. Upregulated lncRNAs were clustered into two groups-early and late upregulated-predicted to be associated with immune activation and the immune response to infection, respectively. Furthermore, an Ingenuity Pathway Analysis revealed canonical pathways and upstream transcription regulators associated with differentially expressed lncRNAs. Several differentially expressed lncRNAs reported elsewhere underwent expressional changes upon M1 or M2 preactivation and subsequent MAC infection. Finally, we showed that expressional change of lncRNAs in MAC-infected BMDMs was mediated by toll-like receptor 2, although there may be other mechanisms that sense MAC infection. We identified differentially expressed lncRNAs in MAC-infected BMDMs, revealing diverse features that imply the distinct roles of these lncRNAs in MAC infection and macrophage polarization.

Abstract LB274: Nuclear transporter HIKESHI is a novel independent poor overall survival for pancreatic cancer; targeting can inhibit pancreatic cancer growth via HSP70 suppression

Abstract HIKESHI, under heat stress conditions, has been reported to transport molecular chaperone HSP70 into the nucleus to prevent cancer cell death from heat stress. The tumor-supporting roles of molecular chaperones in various solid tumors are explored extensively. Yet, the sole targeting of the molecular chaperones could not significantly impact anticancer treatment. Thus, we hypothesized that the specific transporter protein HIKESHI might improve this outcome. However, understanding of the clinical relevance of HIKESHI in pancreatic cancer remained unclear, and few studies addressed the HIKESHI role in cancer and whether targeting can inhibit pancreatic cancer cell viability and the other HSPs. Here, we first investigated the expression pattern of the HIKESHI protein in surgically resected pancreatic cancer samples (n=142) and analyzed the clinical relevance of the HIKESHI. We found that the HIKESHI protein was overexpressed in cancer compared to the adjacent normal tissue, and positive expression of cancer HIKESHI was related to poor overall survival. In addition, HIKESHI-positive expression in cancer was identified as an independent prognostic marker for overall survival. Next, using shRNA, we established the stable HIKESHI-suppressed pancreatic cancer cell line (PANC-1). HIKESHI suppression inhibited cell growth and colony formation ability in vitro. The expression and nuclear translocation of HSP70 and HSP110/105 were suppressed in HIKESHI-suppressed cells compared to the control shRNA cells despite no heat stress conditions. Furthermore, we found that the tumor growth, HIKESHI, and Ki67 expression in HIKESHI-suppressed tumors were inhibited compared to the control shRNA tumors. Next, we performed Cap Analysis of Gene Expression and IP-MS analysis to discover the potential downstream and novel interactive partner of the HIKESHI protein. The extensive mechanism of these analyses is under investigation. Our findings suggest that targeting HIKESHI may be a promising therapeutic strategy against refractory pancreatic cancer. Citation Format: Bilguun Erkhem-Ochir, Takehiko Yokobori, Gendensuren Dorjkhorloo, Haruka Okami, Takaomi Seki, Norifumi Harimoto, Ryo Muranushi, Kouki Hoshino, Kei Hagiwara, Norihiro Ishii, Mariko Tsukagoshi, Kenichiro Araki, Hiroshi Saeki, Ken Shirabe. Nuclear transporter HIKESHI is a novel independent poor overall survival for pancreatic cancer; targeting can inhibit pancreatic cancer growth via HSP70 suppression [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 2 (Late-Breaking, Clinical Trial, and Invited Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(7_Suppl):Abstract nr LB274.

RNAirport: a deep neural network-based database characterizing representative gene models in plants

A 5′-leader, known initially as the 5′-untranslated region, contains multiple isoforms due to alternative splicing (aS) and alternative transcription start site (aTSS). Therefore, a representative 5′-leader is demanded to examine the embedded RNA regulatory elements in controlling translation efficiency. Here, we develop a ranking algorithm and a deep-learning model to annotate representative 5′-leaders for five plant species. We rank the intra-sample and inter-sample frequency of aS-mediated transcript isoforms using the Kruskal–Wallis test-based algorithm and identify the representative aS-5′-leader. To further assign a representative 5′-end, we train the deep-learning model 5′leaderP to learn aTSS-mediated 5′-end distribution patterns from cap-analysis gene expression data. The model accurately predicts the 5′-end, confirmed experimentally in Arabidopsis and rice. The representative 5′-leader-contained gene models and 5′leaderP can be accessed at RNAirport (http://www.rnairport.com/leader5P/). The Stage 1 annotation of 5′-leader records 5′-leader diversity and will pave the way to Ribo-Seq open-reading frame annotation, identical to the project recently initiated by human GENCODE.

A single-cell genomic strategy for alternative transcript start sites identification.

Alternative transcription start sites (TSSs) usage plays a critical role in gene transcription regulation in mammals. However, precisely identifying alternative TSSs remains challenging at the genome-wide level. We report a single-cell genomic technology for alternative TSSs annotation and cell heterogeneity detection. In the method, we utilize Fluidigm C1 system to capture individual cells of interest, SMARTer cDNA synthesis kit to recover full-length cDNAs, then dual priming oligonucleotide system to specifically enrich TSSs for genomic analysis. We apply this method to a genome-wide study of alternative TSSs identification in two different IFN-β stimulated mouse embryonic fibroblasts (MEFs). The data clearly discriminate two IFN-β stimulated MEFs. Moreover, our results indicate 81% expressed genes in these two cell types containing multiple TSSs, which is much higher than previous predictions based on Cap-Analysis Gene Expression (CAGE) (58%) or empirical determination (54%) in various cell types. This indicates that alternative TSSs are more pervasive than expected and implies our strategy could position them at an unprecedented sensitivity. It would be helpful for elucidating their biological insights in future.

NTRK2 expression in gastrointestinal stromal tumors with a special emphasis on the clinicopathological and prognostic impacts

Gastrointestinal stromal tumors (GISTs) are typically characterized by activating mutations of the KIT proto-oncogene receptor tyrosine kinase (KIT) or platelet-derived growth factor receptor alpha (PDGFRA). Recently, the neurotrophic tyrosine receptor kinase (NTRK) fusion was reported in a small subset of wild-type GIST. We examined trk IHC and NTRK gene expressions in GIST. Pan-trk immunohistochemistry (IHC) was positive in 25 (all 16 duodenal and 9 out of 16 small intestinal GISTs) of 139 cases, and all pan-trk positive cases showed diffuse and strong expression of c-kit. Interestingly, all of these cases showed only trkB but not trkA/trkC expression. Cap analysis of gene expression (CAGE) analysis identified increased number of genes whose promoters were activated in pan-trk/trkB positive GISTs. Imbalanced expression of NTRK2, which suggests the presence of NTRK2 fusion, was not observed in any of trkB positive GISTs, despite higher mRNA expression. TrkB expression was found in duodenal GISTs and more than half of small intestinal GISTs, and this subset of cases showed poor prognosis. However, there was not clear difference in clinical outcomes according to the trkB expression status in small intestinal GISTs. These findings may provide a possible hypothesis for trkB overexpression contributing to the tumorigenesis and aggressive clinical outcome in GISTs of duodenal origin.

A budding yeast CAGE dataset comprising two cell types.

The budding yeast Saccharomyces cerevisiae is an excellent model organism for studying chromatin regulation with high-resolution genome-wide analyses. Since newly generated genome-wide data are often compared with publicly available datasets, expanding our dataset repertoire will be beneficial for the field. Information on transcription start sites (TSSs) determined at base pair resolution is essential for elucidating mechanisms of transcription and related chromatin regulation, yet no datasets that cover two different cell types are available. Here, we present a CAGE (cap analysis of gene expression) dataset for a-cells and α-cells grown in defined and rich media. Cell type-specific genes were differentially expressed as expected, ensuring the reliability of the data. Some of the differentially expressed TSSs were medium-specific or detected due to unrecognized chromosome rearrangement. By comparing the CAGE data with a high-resolution nucleosome map, major TSSs were primarily found in +1 nucleosomes, with a peak approximately 30 bp from the promoter-proximal end of the nucleosome. The dataset is available at DDBJ/GEA.

Core promoterome of barley embryo

Precise localization and dissection of gene promoters are key to understanding transcriptional gene regulation and to successful bioengineering applications. The core RNA polymerase II initiation machinery is highly conserved among eukaryotes, leading to a general expectation of equivalent underlying mechanisms. Still, less is known about promoters in the plant kingdom. In this study, we employed cap analysis of gene expression (CAGE) at three embryonic developmental stages in barley to accurately map, annotate, and quantify transcription initiation events. Unsupervised discovery of de novo sequence clusters grouped promoters based on characteristic initiator and position-specific core-promoter motifs. This grouping was complemented by the annotation of transcription factor binding site (TFBS) motifs. Integration with genome-wide epigenomic data sets and gene ontology (GO) enrichment analysis further delineated the chromatin environments and functional roles of genes associated with distinct promoter categories. The TATA-box presence governs all features explored, supporting the general model of two separate genomic regulatory environments. We describe the extent and implications of alternative transcription initiation events, including those that are specific to developmental stages, which can affect the protein sequence or the presence of regions that regulate translation. The generated promoterome dataset provides a valuable genomic resource for enhancing the functional annotation of the barley genome. It also offers insights into the transcriptional regulation of individual genes and presents opportunities for the informed manipulation of promoter architecture, with the aim of enhancing traits of agronomic importance.

OVCH1 Antisense RNA 1 is differentially expressed between non-frail and frail old adults.

While some old adults stay healthy and non-frail up to late in life, others experience multimorbidity and frailty often accompanied by a pro-inflammatory state. The underlying molecular mechanisms for those differences are still obscure. Here, we used gene expression analysis to understand the molecular underpinning between non-frail and frail individuals in old age. Twenty-four adults (50% non-frail and 50% frail) from InCHIANTI study were included. Total RNA extracted from whole blood was analyzed by Cap Analysis of Gene Expression (CAGE). CAGE identified transcription start site (TSS) and active enhancer regions. We identified a set of differentially expressed (DE) TSS and enhancer between non-frail and frail and male and female participants. Several DE TSSs were annotated as lncRNA (XIST and TTTY14) and antisense RNAs (ZFX-AS1 and OVCH1 Antisense RNA 1). The promoter region chr6:366,786,54-366,787,97;+ was DE and overlapping the longevity CDKN1A gene. GWAS-LD enrichment analysis identifies overlapping LD-blocks with the DE regions with reported traits in GWAS catalog (isovolumetric relaxation time and urinary tract infection frequency). Furthermore, we used weighted gene co-expression network analysis (WGCNA) to identify changes of gene expression associated with clinical traits and identify key gene modules. We performed functional enrichment analysis of the gene modules with significant trait/module correlation. One gene module is showing a very distinct pattern in hub genes. Glycogen Phosphorylase L (PYGL) was the top ranked hub gene between non-frail and frail. We predicted transcription factor binding sites (TFBS) and motif activity. TF involved in age-related pathways (e.g., FOXO3 and MYC) shows different expression patterns between non-frail and frail participants. Expanding the study of OVCH1 Antisense RNA 1 and PYGL may help understand the mechanisms leading to loss of homeostasis that ultimately causes frailty.

Age-Related Gene and Protein Expression in Mouse Mandibular Condyle Analyzed by Cap Analysis of Gene Expression and Immunohistochemistry

Introduction: Aging, an inevitable physiological process, leads to morphological and histological degenerative changes in the mandibular condylar cartilage (MCC); however, the molecular mechanism has not yet been elucidated, and little information is available on age-related factors. Therefore, this study was designed to identify age-related factors by investigating the age-related differentially expressed genes (DEGs) and localization of their translated protein expression in the mandibular condyle. Methods: Mandibular condyles were collected from 10- and 50-week-old mice. Total RNA was extracted from the samples and then analyzed using cap analysis of gene expression (CAGE) to identify age-related DEGs. Gene ontology (GO) enrichment analysis was performed to determine which biological processes were most affected by aging in terms of gene expression using Metascape. The mandibular condyle samples were processed for histology to investigate morphological changes caused by aging and for immunohistochemistry to localize the protein expression encoded by age-related genes identified with CAGE. Semi-quantitative immunohistochemistry was performed to assess age-related extracellular matrix (ECM) protein levels in the MCC. The histological sections were also used for Alcian blue histochemistry to detect glycosaminoglycans (GAGs). Results: GO enrichment analysis revealed that the genes related to “extracellular matrix organization,” including Acan, Col1a1, Col1a2, Col2a1, Mmp3, Mmp9, and Mmp13, were most differentially expressed in the aged mandibular condyle. Among these seven genes, Mmp3 was upregulated, and the others were downregulated with aging. Histological examination showed the age-related morphological and histological changes in the MCC. Immunohistochemical investigation showed the localization of matrix metalloproteinases (MMPs)-3, -9, and -13 and their substrate proteins, aggrecan, type I collagen, and type II collagen, in the mandibular condyle at 10 and 50 weeks, indicating different localizations between the young and the aged. In the aged MCC, semi-quantitative immunohistochemistry showed a significant decrease in the aggrecan protein level, and Alcian blue histochemistry showed a decrease in GAGs. Conclusion: MMP-3, MMP-9, and MMP-13 contribute to the remodeling of the ECM of the MCC and subchondral bone during aging by degrading ECM proteins at specific times and sites under the regulation of their production and secretion.

Cap Analysis of Gene Expression Clarifies Transcriptomic Divergence Within Monozygotic Twin Pairs.

Phenotypic variation is the result of gene expression based on complex interaction between genetic and environmental factors. It is well known that genetic and environmental factors influence gene expression, but our understanding of their relative importance remains limited. To obtain a hint for the understanding of their contributions, we took advantage of monozygotic twins, as they share genetic and shared environmental factors but differ in nonshared factors, such as environmental differences and stochastic factors. In this study, we performed cap analysis of gene expression on three pairs of twins and clustered each individual based on their expression profiles of annotated genes. The dendrogram of annotated gene transcripts showed a monophyletic clade for each twin pair. We also analyzed the expression of retrotransposons, such as human endogenous retroviruses (HERVs) and long interspersed nuclear elements (LINEs), given their abundance in the genome. Clustering analyses demonstrated that HERV and LINE expression diverged even within monozygotic twin pairs. Thus, HERVs and LINEs are more susceptible to nonshared factors than annotated genes. Motif analysis of differentially expressed annotated genes suggests that specificity protein/Krüppel-like factor family transcription factors are involved in the expression divergence of annotated gene influenced by nonshared factors. Collectively, our findings suggest that expressions of annotated genes and retrotransposons are differently regulated, and that the expression of retrotransposons is more susceptible to nonshared factors than annotated genes.

Prioritizing genes associated with brain disorders by leveraging enhancer-promoter interactions in diverse neural cells and tissues

BackgroundPrioritizing genes that underlie complex brain disorders poses a considerable challenge. Despite previous studies have found that they shared symptoms and heterogeneity, it remained difficult to systematically identify the risk genes associated with them.MethodsBy using the CAGE (Cap Analysis of Gene Expression) read alignment files for 439 human cell and tissue types (including primary cells, tissues and cell lines) from FANTOM5 project, we predicted enhancer-promoter interactions (EPIs) of 439 cell and tissue types in human, and examined their reliability. Then we evaluated the genetic heritability of 17 diverse brain disorders and behavioral-cognitive phenotypes in each neural cell type, brain region, and developmental stage. Furthermore, we prioritized genes associated with brain disorders and phenotypes by leveraging the EPIs in each neural cell and tissue type, and analyzed their pleiotropy and functionality for different categories of disorders and phenotypes. Finally, we characterized the spatiotemporal expression dynamics of these associated genes in cells and tissues.ResultsWe found that identified EPIs showed activity specificity and network aggregation in cell and tissue types, and enriched TF binding in neural cells played key roles in synaptic plasticity and nerve cell development, i.e., EGR1 and SOX family. We also discovered that most neurological disorders exhibit heritability enrichment in neural stem cells and astrocytes, while psychiatric disorders and behavioral-cognitive phenotypes exhibit enrichment in neurons. Furthermore, our identified genes recapitulated well-known risk genes, which exhibited widespread pleiotropy between psychiatric disorders and behavioral-cognitive phenotypes (i.e., FOXP2), and indicated expression specificity in neural cell types, brain regions, and developmental stages associated with disorders and phenotypes. Importantly, we showed the potential associations of brain disorders with brain regions and developmental stages that have not been well studied.ConclusionsOverall, our study characterized the gene-enhancer regulatory networks and genetic mechanisms in the human neural cells and tissues, and illustrated the value of reanalysis of publicly available genomic datasets.

LINC01137/miR-186-5p/WWOX: a novel axis identified from WWOX-related RNA interactome in bladder cancer.

Introduction: The discovery of non-coding RNA (ncRNA) dates back to the pre-genomics era, but the progress in this field is still dynamic and leverages current post-genomics solutions. WWOX is a global gene expression modulator that is scarcely investigated for its role in regulating cancer-related ncRNAs. In bladder cancer (BLCA), the link between WWOX and ncRNA remains unexplored. The description of AP-2α and AP-2γ transcription factors, known as WWOX-interacting proteins, is more commonplace regarding ncRNA but still merits investigation. Therefore, this in vitro and in silico study aimed to construct an ncRNA-containing network with WWOX/AP-2 and to investigate the most relevant observation in the context of BLCA cell lines and patients. Methods: RT-112, HT-1376, and CAL-29 cell lines were subjected to two stable lentiviral transductions. High-throughput sequencing of cellular variants (deposited in the Gene Expression Omnibus database under the GSE193659 record) enabled the investigation of WWOX/AP-2-dependent differences using various bioinformatics tools (e.g., limma-voom, FactoMineR, multiple Support Vector Machine Recursive Feature Elimination (mSVM-RFE), miRDB, Arena-Idb, ncFANs, RNAhybrid, TargetScan, Protein Annotation Through Evolutionary Relationships (PANTHER), Gene Transcription Regulation Database (GTRD), or Evaluate Cutpoints) and repositories such as The Cancer Genome Atlas (TCGA) and Cancer Cell Line Encyclopedia. The most relevant observations from cap analysis gene expression sequencing (CAGE-seq) were confirmed using real-time PCR, whereas TCGA data were validated using the GSE31684 cohort. Results: The first stage of the whole study justified focusing solely on WWOX rather than on WWOX combined with AP-2α/γ. The most relevant observation of the developed ncRNA-containing network was LINC01137, i.e., long non-coding RNAs (lncRNAs) that unraveled the core network containing UPF1, ZC3H12A, LINC01137, WWOX, and miR-186-5p, the last three being a novel lncRNA/miRNA/mRNA axis. Patients' data confirmed the LINC01137/miR-186-5p/WWOX relationship and provided a set of dependent genes (i.e., KRT18, HES1, VCP, FTH1, IFITM3, RAB34, and CLU). Together with the core network, the gene set was subjected to survival analysis for both TCGA-BLCA and GSE31684 patients, which indicated that the increased expression of WWOX or LINC01137 is favorable, similar to their combination with each other (WWOX↑ and LINC01137↑) or with MIR186 (WWOX↑/LINC01137↑ but MIR186↓). Conclusion: WWOX is implicated in the positive feedback loop with LINC01137 that sponges WWOX-targeting miR-186-5p. This novel WWOX-containing lncRNA/miRNA/mRNA axis should be further investigated to depict its relationships in a broader context, which could contribute to BLCA research and treatment.

Global Analysis of Transcription Start Sites and Enhancers in Endometrial Stromal Cells and Differences Associated with Endometriosis.

Identifying tissue-specific molecular signatures of active regulatory elements is critical to understanding gene regulatory mechanisms. In this study, transcription start sites (TSS) and enhancers were identified using Cap analysis of gene expression (CAGE) across endometrial stromal cell (ESC) samples obtained from women with (n = 4) and without endometriosis (n = 4). ESC TSSs and enhancers were compared to those reported in other tissue and cell types in FANTOM5 and were integrated with RNA-seq and ATAC-seq data from the same samples for regulatory activity and network analyses. CAGE tag count differences between women with and without endometriosis were statistically tested and tags within close proximity to genetic variants associated with endometriosis risk were identified. Over 90% of tag clusters mapping to promoters were observed in cells and tissues in FANTOM5. However, some potential cell-type-specific promoters and enhancers were also observed. Regions of open chromatin identified using ATAC-seq provided further evidence of the active transcriptional regions identified by CAGE. Despite the small sample number, there was evidence of differences associated with endometriosis at 210 consensus clusters, including IGFBP5, CALD1 and OXTR. ESC TSSs were also located within loci associated with endometriosis risk from genome-wide association studies. This study provides novel evidence of transcriptional differences in endometrial stromal cells associated with endometriosis and provides a valuable cell-type specific resource of active TSSs and enhancers in endometrial stromal cells.

Identification and characterization of transcribed enhancers during cerebellar development through enhancer RNA analysis

BackgroundThe development of the brain requires precise coordination of molecular processes across many cell-types. Underpinning these events are gene expression programs which require intricate regulation by non-coding regulatory sequences known as enhancers. In the context of the developing brain, transcribed enhancers (TEs) regulate temporally-specific expression of genes critical for cell identity and differentiation. Transcription of non-coding RNAs at active enhancer sequences, known as enhancer RNAs (eRNAs), is tightly associated with enhancer activity and has been correlated with target gene expression. TEs have been characterized in a multitude of developing tissues, however their regulatory role has yet to be described in the context of embryonic and early postnatal brain development. In this study, eRNA transcription was analyzed to identify TEs active during cerebellar development, as a proxy for the developing brain. Cap Analysis of Gene Expression followed by sequencing (CAGE-seq) was conducted at 12 stages throughout embryonic and early postnatal cerebellar development.ResultsTemporal analysis of eRNA transcription identified clusters of TEs that peak in activity during either embryonic or postnatal times, highlighting their importance for temporally specific developmental events. Functional analysis of putative target genes identified molecular mechanisms under TE regulation revealing that TEs regulate genes involved in biological processes specific to neurons. We validate enhancer activity using in situ hybridization of eRNA expression from TEs predicted to regulate Nfib, a gene critical for cerebellar granule cell differentiation.ConclusionsThe results of this analysis provide a valuable dataset for the identification of cerebellar enhancers and provide insight into the molecular mechanisms critical for brain development under TE regulation. This dataset is shared with the community through an online resource (https://goldowitzlab.shinyapps.io/trans-enh-app/).

Improving the annotation of the cattle genome by annotating transcription start sites in a diverse set of tissues and populations using Cap Analysis Gene Expression sequencing.

Understanding the genomic control of tissue-specific gene expression and regulation can help to inform the application of genomic technologies in farm animal breeding programs. The fine mapping of promoters [transcription start sites (TSS)] and enhancers (divergent amplifying segments of the genome local to TSS) in different populations of cattle across a wide diversity of tissues provides information to locate and understand the genomic drivers of breed- and tissue-specific characteristics. To this aim, we used Cap Analysis Gene Expression (CAGE) sequencing, of 24 different tissues from 3 populations of cattle, to define TSS and their coexpressed short-range enhancers (<1 kb) in the ARS-UCD1.2_Btau5.0.1Y reference genome (1000bulls run9) and analyzed tissue and population specificity of expressed promoters. We identified 51,295 TSS and 2,328 TSS-Enhancer regions shared across the 3 populations (dairy, beef-dairy cross, and Canadian Kinsella composite cattle from 2 individuals, 1 of each sex, per population). Cross-species comparative analysis of CAGE data from 7 other species, including sheep, revealed a set of TSS and TSS-Enhancers that were specific to cattle. The CAGE data set will be combined with other transcriptomic information for the same tissues to create a new high-resolution map of transcript diversity across tissues and populations in cattle for the BovReg project. Here we provide the CAGE data set and annotation tracks for TSS and TSS-Enhancers in the cattle genome. This new annotation information will improve our understanding of the drivers of gene expression and regulation in cattle and help to inform the application of genomic technologies in breeding programs.

Connectome and regulatory hubs of CAGE highly active enhancers

Evidence indicates that enhancers are transcriptionally active. Herein, we investigated transcriptionally active enhancers by using cap analysis of gene expression (CAGE) combined with epigenetic marks and chromatin interactions. We identified CAGE-tag highly active (CHA) enhancers as distant regulatory elements with CAGE-tag ≥ 90th percentile and overlapping with H3K27ac peaks (4.5% of enhancers). CHA enhancers were conserved between mouse and man and were independent from super-enhancers in predicting cell identity with lower P-values. CHA enhancers had increased open chromatin and a higher recruitment of cell-specific transcription factors as well as molecules involved in 3D genome interactions. HiChIP analysis of enhancer-promoter looping indicated that CHA enhancers had a higher density of anchor loops when compared to regular enhancers. A subset of CHA enhancers and promoters characterized by a high density of chromatin loops and forming hub regulatory units were connected to the promoter of immediate early response genes, genes involved in cancer and encoding for transcription factors. Promoter of genes within hub CHA regulatory units were less likely to be paused. CHA enhancers were enriched in gene variants associated with autoimmune disorders and had looping with causal candidate genes as revealed by Mendelian randomization. Hence, CHA enhancers form a dense hierarchical network of chromatin interactions between regulatory elements and genes involved in cell identity and disorders.

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

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. 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. 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. 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.

Long noncoding RNAs are substrates for cytoplasmic capping enzyme.

Cytoplasmic capping returns a cap to specific mRNAs, thus protecting uncapped RNAs from decay. Prior to the identification of cytoplasmic capping, uncapped mRNAs were thought to be degraded. Here, we test whether long noncoding RNAs (lncRNAs) are substrates of the cytoplasmic capping enzyme (cCE). The subcellular localisation of 14 lncRNAs associated with sarcomas were examined in U2OS osteosarcoma cells. We used 5' rapid amplification of cDNA ends (RACE) to assay uncapped forms of these lncRNAs. Inhibiting cytoplasmic capping elevated uncapped forms of selected lncRNAs indicating a plausible role of cCE in targeting them. Analysis of published cap analysis of gene expression (CAGE) data shows increased prevalence of certain 5'-RACE cloned sequences, suggesting that these uncapped lncRNAs are targets of cytoplasmic capping.

An atlas of transcribed human cardiac promoters and enhancers reveals an important role of regulatory elements in heart failure

A deeper knowledge of the dynamic transcriptional activity of promoters and enhancers is needed to improve mechanistic understanding of the pathogenesis of heart failure and heart diseases. In this study, we used cap analysis of gene expression (CAGE) to identify and quantify the activity of transcribed regulatory elements (TREs) in the four cardiac chambers of 21 healthy and ten failing adult human hearts. We identified 17,668 promoters and 14,920 enhancers associated with the expression of 14,519 genes. We showed how these regulatory elements are alternatively transcribed in different heart regions, in healthy versus failing hearts and in ischemic versus non-ischemic heart failure samples. Cardiac-disease-related single-nucleotide polymorphisms (SNPs) appeared to be enriched in TREs, potentially affecting the allele-specific transcription factor binding. To conclude, our open-source heart CAGE atlas will serve the cardiovascular community in improving the understanding of the role of the cardiac gene regulatory networks in cardiovascular disease and therapy. Deviatiiarov, Gams et al. provide an atlas of the actively transcribed regulatory elements (promoters and enhancers) in the human heart by performing cap analysis of gene expression (CAGE) on 21 healthy donor hearts and ten failing hearts from ischemic and non-ischemic cardiomyopathy. The data show how the alternative use of regulatory elements modulates the differential gene expression across different heart chambers, disease states and disease types.