Negative Selection Allows DNA Mismatch Repair-Deficient Mouse Fibroblasts In Vitro to Tolerate High Levels of Somatic Mutations.

Substantial numbers of somatic mutations have been found to accumulate with age in different human tissues. Clonal cellular amplification of some of these mutations can cause cancer and other diseases. However, it is as yet unclear if and to what extent an increased burden of random mutations can affect cellular function without clonal amplification. We tested this in cell culture, which avoids the limitation that an increased mutation burden in vivo typically leads to cancer. We performed single-cell whole-genome sequencing of primary fibroblasts from DNA mismatch repair (MMR) deficient Msh2-/- mice and littermate control animals after long-term passaging. Apart from analyzing somatic mutation burden we analyzed clonality, mutational signatures, and hotspots in the genome, characterizing the complete landscape of somatic mutagenesis in normal and MMR-deficient mouse primary fibroblasts during passaging. While growth rate of Msh2-/- fibroblasts was not significantly different from the controls, the number of de novo single-nucleotide variants (SNVs) increased linearly up until at least 30,000 SNVs per cell, with the frequency of small insertions and deletions (INDELs) plateauing in the Msh2-/- fibroblasts to about 10,000 INDELS per cell. We provide evidence for negative selection and large-scale mutation-driven population changes, including significant clonal expansion of preexisting mutations and widespread cell-strain-specific hotspots. Overall, our results provide evidence that increased somatic mutation burden drives significant cell evolutionary changes in a dynamic cell culture system without significant effects on growth. Since similar selection processes against mutations preventing organ and tissue dysfunction during aging are difficult to envision, these results suggest that increased somatic mutation burden can play a causal role in aging and diseases other than cancer.

Analyzing somatic mutations by single-cell whole-genome sequencing.

Somatic mutations are the cause of cancer and have been implicated in other, noncancerous diseases and aging. While clonally expanded mutations can be studied by deep sequencing of bulk DNA, very few somatic mutations expand clonally, and most are unique to each cell. We describe a detailed protocol for single-cell whole-genome sequencing to discover and analyze somatic mutations in tissues and organs. The protocol comprises single-cell multiple displacement amplification (SCMDA), which ensures efficiency and high fidelity in amplification, and the SCcaller software tool to call single-nucleotide variations and small insertions and deletions from the sequencing data by filtering out amplification artifacts. With SCMDA and SCcaller at its core, this protocol describes a complete procedure for the comprehensive analysis of somatic mutations in a single cell, covering (1) single-cell or nucleus isolation, (2) single-cell or nucleus whole-genome amplification, (3) library preparation and sequencing, and (4) computational analyses, including alignment, variant calling, and mutation burden estimation. Methods are also provided for mutation annotation, hotspot discovery and signature analysis. The protocol takes 12-15 h from single-cell isolation to library preparation and 3-7 d of data processing. Compared with other single-cell amplification methods or single-molecular sequencing, it provides high genomic coverage, high accuracy in single-nucleotide variation and small insertions and deletion calling from the same single-cell genome, and fewer processing steps. SCMDA and SCcaller require basic experience in molecular biology and bioinformatics. The protocol can be utilized for studying mutagenesis and genome mosaicism in normal and diseased human and animal tissues under various conditions.

Expression of retrotransposons contributes to aging in Drosophila.

Retrotransposons are a class of transposable elements capable of self-replication and insertion into new genomic locations. Across species, the mobilization of retrotransposons in somatic cells has been suggested to contribute to the cell and tissue functional decline that occurs during aging. Retrotransposons are broadly expressed across cell types, and de novo insertions have been observed to correlate with tumorigenesis. However, the extent to which new retrotransposon insertions occur during normal aging and their effect on cellular and animal function remains understudied. Here, we use a single nucleus whole genome sequencing approach in Drosophila to directly test whether transposon insertions increase with age in somatic cells. Analyses of nuclei from thoraces and indirect flight muscles using a newly developed pipeline, Retrofind, revealed no significant increase in the number of transposon insertions with age. Despite this, reducing the expression of two different retrotransposons, 412 and Roo, extended lifespan, but did not alter indicators of health such as stress resistance. This suggests a key role for transposon expression and not insertion in regulating longevity. Transcriptomic analyses revealed similar changes to gene expression in 412 and Roo knockdown flies and highlighted changes to genes involved in proteolysis and immune function as potential contributors to the observed changes in longevity. Combined, our data show a clear link between retrotransposon expression and aging.

