Recounting the FANTOM CAGE-Associated Transcriptome

The annotation of GBA1 has been concealed by its protein-coding pseudogene GBAP1

Mutations in GBA1 cause Gaucher disease and are the most important genetic risk factor for Parkinson's disease. However, analysis of transcription at this locus is complicated by its highly homologous pseudogene, GBAP1. We show that >50% of short RNA-sequencing reads mapping to GBA1 also map to GBAP1. Thus, we used long-read RNA sequencing in the human brain, which allowed us to accurately quantify expression from both GBA1 and GBAP1. We discovered significant differences in expression compared to short-read data and identify currently unannotated transcripts of both GBA1 and GBAP1. These included protein-coding transcripts from both genes that were translated in human brain, but without the known lysosomal function-yet accounting for almost a third of transcription. Analyzing brain-specific cell types using long-read and single-nucleus RNA sequencing revealed region-specific variations in transcript expression. Overall, these findings suggest nonlysosomal roles for GBA1 and GBAP1 with implications for our understanding of the role of GBA1 in health and disease.

Loss of heterozygosity impacts MHC expression on the immune microenvironment in CDK12-mutated prostate cancer

BACKGROUND

In prostate cancer (PCa), well-established biomarkers such as MSI status, TMB high, and PDL1 expression serve as reliable indicators for favorable responses to immunotherapy. Recent studies have suggested a potential association between CDK12 mutations and immunotherapy response; however, the precise mechanisms through which CDK12 mutation may influence immune response remain unclear. A plausible explanation for immune evasion in this subset of CDK12-mutated PCa may be reduced MHC expression.

RESULTS

Using genomic data of CDK12-mutated PCa from 48 primary and 10 metastatic public domain samples and a retrospective cohort of 53 low-intermediate risk primary PCa, we investigated how variation in the expression of the MHC genes affected associated downstream pathways. We classified the patients based on gene expression quartiles of MHC-related genes and categorized the tumors into "High" and "Low" expression levels. CDK12-mutated tumors with higher MHC-expressed pathways were associated with the immune system and elevated PD-L1, IDO1, and TIM3 expression. Consistent with an inflamed tumor microenvironment (TME) phenotype, digital cytometric analyses identified increased CD8 + T cells, B cells, γδ T cells, and M1 Macrophages in this group. In contrast, CDK12-mutated tumors with lower MHC expression exhibited features consistent with an immune cold TME phenotype and immunoediting. Significantly, low MHC expression was also associated with chromosome 6 loss of heterozygosity (LOH) affecting the entire HLA gene cluster. These LOH events were observed in both major clonal and minor subclonal populations of tumor cells. In our retrospective study of 53 primary PCa cases from this Institute, we found a 4% (2/53) prevalence of CDK12 mutations, with the confirmation of this defect in one tumor through Sanger sequencing. In keeping with our analysis of public domain data this tumor exhibited low MHC expression at the RNA level. More extensive studies will be required to determine whether reduced HLA expression is generally associated with primary tumors or is a specific feature of CDK12 mutated PCa.

CONCLUSIONS

These data show that analysis of CDK12 alteration, in the context of MHC expression levels, and LOH status may offer improved predictive value for outcomes in this potentially actionable genomic subgroup of PCa. In addition, these findings highlight the need to explore novel therapeutic strategies to enhance MHC expression in CDK12-defective PCa to improve immunotherapy responses.

Identification of androgen response-related lncRNAs in prostate cancer

BACKGROUND

Long noncoding RNAs (lncRNAs) are RNA molecules with over 200 nucleotides that do not code for proteins, but are known to be widely expressed and have key roles in gene regulation and cellular functions. They are also found to be involved in the onset and development of various cancers, including prostate cancer (PCa). Since PCa are commonly driven by androgen regulated signaling, mainly stimulated pathways, identification and determining the influence of lncRNAs in androgen response is useful and necessary. LncRNAs regulated by the androgen receptor (AR) can serve as potential biomarkers for PCa. In the present study, gene expression data analysis were performed to distinguish lncRNAs related to the androgen response pathway.

