RNA Oncological Therapeutics: Intracellular Hairpin RNA Assembly Enables MicroRNA-Triggered Anticancer Functionality.

RNA therapeutics are of global interest because of their versatility in targeting a variety of intracellular and extracellular biomolecules. In that context, long double-stranded RNA (dsRNA) has been studied as an antitumor agent that activates the immune response. However, its performance is constrained by poor cancer selectivity and cell-penetration ability. Here, we designed and synthesized an oncolytic RNA hairpin pair (oHP) that was selectively cytotoxic toward cancer cells expressing abundant oncogenic microRNA-21 (miR-21). Although the structure of each hairpin RNA was thermodynamically metastable, catalytic miR-21 input triggered it to open to generate a long nicked dsRNA. We demonstrated that oHP functioned as a cytotoxic amplifier of information in the presence of miR-21 in various cancer cells and tumor-bearing mice. This work represents the first example of the use of short RNA molecules as build-up-type anticancer agents that are triggered by an oncogenic miRNA.

Aberrant proximal tubule DNA methylation underlies phenotypic changes related to kidney dysfunction in patients with diabetes.

Epigenetic mechanisms are considered to contribute to diabetic nephropathy by maintaining memory of poor glycemic control during the early stages of diabetes. However, DNA methylation changes in the human kidney are poorly characterized, because of the lack of cell type-specific analysis. We examined DNA methylation in proximal tubules purified from diabetic nephropathy patients and identified differentially methylated CpG sites, given the critical role of proximal tubules in the kidney injury. Hypermethylation was observed at CpG sites annotated to genes responsible for proximal tubule functions, including gluconeogenesis, nicotinamide adenine dinucleotide synthesis, transporters of glucose, water, phosphate, and drugs, in diabetic kidneys, while genes involved in oxidative stress and the cytoskeleton exhibited demethylation. Methylation levels of CpG sites annotated to ACTN1, BCAR1, MYH9, UBE4B, AFMID, TRAF2, TXNIP, FOXO3, and HNF4A were correlated with the estimated glomerular filtration rate, while methylation of the CpG site in RUNX1 was associated with interstitial fibrosis and tubular atrophy. Hypermethylation of G6PC and HNF4A was accompanied by decreased expression in diabetic kidneys. Proximal tubule-specific hypomethylation of metabolic genes related to HNF4A observed in control kidneys was compromised in diabetic kidneys, suggesting a role for aberrant DNA methylation in the dedifferentiation process. Multiple genes with aberrant DNA methylation in diabetes overlapped genes with altered expressions in maladaptive proximal tubule cells, including transcription factors PPARA and RREB1. In conclusion, DNA methylation derangement in the proximal tubules of patients with diabetes may drive phenotypic changes, characterized by inflammatory and fibrotic features, along with impaired function in metabolism and transport.

Region-specific DNA hydroxymethylation along the malignant progression of IDH-mutant gliomas.

The majority of low-grade isocitrate dehydrogenase-mutant (IDHmt) gliomas undergo malignant progression (MP), but their underlying mechanism remains unclear. IDHmt gliomas exhibit global DNA methylation, and our previous report suggested that MP could be partly attributed to passive demethylation caused by accelerated cell cycles. However, during MP, there is also active demethylation mediated by ten-eleven translocation, such as DNA hydroxymethylation. Hydroxymethylation is reported to potentially contribute to gene expression regulation, but its role in MP remains under investigation. Therefore, we conducted a comprehensive analysis of hydroxymethylation during MP of IDHmt astrocytoma. Five primary/malignantly progressed IDHmt astrocytoma pairs were analyzed with oxidative bisulfite and the Infinium EPIC methylation array, detecting 5-hydroxymethyl cytosine at over 850,000 locations for region-specific hydroxymethylation assessment. Notably, we observed significant sharing of hydroxymethylated genomic regions during MP across the samples. Hydroxymethylated CpGs were enriched in open sea and intergenic regions (p < 0.001), and genes undergoing hydroxymethylation were significantly associated with cancer-related signaling pathways. RNA sequencing data integration identified 91 genes with significant positive/negative hydroxymethylation-expression correlations. Functional analysis suggested that positively correlated genes are involved in cell-cycle promotion, while negatively correlated ones are associated with antineoplastic functions. Analyses of The Cancer Genome Atlas clinical data on glioma were in line with these findings. Motif-enrichment analysis suggested the potential involvement of the transcription factor KLF4 in hydroxymethylation-based gene regulation. Our findings shed light on the significance of region-specific DNA hydroxymethylation in glioma MP and suggest its potential role in cancer-related gene expression and IDHmt glioma malignancy.

