Mitochondrial DNA mosaicism in normal human somatic cells.

Somatic cells accumulate genomic alterations with age; however, our understanding of mitochondrial DNA (mtDNA) mosaicism remains limited. Here we investigated the genomes of 2,096 clones derived from three cell types across 31 donors, identifying 6,451 mtDNA variants with heteroplasmy levels of ≳0.3%. While the majority of these variants were unique to individual clones, suggesting stochastic acquisition with age, 409 variants (6%) were shared across multiple embryonic lineages, indicating their origin from heteroplasmy in fertilized eggs. The mutational spectrum exhibited replication-strand bias, implicating mtDNA replication as a major mutational process. We evaluated the mtDNA mutation rate (5.0 × 10-8 per base pair) and a turnover frequency of 10-20 per year, which are fundamental components shaping the landscape of mtDNA mosaicism over a lifetime. The expansion of mtDNA-truncating mutations toward homoplasmy was substantially suppressed. Our findings provide comprehensive insights into the origins, dynamics and functional consequences of mtDNA mosaicism in human somatic cells.

CORRIGENDUM: Correction of the ACKNOWLEDGEMENTS (Fund/Grant Support). Epigenetic modulation inhibits epithelial-mesenchymal transition-driven fibrogenesis and enhances characteristics of chemically-derived hepatic progenitors.

[This corrects the article on p. 274 in vol. 106, PMID: 38725803.].

Epigenetic modulation inhibits epithelial-mesenchymal transition-driven fibrogenesis and enhances characteristics of chemically-derived hepatic progenitors.

Purpose: One of the novel cell sources of cell-based liver regenerative medicine is human chemically-derived hepatic progenitors (hCdHs). We previously established this cell by direct hepatocyte reprogramming with a combination of small molecules (hepatocyte growth factor, A83-01, CHIR99021). However, there have been several issues concerning the cell's stability and maintenance, namely the occurrences of epithelial-mesenchymal transition (EMT) that develop fibrotic phenotypes, resulting in the loss of hepatic progenitor characteristics. These hepatic progenitor attributes are thought to be regulated by SOX9, a transcription factor essential for hepatic progenitor cells and cholangiocytes. Methods: To suppress the fibrotic phenotype and improve our long-term hCdHs culture technology, we utilized the epigenetic modulating drugs DNA methyltransferase inhibitor (5-azacytidine) and histone deacetylase inhibitor (sodium butyrate) that have been reported to suppress and revert hepatic fibrosis. To confirm the essential role of SOX9 to our cell, we used clustered regularly interspaced short palindromic repeats-interference (CRISPRi) to repress the SOX9 expression. Results: The treatment of only 5-azacytidine significantly reduces the fibrosis/mesenchymal marker and EMT-related transcription factor expression level in the early passages. Interestingly, this treatment also increased the hepatic progenitor markers expression, even during the reprogramming phase. Then, we confirmed the essential role of SOX9 by repressing the SOX9 expression with CRISPRi which resulted in the downregulation of several essential hepatic progenitor cell markers. Conclusion: These results highlight the capacity of 5-azacytidine to inhibit EMT-driven hepatic fibrosis and the significance of SOX9 on hepatic progenitor cell stemness properties.

Quantitative and qualitative mutational impact of ionizing radiation on normal cells.

The comprehensive genomic impact of ionizing radiation (IR), a carcinogen, on healthy somatic cells remains unclear. Using large-scale whole-genome sequencing (WGS) of clones expanded from irradiated murine and human single cells, we revealed that IR induces a characteristic spectrum of short insertions or deletions (indels) and structural variations (SVs), including balanced inversions, translocations, composite SVs (deletion-insertion, deletion-inversion, and deletion-translocation composites), and complex genomic rearrangements (CGRs), including chromoplexy, chromothripsis, and SV by breakage-fusion-bridge cycles. Our findings suggest that 1 Gy IR exposure causes an average of 2.33 mutational events per Gb genome, comprising 2.15 indels, 0.17 SVs, and 0.01 CGRs, despite a high level of inter-cellular stochasticity. The mutational burden was dependent on total irradiation dose, regardless of dose rate or cell type. The findings were further validated in IR-induced secondary cancers and single cells without clonalization. Overall, our study highlights a comprehensive and clear picture of IR effects on normal mammalian genomes.

