Histone lysine methyltransferases (HMTs), a large class of enzymes that catalyze site-specific methylation of lysine residues on histones and other proteins, play critical roles in controlling transcription, chromatin architecture, and cellular differentiation. However, the genomic landscape and clinical significance of HMTs in breast cancer remain poorly characterized. Here, we conducted a meta-analysis of approximately 50 HMTs in breast cancer and identified associations among recurrent copy number alterations, mutations, gene expression, and clinical outcome. We identified 12 HMTs with the highest frequency of genetic alterations, including 8 with high-level amplification, 2 with putative homozygous deletion, and 2 with somatic mutation. Different subtypes of breast cancer have different patterns of copy number and expression for each HMT gene. In addition, chromosome 1q contains four HMTs that are concurrently or independently amplified or overexpressed in breast cancer. Copy number or mRNA expression of several HMTs was significantly associated with basal-like breast cancer and shorter patient survival. Integrative analysis identified 8 HMTs (SETDB1, SMYD3, ASH1L, SMYD2, WHSC1L1, SUV420H1, SETDB2, and KMT2C) that are dysregulated by genetic alterations, classifying them as candidate therapeutic targets. Together, our findings provide a strong foundation for further mechanistic research and therapeutic options using HMTs to treat breast cancer.

Histone lysine demethylases (KDMs) comprise a large class of enzymes that catalyze site-specific demethylation of lysine residues on histones and other proteins. They play critical roles in controlling transcription, chromatin architecture, and cellular differentiation. However, the genomic landscape and clinical significance of KDMs in breast cancer remain poorly characterized. Here, we conducted a meta-analysis of 24 KDMs in breast cancer and identified associations among recurrent copy number alterations, gene expression, breast cancer subtypes, and clinical outcome. Two KDMs, KDM2A and KDM5B, had the highest frequency of genetic amplification and overexpression. Furthermore, among the 24 KDM genes, KDM2A had the highest correlation between copy number and mRNA expression, and high mRNA levels of KDM2A were significantly associated with shorter survival of breast cancer patients. KDM2A has two isoforms: the long isoform is comprised of a JmjC domain, CXXC-zinc finger, PHD zinc finger, F-box, and the AMN1 protein domain; whereas the short isoform of KDM2A lacks the N-terminal JmjC domain but contains all other motifs. Detailed characterization of KDM2A in breast cancer revealed that the short isoform of KDM2A is more abundant than the long isoform at DNA, mRNA, and protein levels in a subset of breast cancers. Furthermore, our data indicate that the short isoform of KDM2A has oncogenic potential and functions as an oncogenic isoform in a subset of breast cancers. Taken together, our findings suggest that amplification and overexpression of the KDM2A short isoform is critical in breast cancer progression.

Tudor domain–containing proteins (TDRDs), which recognize and bind to methyl-lysine/arginine residues on histones and non-histone proteins, play critical roles in regulating chromatin architecture, transcription, genomic stability, and RNA metabolism. Dysregulation of several TDRDs have been observed in various types of cancer. However, neither the genomic landscape nor clinical significance of TDRDs in breast cancer has been explored comprehensively. Here, we performed an integrated genomic and transcriptomic analysis of 41 TDRD genes in breast cancer (TCGA and METABRIC datasets) and identified associations among recurrent copy number alterations, gene expressions, clinicopathological features, and survival of patients. Among seven TDRDs that had the highest frequency (>10%) of gene amplification, the plant homeodomain finger protein 20-like 1 (PHF20L1) was the most commonly amplified (17.62%) TDRD gene in TCGA breast cancers. Different subtypes of breast cancer had different patterns of copy number and expression for each TDRD. Notably, amplification and overexpression of PHF20L1 were more prevalent in aggressive basal-like and Luminal B subtypes and were significantly associated with shorter survival of breast cancer patients. Furthermore, knockdown of PHF20L1 inhibited cell proliferation in PHF20L1-amplified breast cancer cell lines. PHF20L1 protein contains N-terminal Tudor and C-terminal plant homeodomain domains. Detailed characterization of PHF20L1 in breast cancer revealed that the Tudor domain likely plays a critical role in promoting cancer. Mechanistically, PHF20L1 might participate in regulating DNA methylation by stabilizing DNA methyltransferase 1 (DNMT1) protein in breast cancer. Thus, our results demonstrated the oncogenic potential of PHF20L1 and its association with poor prognostic parameters in breast cancer.

Chromodomain helicase DNA binding proteins (CHDs) are characterized by N‐terminal tandem chromodomains and a central adenosine triphosphate‐dependent helicase domain. CHDs govern the cellular machinery's access to DNA, thereby playing critical roles in various cellular processes including transcription, proliferation, and DNA damage repair. Accumulating evidence demonstrates that mutation and dysregulation of CHDs are implicated in the pathogenesis of developmental disorders and cancer. However, we know little about genomic and transcriptomic alterations and the clinical significance of most CHDs in human cancer. We used TCGA and METABRIC datasets to perform integrated genomic and transcriptomic analyses of nine CHD genes in more than 10 000 primary cancer specimens from 32 tumor types, focusing on breast cancers. We identified associations among recurrent copy number alteration, gene expression, clinicopathological features, and patient survival. We found that CHD7 was the most commonly gained/amplified and mutated, whereas CHD3 was the most deleted across the majority of tumor types, including breast cancer. Overexpression of CHD7 was more prevalent in aggressive subtypes of breast cancer and was significantly correlated with high tumor grade and poor prognosis. CHD7 is required to maintain open, accessible chromatin, thus providing fine‐tuning of transcriptional regulation of certain classes of genes. We found that CHD7 expression was positively correlated with a small subset of classical oncogenes, notably NRAS, in breast cancer. Knockdown of CHD7 inhibits cell proliferation and decreases gene expression of several CHD7 targets, including NRAS, in breast cancer cell lines. Thus, our results demonstrate the oncogenic potential of CHD7 and its association with poor prognostic parameters in human cancer.