Transcriptomic and epigenetic dissection of spinal ependymoma (SP-EPN) identifies clinically relevant subtypes enriched for tumors with and without NF2 mutation.

Ependymomas encompass multiple clinically relevant tumor types based on localization and molecular profiles. Tumors of the methylation class "spinal ependymoma" (SP-EPN) represent the most common intramedullary neoplasms in children and adults. However, their developmental origin is ill-defined, molecular data are scarce, and the potential heterogeneity within SP-EPN remains unexplored. The only known recurrent genetic events in SP-EPN are loss of chromosome 22q and NF2 mutations, but neither types and frequency of these alterations nor their clinical relevance have been described in a large, epigenetically defined series. Transcriptomic (n = 72), epigenetic (n = 225), genetic (n = 134), and clinical data (n = 112) were integrated for a detailed molecular overview on SP-EPN. Additionally, we mapped SP-EPN transcriptomes to developmental atlases of the developing and adult spinal cord to uncover potential developmental origins of these tumors. The integration of transcriptomic ependymoma data with single-cell atlases of the spinal cord revealed that SP-EPN display the highest similarities to mature adult ependymal cells. Unsupervised hierarchical clustering of transcriptomic data together with integrated analysis of methylation profiles identified two molecular SP-EPN subtypes. Subtype A tumors primarily carried previously known germline or sporadic NF2 mutations together with 22q loss (bi-allelic NF2 loss), resulting in decreased NF2 expression. Furthermore, they more often presented as multilocular disease and demonstrated a significantly reduced progression-free survival as compared to SP-EP subtype B. In contrast, subtype B predominantly contained samples without NF2 mutation detected in sequencing together with 22q loss (monoallelic NF2 loss). These tumors showed regular NF2 expression but more extensive global copy number alterations. Based on integrated molecular profiling of a large multi-center cohort, we identified two distinct SP-EPN subtypes with important implications for genetic counseling, patient surveillance, and drug development priorities.

Loss of viral genome with altered immune microenvironment during tumour progression of Epstein-Barr virus-associated gastric carcinoma.

Epstein-Barr virus (EBV) is one of the major drivers of gastric carcinogenesis. EBV infection is established before tumour initiation and is generally maintained throughout tumour development; however, the significance of EBV in tumour maintenance and progression remains to be elucidated. Here, we report eight cases of EBV-associated gastric carcinoma (EBVaGC) with intratumoural heterogenous expression of EBV-encoded small RNA (EBER), a highly expressed latent gene of EBV, and demonstrate clinicopathological characteristics of these rare cases. By performing detailed histological assessment of EBER-positive and -negative components of each case, detection of EBV genome in tumour cells by fluorescence in situ hybridisation, TP73 methylation analysis, whole exome sequencing, and targeted gene panel sequencing, we identified tumours in two patients to be collision tumours of different origins. In the other six patients, some genetic/epigenetic alterations were shared between EBER-positive and -negative components, suggesting that EBV was eliminated from tumour cells during progression. Interestingly, in both tumour types, programmed death ligand 1 and intratumoural infiltration of CD8+ T lymphocytes were lower in EBER-negative than in EBER-positive components, suggesting an immunogenic role of EBV. To the best of our knowledge, this study is the first to demonstrate the detailed histological features and genetic/epigenetic alterations in EBVaGC with heterogenous EBER expression; the loss of EBV may benefit tumour progression and immune evasion and might be clinically important for selecting treatment strategies for such cancers. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

TET1 upregulation drives cancer cell growth through aberrant enhancer hydroxymethylation of HMGA2 in hepatocellular carcinoma.