Single-cell analysis of somatic mutations in human bronchial epithelial cells in relation to aging and smoking.

Although lung cancer risk among smokers is dependent on smoking dose, it remains unknown if this increased risk reflects an increased rate of somatic mutation accumulation in normal lung cells. Here, we applied single-cell whole-genome sequencing of proximal bronchial basal cells from 33 participants aged between 11 and 86 years with smoking histories varying from never-smoking to 116 pack-years. We found an increase in the frequency of single-nucleotide variants and small insertions and deletions with chronological age in never-smokers, with mutation frequencies significantly elevated among smokers. When plotted against smoking pack-years, mutations followed the linear increase in cancer risk until about 23 pack-years, after which no further increase in mutation frequency was observed, pointing toward individual selection for mutation avoidance. Known lung cancer-defined mutation signatures tracked with both age and smoking. No significant enrichment for somatic mutations in lung cancer driver genes was observed.

Single-molecule, quantitative detection of low-abundance somatic mutations by high-throughput sequencing.

Postzygotic somatic mutations have been found associated with human disease, including diseases other than cancer. Most information on somatic mutations has come from studying clonally amplified mutant cells, based on a growth advantage or genetic drift. However, almost all somatic mutations are unique for each cell, and the quantitative analysis of these low-abundance mutations in normal tissues remains a major challenge in biology. Here, we introduce single-molecule mutation sequencing (SMM-seq), a novel approach for quantitative identification of point mutations in normal cells and tissues.

Single-cell analysis of somatic mutation burden in mammary epithelial cells of pathogenic BRCA1/2 mutation carriers.

Inherited germline mutations in the breast cancer gene 1 (BRCA1) or BRCA2 genes (herein BRCA1/2) greatly increase the risk of breast and ovarian cancer, presumably by elevating somatic mutational errors as a consequence of deficient DNA repair. However, this has never been directly demonstrated by a comprehensive analysis of the somatic mutational landscape of primary, noncancer, mammary epithelial cells of women diagnosed with pathogenic BRCA1/2 germline mutations. Here, we used an accurate, single-cell whole-genome sequencing approach to first show that telomerized primary mammary epithelial cells heterozygous for a highly penetrant BRCA1 variant displayed a robustly elevated mutation frequency as compared with their isogenic control cells. We then demonstrated a small but statistically significant increase in mutation frequency in mammary epithelial cells isolated from the breast of BRCA1/2 mutation carriers as compared with those obtained from age-matched controls with no genetically increased risk for breast cancer.

FOXO3a acts to suppress DNA double-strand break-induced mutations.

Genomic instability is one of the hallmarks of aging, and both DNA damage and mutations have been found to accumulate with age in different species. Certain gene families, such as sirtuins and the FoxO family of transcription factors, have been shown to play a role in lifespan extension. However, the mechanism(s) underlying the increased longevity associated with these genes remains largely unknown and may involve the regulation of responses to cellular stressors, such as DNA damage. Here, we report that FOXO3a reduces genomic instability in cultured mouse embryonic fibroblasts (MEFs) treated with agents that induce DNA double-strand breaks (DSBs), that is, clastogens. We show that DSB treatment of both primary human and mouse fibroblasts upregulates FOXO3a expression. FOXO3a ablation in MEFs harboring the mutational reporter gene lacZ resulted in an increase in genome rearrangements after bleomycin treatment; conversely, overexpression of human FOXO3a was found to suppress mutation accumulation in response to bleomycin. We also show that overexpression of FOXO3a in human primary fibroblasts decreases DSB-induced γH2AX foci. Knocking out FOXO3a in mES cells increased the frequency of homologous recombination and non-homologous end-joining events. These results provide the first direct evidence that FOXO3a plays a role in suppressing genome instability, possibly by suppressing genome rearrangements.