RNA methylation, catalysed by a set of RNA methyltransferases (RNMTs), modulates RNA structures, properties, and biological functions. RNMTs are increasingly documented to be dysregulated in various human diseases, particularly developmental disorders and cancer. However, the genomic and transcriptomic alterations of RNMTs, as well as their functional roles in human cancer, are limited. In this study, we utilized an unbiased approach to examine copy number alterations and mutation rates of 58 RNMTs in more than 10,000 clinical samples across 32 human cancer types. We also investigated these alterations and RNMT expression level as they related to clinical features such as tumour subtype, grade, and survival in a large cohort of tumour samples, focusing on breast cancer. Loss-of-function analysis was performed to examine RNMT candidates with important roles in growth and viability of breast cancer cells. We identified a subset of RNMTs, notably TRMT12, NSUN2, TARBP1, and FTSJ3, that were amplified or mutated in a subset of human cancers. Several RNMTs were significantly associated with breast cancer aggressiveness and poor prognosis. Loss-of-function analysis indicated FTSJ3, a 2ʹ-O-Me methyltransferase, as a candidate RNMT with functional roles in promoting cancer growth and survival. A subset of RNMTs, like FTSJ3, represents promising novel targets for anticancer drug discovery. Our findings provide a framework for further study of the functional consequences of RNMT alterations in human cancer and for developing therapies that target cancer-promoting RNMTs in the future.

A wide range of the epigenetic effectors that regulate chromatin modification, gene expression, genomic stability, and DNA repair contain structurally conserved domains called plant homeodomain (PHD) fingers. Alternations of several PHD finger-containing proteins (PHFs) due to genomic amplification, mutations, deletions, and translocations have been linked directly to various types of cancer. However, little is known about the genomic landscape and the clinical significance of PHFs in breast cancer. Hence, we performed a large-scale genomic and transcriptomic analysis of 98 PHF genes in breast cancer using TCGA and METABRIC datasets and correlated the recurrent alterations with clinicopathological features and survival of patients. Different subtypes of breast cancer had different patterns of copy number and expression for each PHF. We identified a subset of PHF genes that was recurrently altered with high prevalence, including PYGO2 (pygopus family PHD finger 2), ZMYND8 (zinc finger, MYND-type containing 8), ASXL1 (additional sex combs like 1) and CHD3 (chromodomain helicase DNA binding protein 3). Copy number increase and overexpression of ZMYND8 were more prevalent in Luminal B subtypes and were significantly associated with shorter survival of breast cancer patients. ZMYND8 was also involved in a positive feedback circuit of the estrogen receptor (ER) pathway, and the expression of ZMYND8 was repressed by the bromodomain and extra terminal (BET) inhibitor in breast cancer. Our findings suggest a promising avenue for future research—to focus on a subset of PHFs to better understand the molecular mechanisms and to identify therapeutic targets in breast cancer.

Human methyltransferase-like (METTL) proteins transfer methyl groups to nucleic acids, proteins, lipids, and other small molecules, subsequently playing important roles in various cellular processes. In this study, we performed integrated genomic, transcriptomic, proteomic, and clinicopathological analyses of 34 METTLs in a large cohort of primary tumor and cell line data. We identified a subset of METTL genes, notably METTL1, METTL7B, and NTMT1, with high frequencies of genomic amplification and/or up-regulation at both the mRNA and protein levels in a spectrum of human cancers. Higher METTL1 expression was associated with high-grade tumors and poor disease prognosis. Loss-of-function analysis in tumor cell lines indicated the biological importance of METTL1, an m7G methyltransferase, in cancer cell growth and survival. Furthermore, functional annotation and pathway analysis of METTL1-associated proteins revealed that, in addition to the METTL1 cofactor WDR4, RNA regulators and DNA packaging complexes may be functionally interconnected with METTL1 in human cancer. Finally, we generated a crystal structure model of the METTL1–WDR4 heterodimeric complex that might aid in understanding the key functional residues. Our results provide new information for further functional study of some METTL alterations in human cancer and might lead to the development of small inhibitors that target cancer-promoting METTLs.

SARS-CoV-2, an emerging coronavirus, has spread rapidly around the world, resulting in over ten million cases and more than half a million deaths as of July 1, 2020. Effective treatments and vaccines for SARS-CoV-2 infection do not currently exist. Previous studies demonstrated that nonstructural protein 16 (nsp16) of coronavirus is an S-adenosyl methionine (SAM)-dependent 2'-O-methyltransferase (2'-O-MTase) that has an important role in viral replication and prevents recognition by the host innate immune system. In the present study, we employed structural analysis, virtual screening, and molecular simulation approaches to identify clinically investigated and approved drugs which can act as promising inhibitors against nsp16 2 0-O-MTase of SARS-CoV-2. Comparative analysis of primary amino acid sequences and crystal structures of seven human CoVs defined the key residues for nsp16 2-O'-MTase functions. Virtual screening and docking analysis ranked the potential inhibitors of nsp16 from more than 4,500 clinically investigated and approved drugs. Furthermore, molecular dynamics simulations were carried out on eight top candidates, including Hesperidin, Rimegepant, Gs-9667, and Sonedenoson, to calculate various structural parameters and understand the dynamic behavior of the drug-protein complexes. Our studies provided the foundation to further test and repurpose these candidate drugs experimentally and/or clinically for COVID-19 treatment.