Ten-eleven translocation 1 (TET1) is an essential methylcytosine dioxygenase of the DNA demethylation pathway. Despite its dysregulation being known to occur in human cancer, the role of TET1 remains poorly understood. In this study, we report that TET1 promotes cell growth in human liver cancer. The transcriptome analysis of 68 clinical liver samples revealed a subgroup of TET1-upregulated hepatocellular carcinoma (HCC), demonstrating hepatoblast-like gene expression signatures. We performed comprehensive cytosine methylation and hydroxymethylation (5-hmC) profiling and found that 5-hmC was aberrantly deposited preferentially in active enhancers. TET1 knockdown in hepatoma cell lines decreased hmC deposition with cell growth suppression. HMGA2 was highly expressed in a TET1high subgroup of HCC, associated with the hyperhydroxymethylation of its intronic region, marked as histone H3K4-monomethylated, where the H3K27-acetylated active enhancer chromatin state induced interactions with its promoter. Collectively, our findings point to a novel type of epigenetic dysregulation, methylcytosine dioxygenase TET1, which promotes cell proliferation via the ectopic enhancer of its oncogenic targets, HMGA2, in hepatoblast-like HCC.

Base-resolution analysis of 5-hydroxymethylcytidine by selective oxidation and reverse transcription arrest.

While 5-hydroxymethylcytidine in RNA (hm5C) is associated with cellular development and differentiation, its distribution and biological function remain largely unexplored because suitable detection methods are lacking. Here, we report a base-resolution sequencing method for hm5C in RNA by applying peroxotungstate-mediated chemical conversion of hm5C to trihydroxylated thymine (thT). Reverse transcription by SuperScript III terminated at the thT site, probably because of its unnatural nucleobase structure producing truncated cDNA. Consequently, base-resolution analysis of the hm5C sites in RNA was achieved with both Sanger sequencing and Illumina sequencing analysis by comparing sequencing data before and after peroxotungstate treatment.

Spatial and temporal expansion of intrahepatic metastasis by molecularly-defined clonality in multiple liver cancers.

Multiple hepatocellular carcinoma (HCC) is divided into two categories: intrahepatic metastasis (IM), which is a true relapse of HCC, and multicentric origin (MO), which is a second primary tumor. Clinical diagnosis of multiple HCC is usually made based on tumor location and/or time to recurrence; however, it is often difficult to distinguish the two types of multiple HCC. Using 41 matched pairs of multiple HCC specimens, we confirmed the accuracy of clinical diagnoses using exome sequence data and investigated the importance of discriminating the type of multiple HCC. Genomic analysis revealed that 18 (43.9%) patients diagnosed as having genomic IM had common mutations in a pair of HCC tumors with the main tumor of these patients being more progressive compared to those with genomic MO. The accuracy of clinical diagnosis based on lobe (Definition 1) and segment (Definition 2) were 68.3% and 78.0%, respectively. Intriguingly, recurrence ≥2 years after initial surgery for 3 patients was IM. The survival of patients with clinical IM was significantly shorter than for those with clinical MO based on both Definition 1 (P = 0.045) and Definition 2 (P = 0.043). However, mean survival was not different between the patients with genomic IM and those with MO (P = 0.364). Taken together, genomic analysis elucidated that liver cancer may spread more extensively and more slowly than previously thought. In addition, distinguishing multiple HCC as IM or MC may have provided biological information but was not of clinical importance with respect to patient prognosis.

Integrated molecular profiling of juvenile myelomonocytic leukemia.