METHODS AND RESULTS

We used publicly available RNA-sequencing and ChIP-seq data to identify lncRNAs that are associated with the androgen response pathway. Using Universal Correlation Coefficient (UCC) and Pearson Correlation Coefficient (PCC) analyses, we found 15 lncRNAs that have (a) highly correlated expression with androgen response genes in PCa and are (b) differentially expressed in the setting of treatment with an androgen agonist as well as antagonist compared to controls. Using publicly available ChIP-seq data, we investigated the role of androgen/AR axis in regulating expression of these lncRNAs. We observed AR binding in the promoter regions of 5 lncRNAs (MIR99AHG, DUBR, DRAIC, PVT1, and COLCA1), showing the direct influence of AR on their expression and highlighting their association with the androgen response pathway.

CONCLUSION

By utilizing publicly available multiomics data and by employing in silico methods, we identified five candidate lncRNAs that are involved in the androgen response pathway. These lncRNAs should be investigated as potential biomarkers for PCa.

Pan-cancer genomic analysis shows hemizygous PTEN loss tumors are associated with immune evasion and poor outcome

In tumors, somatic mutations of the PTEN suppressor gene are associated with advanced disease, chemotherapy resistance, and poor survival. PTEN loss of function may occur by inactivating mutation, by deletion, either affecting one copy (hemizygous loss) leading to reduced gene expression or loss of both copies (homozygous) with expression absent. Various murine models have shown that minor reductions in PTEN protein levels strongly influence tumorigenesis. Most PTEN biomarker assays dichotomize PTEN (i.e. presence vs. absence) ignoring the role of one copy loss. We performed a PTEN copy number analysis of 9793 TCGA cases from 30 different tumor types. There were 419 (4.28%) homozygous and 2484 (25.37%) hemizygous PTEN losses. Hemizygous deletions led to reduced PTEN gene expression, accompanied by increased levels of instability and aneuploidy across tumor genomes. Outcome analysis of the pan-cancer cohort showed that losing one copy of PTEN reduced survival to comparable levels as complete loss, and was associated with transcriptomic changes controlling immune response and the tumor microenvironment. Immune cell abundances were significantly altered for PTEN loss, with changes in head and neck, cervix, stomach, prostate, brain, and colon more evident in hemizygous loss tumors. These data suggest that reduced expression of PTEN in tumors with hemizygous loss leads to tumor progression and influences anticancer immune response pathways.

Pan-Cancer Analysis Reveals Functional Similarity of Three lncRNAs across Multiple Tumors

Long non-coding RNAs (lncRNAs) are emerging as key regulators in many biological processes. The dysregulation of lncRNA expression has been associated with many diseases, including cancer. Mounting evidence suggests lncRNAs to be involved in cancer initiation, progression, and metastasis. Thus, understanding the functional implications of lncRNAs in tumorigenesis can aid in developing novel biomarkers and therapeutic targets. Rich cancer datasets, documenting genomic and transcriptomic alterations together with advancement in bioinformatics tools, have presented an opportunity to perform pan-cancer analyses across different cancer types. This study is aimed at conducting a pan-cancer analysis of lncRNAs by performing differential expression and functional analyses between tumor and non-neoplastic adjacent samples across eight cancer types. Among dysregulated lncRNAs, seven were shared across all cancer types. We focused on three lncRNAs, found to be consistently dysregulated among tumors. It has been observed that these three lncRNAs of interest are interacting with a wide range of genes across different tissues, yet enriching substantially similar biological processes, found to be implicated in cancer progression and proliferation.

Exploratory meta-analysis of hypoxic transcriptomes using a precise transcript reference sequence set

