Common mitochondrial deletions in RNA-Seq: evaluation of bulk, single-cell, and spatial transcriptomic datasets.

Common mitochondrial DNA (mtDNA) deletions are large structural variants in the mitochondrial genome that accumulate in metabolically active tissues with age and have been investigated in various diseases. We applied the Splice-Break2 pipeline (designed for high-throughput quantification of mtDNA deletions) to human RNA-Seq datasets and describe the methodological considerations for evaluating common deletions in bulk, single-cell, and spatial transcriptomics datasets. A robust evaluation of 1570 samples from 14 RNA-Seq studies showed: (i) the abundance of some common deletions detected in PCR-amplified mtDNA correlates with levels observed in RNA-Seq data; (ii) RNA-Seq library preparation method has a strong effect on deletion detection; (iii) deletions had a significant, positive correlation with age in brain and muscle; (iv) deletions were enriched in cortical grey matter, specifically in layers 3 and 5; and (v) brain regions with dopaminergic neurons (i.e., substantia nigra, ventral tegmental area, and caudate nucleus) had remarkable enrichment of common mtDNA deletions.

Large Common Mitochondrial DNA Deletions Are Associated with a Mitochondrial SNP T14798C Near the 3' Breakpoints.

Introduction: Large somatic deletions of mitochondrial DNA (mtDNA) accumulate with aging in metabolically active tissues such as the brain. We have cataloged the breakpoints and frequencies of large mtDNA deletions in the human brain. Methods: We quantified 112 high-frequency mtDNA somatic deletions across four human brain regions with the Splice-Break2 pipeline. In addition, we utilized PLINK/Seq to test the association of mitochondrial genotypes with the abundance of these high-frequency mtDNA deletions. A conservative p value threshold of 5E-08 was used to find the significant loci. Results: One mtDNA SNP (T14798C) was significantly associated with mtDNA deletions in two brain regions, the dorsolateral prefrontal cortex (DLPFC) and the superior temporal gyrus. Since the DLPFC showed the most robust association between T14798C and two deletion breakpoints (7816-14807 and 5462-14807), this association was tested in the DLPFC of a replication sample and validated the first results. Incorporating the C allele at 14,798 bp increased the perfect/imperfect length of the repeat at the 3' breakpoint of the two associated deletions. Conclusion: This is the first study to identify the association of mtDNA SNP with large mtDNA deletions in the human brain. The T14798C allele located in the MT-CYB gene is a common polymorphism that occurs in several mitochondrial haplogroups. We hypothesize that the T14798C association with two deletions occurs by extending the repeat length around the 3' deletion breakpoints. This simple mechanism suggests that mtDNA SNPs can affect the mitochondrial genome structure, especially in brain where high levels of reactive oxygen species lead to deletion accumulation with aging.

Relapse-associated AURKB blunts the glucocorticoid sensitivity of B cell acute lymphoblastic leukemia.

Glucocorticoids (GCs) are used in combination chemotherapies as front-line treatment for B cell acute lymphoblastic leukemia (B-ALL). Although effective, many patients relapse and become resistant to chemotherapy and GCs in particular. Why these patients relapse is not clear. We took a comprehensive, functional genomics approach to identify sources of GC resistance. A genome-wide shRNA screen identified the transcriptional coactivators EHMT2, EHMT1, and CBX3 as important contributors to GC-induced cell death. This complex selectively supports GC-induced expression of genes contributing to cell death. A metaanalysis of gene expression data from B-ALL patient specimens revealed that Aurora kinase B (AURKB), which restrains GC signaling by phosphorylating EHMT1-2, is overexpressed in relapsed B-ALL, suggesting it as a potential contributor to relapse. Inhibition of AURKB enhanced GC-induced expression of cell death genes, resulting in potentiation of GC cytotoxicity in cell lines and relapsed B-ALL patient samples. This function for AURKB is distinct from its canonical role in the cell cycle. These results show the utility of functional genomics in understanding mechanisms of resistance and rapidly identifying combination chemotherapeutics.

A post-translational modification switch controls coactivator function of histone methyltransferases G9a and GLP.

Like many transcription regulators, histone methyltransferases G9a and G9a-like protein (GLP) can act gene-specifically as coregulators, but mechanisms controlling this specificity are mostly unknown. We show that adjacent post-translational methylation and phosphorylation regulate binding of G9a and GLP to heterochromatin protein 1 gamma (HP1γ), formation of a ternary complex with the glucocorticoid receptor (GR) on chromatin, and function of G9a and GLP as coactivators for a subset of GR target genes. HP1γ is recruited by G9a and GLP to GR binding sites associated with genes that require G9a, GLP, and HP1γ for glucocorticoid-stimulated transcription. At the physiological level, G9a and GLP coactivator function is required for glucocorticoid activation of genes that repress cell migration in A549 lung cancer cells. Thus, regulated methylation and phosphorylation serve as a switch controlling G9a and GLP coactivator function, suggesting that this mechanism may be a general paradigm for directing specific transcription factor and coregulator actions on different genes.

