A Genome-Wide Investigation of Effects of Aberrant DNA Methylation on the Usage of Alternative Promoters in Hepatocellular Carcinoma

BackgroundThe alternative usage of promoters provides a way to regulate gene expression, has a significant influence on the transcriptome, and contributes to the cellular transformation of cancer. However, the function of alternative promoters (APs) in hepatocellular carcinoma (HCC) has not been systematically studied yet. In addition, the potential mechanism of regulation to the usage of APs remains unclear. DNA methylation, one of the most aberrant epigenetic modifications in cancers, is known to regulate transcriptional activity. Whether DNA methylation regulates the usage of APs needs to be explored. Here, we aim to investigate the effects of DNA methylation on usage of APs in HCC.MethodsPromoter activities were calculated based on RNA-seq data. Functional enrichment analysis was implemented to conduct GO terms. Correlation tests were used to detect the correlation between promoter activity and methylation status. The LASSO regression model was used to generate a diagnostic model. Kaplan–Meier analysis was used to compare the overall survival between high and low methylation groups. RNA-seq and whole-genome bisulfite sequencing (WGBS) in HCC samples were performed to validate the correlation of promoter activity and methylation.ResultsWe identified 855 APs in total, which could be well used to distinguish cancer from normal samples. The correlation of promoter activity and DNA methylation in APs was observed, and the APs with negative correlation were defined as methylation-regulated APs (mrAPs). Six mrAPs were identified to generate a diagnostic model with good performance (AUC = 0.97). Notably, the majority of mrAPs had CpG sites that could be used to predict clinical outcomes by methylation status. Finally, we verified 85.6% of promoter activity variation and 92.3% of methylation changes in our paired RNA-seq and WGBS samples, respectively. The negative correlation between promoter activity and methylation status was further confirmed in our HCC samples.ConclusionThe aberrant methylation status plays a critical role in the precision usage of APs in HCC, which sheds light on the mechanism of cancer development and provides a new insight into cancer screening and treatment.

SMAD4–201 transcript as a putative biomarker in colorectal cancer

Background Transcripts with alternative 5′-untranslated regions (UTRs) result from the activity of alternative promoters and they can determine gene expression by influencing its stability and translational efficiency, thus executing complex regulation of developmental, physiological and pathological processes. Transcriptional regulation of human SMAD4, a key tumor suppressor deregulated in most gastrointestinal cancers, entails four alternative promoters. These promoters and alternative transcripts they generate remain unexplored as contributors to the SMAD4 deregulation in cancer. The aim of this study was to investigate the relative abundance of the transcript SMAD4–201 in colorectal cell lines and tissues in order to establish if its fluctuations may be associated with colorectal cancer (CRC). Methods Relative abundance of SMAD4–201 in total SMAD4 mRNA was analyzed using quantitative PCR in a set of permanent human colon cell lines and tumor and corresponding healthy tissue samples from patients with CRC. Results The relative abundance of SMAD4–201 in analyzed cell lines varied between 16 and 47%. A similar relative abundance of SMAD4–201 transcript was found in the majority of analyzed human tumor tissue samples, and it was averagely 20% lower in non-malignant in comparison to malignant tissue samples (p = 0.001). Transcript SMAD4–202 was not detectable in any of the analyzed samples, so the observed fluctuations in the composition of SMAD4 transcripts can be attributed to transcripts other than SMAD4–201 and SMAD4–202. Conclusion The expression profile of SMAD4–201 in human tumor and non-tumor tissue samples may indicate the translational potential of this molecule in CRC, but further research is needed to clarify its usability as a potential biomarker for early diagnosis.

Epigenetics of Mitochondria-Associated Genes in Striated Muscle

Striated muscle has especially large energy demands. We identified 97 genes preferentially expressed in skeletal muscle and heart, but not in aorta, and found significant enrichment for mitochondrial associations among them. We compared the epigenomic and transcriptomic profiles of the 27 genes associated with striated muscle and mitochondria. Many showed strong correlations between their tissue-specific transcription levels, and their tissue-specific promoter, enhancer, or open chromatin as well as their DNA hypomethylation. Their striated muscle-specific enhancer chromatin was inside, upstream, or downstream of the gene, throughout much of the gene as a super-enhancer (CKMT2, SLC25A4, and ACO2), or even overlapping a neighboring gene (COX6A2, COX7A1, and COQ10A). Surprisingly, the 3′ end of the 1.38 Mb PRKN (PARK2) gene (involved in mitophagy and linked to juvenile Parkinson’s disease) displayed skeletal muscle/myoblast-specific enhancer chromatin, a myoblast-specific antisense RNA, as well as brain-specific enhancer chromatin. We also found novel tissue-specific RNAs in brain and embryonic stem cells within PPARGC1A (PGC-1α), which encodes a master transcriptional coregulator for mitochondrial formation and metabolism. The tissue specificity of this gene’s four alternative promoters, including a muscle-associated promoter, correlated with nearby enhancer chromatin and open chromatin. Our in-depth epigenetic examination of these genes revealed previously undescribed tissue-specific enhancer chromatin, intragenic promoters, regions of DNA hypomethylation, and intragenic noncoding RNAs that give new insights into transcription control for this medically important set of genes.