Juvenile myelomonocytic leukemia (JMML), a rare and aggressive myelodysplastic/myeloproliferative neoplasm that occurs in infants and during early childhood, is characterized by excessive myelomonocytic cell proliferation. More than 80% of patients harbor germ line and somatic mutations in RAS pathway genes (eg, PTPN11, NF1, NRAS, KRAS, and CBL), and previous studies have identified several biomarkers associated with poor prognosis. However, the molecular pathogenesis of 10% to 20% of patients and the relationships among these biomarkers have not been well defined. To address these issues, we performed an integrated molecular analysis of samples from 150 JMML patients. RNA-sequencing identified ALK/ROS1 tyrosine kinase fusions (DCTN1-ALK, RANBP2-ALK, and TBL1XR1-ROS1) in 3 of 16 patients (18%) who lacked canonical RAS pathway mutations. Crizotinib, an ALK/ROS1 inhibitor, markedly suppressed ALK/ROS1 fusion-positive JMML cell proliferation in vitro. Therefore, we administered crizotinib to a chemotherapy-resistant patient with the RANBP2-ALK fusion who subsequently achieved complete molecular remission. In addition, crizotinib also suppressed proliferation of JMML cells with canonical RAS pathway mutations. Genome-wide methylation analysis identified a hypermethylation profile resembling that of acute myeloid leukemia (AML), which correlated significantly with genetic markers with poor outcomes such as PTPN11/NF1 gene mutations, 2 or more genetic mutations, an AML-type expression profile, and LIN28B expression. In summary, we identified recurrent activated ALK/ROS1 fusions in JMML patients without canonical RAS pathway gene mutations and revealed the relationships among biomarkers for JMML. Crizotinib is a promising candidate drug for the treatment of JMML, particularly in patients with ALK/ROS1 fusions.

Correction: Integrated Multiregional Analysis Proposing a New Model of Colorectal Cancer Evolution.

[This corrects the article DOI: 10.1371/journal.pgen.1005778.].

Differential regulation of CpG island methylation within divergent and unidirectional promoters in colorectal cancer.

The silencing of tumor suppressor genes by promoter CpG island (CGI) methylation is an important cause of oncogenesis. Silencing of MLH1 and BRCA1, two examples of oncogenic events, results from promoter CGI methylation. Interestingly, both MLH1 and BRCA1 have a divergent promoter, from which another gene on the opposite strand is also transcribed. Although studies have shown that divergent transcription is an important factor in transcriptional regulation, little is known about its implication in aberrant promoter methylation in cancer. In this study, we analyzed the methylation status of CGI in divergent promoters using a recently enriched transcriptome database. We measured the extent of CGI methylation in 119 colorectal cancer (CRC) clinical samples (65 microsatellite instability high [MSI-H] CRC with CGI methylator phenotype, 28 MSI-H CRC without CGI methylator phenotype and 26 microsatellite stable CRC) and 21 normal colorectal tissues using Infinium MethylationEPIC BeadChip. We found that CGI within divergent promoters are less frequently methylated than CGI within unidirectional promoters in normal cells. In the genome of CRC cells, CGI within unidirectional promoters are more vulnerable to aberrant methylation than CGI within divergent promoters. In addition, we identified three DNA sequence motifs that correlate with methylated CGI. We also showed that methylated CGI are associated with genes whose expression is low in normal cells. Thus, we here provide fundamental observations regarding the methylation of divergent promoters that are essential for the understanding of carcinogenesis and development of cancer prevention strategies.

Integrated Multiregional Analysis Proposing a New Model of Colorectal Cancer Evolution.

Understanding intratumor heterogeneity is clinically important because it could cause therapeutic failure by fostering evolutionary adaptation. To this end, we profiled the genome and epigenome in multiple regions within each of nine colorectal tumors. Extensive intertumor heterogeneity is observed, from which we inferred the evolutionary history of the tumors. First, clonally shared alterations appeared, in which C>T transitions at CpG site and CpG island hypermethylation were relatively enriched. Correlation between mutation counts and patients' ages suggests that the early-acquired alterations resulted from aging. In the late phase, a parental clone was branched into numerous subclones. Known driver alterations were observed frequently in the early-acquired alterations, but rarely in the late-acquired alterations. Consistently, our computational simulation of the branching evolution suggests that extensive intratumor heterogeneity could be generated by neutral evolution. Collectively, we propose a new model of colorectal cancer evolution, which is useful for understanding and confronting this heterogeneous disease.

Aberrant DNA methylation of pregnane X receptor underlies metabolic gene alterations in the diabetic kidney.

