A Pan-cancer Transcriptome Analysis Reveals Pervasive Regulation through Alternative Promoters

Pan-cancer transcriptome analysis reveals widespread regulation through alternative tandem transcription initiation

Abnormal transcription initiation from alternative first exon has been reported to promote tumorigenesis. However, the prevalence and impact of gene expression regulation mediated by alternative tandem transcription initiation were mostly unknown in cancer. Here, we developed a robust computational method to analyze alternative tandem transcription start site (TSS) usage from standard RNA sequencing data. Applying this method to pan-cancer RNA sequencing datasets, we observed widespread dysregulation of tandem TSS usage in tumors, many of which were independent of changes in overall expression level or alternative first exon usage. We showed that the dynamics of tandem TSS usage was associated with epigenomic modulation. We found that significant 5' untranslated region shortening of gene TIMM13 contributed to increased protein production, and up-regulation of TIMM13 by CRISPR-mediated transcriptional activation promoted proliferation and migration of lung cancer cells. Our findings suggest that dysregulated tandem TSS usage represents an addtional layer of cancer-associated transcriptome alterations.

Integrative analysis of ultra-deep RNA-seq reveals alternative promoter usage as a mechanism of activating oncogenic programmes during prostate cancer progression

Transcription factor (TF) proteins regulate gene activity by binding to regulatory regions, most importantly at gene promoters. Many genes have alternative promoters (APs) bound by distinct TFs. The role of differential TF activity at APs during tumour development is poorly understood. Here we show, using deep RNA sequencing in 274 biopsies of benign prostate tissue, localized prostate tumours and metastatic castration-resistant prostate cancer, that AP usage increases as tumours progress and APs are responsible for a disproportionate amount of tumour transcriptional activity. Expression of the androgen receptor (AR), the key driver of prostate tumour activity, is correlated with elevated AP usage. We identified AR, FOXA1 and MYC as potential drivers of AP activation. DNA methylation is a likely mechanism for AP activation during tumour progression and lineage plasticity. Our data suggest that prostate tumours activate APs to magnify the transcriptional impact of tumour drivers, including AR and MYC.

TAD hierarchy restricts poised LTR activation and loss of TAD hierarchy promotes LTR co-option in cancer

Transposable elements (TEs) are abundant in the human genome, and they provide the sources for genetic and functional diversity. The regulation of TEs expression and their functional consequences in physiological conditions and cancer development remain to be fully elucidated. Previous studies suggested TEs are repressed by DNA methylation and chromatin modifications. The effect of 3D chromatin topology on TE regulation remains elusive. Here, by integrating transcriptome and 3D genome architecture studies, we showed that haploinsufficient loss of NIPBL selectively activates alternative promoters at the long terminal repeats (LTRs) of the TE subclasses. This activation occurs through the reorganization of topologically associating domain (TAD) hierarchical structures and recruitment of proximal enhancers. These observations indicate that TAD hierarchy restricts transcriptional activation of LTRs that already possess open chromatin features. In cancer, perturbation of the hierarchical chromatin topology can lead to co-option of LTRs as functional alternative promoters in a context-dependent manner and drive aberrant transcriptional activation of novel oncogenes and other divergent transcripts. These data uncovered a new layer of regulatory mechanism of TE expression beyond DNA and chromatin modification in human genome. They also posit the TAD hierarchy dysregulation as a novel mechanism for alternative promoter-mediated oncogene activation and transcriptional diversity in cancer, which may be exploited therapeutically.

Machine Learning-Based Integrated Multiomics Characterization of Colorectal Cancer Reveals Distinctive Metabolic Signatures

The metabolic signature identification of colorectal cancer is critical for its early diagnosis and therapeutic approaches that will significantly block cancer progression and improve patient survival. Here, we combined an untargeted metabolic analysis strategy based on internal extractive electrospray ionization mass spectrometry and the machine learning approach to analyze metabolites in 173 pairs of cancer samples and matched normal tissue samples to build robust metabolic signature models for diagnostic purposes. Screening and independent validation of metabolic signatures from colorectal cancers via machine learning methods (Logistic Regression_L1 for feature selection and eXtreme Gradient Boosting for classification) was performed to generate a panel of seven signatures with good diagnostic performance (the accuracy of 87.74%, sensitivity of 85.82%, and specificity of 89.66%). Moreover, seven signatures were evaluated according to their ability to distinguish between cancer and normal tissues, with the metabolic molecule PC (30:0) showing good diagnostic performance. In addition, genes associated with PC (30:0) were identified by multiomics analysis (combining metabolic data with transcriptomic data analysis) and our results showed that PC (30:0) could promote the proliferation of colorectal cancer cell SW480, revealing the correlation between genetic changes and metabolic dysregulation in cancer. Overall, our results reveal potential determinants affecting metabolite dysregulation, paving the way for a mechanistic understanding of altered tissue metabolites in colorectal cancer and design interventions for manipulating the levels of circulating metabolites.

A Bayesian Change Point Model for Dynamic Alternative Transcription Start Site Usage During Cellular Differentiation

ABSTRACT An alternative transcription start site (ATSS) is a major driving force for increasing the complexity of transcripts in human tissues. As a transcriptional regulatory mechanism, ATSS has biological significance. Many studies have confirmed that ATSS plays an important role in diseases and cell development and differentiation. However, exploration of its dynamic mechanisms remains insufficient. Identifying ATSS change points during cell differentiation is critical for elucidating potential dynamic mechanisms. For relative ATSS usage as percentage data, the existing methods lack sensitivity to detect the change point for ATSS longitudinal data. In addition, some methods have strict requirements for data distribution and cannot be applied to deal with this problem. In this study, the Bayesian change point detection model was first constructed using reparameterization techniques for two parameters of a beta distribution for the percentage data type, and the posterior distributions of parameters and change points were obtained using Markov Chain Monte Carlo (MCMC) sampling. With comprehensive simulation studies, the performance of the Bayesian change point detection model is found to be consistently powerful and robust across most scenarios with different sample sizes and beta distributions. Second, differential ATSS events in the real data, whose change points were identified using our method, were clustered according to their change points. Last, for each change point, pathway and transcription factor motif analyses were performed on its differential ATSS events. The results of our analyses demonstrated the effectiveness of the Bayesian change point detection model and provided biological insights into cell differentiation.

Atherosclerosis-related biomarker PABPC1 predicts pan-cancer events

BACKGROUND

Atherosclerosis (AS) and tumours are the leading causes of death worldwide and share common risk factors, detection methods and molecular markers. Therefore, searching for serum markers shared by AS and tumours is beneficial to the early diagnosis of patients.

METHODS

The sera of 23 patients with AS-related transient ischaemic attack were screened by serological identification of antigens through recombinant cDNA expression cloning (SEREX), and cDNA clones were identified. Pathway function enrichment analysis was performed on cDNA clones to identify their biological pathways and determine whether they were related to AS or tumours. Subsequently, gene-gene and protein-protein interactions were performed and AS-associated markers would be discovered. The expression of AS biomarkers in human normal organs and pan-cancer tumour tissues were explored. Then, immune infiltration level and tumour mutation burden of various immune cells were evaluated. Survival curves analysis could show the expression of AS markers in pan-cancer.

RESULTS

AS-related sera were screened by SEREX, and 83 cDNA clones with high homology were obtained. Through functional enrichment analysis, it was found that their functions were closely related to AS and tumour functions. After multiple biological information interaction screening and the external cohort validating, poly(A) binding protein cytoplasmic 1 (PABPC1) was found to be a potential AS biomarker. To assess whether PABPC1 was related to pan-cancer, its expression in different tumour pathological stages and ages was screened. Since AS-associated proteins were closely related to cancer immune infiltration, we investigated and found that PABPC1 had the same role in pan-cancer. Finally, analysis of Kaplan-Meier survival curves revealed that high PABPC1 expression in pan-cancer was associated with high risk of death.

CONCLUSIONS

Through the findings of SEREX and bioinformatics pan-cancer analysis, we concluded that PABPC1 might serve as a potential biomarker for the prediction and diagnosis of AS and pan-cancer.