Gene expression studies are intrinsically biased, with many studies influenced by concomitant information such as gene-disease associations. This limitation can be overcome using a data-driven analysis approach without relying on ancillary information. The FANTOM CAGE-Associated Transcriptome project provides a comprehensive meta-assembly of the human transcriptome using coding and noncoding genes. Hypoxia strongly influences gene expression; in addition, noncoding RNA (ncRNA) metabolism is down-regulated in response to hypoxic stimuli. We evaluated the differential response of various transcripts to hypoxia by determining their hypoxia responsiveness scores. Enrichment analysis revealed that several genes associated with ncRNA metabolism, particularly those involved in ribosomal RNA processing, were down-regulated in response to hypoxia. Previously published information from the FANTOM CAGE-Associated Transcriptome project was suitable for meta-analysis of the transcriptome sequencing data from both coding and ncRNAs and to evaluate the hypoxia responsiveness of target transcripts and relationship between sense-antisense transcripts from the same locus. Our results may facilitate functional annotation of various transcripts including ncRNAs, allowing for both sense and antisense and coding and noncoding evaluations.

De novo annotation of lncRNA HOTAIR transcripts by long-read RNA capture-seq reveals a differentiation-driven isoform switch

Background

LncRNAs are tissue-specific and emerge as important regulators of various biological processes and as disease biomarkers. HOTAIR is a well-established pro-oncogenic lncRNA which has been attributed a variety of functions in cancer and native contexts. However, a lack of an exhaustive, cell type-specific annotation questions whether HOTAIR functions are supported by the expression of multiple isoforms.

Results

Using a capture long-read sequencing approach, we characterize HOTAIR isoforms expressed in human primary adipose stem cells. We find HOTAIR isoforms population displays varied splicing patterns, frequently leading to the exclusion or truncation of canonical LSD1 and PRC2 binding domains. We identify a highly cell type-specific HOTAIR isoform pool regulated by distinct promoter usage, and uncover a shift in the HOTAIR TSS usage that modulates the balance of HOTAIR isoforms at differentiation onset.

Conclusion

Our results highlight the complexity and cell type-specificity of HOTAIR isoforms and open perspectives on functional implications of these variants and their balance to key cellular processes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08887-w.

RNA-sequencing highlights differential regulated pathways involved in cell cycle and inflammation in orbitofacial neurofibromas

Although most commonly benign, neurofibromas (NFs) can have devastating functional and cosmetic effects in addition to the possibility of malignant transformation. Orbitofacial NFs, in particular, may cause progressive, disfiguring tumors of the lid, brow, temple, face, and orbit, and clinical evidence suggests that they may have increased local aggressiveness compared to NFs developing at other sites. The purpose of this study was to identify biological differences between orbitofacial NFs and those occurring at other anatomic sites. We performed RNA-sequencing in orbitofacial (n = 10) and non-orbitofacial (n = 9) NFs. Differential gene expression analysis demonstrated that a variety of gene sets including genes involved in cell proliferation, interferon, and immune-related pathways were enriched in orbitofacial NF. Comparisons with publicly available databases of various Schwann cell tumors and malignant peripheral nerve sheath tumor (MPNST) revealed a significant overlap of differentially expressed genes between orbitofacial versus non-orbitofacial NF and plexiform NF versus MPNST. In summary, we identified gene expression differences between orbitofacial NF and NFs occurring at other locations. Further investigation may be warranted, given that orbitofacial NF are notoriously difficult to treat and associated with disproportionate morbidity.

Molecular Insight Into the Therapeutic Potential of Long Non-coding RNA-Associated Competing Endogenous RNA Axes in Alzheimer's Disease: A Systematic Scoping Review