Epigenetic abnormalities have been suggested to mediate metabolic memory observed in diabetic complications. We have shown that epigenetic alterations may induce persistent phenotypic changes in the proximal tubules of the diabetic kidneys. In this study, we show that pregnane X receptor (PXR), a xenobiotic nuclear receptor, is epigenetically altered and upregulated and may have a possible function in the diabetic kidney. PXR has been shown to play a critical role in metabolic changes in obesity and diabetes; however, its distribution and function in the kidney are unknown. In the normal kidney, Pxr was selectively expressed in the proximal tubular cells with demethylation in the promoter DNA. In db/db mice, significant increases in Pxr mRNA, further demethylation of DNA, and stimulatory histone marks in the promoter were observed. Epigenetic changes are likely to play a causative role in PXR induction, since a DNA methyltransferase inhibitor increased PXR mRNA in cultured human proximal tubular cells. Administration of a PXR agonist increased mRNA levels of solute carrier organic anion transporter family member 2B1 ( Slco2b1), a xenobiotic transporter; response gene to complement 32 ( Rgc32), a molecule known to exert fibrotic effects in the kidney; and phosphoenolpyruvate carboxykinase 1 ( Pck1), a gluconeogenic enzyme in the kidney. The expressions of these genes were inhibited by PXR small interfering RNA in cultured proximal tubular cells. Increased mRNA levels of Slco2b1, Rgc32, and Pck1 were also observed in the kidney of db/db mice. These data indicate that PXR is upregulated in the diabetic kidney with aberrant epigenetic modifications and may modulate the course of diabetic kidney disease through the activation of these genes.

Transduced caudal-type homeobox (CDX) 2/CDX1 can induce growth inhibition on CDX-deficient gastric cancer by rapid intestinal differentiation.

Intestinal metaplasia induced by ectopic expression of caudal-type homeobox (CDX)2 and/or CDX1 (CDX) is frequently observed around gastric cancer (GC). Abnormal expression of CDX is also observed in GC and suggests that inappropriate gastrointestinal differentiation plays essential roles in gastric tumorigenesis, but their roles on tumorigenesis remain unelucidated. Publicly available databases show that GC patients with higher CDX expression have significantly better clinical outcomes. We introduced CDX2 and CDX1 genes separately into GC-originated MKN7 and TMK1 cells deficient in CDX. Marked suppression of cell growth and dramatic morphological change into spindle-shaped flat form were observed along with induction of intestinal marker genes. G0-G1 growth arrest was accompanied by changed expression of cell cycle-related genes but not with apoptosis or senescence. Microarray analyses additionally showed decreased expression of gastric marker genes and increased expression of stemness-associated genes. Hierarchical clustering of 111 GC tissues and 21 non-cancerous gastric tissues by selected 18 signature genes based on our transcriptome analyses clearly categorized the 132 tissues into non-cancer, "CDX signature"-positive GC, and "CDX signature"-negative GC. Gene set enrichment analysis indicated that "CDX signature"-positive GC has lower malignant features. Immunohistochemistry of 89 GC specimens showed that 50.6% were CDX2-deficient, 66.3% were CDX1-deficient, and 44.9% were concomitant CDX2/CDX1-deficient, suggesting that potentially targetable GC cases by induced intestinal differentiation are quite common. In conclusion, exogenous expression of CDX2/CDX1 can lead to efficient growth inhibition of CDX-deficient GC cells. It is based on rapidly induced intestinal differentiation, which may be a future therapeutic strategy.

Distinct molecular profile of diffuse cerebellar gliomas.