Genetic variation is a key determinant of chromatin accessibility and drives differences in the regulatory landscape of C57BL/6J and 129S1/SvImJ mice

Most common genetic variants associated with disease are located in non-coding regions of the genome. One mechanism by which they function is through altering transcription factor (TF) binding. In this study, we explore how genetic variation is connected to differences in the regulatory landscape of livers from C57BL/6J and 129S1/SvImJ mice fed either chow or a high-fat diet. To identify sites where regulatory variation affects TF binding and nearby gene expression, we employed an integrative analysis of H3K27ac ChIP-seq (active enhancers), ATAC-seq (chromatin accessibility) and RNA-seq (gene expression). We show that, across all these assays, the genetically driven (i.e. strain-specific) differences in the regulatory landscape are more pronounced than those modified by diet. Most notably, our analysis revealed that differentially accessible regions (DARs, N = 29635, FDR < 0.01 and fold change > 50%) are almost always strain-specific and enriched with genetic variation. Moreover, proximal DARs are highly correlated with differentially expressed genes. We also show that TF binding is affected by genetic variation, which we validate experimentally using ChIP-seq for TCF7L2 and CTCF. This study provides detailed insights into how non-coding genetic variation alters the gene regulatory landscape, and demonstrates how this can be used to study the regulatory variation influencing TF binding.

Identification of an alternative short ARID5B isoform associated with B-ALL survival

Utilizing RNA sequence (RNA-Seq) splice junction data from a cohort of 1841 B-cell acute lymphoblastic leukemia (B-ALL) patients we define transcriptionally distinct isoforms of ARID5B, a risk-associated gene identified in genome wide association studies (GWAS), which associate with disease survival. Short (S) and long (L) ARID5B transcripts, which differ in an encoded BAH-like chromatin interaction domain, show remarkable correlation to the isoform splicing pattern. Testing of the ARID5B proximal promoter of the S & L isoforms indicated that both are functionally independent in luciferase reporter assays. Increased short isoform expression is associated with decreased event-free and overall survival. The abundance of short and long transcripts strongly correlates to B-ALL prognostic stratification, where B-ALL subtypes with poor outcomes express a higher proportion of the S-isoform. These data demonstrate that the analysis of independent promoters and alternative splicing events are essential for improved risk stratification and a more complete understanding of disease pathology.

hsa_circ_0087100/hsa-miR-6743-5p affects Th1 cell differentiation by regulating STAT1 in diabetic retinopathy

Objective: To elucidate the role of the competitive endogenous RNA (ceRNA) network in immune infiltration of diabetic retinopathy (DR). Methods: We obtained differentially expressed (DE) circRNAs, miRNAs and mRNAs from the Gene Expression Omnibus database. Then, we identified immune infiltration by CIBERSORT and single-sample gene set enrichment analysis and discovered co-expression genes by weighted gene co-expression network analysis. Furthermore, STAT1-mediated Th1 differentiation was determined in DR cell models, DR patients and DR mouse models. Results: hsa_circ_0087100/hsa-miR-6743-5p/STAT1 was involved in immune infiltration of Th1 cells. Aberrant expression of the ceRNA network and STAT1-mediated Th1 differentiation was thus verified in vitro and in vivo. Conclusion: hsa_circ_0087100/hsa-miR-6743-5p/STAT1 may affect Th1 cell differentiation in DR.

(Alternative) transcription start sites as regulators of RNA processing

Alternative transcription start site usage (ATSS) is a widespread regulatory strategy that enables genes to choose between multiple genomic loci for initiating transcription. This mechanism is tightly controlled during development and is often altered in disease states. In this review, we examine the growing evidence highlighting a role for transcription start sites (TSSs) in the regulation of mRNA isoform selection during and after transcription. We discuss how the choice of transcription initiation sites influences RNA processing and the importance of this crosstalk for cell identity and organism function. We also speculate on possible mechanisms underlying the integration of transcriptional and post-transcriptional processes.

Quantitative profiling of human translation initiation reveals regulatory elements that potently affect endogenous and therapeutically modified mRNAs

mRNA therapeutics offer a potentially universal strategy for the efficient development and delivery of therapeutic proteins. Current mRNA vaccines include chemically modified nucleotides to reduce cellular immunogenicity. Here, we develop an efficient, high-throughput method to measure human translation initiation on therapeutically modified as well as endogenous RNAs. Using systems-level biochemistry, we quantify ribosome recruitment to tens of thousands of human 5' untranslated regions and identify sequences that mediate 250-fold effects. We observe widespread effects of coding sequences on translation initiation and identify small regulatory elements of 3-6 nucleotides that are sufficient to potently affect translational output. Incorporation of N1-methylpseudouridine (m1Ψ) selectively enhances translation by specific 5' UTRs that we demonstrate surpass those of current mRNA vaccines. Our approach is broadly applicable to dissect mechanisms of human translation initiation and engineer more potent therapeutic mRNAs.

Highlights

Measurement of >30,000 human 5' UTRs reveals a 250-fold range of translation outputSystematic mutagenesis demonstrates the causality of short (3-6nt) regulatory elementsN1-methylpseudouridine alters translation initiation in a sequence-specific mannerOptimal modified 5' UTRs outperform those in the current class of mRNA vaccines.

A zebrafish NLRX1 isoform downregulates fish IFN responses by targeting the adaptor STING

In mammals, NLRX1 is a unique member of the nucleotide-binding domain and leucine-rich repeat (NLR) family showing an ability to negatively regulate IFN antiviral immunity. Intron-containing genes, including NLRX1, have more than one transcript due to alternative splicing; however, little is known about the function of its splicing variants. Here, we identified a transcript variant of NLRX1 in zebrafish (Danio rerio), termed NLRX1-tv4, as a negative regulator of fish IFN response. Zebrafish NLRX1-tv4 was slightly induced by viral infection, with an expression pattern similar to the full-length NLRX1. Despite the lack of an N-terminal domain that exists in the full-length NLRX1, overexpression of NLRX1-tv4 still impaired fish IFN antiviral response and promoted viral replication in fish cells, similar to the full-length NLRX1. Mechanistically, NLRX1-tv4 targeted STING for proteasome-dependent protein degradation by recruiting an E3 ubiquitin ligase RNF5 to drive the K48-linked ubiquitination, eventually downregulating the IFN antiviral response. Mapping of NLRX1-tv4 domains showed that its N-terminal and C-terminal regions exhibited a similar potential to inhibit STING-mediated IFN antiviral response. Our findings reveal that like the full-length NLRX1, zebrafish NLRX-tv4 functions as an inhibitor to shape fish IFN antiviral response.IMPORTANCEIn this study, we demonstrate that a transcript variant of zebrafish NLRX1, termed NLRX1-tv4, downregulates fish IFN response and promotes virus replication by targeting STING for protein degradation and impairing the interaction of STING and TBK1 and that its N- and C-terminus exhibit a similar inhibitory potential. Our results are helpful in clarifying the current contradictory understanding of structure and function of vertebrate NLRX1s.

A type 1 immunity-restricted promoter of the IL-33 receptor gene directs antiviral T-cell responses

The pleiotropic alarmin interleukin-33 (IL-33) drives type 1, type 2 and regulatory T-cell responses via its receptor ST2. Subset-specific differences in ST2 expression intensity and dynamics suggest that transcriptional regulation is key in orchestrating the context-dependent activity of IL-33-ST2 signaling in T-cell immunity. Here, we identify a previously unrecognized alternative promoter in mice and humans that is located far upstream of the curated ST2-coding gene and drives ST2 expression in type 1 immunity. Mice lacking this promoter exhibit a selective loss of ST2 expression in type 1- but not type 2-biased T cells, resulting in impaired expansion of cytotoxic T cells (CTLs) and T-helper 1 cells upon viral infection. T-cell-intrinsic IL-33 signaling via type 1 promoter-driven ST2 is critical to generate a clonally diverse population of antiviral short-lived effector CTLs. Thus, lineage-specific alternative promoter usage directs alarmin responsiveness in T-cell subsets and offers opportunities for immune cell-specific targeting of the IL-33-ST2 axis in infections and inflammatory diseases.