A direct comparison of interphase FISH versus low-coverage single cell sequencing to detect aneuploidy reveals respective strengths and weaknesses.

Aneuploidy has been reported to occur at remarkably high levels in normal somatic tissues using Fluorescence In Situ Hybridization (FISH). Recently, these reports were contradicted by single-cell low-coverage whole genome sequencing (scL-WGS) analyses, which showed aneuploidy frequencies at least an order of magnitude lower. To explain these seemingly contradictory findings, we used both techniques to analyze artificially generated mock aneuploid cells and cells with natural random aneuploidy. Our data indicate that while FISH tended to over-report aneuploidies, a modified 2-probe approach can accurately detect low levels of aneuploidy. Further, scL-WGS tends to underestimate aneuploidy levels, especially in a polyploid background.

Single-cell whole-genome sequencing reveals the functional landscape of somatic mutations in B lymphocytes across the human lifespan.

Accumulation of mutations in somatic cells has been implicated as a cause of aging since the 1950s. However, attempts to establish a causal relationship between somatic mutations and aging have been constrained by the lack of methods to directly identify mutational events in primary human tissues. Here we provide genome-wide mutation frequencies and spectra of human B lymphocytes from healthy individuals across the entire human lifespan using a highly accurate single-cell whole-genome sequencing method. We found that the number of somatic mutations increases from <500 per cell in newborns to >3,000 per cell in centenarians. We discovered mutational hotspot regions, some of which, as expected, were located at Ig genes associated with somatic hypermutation (SHM). B cell-specific mutation signatures associated with development, aging, or SHM were found. The SHM signature strongly correlated with the signature found in human B cell tumors, indicating that potential cancer-causing events are already present even in B cells of healthy individuals. We also identified multiple mutations in sequence features relevant to cellular function (i.e., transcribed genes and gene regulatory regions). Such mutations increased significantly during aging, but only at approximately one-half the rate of the genome average, indicating selection against mutations that impact B cell function. This full characterization of the landscape of somatic mutations in human B lymphocytes indicates that spontaneous somatic mutations accumulating with age can be deleterious and may contribute to both the increased risk for leukemia and the functional decline of B lymphocytes in the elderly.

Bleomycin-induced genome structural variations in normal, non-tumor cells.

Many anticancer drugs are genotoxic agents inducing DNA breaks in actively proliferating cancer cells. However, these same drugs also induce mutations, mostly genome structural variations (GSVs). The detection of GSVs in normal cells and tissues is a major challenge due to the very low abundance of these mutations, which are essentially only detectable in clonal outgrowths, such as tumors. Previously we developed Structural Variant Search (SVS) - an NGS-based assay for the quantitative detection of somatic GSVs in normal cells. Using an improved version of SVS we now demonstrate that the same dose of the anti-cancer drug bleomycin induces about 5 times more somatic GSVs in quiescent primary human fibroblasts than in proliferating cells. GVS induction in non-dividing, normal cells was subsequently confirmed in vivo by demonstrating that a single dose of bleomycin leads to a significant increase of GSV frequency in mouse liver and heart, two postmitotic tissues. Our findings suggest that normal non-cycling differentiated cells may serve as a reservoir of iatrogenically induced mutations. These results provide more insight into the possible molecular mechanisms that underlie late-life morbidities in cancer survivors exposed to chemotherapy.

Global, integrated analysis of methylomes and transcriptomes from laser capture microdissected bronchial and alveolar cells in human lung.