Recent studies have demonstrated that tumor-driving alterations are often different among gliomas that originated from different brain regions and have underscored the importance of analyzing molecular characteristics of gliomas stratified by brain region. Therefore, to elucidate molecular characteristics of diffuse cerebellar gliomas (DCGs), 27 adult, mostly glioblastoma cases were analyzed. Comprehensive analysis using whole-exome sequencing, RNA sequencing, and Infinium methylation array (n = 17) demonstrated their distinct molecular profile compared to gliomas in other brain regions. Frequent mutations in chromatin-modifier genes were identified including, noticeably, a truncating mutation in SETD2 (n = 4), which resulted in loss of H3K36 trimethylation and was mutually exclusive with H3F3A K27M mutation (n = 3), suggesting that epigenetic dysregulation may lead to DCG tumorigenesis. Alterations that cause loss of p53 function including TP53 mutation (n = 9), PPM1D mutation (n = 2), and a novel type of PPM1D fusion (n = 1), were also frequent. On the other hand, mutations and copy number changes commonly observed in cerebral gliomas were infrequent. DNA methylation profile analysis demonstrated that all DCGs except for those with H3F3A mutations were categorized in the "RTK I (PDGFRA)" group, and those DCGs had a gene expression signature that was highly associated with PDGFRA. Furthermore, compared with the data of 315 gliomas derived from different brain regions, promoter methylation of transcription factors genes associated with glial development showed a characteristic pattern presumably reflecting their tumor origin. Notably, SOX10, a key transcription factor associated with oligodendroglial differentiation and PDGFRA regulation, was up-regulated in both DCG and H3 K27M-mutant diffuse midline glioma, suggesting their developmental and biological commonality. In contrast, SOX10 was silenced by promoter methylation in most cerebral gliomas. These findings may suggest potential tailored targeted therapy for gliomas according to their brain region, in addition to providing molecular clues to identify the region-related cellular origin of DCGs.

Base-Resolution Analysis of 5-Hydroxymethylcytosine by One-Pot Bisulfite-Free Chemical Conversion with Peroxotungstate.

5-Hydroxymethylcytosine (hmC) is an essential intermediate in the active DNA demethylation pathway. Here we report a new base-resolution method for measuring hmC by combining peroxotungstate-mediated oxidation and sequencing analysis. We reveal that an oxidized product of hmC, trihydroxylated thymine (thT), tolerated the incorporation of dATP as a substrate in the process of DNA polymerase elongation. By comparing the results of Sanger sequencing before and after the oxidation, we observed that hmC sites on single-stranded DNAs could be discriminated from unmethylated cytosines. We found that a thermal cycle condition during peroxotungstate treatment enhanced the oxidation reaction of hmC in double-stranded DNA. Furthermore, Illumina sequencing analysis of hmC-containing synthetic genome fragments enabled us to identify simultaneously the positions of hmC in base resolution. This bisulfite-free simple hmC detection technique could facilitate the acquisition of epigenomic information.

Prognostic significance of CpG island methylator phenotype in surgically resected small cell lung carcinoma.

Methylation is closely involved in the development of various carcinomas. However, few datasets are available for small cell lung cancer (SCLC) due to the scarcity of fresh tumor samples. The aim of the present study is to clarify relationships between clinicopathological features and results of the comprehensive genome-wide methylation profile of SCLC. We investigated the genome-wide DNA methylation status of 28 tumor and 13 normal lung tissues, and gene expression profiling of 25 SCLC tissues. Following unsupervised hierarchical clustering and non-negative matrix factorization, gene ontology analysis was performed. Clustering of SCLC led to the important identification of a CpG island methylator phenotype (CIMP) of the tumor, with a significantly poorer prognosis (P = 0.002). Multivariate analyses revealed that postoperative chemotherapy and non-CIMP were significantly good prognostic factors. Ontology analyses suggested that the extrinsic apoptosis pathway was suppressed, including TNFRSF1A, TNFRSF10A and TRADD in CIMP tumors. Here we revealed that CIMP was an important prognostic factor for resected SCLC. Delineation of this phenotype may also be useful for the development of novel apoptosis-related chemotherapeutic agents for treatment of the aggressive tumor.

Hes1 promotes blast crisis in chronic myelogenous leukemia through MMP-9 upregulation in leukemic cells.