FLIBase: a comprehensive repository of full-length isoforms across human cancers and tissues

Regulatory processes at the RNA transcript level play a crucial role in generating transcriptome diversity and proteome composition in human cells, impacting both physiological and pathological states. This study introduces FLIBase (www.FLIBase.org), a specialized database that focuses on annotating full-length isoforms using long-read sequencing techniques. We collected and integrated long-read (351 samples) and short-read (12 469 samples) RNA sequencing data from diverse normal and cancerous human tissues and cells. The current version of FLIBase comprises a total of 983 789 full-length spliced isoforms, identified through long-read sequences and verified using short-read exon-exon splice junctions. Of these, 188 248 isoforms have been annotated, while 795 541 isoforms remain unannotated. By overcoming the limitations of short-read RNA sequencing methods, FLIBase provides an accurate and comprehensive representation of full-length transcripts. These comprehensive annotations empower researchers to undertake various downstream analyses and investigations. Importantly, FLIBase exhibits a significant advantage in identifying a substantial number of previously unannotated isoforms and tumor-specific RNA transcripts. These tumor-specific RNA transcripts have the potential to serve as a source of immunogenic recurrent neoantigens. This remarkable discovery holds tremendous promise for advancing the development of tailored RNA-based diagnostic and therapeutic strategies for various types of human cancer.

Quantitative analysis of transcription start site selection reveals control by DNA sequence, RNA polymerase II activity and NTP levels

Transcription start site (TSS) selection is a key step in gene expression and occurs at many promoter positions over a wide range of efficiencies. Here we develop a massively parallel reporter assay to quantitatively dissect contributions of promoter sequence, nucleoside triphosphate substrate levels and RNA polymerase II (Pol II) activity to TSS selection by 'promoter scanning' in Saccharomyces cerevisiae (Pol II MAssively Systematic Transcript End Readout, 'Pol II MASTER'). Using Pol II MASTER, we measure the efficiency of Pol II initiation at 1,000,000 individual TSS sequences in a defined promoter context. Pol II MASTER confirms proposed critical qualities of S. cerevisiae TSS -8, -1 and +1 positions, quantitatively, in a controlled promoter context. Pol II MASTER extends quantitative analysis to surrounding sequences and determines that they tune initiation over a wide range of efficiencies. These results enabled the development of a predictive model for initiation efficiency based on sequence. We show that genetic perturbation of Pol II catalytic activity alters initiation efficiency mostly independently of TSS sequence, but selectively modulates preference for the initiating nucleotide. Intriguingly, we find that Pol II initiation efficiency is directly sensitive to guanosine-5'-triphosphate levels at the first five transcript positions and to cytosine-5'-triphosphate and uridine-5'-triphosphate levels at the second position genome wide. These results suggest individual nucleoside triphosphate levels can have transcript-specific effects on initiation, representing a cryptic layer of potential regulation at the level of Pol II biochemical properties. The results establish Pol II MASTER as a method for quantitative dissection of transcription initiation in eukaryotes.

Pan-cancer analysis of promoter activity quantitative trait loci

Altered promoter activity has been generally observed in diverse biological processes, including tumorigenesis. Accumulating evidence suggests that employing a quantitative trait locus mapping approach is effective in comprehending the genetic basis of promoter activity. By utilizing genotype data from The Cancer Genome Atlas and calculating corresponding promoter activity values using proActiv, we systematically evaluated the impact of genetic variants on promoter activity and identified >1.0 million promoter activity quantitative trait loci (paQTLs) as both cis- and trans-acting. Additionally, leveraging data from the genome-wide association study (GWAS) catalog, we discovered >1.3 million paQTLs that overlap with known GWAS linkage disequilibrium regions. Remarkably, ∼9324 paQTLs exhibited significant associations with patient prognosis. Moreover, investigating the impact of promoter activity on >1000 imputed antitumor therapy responses among pan-cancer patients revealed >43 000 million significant associations. Furthermore, ∼25 000 significant associations were identified between promoter activity and immune cell abundance. Finally, a user-friendly data portal, Pancan-paQTL (https://www.hbpding.com/PancanPaQTL/), was constructed for users to browse, search and download data of interest. Pancan-paQTL serves as a comprehensive multidimensional database, enabling functional and clinical investigations into genetic variants associated with promoter activity, drug responses and immune infiltration across multiple cancer types.

Regulation of alternative splicing and polyadenylation in neurons

Cell-type-specific gene expression is a fundamental feature of multicellular organisms and is achieved by combinations of regulatory strategies. Although cell-restricted transcription is perhaps the most widely studied mechanism, co-transcriptional and post-transcriptional processes are also central to the spatiotemporal control of gene functions. One general category of expression control involves the generation of multiple transcript isoforms from an individual gene, whose balance and cell specificity are frequently tightly regulated via diverse strategies. The nervous system makes particularly extensive use of cell-specific isoforms, specializing the neural function of genes that are expressed more broadly. Here, we review regulatory strategies and RNA-binding proteins that direct neural-specific isoform processing. These include various classes of alternative splicing and alternative polyadenylation events, both of which broadly diversify the neural transcriptome. Importantly, global alterations of splicing and alternative polyadenylation are characteristic of many neural pathologies, and recent genetic studies demonstrate how misregulation of individual neural isoforms can directly cause mutant phenotypes.

Transcript CD81-215 may be a long noncoding RNA of stromal origin with tumor-promoting role in colon cancer

The role of tetraspanin CD81 in malignant transformation is best studied in colorectal cancer, and it appears that other transcripts beside the fully coding mRNA may also be dysregulated in malignant cells. Recent data from a comprehensive pan-cancer transcriptome analysis demonstrated differential activity of two alternative CD81 gene promoters in malignant versus nonmalignant gut mucosa. The promoter active in gut mucosa gives rise to transcripts CD81-203 and CD81-213, while the promoter active in colon and rectal cancer gives rise to transcripts CD81-205 and CD81-215. Our study aimed to explore the biomarker potential of the transcripts from the alternative CD81 gene promoters in colon cancer, as well as to investigate their structure and potential function using in silico tools. The analysis of the transcripts' expression in several colon cell lines cultivated in 2D and 3D and a set of colon cancer and healthy gut mucosa samples by qPCR and RNA sequencing suggested their low expression and stromal origin. Expression patterns in tumor and nontumor tissue along with in silico data suppose that the transcript CD81-215 may be a noncoding RNA of stromal origin with possible involvement in signaling related to malignant transformation.

Intra-promoter switch of transcription initiation sites in proliferation signaling-dependent RNA metabolism

Global changes in transcriptional regulation and RNA metabolism are crucial features of cancer development. However, little is known about the role of the core promoter in defining transcript identity and post-transcriptional fates, a potentially crucial layer of transcriptional regulation in cancer. In this study, we use CAGE-seq analysis to uncover widespread use of dual-initiation promoters in which non-canonical, first-base-cytosine (C) transcription initiation occurs alongside first-base-purine initiation across 59 human cancers and healthy tissues. C-initiation is often followed by a 5' terminal oligopyrimidine (5'TOP) sequence, dramatically increasing the range of genes potentially subjected to 5'TOP-associated post-transcriptional regulation. We show selective, dynamic switching between purine and C-initiation site usage, indicating transcription initiation-level regulation in cancers. We additionally detail global metabolic changes in C-initiation transcripts that mark differentiation status, proliferative capacity, radiosensitivity, and response to irradiation and to PI3K-Akt-mTOR and DNA damage pathway-targeted radiosensitization therapies in colorectal cancer organoids and cancer cell lines and tissues.

Co-expression of PADI isoforms during progenitor differentiation enables functional diversity