Gene regulatory analysis of highly diverse human tissues in vivo is essentially constrained by the challenge of performing genome-wide, integrated epigenetic and transcriptomic analysis in small selected groups of specific cell types. Here we performed genome-wide bisulfite sequencing and RNA-seq from the same small groups of bronchial and alveolar cells isolated by laser capture microdissection from flash-frozen lung tissue of 12 donors and their peripheral blood T cells. Methylation and transcriptome patterns differed between alveolar and bronchial cells, while each of these epithelia showed more differences from mesodermally-derived T cells. Differentially methylated regions (DMRs) between alveolar and bronchial cells tended to locate at regulatory regions affecting promoters of 4,350 genes. A large number of pathways enriched for these DMRs including GTPase signal transduction, cell death, and skeletal muscle. Similar patterns of transcriptome differences were observed: 4,108 differentially expressed genes (DEGs) enriched in GTPase signal transduction, inflammation, cilium assembly, and others. Prioritizing using DMR-DEG regulatory network, we highlighted genes, e.g., ETS1, PPARG, and RXRG, at prominent alveolar vs. bronchial cell discriminant nodes. Our results show that multi-omic analysis of small, highly specific cells is feasible and yields unique physiologic loci distinguishing human lung cell types in situ.

The progeroid phenotype of Ku80 deficiency is dominant over DNA-PKCS deficiency.

Ku80 and DNA-PKCS are both involved in the repair of double strand DNA breaks via the nonhomologous end joining (NHEJ) pathway. While ku80-/- mice exhibit a severely reduced lifespan and size, this phenotype is less pronounced in dna-pkcs-/- mice. However, these observations are based on independent studies with varying genetic backgrounds. Here, we generated ku80-/-, dna-pkcs-/- and double knock out mice in a C57Bl6/J*FVB F1 hybrid background and compared their lifespan, end of life pathology and mutation frequency in liver and spleen using a lacZ reporter. Our data confirm that inactivation of Ku80 and DNA-PKCS causes reduced lifespan and bodyweights, which is most severe in ku80-/- mice. All mutant mice exhibited a strong increase in lymphoma incidence as well as other aging-related pathology (skin epidermal and adnexal atrophy, trabacular bone reduction, kidney tubular anisokaryosis, and cortical and medullar atrophy) and severe lymphoid depletion. LacZ mutation frequency analysis did not show strong differences in mutation frequencies between knock out and wild type mice. The ku80-/- mice had the most severe phenotype and the Ku80-mutation was dominant over the DNA-PKCS-mutation. Presumably, the more severe degenerative effect of Ku80 inactivation on lifespan compared to DNA-PKCS inactivation is caused by additional functions of Ku80 or activity of free Ku70 since both Ku80 and DNA-PKCS are essential for NHEJ.

Deletion of individual Ku subunits in mice causes an NHEJ-independent phenotype potentially by altering apurinic/apyrimidinic site repair.

Ku70 and Ku80 form a heterodimer called Ku that forms a holoenzyme with DNA dependent-protein kinase catalytic subunit (DNA-PKCS) to repair DNA double strand breaks (DSBs) through the nonhomologous end joining (NHEJ) pathway. As expected mutating these genes in mice caused a similar DSB repair-defective phenotype. However, ku70(-/-) cells and ku80(-/-) cells also appeared to have a defect in base excision repair (BER). BER corrects base lesions, apurinic/apyrimidinic (AP) sites and single stand breaks (SSBs) utilizing a variety of proteins including glycosylases, AP endonuclease 1 (APE1) and DNA Polymerase ß (Pol ß). In addition, deleting Ku70 was not equivalent to deleting Ku80 in cells and mice. Therefore, we hypothesized that free Ku70 (not bound to Ku80) and/or free Ku80 (not bound to Ku70) possessed activity that influenced BER. To further test this hypothesis we performed two general sets of experiments. The first set showed that deleting either Ku70 or Ku80 caused an NHEJ-independent defect. We found ku80(-/-) mice had a shorter life span than dna-pkcs(-/-) mice demonstrating a phenotype that was greater than deleting the holoenzyme. We also found Ku70-deletion induced a p53 response that reduced the level of small mutations in the brain suggesting defective BER. We further confirmed that Ku80-deletion impaired BER via a mechanism that was not epistatic to Pol ß. The second set of experiments showed that free Ku70 and free Ku80 could influence BER. We observed that deletion of either Ku70 or Ku80, but not both, increased sensitivity of cells to CRT0044876 (CRT), an agent that interferes with APE1. In addition, free Ku70 and free Ku80 bound to AP sites and in the case of Ku70 inhibited APE1 activity. These observations support a novel role for free Ku70 and free Ku80 in altering BER.