Automated identification of cell-type–specific genes and alternative promoters

Background: Identifying key transcriptional features, such as genes or transcripts, involved in cellular differentiation remains a challenging problem. Current methods for identifying key transcriptional features predominantly rely on pairwise comparisons among different cell types. These methods also identify long lists of differentially expressed transcriptional features. Combining the results from many such pairwise comparisons to find the transcriptional features specific only to one cell type is not straightforward. Thus, one must have a principled method for amalgamating pairwise cell type comparisons that makes full use of prior knowledge about the developmental relationships between cell types. Method: We developed Cell Lineage Analysis (CLA), a computational method which identifies transcriptional features with expression patterns that discriminate cell types, incorporating Cell Ontology knowledge on the relationship between different cell types. CLA uses random forest classification with a stratified bootstrap to increase the accuracy of binary classifiers when each cell type have a different number of samples. Regularized random forest results in a classifier that selects few but important transcriptional features. For each cell type pair, CLA runs multiple instances of regularized random forest and reports the transcriptional features consistently selected. CLA not only discriminates individual cell types but can also discriminate lineages of cell types related in the developmental hierarchy. Results: We applied CLA to Functional Annotation of the Mammalian Genome 5 (FANTOM5) data and identified discriminative transcription factor and long non-coding RNA (lncRNA) genes for 71 human cell types. With capped analysis of gene expression (CAGE) data, CLA identified individual cell-type–specific alternative promoters for cell surface markers. Compared to random forest with a standard bootstrap approach, CLA's stratified bootstrap approach improved the accuracy of gene expression classification models for more than 95% of 2060 cell type pairs examined. Applied on 10X Genomics single-cell RNA-seq data for CD14+ monocytes and FCGR3A+ monocytes, CLA selected only 13 discriminative genes. These genes included the top 9 out of 370 significantly differentially expressed genes obtained from conventional differential expression analysis methods. Discussion: Our CLA method combines tools to simplify the interpretation of transcriptome datasets from many cell types. It automates the identification of the most differentially expressed genes for each cell type pairs CLA's lineage score allows easy identification of the best transcriptional markers for each cell type and lineage in both bulk and single-cell transcriptomic data. Availability: CLA is available at https://cla.hoffmanlab.org. We deposited the version of the CLA source with which we ran our experiments at https://doi.org/10.5281/zenodo.3630670. We deposited other analysis code and results at https://doi.org/10.5281/zenodo.5735636.

Concordant Androgen-Regulated Expression of Divergent Rhox5 Promoters in Sertoli Cells

Concordant transcriptional regulation can generate multiple gene products that collaborate to achieve a common goal. Here we report a case of concordant transcriptional regulation that instead drives a single protein to be produced in the same cell type from divergent promoters. This gene product—the RHOX5 homeobox transcription factor—is translated from 2 different mRNAs with different 5′ untranslated regions (UTRs) transcribed from alternative promoters. Despite the fact that these 2 promoters—the proximal promoter (Pp) and the distal promoter (Pd)—exhibit different patterns of tissue-specific activity, share no obvious sequence identity, and depend on distinct transcription factors for expression, they exhibit a remarkably similar expression pattern in the testes. In particular, both depend on androgen signaling for expression in the testes, where they are specifically expressed in Sertoli cells and have a similar stage-specific expression pattern during the seminiferous epithelial cycle. We report evidence for 3 mechanisms that collaborate to drive concordant Pp/Pd expression. First, both promoters have an intrinsic ability to respond to androgen receptor and androgen. Second, the Pp acts as an enhancer to promote androgen-dependent transcription from the Pd. Third, Pd transcription is positively autoregulated by the RHOX5 protein, which is first produced developmentally from the Pp. Together, our data support a model in which the Rhox5 homeobox gene evolved multiple mechanisms to activate both of its promoters in Sertoli cells to produce Rhox5 in an androgen-dependent manner during different phases of spermatogenesis.

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