Protein isoforms, generated through alternative splicing or promoter usage, contribute to tissue function. Here, we characterize the expression of predicted Padi3α and Padi3β isoforms in hair follicles and describe expression of Padi2β, a hitherto unknown PADI2 isoform, in the oligodendrocyte lineage. Padi2β transcription is initiated from a downstream intronic promoter, generating an N-terminally truncated, unstable, PADI2β. By contrast to the established role of the canonical PADI2 (PADI2α) (Falcao et al. 2019 Cell Rep. 27, 1090-1102.e10. (doi:10.1016/j.celrep.2019.03.108)), PADI2β inhibits oligodendrocyte differentiation, suggesting that PADI2 isoforms exert opposing effects on oligodendrocyte lineage progression. We localize Padi3α and Padi3β to developing hair follicles and find that both transcripts are expressed at low levels in progenitor cells, only to increase in expression concomitant with differentiation. When expressed in vitro, PADI3α and PADI3β are enriched in the cytoplasm and precipitate together. Whereas PADI3β protein stability is low and PADI3β fails to induce protein citrullination, we find that the enzymatic activity and protein stability of PADI3α is reduced in the presence of PADI3β. We propose that PADI3β modulates PADI3α activity by direct binding and heterodimer formation. Here, we establish expression and function of Padi2 and Padi3 isoforms, expanding on the mechanisms in place to regulate citrullination in complex tissues. This article is part of the Theo Murphy meeting issue 'The virtues and vices of protein citrullination'.

Transposable element-initiated enhancer-like elements generate the subgenome-biased spike specificity of polyploid wheat

Transposable elements (TEs) comprise ~85% of the common wheat genome, which are highly diverse among subgenomes, possibly contribute to polyploid plasticity, but the causality is only assumed. Here, by integrating data from gene expression cap analysis and epigenome profiling via hidden Markov model in common wheat, we detect a large proportion of enhancer-like elements (ELEs) derived from TEs producing nascent noncoding transcripts, namely ELE-RNAs, which are well indicative of the regulatory activity of ELEs. Quantifying ELE-RNA transcriptome across typical developmental stages reveals that TE-initiated ELE-RNAs are mainly from RLG_famc7.3 specifically expanded in subgenome A. Acquisition of spike-specific transcription factor binding likely confers spike-specific expression of RLG_famc7.3-initiated ELE-RNAs. Knockdown of RLG_famc7.3-initiated ELE-RNAs resulted in global downregulation of spike-specific genes and abnormal spike development. These findings link TE expansion to regulatory specificity and polyploid developmental plasticity, highlighting the functional impact of TE-driven regulatory innovation on polyploid evolution.

Epigenetic Landscape and Therapeutic Implication of Gene Isoforms of Doublecortin-Like Kinase 1 for Cancer Stem Cells

While significant strides have been made in understanding cancer biology, the enhancement in patient survival is limited, underscoring the urgency for innovative strategies. Epigenetic modifications characterized by hereditary shifts in gene expression without changes to the DNA sequence play a critical role in producing alternative gene isoforms. When these processes go awry, they influence cancer onset, growth, spread, and cancer stemness. In this review, we delve into the epigenetic and isoform nuances of the protein kinase, doublecortin-like kinase 1 (DCLK1). Recognized as a hallmark of tumor stemness, DCLK1 plays a pivotal role in tumorigenesis, and DCLK1 isoforms, shaped by alternative promoter usage and splicing, can reveal potential therapeutic touchpoints. Our discussion centers on recent findings pertaining to the specific functions of DCLK1 isoforms and the prevailing understanding of its epigenetic regulation via its two distinct promoters. It is noteworthy that all DCLK1 isoforms retain their kinase domain, suggesting that their unique functionalities arise from non-kinase mechanisms. Consequently, our research has pivoted to drugs that specifically influence the epigenetic generation of these DCLK1 isoforms. We posit that a combined therapeutic approach, harnessing both the epigenetic regulators of specific DCLK1 isoforms and DCLK1-targeted drugs, may prove more effective than therapies that solely target DCLK1.

CamoTSS: analysis of alternative transcription start sites for cellular phenotypes and regulatory patterns from 5' scRNA-seq data

Five-prime single-cell RNA-seq (scRNA-seq) has been widely employed to profile cellular transcriptomes, however, its power of analysing transcription start sites (TSS) has not been fully utilised. Here, we present a computational method suite, CamoTSS, to precisely identify TSS and quantify its expression by leveraging the cDNA on read 1, which enables effective detection of alternative TSS usage. With various experimental data sets, we have demonstrated that CamoTSS can accurately identify TSS and the detected alternative TSS usages showed strong specificity in different biological processes, including cell types across human organs, the development of human thymus, and cancer conditions. As evidenced in nasopharyngeal cancer, alternative TSS usage can also reveal regulatory patterns including systematic TSS dysregulations.

Alternative Transcripts Diversify Genome Function for Phenome Relevance to Health and Diseases

Manipulation using alternative exon splicing (AES), alternative transcription start (ATS), and alternative polyadenylation (APA) sites are key to transcript diversity underlying health and disease. All three are pervasive in organisms, present in at least 50% of human protein-coding genes. In fact, ATS and APA site use has the highest impact on protein identity, with their ability to alter which first and last exons are utilized as well as impacting stability and translation efficiency. These RNA variants have been shown to be highly specific, both in tissue type and stage, with demonstrated importance to cell proliferation, differentiation and the transition from fetal to adult cells. While alternative exon splicing has a limited effect on protein identity, its ubiquity highlights the importance of these minor alterations, which can alter other features such as localization. The three processes are also highly interwoven, with overlapping, complementary, and competing factors, RNA polymerase II and its CTD (C-terminal domain) chief among them. Their role in development means dysregulation leads to a wide variety of disorders and cancers, with some forms of disease disproportionately affected by specific mechanisms (AES, ATS, or APA). Challenges associated with the genome-wide profiling of RNA variants and their potential solutions are also discussed in this review.

Cancer-Immunity Marker RNA Expression Levels across Gynecologic Cancers: Implications for Immunotherapy

Our objective was to characterize cancer-immunity marker expression in gynecologic cancers and compare immune landscapes between gynecologic tumor subtypes and with nongynecologic solid tumors. RNA expression levels of 51 cancer-immunity markers were analyzed in patients with gynecologic cancers versus nongynecologic cancers, and normalized to a reference population of 735 control cancers, ranked from 0 to 100, and categorized as low (0-24), moderate (25-74), or high (75-100) percentile rank. Of the 72 patients studied, 43 (60%) had ovarian, 24 (33%) uterine, and 5 (7%) cervical cancer. No two immune profiles were identical according to expression rank (0-100) or rank level (low, moderate, or high). Patients with cervical cancer had significantly higher expression level ranks of immune activating, proinflammatory, tumor-infiltrating lymphocyte markers, and checkpoints than patients with uterine or ovarian cancer (P < 0.001 for all comparisons). However, there were no significant differences in immune marker expression between uterine and ovarian cancers. Tumors with PD-L1 tumor proportional score (TPS) ≥1% versus 0% had significantly higher expression levels of proinflammatory markers (58 vs. 49%, P = 0.0004). Compared to patients with nongynecologic cancers, more patients with gynecologic cancers express high levels of IDO-1 (44 vs. 13%, P < 0.001), LAG3 (35 vs. 21%, P = 0.008), and IL10 (31 vs. 15%, P = 0.002.) Patients with gynecologic cancers have complex and heterogeneous immune landscapes that are distinct from patient to patient and from other solid tumors. High levels of IDO1 and LAG3 suggest that clinical trials with IDO1 inhibitors or LAG3 inhibitors, respectively, may be warranted in gynecologic cancers.

Novel pyroptosis-related lncRNAs and ceRNAs predict osteosarcoma prognosis and indicate immune microenvironment signatures

Objective

To study pyroptosis-related biomarkers that are associated with the prognosis and immune microenvironment characteristics of osteosarcoma (OS). The goal is to establish a foundation for the prognosis and treatment of OS.

Methods

We retrieved transcriptome and clinical data from The Cancer Genome Atlas (TCGA) database for 88 OS patients. Using this data, we constructed a prognostic model to identify pyroptosis-related genes (PRGs) associated with OS prognosis. To further explore the biological function of these PRGs, we performed enrichment analysis. To identify pyroptosis-related long non-coding RNAs (PRLncs) associated with the prognosis of OS, we performed co-expression analysis. Subsequently, a risk prognostic model was constructed using these PRLncs to generate a risk score, termed as PRLncs-score, thereby obtaining PRLncs associated with the prognosis of OS. The accuracy of the prognostic model was verified through survival analysis, risk curve, independent prognostic analysis, receiver operating characteristic (ROC) curve, difference analysis between high- and low-risk groups, and clinical correlation analysis. And to determine whether PRLncs-score is independent prognostic factor for OS. In addition, we further conducted external and internal validation for the risk prognosis model. Further analyses of immune cell infiltration and tumor microenvironment were performed. A pyroptosis-related competitive endogenous RNA (PRceRNA) network was constructed to obtain PRceRNAs associated with the prognosis of OS and performed gene set enrichment analysis (GSEA) on PRceRNA genes.