High levels of HES1 expression are frequently found in BCR-ABL(+) chronic myelogenous leukemia in blast crisis (CML-BC). In mouse bone marrow transplantation (BMT) models, co-expression of BCR-ABL and Hes1 induces CML-BC-like disease; however, the underlying mechanism remained elusive. Here, based on gene expression analysis, we show that MMP-9 is upregulated by Hes1 in common myeloid progenitors (CMPs). Analysis of promoter activity demonstrated that Hes1 upregulated MMP-9 by activating NF-κB. Analysis of 20 samples from CML-BC patients showed that MMP-9 was highly expressed in three, with two exhibiting high levels of HES1 expression. Interestingly, MMP-9 deficiency impaired the cobblestone area-forming ability of CMPs expressing BCR-ABL and Hes1 that were in conjunction with a stromal cell layer. In addition, CMPs expressing BCR-ABL and Hes1 secreted MMP-9, promoting the release of soluble Kit-ligand (sKitL) from stromal cells, thereby enhancing proliferation of the leukemic cells. In accordance, mice transplanted with CMPs expressing BCR-ABL and Hes1 exhibited high levels of sKitL as well as MMP-9 in the serum. Importantly, MMP-9 deficiency impaired the development of CML-BC-like disease induced by BCR-ABL and Hes1 in mouse BMT models. The present results suggest that Hes1 promotes the development of CML-BC, partly through MMP-9 upregulation in leukemic cells.

Whole-exome sequencing of human pancreatic cancers and characterization of genomic instability caused by MLH1 haploinsufficiency and complete deficiency.

Whole-exome sequencing (Exome-seq) has been successfully applied in several recent studies. We here sequenced the exomes of 15 pancreatic tumor cell lines and their matched normal samples. We captured 162,073 exons of 16,954 genes and sequenced the targeted regions to a mean coverage of 56-fold. This study identified a total of 1517 somatic mutations and validated 934 mutations by transcriptome sequencing. We detected recurrent mutations in 56 genes. Among them, 41 have not been described. The mutation rates varied widely among cell lines. The diversity of the mutation rates was significantly correlated with the distinct MLH1 copy-number status. Exome-seq revealed intensive genomic instability in a cell line with MLH1 homozygous deletion, indicated by a dramatically elevated rate of somatic substitutions, small insertions/deletions (indels), as well as indels in microsatellites. Notably, we found that MLH1 expression was decreased by nearly half in cell lines with an allelic loss of MLH1. While these cell lines were negative in conventional microsatellite instability assay, they showed a 10.5-fold increase in the rate of somatic indels, e.g., truncating indels in TP53 and TGFBR2, indicating MLH1 haploinsufficiency in the correction of DNA indel errors. We further analyzed the exomes of 15 renal cell carcinomas and confirmed MLH1 haploinsufficiency. We observed a much higher rate of indel mutations in the affected cases and identified recurrent truncating indels in several cancer genes such as VHL, PBRM1, and JARID1C. Together, our data suggest that MLH1 hemizygous deletion, through increasing the rate of indel mutations, could drive the development and progression of sporadic cancers.

BRCC3 mutations in myeloid neoplasms.

Next generation sequencing technologies have provided insights into the molecular heterogeneity of various myeloid neoplasms, revealing previously unknown somatic genetic events. In our cohort of 1444 cases analyzed by next generation sequencing, somatic mutations in the gene BRCA1-BRCA2-containing complex 3 (BRCC3) were identified in 28 cases (1.9%). BRCC3 is a member of the JAMM/MPN+ family of zinc metalloproteases capable of cleaving Lys-63 linked polyubiquitin chains, and is implicated in DNA repair. The mutations were located throughout its coding region. The average variant allelic frequency of BRCC3 mutations was 30.1%, and by a serial sample analysis at two different time points a BRCC3 mutation was already identified in the initial stage of a myelodysplastic syndrome. BRCC3 mutations commonly occurred in nonsense (n=12), frameshift (n=4), and splice site (n=5) configurations. Due to the marginal male dominance (odds ratio; 2.00, 0.84-4.73) of BRCC3 mutations, the majority of mutations (n=23; 82%) were hemizygous. Phenotypically, BRCC3 mutations were frequently observed in myelodysplastic syndromes and myelodysplastic/myeloproliferative neoplasms and associated with -Y abnormality (odds ratio; 3.70, 1.25-11.0). Clinically, BRCC3 mutations were also related to higher age (P=0.01), although prognosis was not affected. Knockdown of Brcc3 gene expression in murine bone marrow lineage negative, Sca1 positive, c-kit positive cells resulted in 2-fold more colony formation and modest differentiation defect. Thus, BRCC3 likely plays a role as tumor-associated gene in myelodysplastic syndromes and myelodysplastic/myeloproliferative neoplasms.