Alzheimer’s disease (AD) is a heterogeneous degenerative brain disorder with a rising prevalence worldwide. The two hallmarks that characterize the AD pathophysiology are amyloid plaques, generated via aggregated amyloid β, and neurofibrillary tangle, generated via accumulated phosphorylated tau. At the post-transcriptional and transcriptional levels, the regulatory functions of non-coding RNAs, in particular long non-coding RNAs (lncRNAs), have been ascertained in gene expressions. It is noteworthy that a number of lncRNAs feature a prevalent role in their potential of regulating gene expression through modulation of microRNAs via a process called the mechanism of competing endogenous RNA (ceRNA). Given the multifactorial nature of ceRNA interaction networks, they might be advantageous in complex disorders (e.g., AD) investigations at the therapeutic targets level. We carried out scoping review in this research to analyze validated loops of ceRNA in AD and focus on ceRNA axes associated with lncRNA. This scoping review was performed according to a six-stage methodology structure and PRISMA guideline. A systematic search of seven databases was conducted to find eligible articles prior to July 2021. Two reviewers independently performed publications screening and data extraction, and quantitative and qualitative analyses were conducted. Fourteen articles were identified that fulfill the inclusion criteria. Studies with different designs reported nine lncRNAs that were experimentally validated to act as ceRNA in AD in human-related studies, including BACE1-AS, SNHG1, RPPH1, NEAT1, LINC00094, SOX21-AS1, LINC00507, MAGI2-AS3, and LINC01311. The BACE1-AS/BACE1 was the most frequent ceRNA pair. Among miRNAs, miR-107 played a key role by regulating three different loops. Understanding the various aspects of this regulatory mechanism can help elucidate the unknown etiology of AD and provide new molecular targets for use in therapeutic and clinical applications.

recount3: summaries and queries for large-scale RNA-seq expression and splicing

We present recount3, a resource consisting of over 750,000 publicly available human and mouse RNA sequencing (RNA-seq) samples uniformly processed by our new Monorail analysis pipeline. To facilitate access to the data, we provide the recount3 and snapcount R/Bioconductor packages as well as complementary web resources. Using these tools, data can be downloaded as study-level summaries or queried for specific exon-exon junctions, genes, samples, or other features. Monorail can be used to process local and/or private data, allowing results to be directly compared to any study in recount3. Taken together, our tools help biologists maximize the utility of publicly available RNA-seq data, especially to improve their understanding of newly collected data. recount3 is available from http://rna.recount.bio .

Transcriptional landscape of PTEN loss in primary prostate cancer

BACKGROUND

PTEN is the most frequently lost tumor suppressor in primary prostate cancer (PCa) and its loss is associated with aggressive disease. However, the transcriptional changes associated with PTEN loss in PCa have not been described in detail. In this study, we highlight the transcriptional changes associated with PTEN loss in PCa.

METHODS

Using a meta-analysis approach, we leveraged two large PCa cohorts with experimentally validated PTEN and ERG status by Immunohistochemistry (IHC), to derive a transcriptomic signature of PTEN loss, while also accounting for potential confounders due to ERG rearrangements. This signature was expanded to lncRNAs using the TCGA quantifications from the FC-R2 expression atlas.

RESULTS

The signatures indicate a strong activation of both innate and adaptive immune systems upon PTEN loss, as well as an expected activation of cell-cycle genes. Moreover, we made use of our recently developed FC-R2 expression atlas to expand this signature to include many non-coding RNAs recently annotated by the FANTOM consortium. Highlighting potential novel lncRNAs associated with PTEN loss and PCa progression.

CONCLUSION

We created a PCa specific signature of the transcriptional landscape of PTEN loss that comprises both the coding and an extensive non-coding counterpart, highlighting potential new players in PCa progression. We also show that contrary to what is observed in other cancers, PTEN loss in PCa leads to increased activation of the immune system. These findings can help the development of new biomarkers and help guide therapy choices.

Epigenetic transcripts of LINC01311 and hsa-miR-146a-5p regulate neural development in a cellular model of Alzheimer's disease