Results

We obtained five PRGs (CHMP4C, BAK1, GSDMA, CASP1, and CASP6) that predicted OS prognosis and seven PRLncs (AC090559.1, AP003119.2, CARD8-AS1, AL390728.4, SATB2-AS1, AL133215.2, and AC009495.3) and one PRceRNA (CARD8-AS1-hsa-miR-21-5p-IL1B) that predicted OS prognosis and indicated characteristics of the OS immune microenvironment. The PRLncs-score, in combination with other clinical features, was established as an independent prognostic factor for OS patients. Subsequent scrutiny of the tumor microenvironment and immune infiltration indicated that patients with low-PRLncs-scores were associated with reduced metastatic risk, improved survival rates, heightened levels of immune cells and stroma, and increased immune activity compared to those with high-PRLncs-scores.

Conclusion

The study's findings offer insight into the prognosis of OS and its immune microenvironment, and hold promise for improving early diagnosis and immunotherapy.

Causal role of a promoter polymorphism in natural variation of the Arabidopsis floral repressor gene FLC

Noncoding polymorphism frequently associates with phenotypic variation, but causation and mechanism are rarely established. Noncoding single-nucleotide polymorphisms (SNPs) characterize the major haplotypes of the Arabidopsis thaliana floral repressor gene FLOWERING LOCUS C (FLC). This noncoding polymorphism generates a range of FLC expression levels, determining the requirement for and the response to winter cold. The major adaptive determinant of these FLC haplotypes was shown to be the autumnal levels of FLC expression. Here, we investigate how noncoding SNPs influence FLC transcriptional output. We identify an upstream transcription start site (uTSS) cluster at FLC, whose usage is increased by an A variant at the promoter SNP-230. This variant is present in relatively few Arabidopsis accessions, with the majority containing G at this site. We demonstrate a causal role for the A variant at -230 in reduced FLC transcriptional output. The G variant upregulates FLC expression redundantly with the major transcriptional activator FRIGIDA (FRI). We demonstrate an additive interaction of SNP-230 with an intronic SNP+259, which also differentially influences uTSS usage. Combinatorial interactions between noncoding SNPs and transcriptional activators thus generate quantitative variation in FLC transcription that has facilitated the adaptation of Arabidopsis accessions to distinct climates.

CanISO: a database of genomic and transcriptomic variations in domestic dog (Canis lupus familiaris)

BACKGROUND

The domestic dog, Canis lupus familiaris, is a companion animal for humans as well as an animal model in cancer research due to similar spontaneous occurrence of cancers as humans. Despite the social and biological importance of dogs, the catalogue of genomic variations and transcripts for dogs is relatively incomplete.

RESULTS

We developed CanISO, a new database to hold a large collection of transcriptome profiles and genomic variations for domestic dogs. CanISO provides 87,692 novel transcript isoforms and 60,992 known isoforms from whole transcriptome sequencing of canine tumors (N = 157) and their matched normal tissues (N = 64). CanISO also provides genomic variation information for 210,444 unique germline single nucleotide polymorphisms (SNPs) from the whole exome sequencing of 183 dogs, with a query system that searches gene- and transcript-level information as well as covered SNPs. Transcriptome profiles can be compared with corresponding human transcript isoforms at a tissue level, or between sample groups to identify tumor-specific gene expression and alternative splicing patterns.

CONCLUSIONS

CanISO is expected to increase understanding of the dog genome and transcriptome, as well as its functional associations with humans, such as shared/distinct mechanisms of cancer. CanISO is publicly available at https://www.kobic.re.kr/caniso/ .

SPAG5 Expression Predicts Poor Prognosis and is Associated With Adverse Immune Infiltration in Lung Adenocarcinomas

Background

Sperm-associated antigen 5 (SPAG5) has been identified as a novel driver oncogene involved in multiple cancers; however, its role in lung adenocarcinoma (LUAD) needs further investigation. Our study aims to elucidate the potential significance of SPAG5 in LUAD prognosis and its implications for the efficacy of immunotherapy.

Methods

In this study, we used bioinformatics analysis and tissue microarray (TMA) staining to examine the potential role of SPAG5 in LUAD survival and response to immunotherapy. We used the Oncomine, TIMER2.0, Gene Expression Profiling Interactive Analysis (GEPIA), Sangerbox, PredicScan, and Kaplan-Meier Plotter databases to examine the expression and prognostic role of SPAG5 in the LUAD of The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO), and other databases. We also used Cancer Single-cell State Atlas (CancerSEA) and Tumor Immune Estimation Resource (TIMER2.0) to analyze the association of SPAG5 with malignant phenotype and tumor immune microenvironment. Furthermore, Immune Cell Abundance Identifier (ImmuCellAI) analysis of TCGA sequencing data was used to predict the role of SPAG5 in determining the response to immune checkpoint blockade (ICB) treatment in LUAD. Co-expression analysis of programmed death-ligand 1 (PD-L1) and SPAG5 was performed using LUAD TMA immunohistochemistry (IHC) analysis.

Results

Our findings indicate that SPAG5 is overexpressed in LUAD and is positively correlated with advanced clinical stage, poor overall survival, relapse-free survival, and progression-free survival outcomes. SPAG5 may be involved in regulating the cell cycle, proliferation, invasion, DNA damage and repair, and tumor immunosuppression. Furthermore, TMA IHC analysis showed a positive correlation between PD-L1 expression in LUAD and SPAG5 which suggests that SPAG5 may serve as a potential predictor of response to ICB therapy in LUAD.

Conclusions

Our results highlight the role of SAPG5 in promoting a tumor malignancy phenotype and immunosuppression in LUAD and suggest that SPAG5 may serve as a potential response marker for ICB therapy.

Lola-I is a promoter pioneer factor that establishes de novo Pol II pausing during development

While the accessibility of enhancers is dynamically regulated during development, promoters tend to be constitutively accessible and poised for activation by paused Pol II. By studying Lola-I, a Drosophila zinc finger transcription factor, we show here that the promoter state can also be subject to developmental regulation independently of gene activation. Lola-I is ubiquitously expressed at the end of embryogenesis and causes its target promoters to become accessible and acquire paused Pol II throughout the embryo. This promoter transition is required but not sufficient for tissue-specific target gene activation. Lola-I mediates this function by depleting promoter nucleosomes, similar to the action of pioneer factors at enhancers. These results uncover a level of regulation for promoters that is normally found at enhancers and reveal a mechanism for the de novo establishment of paused Pol II at promoters.

Reverse-ChIP Techniques for Identifying Locus-Specific Proteomes: A Key Tool in Unlocking the Cancer Regulome

A phenotypic hallmark of cancer is aberrant transcriptional regulation. Transcriptional regulation is controlled by a complicated array of molecular factors, including the presence of transcription factors, the deposition of histone post-translational modifications, and long-range DNA interactions. Determining the molecular identity and function of these various factors is necessary to understand specific aspects of cancer biology and reveal potential therapeutic targets. Regulation of the genome by specific factors is typically studied using chromatin immunoprecipitation followed by sequencing (ChIP-Seq) that identifies genome-wide binding interactions through the use of factor-specific antibodies. A long-standing goal in many laboratories has been the development of a 'reverse-ChIP' approach to identify unknown binding partners at loci of interest. A variety of strategies have been employed to enable the selective biochemical purification of sequence-defined chromatin regions, including single-copy loci, and the subsequent analytical detection of associated proteins. This review covers mass spectrometry techniques that enable quantitative proteomics before providing a survey of approaches toward the development of strategies for the purification of sequence-specific chromatin as a 'reverse-ChIP' technique. A fully realized reverse-ChIP technique holds great potential for identifying cancer-specific targets and the development of personalized therapeutic regimens.