p57Kip2 imposes the reserve stem cell state of gastric chief cells.

Adult stem cells constantly react to local changes to ensure tissue homeostasis. In the main body of the stomach, chief cells produce digestive enzymes; however, upon injury, they undergo rapid proliferation for prompt tissue regeneration. Here, we identified p57Kip2 (p57) as a molecular switch for the reserve stem cell state of chief cells in mice. During homeostasis, p57 is constantly expressed in chief cells but rapidly diminishes after injury, followed by robust proliferation. Both single-cell RNA sequencing and dox-induced lineage tracing confirmed the sequential loss of p57 and activation of proliferation within the chief cell lineage. In corpus organoids, p57 overexpression induced a long-term reserve stem cell state, accompanied by altered niche requirements and a mature chief cell/secretory phenotype. Following the constitutive expression of p57 in vivo, chief cells showed an impaired injury response. Thus, p57 is a gatekeeper that imposes the reserve stem cell state of chief cells in homeostasis.

Patient-derived organoids as a preclinical platform for precision medicine in colorectal cancer.

Patient-derived organoids are being considered as models that can help guide personalized therapy through in vitro anticancer drug response evaluation. However, attempts to quantify in vitro drug responses in organoids and compare them with responses in matched patients remain inadequate. In this study, we investigated whether drug responses of organoids correlate with clinical responses of matched patients and disease progression of patients. Organoids were established from 54 patients with colorectal cancer who (except for one patient) did not receive any form of therapy before, and tumor organoids were assessed through whole-exome sequencing. For comparisons of in vitro drug responses in matched patients, we developed an 'organoid score' based on the variable anticancer treatment responses observed in organoids. Very interestingly, a higher organoid score was significantly correlated with a lower tumor regression rate after the standard-of-care treatment in matched patients. Additionally, we confirmed that patients with a higher organoid score (≥ 2.5) had poorer progression-free survival compared with those with a lower organoid score (< 2.5). Furthermore, to assess potential drug repurposing using an FDA-approved drug library, ten tumor organoids derived from patients with disease progression were applied to a simulation platform. Taken together, organoids and organoid scores can facilitate the prediction of anticancer therapy efficacy, and they can be used as a simulation model to determine the next therapeutic options through drug screening. Organoids will be an attractive platform to enable the implementation of personalized therapy for colorectal cancer patients.

Clonal dynamics in early human embryogenesis inferred from somatic mutation.

Cellular dynamics and fate decision in early human embryogenesis remain largely unknown owing to the challenges of performing studies in human embryos1. Here, we explored whole-genomes of 334 single-cell colonies and targeted deep sequences of 379 bulk tissues obtained from various anatomical locations of seven recently deceased adult human donors. Using somatic mutations as an intrinsic barcode, we reconstructed early cellular phylogenies that demonstrate (1) an endogenous mutational rate that is higher in the first cell division but decreases to approximately one per cell per cell division later in life; (2) universal unequal contribution of early cells to embryo proper, resulting from early cellular bottlenecks that stochastically set aside epiblast cells within the embryo; (3) examples of varying degrees of early clonal imbalances between tissues on the left and right sides of the body, different germ layers and specific anatomical parts and organs; (4) emergence of a few ancestral cells that will substantially contribute to adult cell pools in blood and liver; and (5) presence of mitochondrial DNA heteroplasmy in the fertilized egg. Our approach also provides insights into the age-related mutational processes and loss of sex chromosomes in normal somatic cells. In sum, this study provides a foundation for future studies to complete cellular phylogenies in human embryogenesis.

Low-Level Brain Somatic Mutations Are Implicated in Schizophrenia.