Emerging evidence has shown that Long noncoding RNAs (LncRNAs) are aberrantly expressed and functionally involved in the development of neurodegenerative disorders. In this work, we investigated the regulatory effects of lncRNA of LINC01311 and its competing endogenous RNA target of hsa-miR-146a-5p in a cellular model of Alzheimer's disease (AD). SH-SY5Y cells were treated with synthetic Βeta-Amyloid Peptide (1-42) (AB1-42) in vitro to induce AD-like neural injuries. Expressions of LINC01311 and hsa-miR-146a-5p were monitored by qRT-PCR. LINC01311 was upregulated and hsa-miR-146a-5p downregulated to examine their functional regulations on AB1-42-induced apoptosis, proliferation slowdown, autophagy, and amyloid precursor protein (APP) accumulations. Hsa-miR-146a-5p was also overexpressed in LINC01311-upregulated SH-SY5Y cells to examine their correlated regulations on AB1-42-induced neural injuries. LINC01311 was downregulated whereas hsa-miR-146a-5p upregulated in AB1-42 treated SH-SY5Y cells. LINC01311 upregulation and hsa-miR-146a-5p downregulation protected AB1-42-induced apoptosis, proliferation slowdown, autophagy, and APP accumulations in SH-SY5Y cells. Hsa-miR-146a-5p overexpression reversed the protection of LINC01311 on AB1-42-induced neural injuries. Our work demonstrated that the epigenetic axis of LINC01311/hsa-miR-146a-5p was involved in the functional regulation of human-lineage neurons in a cellular model of AD, thus suggesting a clinical potential of exploring epigenetic network for treating AD patients.

LncRNAs and Available Databases

Long noncoding RNAs (lncRNAs) are involved in many regulatory mechanisms in practically every step of the RNA cycle, from transcription to RNA stability and translation. They are a highly heterogeneous class of molecules in terms of site of production, interaction networks, and functions. More and more databases are available on the web with the aim to make public information about lncRNA accessible to the scientific community. Here we review the most interesting resources with the purpose to organize a compendium of useful tools to interrogate before studying a lncRNA of interest.

Global Analysis of Transcription Start Sites in the New Ovine Reference Genome (Oar rambouillet v1.0)

The overall aim of the Ovine FAANG project is to provide a comprehensive annotation of the new highly contiguous sheep reference genome sequence (Oar rambouillet v1.0). Mapping of transcription start sites (TSS) is a key first step in understanding transcript regulation and diversity. Using 56 tissue samples collected from the reference ewe Benz2616, we have performed a global analysis of TSS and TSS-Enhancer clusters using Cap Analysis Gene Expression (CAGE) sequencing. CAGE measures RNA expression by 5' cap-trapping and has been specifically designed to allow the characterization of TSS within promoters to single-nucleotide resolution. We have adapted an analysis pipeline that uses TagDust2 for clean-up and trimming, Bowtie2 for mapping, CAGEfightR for clustering, and the Integrative Genomics Viewer (IGV) for visualization. Mapping of CAGE tags indicated that the expression levels of CAGE tag clusters varied across tissues. Expression profiles across tissues were validated using corresponding polyA+ mRNA-Seq data from the same samples. After removal of CAGE tags with <10 read counts, 39.3% of TSS overlapped with 5' ends of 31,113 transcripts that had been previously annotated by NCBI (out of a total of 56,308 from the NCBI annotation). For 25,195 of the transcripts, previously annotated by NCBI, no TSS meeting stringent criteria were identified. A further 14.7% of TSS mapped to within 50 bp of annotated promoter regions. Intersecting these predicted TSS regions with annotated promoter regions (±50 bp) revealed 46% of the predicted TSS were "novel" and previously un-annotated. Using whole-genome bisulfite sequencing data from the same tissues, we were able to determine that a proportion of these "novel" TSS were hypo-methylated (32.2%) indicating that they are likely to be reproducible rather than "noise". This global analysis of TSS in sheep will significantly enhance the annotation of gene models in the new ovine reference assembly. Our analyses provide one of the highest resolution annotations of transcript regulation and diversity in a livestock species to date.

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.