Multi-landmark alignment of genomic signals reveals conserved expression patterns across transcription start sites

The prevalent one-dimensional alignment of genomic signals to a reference landmark is a cornerstone of current methods to study transcription and its DNA-dependent processes but it is prone to mask potential relations among multiple DNA elements. We developed a systematic approach to align genomic signals to multiple locations simultaneously by expanding the dimensionality of the genomic-coordinate space. We analyzed transcription in human and uncovered a complex dependence on the relative position of neighboring transcription start sites (TSSs) that is consistently conserved among cell types. The dependence ranges from enhancement to suppression of transcription depending on the relative distances to the TSSs, their intragenic position, and the transcriptional activity of the gene. Our results reveal a conserved hierarchy of alternative TSS usage within a previously unrecognized level of genomic organization and provide a general methodology to analyze complex functional relationships among multiple types of DNA elements.

Quantitative Analysis of Isoform Switching in Cancer

Over the past 8 years, multiple studies examined the phenomenon of isoform switching in human cancers and discovered that isoform switching is widespread, with hundreds to thousands of such events per cancer type. Although all of these studies used slightly different definitions of isoform switching, which in part led to a rather poor overlap of their results, they all leveraged transcript usage, a proportion of the transcript's expression in the total expression level of the parent gene, to detect isoform switching. However, how changes in transcript usage correlate with changes in transcript expression is not sufficiently explored. In this article, we adopt the most common definition of isoform switching and use a state-of-the-art tool for the analysis of differential transcript usage, SatuRn, to detect isoform switching events in 12 cancer types. We analyze the detected events in terms of changes in transcript usage and the relationship between transcript usage and transcript expression on a global scale. The results of our analysis suggest that the relationship between changes in transcript usage and changes in transcript expression is far from straightforward, and that such quantitative information can be effectively used for prioritizing isoform switching events for downstream analyses.

Integrative genomic analyses of promoter G-quadruplexes reveal their selective constraint and association with gene activation

G-quadruplexes (G4s) regulate DNA replication and gene transcription, and are enriched in promoters without fully appreciated functional relevance. Here we show high selection pressure on putative G4 (pG4) forming sequences in promoters through investigating genetic and genomic data. Analyses of 76,156 whole-genome sequences reveal that G-tracts and connecting loops in promoter pG4s display lower or higher allele frequencies, respectively, than pG4-flanking regions, and central guanines (Gs) in G-tracts show higher selection pressure than other Gs. Additionally, pG4-promoters produce over 72.4% of transcripts, and promoter G4-containing genes are expressed at relatively high levels. Most genes repressed by TMPyP4, a G4-ligand, regulate epigenetic processes, and promoter G4s are enriched with gene activation histone marks, chromatin remodeler and transcription factor binding sites. Consistently, cis-expression quantitative trait loci (cis-eQTLs) are enriched in promoter pG4s and their G-tracts. Overall, our study demonstrates selective constraint of promoter G4s and reinforces their stimulative role in gene expression.

Alternative promoters in CpG depleted regions are prevalently associated with epigenetic misregulation of liver cancer transcriptomes

Transcriptional regulation is commonly governed by alternative promoters. However, the regulatory architecture in alternative and reference promoters, and how they differ, remains elusive. In 100 CAGE-seq libraries from hepatocellular carcinoma patients, here we annotate 4083 alternative promoters in 2926 multi-promoter genes, which are largely undetected in normal livers. These genes are enriched in oncogenic processes and predominantly show association with overall survival. Alternative promoters are narrow nucleosome depleted regions, CpG island depleted, and enriched for tissue-specific transcription factors. Globally tumors lose DNA methylation. We show hierarchical retention of intragenic DNA methylation with CG-poor regions rapidly losing methylation, while CG-rich regions retain it, a process mediated by differential SETD2, H3K36me3, DNMT3B, and TET1 binding. This mechanism is validated in SETD2 knockdown cells and SETD2-mutated patients. Selective DNA methylation loss in CG-poor regions makes the chromatin accessible for alternative transcription. We show alternative promoters can control tumor transcriptomes and their regulatory architecture.

The landscape of cryptic antisense transcription in human cancers reveals an oncogenic noncoding RNA in lung cancer

Cryptic transcription initiation has been previously linked to activation of oncogenic transcripts. However, the prevalence and impact of cryptic antisense transcription from the opposite strand of protein-coding genes were mostly unknown in cancer. Applying a robust computational pipeline to publicly available transcriptome and epigenome datasets, we identified hundreds of previously unannotated cryptic antisense polyadenylated transcripts (CAPTs) that were enriched in tumor samples. We showed that the activation of cryptic antisense transcription was associated with increased chromatin accessibility and active histone marks. Accordingly, we found that many of the antisense transcripts were inducible by treatment of epigenetic drugs. Moreover, CRISPR-mediated epigenetic editing assays revealed that transcription of a noncoding RNA LRRK1-CAPT promoted LUSC cell proliferation, suggesting its oncogenic role. Our findings largely expand our understanding of cancer-associated transcription events, which may facilitate the development of novel strategies for cancer diagnosis and treatment.

Chromatin and noncoding RNA-mediated mechanisms of gastric tumorigenesis

Gastric cancer (GC) is one of the most common and deadly cancers in the world. It is a multifactorial disease highly influenced by environmental factors, which include radiation, smoking, diet, and infectious pathogens. Accumulating evidence suggests that epigenetic regulators are frequently altered in GC, playing critical roles in gastric tumorigenesis. Epigenetic regulation involves DNA methylation, histone modification, and noncoding RNAs. While it is known that environmental factors cause widespread alterations in DNA methylation, promoting carcinogenesis, the chromatin- and noncoding RNA-mediated mechanisms of gastric tumorigenesis are still poorly understood. In this review, we focus on discussing recent discoveries addressing the roles of histone modifiers and noncoding RNAs and the mechanisms of their interactions in gastric tumorigenesis. A better understanding of epigenetic regulation would likely facilitate the development of novel therapeutic approaches targeting specific epigenetic regulators in GC.

Integrative analysis of Iso-Seq and RNA-seq reveals dynamic changes of alternative promoter, alternative splicing and alternative polyadenylation during Angiotensin II-induced senescence in rat primary aortic endothelial cells

In eukaryotes, alternative promoter (AP), alternative splicing (AS), and alternative polyadenylation (APA) are three crucial regulatory mechanisms that modulate message RNA (mRNA) diversity. Although AP, AS and APA are involved in diverse biological processess, whether they have dynamic changes in Angiotensin II (Ang II) induced senescence in rat primary aortic endothelial cells (RAECs), an important cellular model for studying cardiovascular disease, remains unclear. Here we integrated both PacBio single-molecule long-read isoform sequencing (Iso-Seq) and Illumina short-read RNA sequencing (RNA-seq) to analyze the changes of AP, AS and APA in Ang II-induced senescent RAECs. Iso-Seq generated 36,278 isoforms from 10,145 gene loci and 65.81% of these isoforms are novel, which were further cross-validated by public data obtained by other techonologies such as CAGE, PolyA-Seq and 3'READS. APA contributed most to novel isoforms, followed by AS and AP. Further investigation showed that AP, AS and APA could all contribute to the regulation of isoform, but AS has more dynamic changes compared to AP and APA upon Ang II stimulation. Genes undergoing AP, AS and APA in Ang II-treated cells are enriched in various pathways related to aging or senescence, suggesting that these molecular changes are involved in functional alterations during Ang II-induced senescence. Together, the present study largely improved the annotation of rat genome and revealed gene expression changes at isoform level, extending the understanding of the complexity of gene regulation in Ang II-treated RAECs, and also provided novel clues for discovering the regulatory mechanism undelying Ang II caused vascular senescence and diseases.

An anoikis-related gene signature predicts prognosis and reveals immune infiltration in hepatocellular carcinoma

Background

Hepatocellular carcinoma (HCC) is a global health burden with poor prognosis. Anoikis, a novel programmed cell death, has a close interaction with metastasis and progression of cancer. In this study, we aimed to construct a novel bioinformatics model for evaluating the prognosis of HCC based on anoikis-related gene signatures as well as exploring the potential mechanisms.