Selective coregulator function and restriction of steroid receptor chromatin occupancy by Hic-5.

Steroid receptors (SRs) bind specific DNA regulatory sequences, thereby activating and repressing gene expression. We previously showed that transcriptional coregulator Hic-5 facilitates glucocorticoid regulation of some genes but blocks glucocorticoid regulation of others. Here, in a genome-wide analysis, Hic-5 depletion dramatically increased the global number of sites occupied by glucocorticoid receptor (GR) α (the major GR isoform), and many binding sites blocked by Hic-5 were associated with genes for which Hic-5 also blocked glucocorticoid-regulated expression. Hic-5 had similar effects on GRγ (a splice variant of GRα) and estrogen receptor α (ERα), facilitating hormonal regulation of some genes and blocking hormonal regulation of others. As with GRα, Hic-5 blocking of hormonal gene regulation mediated by GRγ and ERα was associated with blocking of GRγ and ERα occupancy at nearby sites. Hic-5 supported hormonal regulation of many more genes for GRα than for GRγ or ERα and thus exhibited selective coregulator functions for different SRs. In contrast, the number of Hic-5-blocked genes was similar for all 3 SRs. In addition to classic coregulator activity, Hic-5 influences the genomic occupancy of multiple SRs and thereby blocks some aspects of hormonal regulation. Thus, Hic-5, because of its tissue-specific expression, could contribute to tissue-specific genomic occupancy and gene regulation by SRs.

Role of distinct surfaces of the G9a ankyrin repeat domain in histone and DNA methylation during embryonic stem cell self-renewal and differentiation.

BACKGROUND: Epigenetic modifications such as histone and DNA methylation are essential for silencing pluripotency genes during embryonic stem cell (ESC) differentiation. G9a is the major histone H3 Lys9 (H3K9) methyltransferase in euchromatin and is required for the de novo DNA methylation of the key regulator of pluripotency Oct3/4 during ESC differentiation. Surprisingly, the catalytic activity of G9a is not required for its role in de novo DNA methylation and the precise molecular mechanisms of G9a in this process are poorly understood. It has been suggested that the G9a ankyrin repeat domain, which can interact with both H3K9me2 and the DNA methyltransferase DNMT3A, could facilitate de novo DNA methylation by bridging the interaction between DNMT3A and H3K9me2-marked chromatin. RESULTS: Here, we demonstrate that the G9a ankyrin domain H3K9me2-binding function is not required for the de novo DNA methylation of Oct3/4 during ESC differentiation. Moreover, we show that the interaction between the G9a ankyrin domain and DNMT3A is not sufficient to ensure efficient de novo DNA methylation. More importantly, we characterize a specific residue of the G9a ankyrin domain (Asp905) that is critical for both maintaining cellular H3K9me2 levels in undifferentiated ESCs and for the establishment of de novo DNA methylation during differentiation. CONCLUSIONS: These results represent an exciting breakthrough, which reveals 1) an unexpected critical biological function of the G9a ankyrin domain in global histone H3K9 methylation and 2) valuable insights into the molecular mechanisms and interaction surfaces through which G9a regulates de novo DNA methylation of Oct3/4 during ESC differentiation.

G9a functions as a molecular scaffold for assembly of transcriptional coactivators on a subset of glucocorticoid receptor target genes.

Histone H3 lysine-9 methyltransferase G9a/EHMT2/KMT1C is a key corepressor of gene expression. However, activation of a limited number of genes by G9a (independent of its catalytic activity) has also been observed, although the precise molecular mechanisms are unknown. By using RNAi in combination with gene expression microarray analysis, we found that G9a functions as a positive and a negative transcriptional coregulator for discrete subsets of genes that are regulated by the hormone-activated Glucocorticoid Receptor (GR). G9a was recruited to GR-binding sites (but not to the gene body) of its target genes and interacted with GR, suggesting recruitment of G9a by GR. In contrast to its corepressor function, positive regulation of gene expression by G9a involved G9a-mediated enhanced recruitment of coactivators CARM1 and p300 to GR target genes. Further supporting a role for G9a as a molecular scaffold for its coactivator function, the G9a-specific methyltransferase inhibitor UNC0646 did not affect G9a coactivator function but selectively decreased G9a corepressor function for endogenous target genes. Overall, G9a functioned as a coactivator for hormone-activated genes and as a corepressor in support of hormone-induced gene repression, suggesting that the positive or negative actions of G9a are determined by the gene-specific regulatory environment and chromatin architecture. These findings indicate distinct mechanisms of G9a coactivator vs. corepressor functions in transcriptional regulation and provide insight into the molecular mechanisms of G9a coactivator function. Our results also suggest a physiological role of G9a in fine tuning the set of genes that respond to glucocorticoids.