Myc-dependent genome instability and lifespan in Drosophila.

The Myc family of transcription factors are key regulators of cell growth and proliferation that are dysregulated in a large number of human cancers. When overexpressed, Myc family proteins also cause genomic instability, a hallmark of both transformed and aging cells. Using an in vivo lacZ mutation reporter, we show that overexpression of Myc in Drosophila increases the frequency of large genome rearrangements associated with erroneous repair of DNA double-strand breaks (DSBs). In addition, we find that overexpression of Myc shortens adult lifespan and, conversely, that Myc haploinsufficiency reduces mutation load and extends lifespan. Our data provide the first evidence that Myc may act as a pro-aging factor, possibly through its ability to greatly increase genome instability.

Loss of the bloom syndrome helicase increases DNA ligase 4-independent genome rearrangements and tumorigenesis in aging Drosophila.

BACKGROUND: The BLM DNA helicase plays a vital role in maintaining genome stability. Mutations in BLM cause Bloom syndrome, a rare disorder associated with cancer predisposition and premature aging. Humans and mice with blm mutations have increased frequencies of spontaneous mutagenesis, but the molecular basis of this increase is not well understood. In addition, the effect of aging on spontaneous mutagenesis in blm mutants has not been characterized. To address this, we used a lacZ reporter system in wild-type and several mutant strains of Drosophila melanogaster to analyze mechanisms of mutagenesis throughout their lifespan. RESULTS: Our data show that Drosophila lacking BLM have an elevated frequency of spontaneous genome rearrangements that increases with age. Although in normal flies most genome rearrangements occur through DNA ligase 4-dependent classical end joining, most rearrangements that accumulate during aging in blm mutants do not require DNA ligase 4, suggesting the influence of an alternative end-joining mechanism. Adult blm mutants also display reduced lifespan and ligase 4-independent enhanced tumorigenesis in mitotically active tissues. CONCLUSIONS: These results suggest that Drosophila BLM suppresses error-prone alternative end-joining repair of DNA double-strand breaks that can result in genome instability and tumor formation during aging. In addition, since loss of BLM significantly affects lifespan and tumorigenesis, the data provide a link between error-prone end joining, genome rearrangements, and tumor formation in a model metazoan.

Pleiotropic impact of constitutive fosB inactivation on nicotine-induced behavioral alterations and stress-related traits in mice.

Multiple genes are thought to influence both susceptibility to nicotine dependence and its comorbid behavioral traits in humans. However, which specific genes contribute to this pleiotropic effect is poorly understood. Previous rodent studies have shown that many addictive substances and stressful stimuli increase the expression of the transcription factor FosB in limbic and associated regions and that the protein products of fosB contribute to certain behavioral effects of cocaine and morphine. However, the role of this gene in nicotine-regulated behaviors and dependence-related behavioral traits is unknown. We tested the hypothesis that a constitutive level of FosB affects nicotine-regulated behaviors and comorbid behavioral traits using constitutive fosB knockout (KO) mice. Following repeated or prolonged nicotine administration, but not a single acute administration, KO mice were impaired in conditioned place preference, oral nicotine intake and motor suppression. In wild-type mice, repeated nicotine injections, but not a single acute injection, increased the expression of FosB and its truncated variant DeltaFosB in the targets but not at the origins of the mesolimbic and nigrostriatal dopamine pathways; no detectable level of FosB/DeltaFosB was found in KO mice. In tasks designed to assess behavioral traits, KO mice exhibited more pronounced behavioral abnormalities when stress levels were high than when they were minimized. Our results suggest that the constitutive absence of fosB has a pleiotropic influence on the behavioral effects of repeated or prolonged nicotine administration and on stress-related behavioral traits in mice.

Monoamine oxidase A knockout mice exhibit impaired nicotine preference but normal responses to novel stimuli.