Materials and methods

We downloaded the RNA expression profiles and clinical data of liver hepatocellular carcinoma from TCGA database, ICGC database and GEO database. DEG analysis was performed using TCGA and verified in the GEO database. The anoikis-related risk score was developed via univariate Cox regression, LASSO Cox regression and multivariate Cox regression, which was then used to categorize patients into high- and low-risk groups. Then GO and KEGG enrichment analyses were performed to investigate the function between the two groups. CIBERSORT was used for determining the fractions of 22 immune cell types, while the ssGSEA analyses was used to estimate the differential immune cell infiltrations and related pathways. The "pRRophetic" R package was applied to predict the sensitivity of administering chemotherapeutic and targeted drugs.

Results

A total of 49 anoikis-related DEGs in HCC were detected and 3 genes (EZH2, KIF18A and NQO1) were selected out to build a prognostic model. Furthermore, GO and KEGG functional enrichment analyses indicated that the difference in overall survival between risk groups was closely related to cell cycle pathway. Notably, further analyses found the frequency of tumor mutations, immune infiltration level and expression of immune checkpoints were significantly different between the two risk groups, and the results of the immunotherapy cohort showed that patients in the high-risk group have a better immune response. Additionally, the high-risk group was found to have higher sensitivity to 5-fluorouracil, doxorubicin and gemcitabine.

Conclusion

The novel signature of 3 anoikis-related genes (EZH2, KIF18A and NQO1) can predict the prognosis of patients with HCC, and provide a revealing insight into personalized treatments in HCC.

A Guide to MethylationToActivity: A Deep Learning Framework That Reveals Promoter Activity Landscapes from DNA Methylomes in Individual Tumors

Genome-wide DNA methylomes have contributed greatly to tumor detection and subclassification. However, interpreting the biological impact of the DNA methylome at the individual gene level remains a challenge. MethylationToActivity (M2A) is a pipeline that uses convolutional neural networks to infer H3K4me3 and H3K27ac enrichment from DNA methylomes and thus infer promoter activity. It was shown to be highly accurate and robust in revealing promoter activity landscapes in various pediatric and adult cancers. The following will present a user-friendly guide through the model pipeline.

SMYD3 associates with the NuRD (MTA1/2) complex to regulate transcription and promote proliferation and invasiveness in hepatocellular carcinoma cells

BACKGROUND

SMYD3, a member of the SET and MYND domain-containing (SMYD) family, is a histone methyltransferase (HMT) and transcription factor that plays an important role in transcriptional regulation in human carcinogenesis.

RESULTS

Using affinity purification and mass spectrometry assays to identify SMYD3-associated proteins in hepatocellular carcinoma (HCC) cells, we found several previously undiscovered SMYD3-interacting proteins, including the NuRD (MTA1/2) complex, the METTL family, and the CRL4B complex. Transcriptomic analysis of the consequences of knocking down SMYD3, MTA1, or MTA2 in HCC cells showed that SMYD3/NuRD complex targets a cohort of genes, some of which are critically involved in cell growth and migration. qChIP analyses showed that SMYD3 knockdown led to a significant reduction in the binding of MTA1 or MTA2 to the promoters of IGFBP4 and led to a significant decrease in H4K20me3 and a marked increase in H4Ac at the IGFBP4 promoter. In addition, we demonstrated that SMYD3 promotes cell proliferation, invasion, and tumorigenesis in vivo and in vitro and found that its expression is markedly upregulated in human liver cancer. Knockdown of MTA1 or MTA2 had the same effect as knockdown of SMYD3 on proliferation and invasion of hepatocellular carcinoma cells. Catalytic mutant SMYD3 could not rescue the phenotypic effects caused by knockdown of SMYD3. Inhibitors of SMYD3 effectively inhibited the proliferation and invasiveness of HCC cells.

CONCLUSIONS

These findings revealed that SMYD3 could transcriptionally repress a cohort of target genes expression by associating with the NuRD (MTA1/2) complex, thereby promoting the proliferation and invasiveness of HCC cells. Our results support the case for pursuing SMYD3 as a practical prognostic marker or therapeutic target against HCC.

Classification and characterization of alternative promoters in 26 lung adenocarcinoma cell lines

BACKGROUND

Genome-wide landscape of alternative promoter use remains unknown. We determined expression profiles of promoters in 26 lung adenocarcinoma cell lines using the transcriptional start site-sequencing data and proposed an index 'canonical promoter usage' to quantify the diversity of alternative promoter usage.

METHODS

Transcriptional start site-sequencing and other datasets were obtained from the DataBase of Transcriptional Start Sites. Transcriptional start site-sequencing read clusters were mapped onto RefGene to determine the promoters. Commonly used promoters were designated as canonical promoters. The sequence logos, CpG islands, DNA methylation and histone modifications of canonical and non-canonical promoters were examined. Canonical promoter usage was calculated by dividing 'read counts of a canonical promoter' by 'read counts of all the units of promoters' on each gene. The expressed genes were subjected to hierarchical clustering according to their canonical promoter usage.

RESULTS

Among 104 455 promoters for 14 297 genes, 8659 canonical and 68 197 non-canonical promoters were identified. Corresponding to higher expression, canonical promoters showed core promoter sequences, higher CpG island positivity, less DNA methylation and higher transcription-promoting histone modifications. Gene ontology enrichment analysis revealed that the clusters with lower canonical promoter usage were related to signalling pathways, whereas clusters of tightly regulated genes with higher canonical promoter usage were related to housekeeping genes.

CONCLUSION

Canonical promoters were regulated by conventional transcriptional machinery, while non-canonical promoters would be targets of 'leaky' expression. Further investigation is warranted to analyse the correlation between alternative promoter usage and biological characteristics contributing to carcinogenesis.

An integrated analysis of the cancer genome atlas data discovers a hierarchical association structure across thirty three cancer types

Cancer cells harbor molecular alterations at all levels of information processing. Genomic/epigenomic and transcriptomic alterations are inter-related between genes, within and across cancer types and may affect clinical phenotypes. Despite the abundant prior studies of integrating cancer multi-omics data, none of them organizes these associations in a hierarchical structure and validates the discoveries in extensive external data. We infer this Integrated Hierarchical Association Structure (IHAS) from the complete data of The Cancer Genome Atlas (TCGA) and compile a compendium of cancer multi-omics associations. Intriguingly, diverse alterations on genomes/epigenomes from multiple cancer types impact transcriptions of 18 Gene Groups. Half of them are further reduced to three Meta Gene Groups enriched with (1) immune and inflammatory responses, (2) embryonic development and neurogenesis, (3) cell cycle process and DNA repair. Over 80% of the clinical/molecular phenotypes reported in TCGA are aligned with the combinatorial expressions of Meta Gene Groups, Gene Groups, and other IHAS subunits. Furthermore, IHAS derived from TCGA is validated in more than 300 external datasets including multi-omics measurements and cellular responses upon drug treatments and gene perturbations in tumors, cancer cell lines, and normal tissues. To sum up, IHAS stratifies patients in terms of molecular signatures of its subunits, selects targeted genes or drugs for precision cancer therapy, and demonstrates that associations between survival times and transcriptional biomarkers may vary with cancer types. These rich information is critical for diagnosis and treatments of cancers.

Mediators of Capillary-to-Venule Conversion in the Chronic Inflammatory Skin Disease Psoriasis

Psoriasis is a chronic inflammatory skin disease characterized by epidermal hyperplasia and hyperkeratosis, immune cell infiltration and vascular remodeling. Despite the emerging recognition of vascular normalization as a potential strategy for managing psoriasis, an in-depth delineation of the remodeled dermal vasculature has been missing. In this study, we exploited 5' single-cell RNA sequencing to investigate the transcriptomic alterations in different subpopulations of blood vascular and lymphatic endothelial cells directly isolated from psoriatic and healthy human skin. Individual subtypes of endothelial cells underwent specific molecular repatterning associated with cell adhesion and extracellular matrix organization. Blood capillaries, in particular, showed upregulation of the melanoma cell adhesion molecule as well as its binding partners and adopted postcapillary venule‒like characteristics during chronic inflammation that are more permissive to leukocyte transmigration. We also identified psoriasis-specific interactions between cis-regulatory enhancers and promoters for each endothelial cell subtype, revealing the dysregulated gene regulatory networks in psoriasis. Together, our results provide more insights into the specific transcriptional responses and epigenetic signatures of endothelial cells lining different vessel compartments in chronic skin inflammation.