BACKGROUND: Somatic mutations arising from the brain have recently emerged as significant contributors to neurodevelopmental disorders, including childhood intractable epilepsy and cortical malformations. However, whether brain somatic mutations are implicated in schizophrenia (SCZ) is not well established. METHODS: We performed deep whole exome sequencing (average read depth > 550×) of matched dorsolateral prefrontal cortex and peripheral tissues from 27 patients with SCZ and 31 age-matched control individuals, followed by comprehensive and strict analysis of somatic mutations, including mutagenesis signature, substitution patterns, and involved pathways. In particular, we explored the impact of deleterious mutations in GRIN2B through primary neural culture. RESULTS: We identified an average of 4.9 and 5.6 somatic mutations per exome per brain in patients with SCZ and control individuals, respectively. These mutations presented with average variant allele frequencies of 8.0% in patients with SCZ and 7.6% in control individuals. Although mutational profiles, such as the number and type of mutations, showed no significant difference between patients with SCZ and control individuals, somatic mutations in SCZ brains were significantly enriched for SCZ-related pathways, including dopamine receptor, glutamate receptor, and long-term potentiation pathways. Furthermore, we showed that brain somatic mutations in GRIN2B (encoding glutamate ionotropic NMDA receptor subunit 2B), which were found in two patients with SCZ, disrupted the location of GRIN2B across the surface of dendrites among primary cultured neurons. CONCLUSIONS: Taken together, this study shows that brain somatic mutations are associated with the pathogenesis of SCZ.

The genome-wide landscape of C:G > T:A polymorphism at the CpG contexts in the human population.

BACKGROUND: The C:G > T:A substitution at the CpG dinucleotide contexts is the most frequent substitution type in genome evolution. The mutational process is obviously ongoing in the human germline; however, its impact on common and rare genomic polymorphisms has not been comprehensively investigated yet. Here we observed the landscape and dynamics of C:G > T:A substitutions from population-scale human genome sequencing datasets including ~ 4300 whole-genomes from the 1000 Genomes and the pan-cancer analysis of whole genomes (PCAWG) Project and ~ 60,000 whole-exomes from the Exome Aggregation Consortium (ExAC) database. RESULTS: Of the 28,084,558 CpG sites in the human reference genome, 26.0% show C:G > T:A substitution in the dataset. Remarkably, CpGs in CpG islands (CGIs) have a much lower frequency of such mutations (5.6%). Interestingly, the mutation frequency of CGIs is not uniform with a significantly higher C:G > T:A substitution rate for intragenic CGIs compared to other types. For non-CGI CpGs, the mutation rate was positively correlated with the distance from the nearest CGI up to 2 kb. Finally, we found the impact of negative selection for coding CpG mutations resulting in amino acid change. CONCLUSIONS: This study provides the first unbiased rate of C:G > T:A substitution at the CpG dinucleotide contexts, using population-scale human genome sequencing data. Our findings provide insights into the dynamics of the mutation acquisition in the human genome.

Tracing Oncogene Rearrangements in the Mutational History of Lung Adenocarcinoma.

Mutational processes giving rise to lung adenocarcinomas (LADCs) in non-smokers remain elusive. We analyzed 138 LADC whole genomes, including 83 cases with minimal contribution of smoking-associated mutational signature. Genomic rearrangements were not correlated with smoking-associated mutations and frequently served as driver events of smoking-signature-low LADCs. Complex genomic rearrangements, including chromothripsis and chromoplexy, generated 74% of known fusion oncogenes, including EML4-ALK, CD74-ROS1, and KIF5B-RET. Unlike other collateral rearrangements, these fusion-oncogene-associated rearrangements were frequently copy-number-balanced, representing a genomic signature of early oncogenesis. Analysis of mutation timing revealed that fusions and point mutations of canonical oncogenes were often acquired in the early decades of life. During a long latency, cancer-related genes were disrupted or amplified by complex rearrangements. The genomic landscape was different between subgroups-EGFR-mutant LADCs had frequent whole-genome duplications with p53 mutations, whereas fusion-oncogene-driven LADCs had frequent SETD2 mutations. Our study highlights LADC oncogenesis driven by endogenous mutational processes.

Identification of a quadruple mutation that confers tenofovir resistance in chronic hepatitis B patients.