Nicotine is thought to act on brain monoamine systems that normally mediate diverse motivational behaviors. How monoamine-related genes contribute to behavioral traits (e.g. responses to novel stimuli) comorbid with the susceptibility to nicotine addiction is still poorly understood. We examined the impact of constitutive monoamine oxidase A (MAOA) deficiency in mice on nicotine reward and responses to novel stimuli. Age-matched, male Maoa-knockout (KO) mice and wild-type (WT) littermates were tested for nicotine-induced conditioned place preference (CPP); voluntary oral nicotine preference/intake; spontaneous locomotor activity in a novel, inescapable open field; and novelty place preference. Nicotine preference in WT mice was reduced in Maoa-KO mice in the CPP and oral preference/intake tests. Control experiments showed that these phenotypes were not due to abnormalities in nicotine metabolism, fluid intake or response to taste. In contrast, Maoa-KO mice were normal in their behavioral response to a novel, inescapable open field and in their preference for a novel place. The observed phenotypes suggest that a constitutive deficiency of MAOA reduces the rewarding effects of nicotine without altering behavioral responses to novel stimuli in mice. Constitutive MAOA activity levels are likely to contribute to the vulnerability or resiliency to nicotine addiction by altering the rewarding effects of nicotine.

DARPP-32 phosphorylation opposes the behavioral effects of nicotine.

BACKGROUND: The addictive properties of nicotine are mediated via dopaminergic pathways and their post-synaptic neurons in the striatum. Because post-synaptic neurons within the striatum contain high levels of the dopamine- and cAMP-regulated phosphoprotein of 32 kDa (DARPP-32), we hypothesized that DARPP-32 may functionally contribute to the behavioral effects of nicotine. METHODS: We examined the behavioral effects of nicotine and the phosphorylation state of DARPP-32 in wild-type (WT) and DARPP-32 knockout (KO) mice. In one experiment, we assessed voluntary nicotine intake (0-50 microg/ml) of WT and KO mice in a two-bottle choice paradigm. In a separate experiment, the motor-depressant effects of acute and repeated nicotine injections (0-.8 mg/kg, subcutaneously [SC]) were assessed. The phosphorylation of DARPP-32 at threonine34 and threonine75 were examined using Western blotting. RESULTS: A heightened responsiveness to nicotine was seen in KO mice when compared with WT mice in oral intake and motor depression. The enhanced responsiveness in KO mice was not due to alterations in taste sensations, fluid intake, or blood nicotine or cotinine levels. Systemic injections of nicotine resulted in increased striatal DARPP-32 phosphorylation at threonine34 and threonine75. CONCLUSIONS: DARPP-32 opposes the behavioral effects of nicotine possibly via concurrent phosphorylation at the two threonine sites.

Chronic treatment with atypical neuroleptics induces striosomal FosB/DeltaFosB expression in rats.

BACKGROUND: Studies have shown that neuroleptics regulate expression of the transcription factor FosB/DeltaFosB in the striatum, including the accumbens and caudate-putamen; however, the striatum is also divided into another structural dimension, the striosome and matrix compartments. The precise distribution of FosB/DeltaFosB induced by chronic neuroleptics in these striatal compartments is poorly understood. METHODS: Rats received either single acute injections or chronic injections of clozapine (0 or 20 mg/kg, intraperitoneally [IP]), olanzapine (0 or 5 mg/kg, IP), or haloperidol (0 or 1.5 mg/kg, IP) for 25 days. The levels and compartmental distribution of FosB/DeltaFosB were examined. RESULTS: Chronic clozapine induced clustered FosB/DeltaFosB expression within striosomes of the caudate-putamen. This pattern was due to increased levels of FosB/DeltaFosB in striosomes within the ventrolateral caudate-putamen and reduced levels of basal FosB/DeltaFosB in the matrix in the entire caudate-putamen. In contrast, chronic haloperidol increased FosB/DeltaFosB equally within the matrix and striosomes throughout the entire caudate-putamen. Chronic olanzapine induced an intermediate pattern. CONCLUSIONS: The relative absence of FosB/DeltaFosB expression in the matrix correlates with the lack of parkinsonism of atypical neuroleptics. Expression of FosB/DeltaFosB in the matrix may contribute to parkinsonism of typical neuroleptics.