Lysine methyltransferase G9a is not required for DNMT3A/3B anchoring to methylated nucleosomes and maintenance of DNA methylation in somatic cells.

BACKGROUND: DNA methylation, histone modifications and nucleosome occupancy act in concert for regulation of gene expression patterns in mammalian cells. Recently, G9a, a H3K9 methyltransferase, has been shown to play a role in establishment of DNA methylation at embryonic gene targets in ES cells through recruitment of de novo DNMT3A/3B enzymes. However, whether G9a plays a similar role in maintenance of DNA methylation in somatic cells is still unclear. RESULTS: Here we show that G9a is not essential for maintenance of DNA methylation in somatic cells. Knockdown of G9a has no measurable effect on DNA methylation levels at G9a-target loci. DNMT3A/3B remain stably anchored to nucleosomes containing methylated DNA even in the absence of G9a, ensuring faithful propagation of methylated states in cooperation with DNMT1 through somatic divisions. Moreover, G9a also associates with nucleosomes in a DNMT3A/3B and DNA methylation-independent manner. However, G9a knockdown synergizes with pharmacologic inhibition of DNMTs resulting in increased hypomethylation and inhibition of cell proliferation. CONCLUSIONS: Taken together, these data suggest that G9a is not involved in maintenance of DNA methylation in somatic cells but might play a role in re-initiation of de novo methylation after treatment with hypomethylating drugs, thus serving as a potential target for combinatorial treatments strategies involving DNMTs inhibitors.

A distinct mechanism for coactivator versus corepressor function by histone methyltransferase G9a in transcriptional regulation.

Histone methyltransferase G9a has been understood primarily as a corepressor of gene expression, but we showed previously that G9a positively regulates nuclear receptor-mediated transcription in reporter gene assays. Here, we show that endogenous G9a contributes to the estradiol (E(2))-dependent induction of some endogenous target genes of estrogen receptor (ER)α in MCF-7 breast cancer cells while simultaneously limiting the E(2)-induced expression of other ERα target genes. Thus, G9a has a dual and selective role as a coregulator for ERα target genes. The ERα binding regions associated with the pS2 gene, which requires G9a for E(2)-induced expression, are transiently occupied by G9a at 15 min after beginning E(2) treatment, suggesting that G9a coactivator function is by direct interaction with ERα target genes. Transient reporter gene assays with deletion mutants of G9a demonstrated that domains previously associated with the corepressor functions of G9a (C-terminal methyltransferase domain, ankyrin repeat domain, and cysteine-rich domain) were unnecessary for G9a coactivator function in ERα-mediated transcription. In contrast, the N-terminal domain of G9a was necessary and sufficient for enhancement of ERα-mediated transcription and for E(2)-induced occupancy of G9a on ERα binding sites associated with endogenous target genes of ERα. In addition to a previously identified activation domain, this region contains a previously uncharacterized ligand-dependent ERα binding function, indicating how G9a is recruited to the target genes. Therefore, the coactivator and corepressor functions of G9a involve different G9a domains and different molecular mechanisms.

A coactivator role of CARM1 in the dysregulation of ß-catenin activity in colorectal cancer cell growth and gene expression.

Aberrant activation of Wnt/ß-catenin signaling, resulting in the expression of Wnt-regulated oncogenes, is recognized as a critical factor in the etiology of colorectal cancer. Occupancy of ß-catenin at promoters of Wnt target genes drives transcription, but the mechanism of ß-catenin action remains poorly understood. Here, we show that CARM1 (coactivator-associated arginine methyltransferase 1) interacts with ß-catenin and positively modulates ß-catenin-mediated gene expression. In colorectal cancer cells with constitutively high Wnt/ß-catenin activity, depletion of CARM1 inhibits expression of endogenous Wnt/ß-catenin target genes and suppresses clonal survival and anchorage-independent growth. We also identified a colorectal cancer cell line (RKO) with a low basal level of ß-catenin, which is dramatically elevated by treatment with Wnt3a. Wnt3a also increased the expression of a subset of endogenous Wnt target genes, and CARM1 was required for the Wnt-induced expression of these target genes and the accompanying dimethylation of arginine 17 of histone H3. Depletion of ß-catenin from RKO cells diminished the Wnt-induced occupancy of CARM1 on a Wnt target gene, indicating that CARM1 is recruited to Wnt target genes through its interaction with ß-catenin and contributes to transcriptional activation by mediating events (including histone H3 methylation) that are downstream from the actions of ß-catenin. Therefore, CARM1 is an important positive modulator of Wnt/ß-catenin transcription and neoplastic transformation, and may thereby represent a novel target for therapeutic intervention in cancers involving aberrantly activated Wnt/ß-catenin signaling.