A continuous increase in the prevalence of heart failure and the lack of adequate therapy highlight poor understanding of the underlying genetic regulatory mechanisms involved in heart failure pathogenesis. Growing evidence has demonstrated a significant contribution of non-coding genome regulatory elements towards transcriptomic changes in heart disease. Thus, there is a pressing need for a comprehensive resource of the human cardiac regulatory network in healthy and failing states. We applied cap analysis of gene expression sequencing to directly measure the expression of RNA associated with enhancers and promoters. Based on this data, we constructed the atlas of transcribed cardiac regulatory elements from 21 healthy and 10 failing (ischemic and non-ischemic cardiomyopathy) human hearts. In total, we have sequenced 109 samples from the left and right atria and ventricles, identifying 17,668 promoters and 14,920 enhancers associated with 14,519 genes. Leveraging this atlas, we provide insights into functional and structural regulatory changes between healthy and failing hearts. Healthy atria and ventricles had distinct pathway enrichment and transcription factor binding patterns, significantly remodeled by heart failure. Using the advantages of deep sequencing that allow effective analysis of cis-regulatory elements-derived RNA, we found that heart failure is associated with the expression of transcripts derived from alternative promoters and a specific set of transcribed enhancers. Furthermore, we identified a high prevalence of single nucleotide polymorphisms associated with cardiovascular diseases within the regulatory regions highlighting their importance in disease pathogenesis. This open-source atlas will serve the cardiovascular community to improve understanding cardiac regulatory network and facilitate the development of novel therapeutics.

Transposable elements and their KZFP controllers are drivers of transcriptional innovation in the developing human brain

Transposable elements (TEs) account for more than 50% of the human genome and many have been co-opted throughout evolution to provide regulatory functions for gene expression networks. Several lines of evidence suggest that these networks are fine-tuned by the largest family of TE controllers, the KRAB-containing zinc finger proteins (KZFPs). One tissue permissive for TE transcriptional activation (termed “transposcription”) is the adult human brain, however comprehensive studies on the extent of this process and its potential contribution to human brain development are lacking. To elucidate the spatiotemporal transposcriptome of the developing human brain, we have analyzed two independent RNA-seq data sets encompassing 16 brain regions from eight weeks postconception into adulthood. We reveal a distinct KZFP:TE transcriptional profile defining the late prenatal to early postnatal transition, and the spatiotemporal and cell type–specific activation of TE-derived alternative promoters driving the expression of neurogenesis-associated genes. Long-read sequencing confirmed these TE-driven isoforms as significant contributors to neurogenic transcripts. We also show experimentally that a co-opted antisense L2 element drives temporal protein relocalization away from the endoplasmic reticulum, suggestive of novel TE dependent protein function in primate evolution. This work highlights the widespread dynamic nature of the spatiotemporal KZFP:TE transcriptome and its importance throughout TE mediated genome innovation and neurotypical human brain development. To facilitate interactive exploration of these spatiotemporal gene and TE expression dynamics, we provide the “Brain TExplorer” web application freely accessible for the community.

Differential regulation of alternative promoters emerges from unified kinetics of enhancer-promoter interaction

Many eukaryotic genes contain alternative promoters with distinct expression patterns. How these promoters are differentially regulated remains elusive. Here, we apply single-molecule imaging to quantify the transcriptional regulation of two alternative promoters (P1 and P2) of the Bicoid (Bcd) target gene hunchback in syncytial blastoderm Drosophila embryos. Contrary to the previous notion that Bcd only activates P2, we find that Bcd activates both promoters via the same two enhancers. P1 activation is less frequent and requires binding of more Bcd molecules than P2 activation. Using a theoretical model to relate promoter activity to enhancer states, we show that the two promoters follow common transcription kinetics driven by sequential Bcd binding at the two enhancers. Bcd binding at either enhancer primarily activates P2, while P1 activation relies more on Bcd binding at both enhancers. These results provide a quantitative framework for understanding the dynamics of complex eukaryotic gene regulation.

The Ends Dictate the Means: Promoter Switching in Herpesvirus Gene Expression

Herpesvirus gene expression is dynamic and complex, with distinct complements of viral genes expressed at specific times in different infection contexts. These complex patterns of viral gene expression arise in part from the integration of multiple cellular and viral signals that affect the transcription of viral genes. The use of alternative promoters provides an increased level of control, allowing different promoters to direct the transcription of the same gene in response to distinct temporal and contextual cues. While once considered rare, herpesvirus alternative promoter usage was recently found to be far more pervasive and impactful than previously thought. Here we review several examples of promoter switching in herpesviruses and discuss the functional consequences on the transcriptional and post-transcriptional regulation of viral gene expression. Expected final online publication date for the Annual Review of Virology, Volume 8 is September 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.