Expression profile and prognostic values of SMC family members in HCC

OBJECTIVE

The structural maintenance of chromosome (SMC) gene family, including 6 proteins, is involved in a wide range of biological functions in different human cancers. Nevertheless, there is little research on the expression patterns, potential functions and prognostic value of SMC genes in hepatocellular carcinoma (HCC). Based on publicly available databases and integrative bioinformatics analysis, we tried to determine the value of SMC gene expression in predicting the risk of developing HCC.

METHODS

The expression and copy number variations data of SMC family members were obtained from TCGA (The Cancer Genome Atlas). We identified the prognostic values of SMC family members and their clinical features. GSEA (Gene Set Enrichment Analysis) was conducted to detect the mechanism underlying the involvement of SMC family members in liver cancer. We used Tumor Immune Estimation Resource database to explore the associations between TIICs (Tumor Immune Infiltrating Cells) and the SMC family members.

RESULTS

Our analysis proved that downregulation of SMC family members was common modification in HCC patients. In HCC, the expression of SMC1A, SMC2, SMC3, SMC4, SMC6 were upregulated. Upregulation of SMC2, SMC3, and SMC4, along with the clinical stage of HCC, were associated with a poor prognosis according to the results of univariate and multivariate Cox proportional hazards regression analysis. SMC2, SMC3, and SMC4 are also related to tumor purity and immune infiltration levels of HCC. The GSEA results proved that SMC family members take part in numerous biological processes underlying tumorigenesis.

CONCLUSION

In this study, we comprehensively analyzed the expression of SMC family members in patients with HCC. This can provide insights for further investigation of the SMC members as potential therapeutic targets in HCC and suggest that the use of SMC inhibitor targeting SMC2, SMC3, and SMC4 can be a practical strategy for the therapy of HCC.

Annotating TSSs in Multiple Cell Types Based on DNA Sequence and RNA-seq Data via DeeReCT-TSS

The accurate annotation of transcription start sites (TSSs) and their usage are critical for the mechanistic understanding of gene regulation in different biological contexts. To fulfill this, specific high-throughput experimental technologies have been developed to capture TSSs in a genome-wide manner, and various computational tools have also been developed for in silico prediction of TSSs solely based on genomic sequences. Most of these computational tools cast the problem as a binary classification task on a balanced dataset, thus resulting in drastic false positive predictions when applied on the genome scale. Here, we present DeeReCT-TSS, a deep learning-based method that is capable of identifying TSSs across the whole genome based on both DNA sequence and conventional RNA sequencing data. We show that by effectively incorporating these two sources of information, DeeReCT-TSS significantly outperforms other solely sequence-based methods on the precise annotation of TSSs used in different cell types. Furthermore, we develop a meta-learning-based extension for simultaneous TSS annotations on 10 cell types, which enables the identification of cell type-specific TSSs. Finally, we demonstrate the high precision of DeeReCT-TSS on two independent datasets by correlating our predicted TSSs with experimentally defined TSS chromatin states. The source code for DeeReCT-TSS is available at https://github.com/JoshuaChou2018/DeeReCT-TSS_release and https://ngdc.cncb.ac.cn/biocode/tools/BT007316.

Signature of gene expression profile of liver sinusoidal endothelial cells in nonalcoholic steatohepatitis

Background: There has been emerging evidence that liver sinusoidal endothelial cells (LSECs) play an important role in the pathogenesis of nonalcoholic steatohepatitis (NASH). This study aims to figure out the signature of the gene expression profile of LSECs in NASH and to explore potential biomarkers related to damaged LSECs in NASH.

Methods and materials: Animal experiments were performed to demonstrate the significant structural damage of LSECs in the NASH model. To further understand the functional changes of these damaged LSECs in NASH, we used the public GEO database that contained microarray data for the gene expression of LSECs in NASH and normal mouse liver. Differentially expressed genes (DEGs) were analyzed, and further Gene Ontology (GO) enrichment analysis was performed to understand the functional changes. The hub genes were then identified and validated via external GEO databases.

Results: There was significant structural damage to LSECs in the NASH model, accompanied by remarkable functional changes of LSECs with 174 DEGs (156 upregulated and 18 downregulated genes). The functions of these DEGs were mainly enriched in the inflammatory reactions and immune responses. Nine specifically expressed hub genes were identified. Among them, CCL4 and ITGAX showed the most significant correlation with NASH, with AUROC of 0.77 and 0.86, respectively. The protein–protein interaction network, mRNA–miRNA interaction network, and ceRNA network were further predicted.

Conclusion: LSECs show significant structural damage and functional changes in NASH. The LSEC-related DEGs, such as CCL4 and ITGAX, might be promising biomarkers as well as potential treatment targets for NASH.

TTC7B Is a Novel Prognostic-Related Biomarker in Glioma Correlating with Immune Infiltrates and Response to Oxidative Stress by Temozolomide

Background

Gliomas are one of the most prevalent malignant brain tumors. Hence, identifying biological markers for glioma is imperative. TTC7B (Tetratricopeptide Repeat Domain 7B) is a gene whose role in cancer in currently identified. To this end, we examined the TTC7B expression as well as its prognostic significance, biological roles, and immune system impacts in patients with glioma.

Methods

We evaluated the function of TTC7B in GBM and LGG through the published CGGA (Chinese Glioma Genome Atlas) and TCGA (The Cancer Genome Atlas) databases. CIBERSORT and TIMER were used to analyze the link between TTC7B and immune cells, while R was used for statistical analysis. In addition, Transwell analysis, including migration and invasion assays, was performed to identify the relationship between TTC7B and temozolomide.

Results

Low expression of TTC7B was observed in GBM and LGG. 1p/19q codeletion, IDH mutation, chemotherapy, and grade were found to have a significant correlation with TTC7B. Besides, low TTC7B expression was linked with low overall survival (OS) in both GBM and LGG. In the Cox analysis, TTC7B was found to independently function as a risk element for OS of patients with glioma. Furthermore, CIBERSORT analysis demonstrated a positive link between TTC7B and multiple immune cells, especially activated NK cells. Transwell analysis, including migration and invasion assays, revealed that temozolomide reduced the migration and invasion capacity of glioma cells and increased the expression of TTC7B.

Conclusion

In all, TTC7B could serve as a promising prognostic indicator of LGG and GBM, and is closely associated with immune infiltration and response to oxidative stress by temozolomide.

Integrated bioinformatics analysis for novel miRNAs markers and ceRNA network in diabetic retinopathy

In order to seek a more outstanding diagnosis and treatment of diabetic retinopathy (DR), we predicted the miRNA biomarkers of DR and explored the pathological mechanism of DR through bioinformatics analysis.

Method: Based on public omics data and databases, we investigated ncRNA (non-coding RNA) functions based on the ceRNA hypothesis.

Result: Among differentially expressed miRNAs (DE-miRNAs), hsa-miR-1179, -4797-3p and -665 may be diagnosis biomarkers of DR. Functional enrichment analysis revealed differentially expressed mRNAs (DE-mRNAs) enriched in mitochondrial transport, cellular respiration and energy derivation. 18 tissue/organ-specific expressed genes, 10 hub genes and gene cluster modules were identified. The ceRNA networks lncRNA FBXL19-AS1/miR-378f/MRPL39 and lncRNA UBL7-AS1/miR-378f/MRPL39 might be potential RNA regulatory pathways in DR.

Conclusion: Differentially expressed hsa-miR-1179, -4797-3p and -665 can be used as powerful markers for DR diagnosis, and the ceRNA network: lncRNA FBXL19-AS1/UBL7-AS1-miR-378f-MRPL39 may represent an important regulatory role in DR progression.