BACKGROUND & AIMS: Tenofovir disoproxil fumarate (TDF) is one the most potent nucleot(s)ide analogues for treating chronic hepatitis B virus (HBV) infection. Phenotypic resistance caused by genotypic resistance to TDF has not been reported. This study aimed to characterize HBV mutations that confer tenofovir resistance. METHODS: Two patients with viral breakthrough during treatment with TDF-containing regimens were prospectively enrolled. The gene encoding HBV reverse transcriptase was sequenced. Eleven HBV clones harboring a series of mutations in the reverse transcriptase gene were constructed by site-directed mutagenesis. Drug susceptibility of each clone was determined by Southern blot analysis and real-time PCR. The relative frequency of mutants was evaluated by ultra-deep sequencing and clonal analysis. RESULTS: Five mutations (rtS106C [C], rtH126Y [Y], rtD134E [E], rtM204I/V, and rtL269I [I]) were commonly found in viral isolates from 2 patients. The novel mutations C, Y, and E were associated with drug resistance. In assays for drug susceptibility, the IC50 value for wild-type HBV was 3.8 ±â€¯0.6 µM, whereas the IC50 values for CYE and CYEI mutants were 14.1 ±â€¯1.8 and 58.1 ±â€¯0.9 µM, respectively. The IC90 value for wild-type HBV was 30 ±â€¯0.5 µM, whereas the IC90 values for CYE and CYEI mutants were 185 ±â€¯0.5 and 790 ±â€¯0.2 µM, respectively. Both tenofovir-resistant mutants and wild-type HBV had similar susceptibility to the capsid assembly modulator NVR 3-778 (IC50 <0.4 µM vs. IC50 = 0.4 µM, respectively). CONCLUSIONS: Our study reveals that the quadruple (CYEI) mutation increases the amount of tenofovir required to inhibit HBV by 15.3-fold in IC50 and 26.3-fold in IC90. These results demonstrate that tenofovir-resistant HBV mutants can emerge, although the genetic barrier is high. LAY SUMMARY: Tenofovir is the most potent nucleotide analogue for the treatment of chronic hepatitis B virus infection and there has been no hepatitis B virus mutation that confers >10-fold resistance to tenofovir up to 8 years. Herein, we identified, for the first time, a quadruple mutation that conferred 15.3-fold (IC50) and 26.3-fold (IC90) resistance to tenofovir in 2 patients who experienced viral breakthrough during tenofovir treatment.

Clonal History and Genetic Predictors of Transformation Into Small-Cell Carcinomas From Lung Adenocarcinomas.

Purpose Histologic transformation of EGFR mutant lung adenocarcinoma (LADC) into small-cell lung cancer (SCLC) has been described as one of the major resistant mechanisms for epidermal growth factor receptor (EGFR) tyrosine kinase inhibitors (TKIs). However, the molecular pathogenesis is still unclear. Methods We investigated 21 patients with advanced EGFR-mutant LADCs that were transformed into EGFR TKI-resistant SCLCs. Among them, whole genome sequencing was applied for nine tumors acquired at various time points from four patients to reconstruct their clonal evolutionary history and to detect genetic predictors for small-cell transformation. The findings were validated by immunohistochemistry in 210 lung cancer tissues. Results We identified that EGFR TKI-resistant LADCs and SCLCs share a common clonal origin and undergo branched evolutionary trajectories. The clonal divergence of SCLC ancestors from the LADC cells occurred before the first EGFR TKI treatments, and the complete inactivation of both RB1 and TP53 were observed from the early LADC stages in sequenced tumors. We extended the findings by immunohistochemistry in the early-stage LADC tissues of 75 patients treated with EGFR TKIs; inactivation of both Rb and p53 was strikingly more frequent in the small-cell-transformed group than in the nontransformed group (82% v 3%; odds ratio, 131; 95% CI, 19.9 to 859). Among patients registered in a predefined cohort (n = 65), an EGFR mutant LADC that harbored completely inactivated Rb and p53 had a 43× greater risk of small-cell transformation (relative risk, 42.8; 95% CI, 5.88 to 311). Branch-specific mutational signature analysis revealed that apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like (APOBEC)-induced hypermutation was frequent in the branches toward small-cell transformation. Conclusion EGFR TKI-resistant SCLCs are branched out early from the LADC clones that harbor completely inactivated RB1 and TP53. The evaluation of RB1 and TP53 status in EGFR TKI-treated LADCs is informative in predicting small-cell transformation.