Contrasting somatic mutation patterns in aging human neurons and oligodendrocytes.

Characterizing somatic mutations in the brain is important for disentangling the complex mechanisms of aging, yet little is known about mutational patterns in different brain cell types. Here, we performed whole-genome sequencing (WGS) of 86 single oligodendrocytes, 20 mixed glia, and 56 single neurons from neurotypical individuals spanning 0.4-104 years of age and identified >92,000 somatic single-nucleotide variants (sSNVs) and small insertions/deletions (indels). Although both cell types accumulate somatic mutations linearly with age, oligodendrocytes accumulated sSNVs 81% faster than neurons and indels 28% slower than neurons. Correlation of mutations with single-nucleus RNA profiles and chromatin accessibility from the same brains revealed that oligodendrocyte mutations are enriched in inactive genomic regions and are distributed across the genome similarly to mutations in brain cancers. In contrast, neuronal mutations are enriched in open, transcriptionally active chromatin. These stark differences suggest an assortment of active mutagenic processes in oligodendrocytes and neurons.

ERα-associated translocations underlie oncogene amplifications in breast cancer.

Focal copy-number amplification is an oncogenic event. Although recent studies have revealed the complex structure1-3 and the evolutionary trajectories4 of oncogene amplicons, their origin remains poorly understood. Here we show that focal amplifications in breast cancer frequently derive from a mechanism-which we term translocation-bridge amplification-involving inter-chromosomal translocations that lead to dicentric chromosome bridge formation and breakage. In 780 breast cancer genomes, we observe that focal amplifications are frequently connected to each other by inter-chromosomal translocations at their boundaries. Subsequent analysis indicates the following model: the oncogene neighbourhood is translocated in G1 creating a dicentric chromosome, the dicentric chromosome is replicated, and as dicentric sister chromosomes segregate during mitosis, a chromosome bridge is formed and then broken, with fragments often being circularized in extrachromosomal DNAs. This model explains the amplifications of key oncogenes, including ERBB2 and CCND1. Recurrent amplification boundaries and rearrangement hotspots correlate with oestrogen receptor binding in breast cancer cells. Experimentally, oestrogen treatment induces DNA double-strand breaks in the oestrogen receptor target regions that are repaired by translocations, suggesting a role of oestrogen in generating the initial translocations. A pan-cancer analysis reveals tissue-specific biases in mechanisms initiating focal amplifications, with the breakage-fusion-bridge cycle prevalent in some and the translocation-bridge amplification in others, probably owing to the different timing of DNA break repair. Our results identify a common mode of oncogene amplification and propose oestrogen as its mechanistic origin in breast cancer.

DOCK2 Sets the Threshold for Entry into the Virtual Memory CD8+ T Cell Compartment by Negatively Regulating Tonic TCR Triggering.

The control of cytoskeletal dynamics by dedicator of cytokinesis 2 (DOCK2), a hematopoietic cell-specific actin effector protein, has been implicated in TCR signaling and T cell migration. Biallelic mutations in Dock2 have been identified in patients with a recessive form of combined immunodeficiency with defects in T, B, and NK cell activation. Surprisingly, we show in this study that certain immune functions of CD8+ T cells are enhanced in the absence of DOCK2. Dock2-deficient mice have a pronounced expansion of their memory T cell compartment. Bone marrow chimera and adoptive transfer studies indicate that these memory T cells develop in a cell-intrinsic manner following thymic egress. Transcriptional profiling, TCR repertoire analyses, and cell surface marker expression indicate that Dock2-deficient naive CD8+ T cells directly convert into virtual memory cells without clonal effector T cell expansion. This direct conversion to memory is associated with a selective increase in TCR sensitivity to self-peptide MHC in vivo and an enhanced response to weak agonist peptides ex vivo. In contrast, the response to strong agonist peptides remains unaltered in Dock2-deficient T cells. Collectively, these findings suggest that the regulation of the actin dynamics by DOCK2 enhances the threshold for entry into the virtual memory compartment by negatively regulating tonic TCR triggering in response to weak agonists.

Polygenic risk scores for autoimmune related diseases are significantly different in cancer exceptional responders.

A small number of cancer patients respond exceptionally well to therapies and survive significantly longer than patients with similar diagnoses. Profiling the germline genetic backgrounds of exceptional responder (ER) patients, with extreme survival times, can yield insights into the germline polymorphisms that influence response to therapy. As ERs showed a high incidence in autoimmune diseases, we hypothesized the differences in autoimmune disease risk could reflect the immune background of ERs and contribute to better cancer treatment responses. We analyzed the germline variants of 51 ERs using polygenic risk score (PRS) analysis. Compared to typical cancer patients, the ERs had significantly elevated PRSs for several autoimmune-related diseases: type 1 diabetes, hypothyroidism, and psoriasis. This indicates that an increased genetic predisposition towards these autoimmune diseases is more prevalent among the ERs. In contrast, ERs had significantly lower PRSs for developing inflammatory bowel disease. The left-skew of type 1 diabetes score was significant for exceptional responders. Variants on genes involved in the T1D PRS model associated with cancer drug response are more likely to co-occur with other variants among ERs. In conclusion, ERs exhibited different risks for autoimmune diseases compared to typical cancer patients, which suggests that changes in a patient's immune set point or immune surveillance specificity could be a potential mechanistic link to their exceptional response. These findings expand upon previous research on immune checkpoint inhibitor-treated patients to include those who received chemotherapy or radiotherapy.

High-resolution detection of copy number alterations in single cells with HiScanner.

Improvements in single-cell whole-genome sequencing (scWGS) assays have enabled detailed characterization of somatic copy number alterations (CNAs) at the single-cell level. Yet, current computational methods are mostly designed for detecting chromosome-scale changes in cancer samples with low sequencing coverage. Here, we introduce HiScanner (High-resolution Single-Cell Allelic copy Number callER), which combines read depth, B-allele frequency, and haplotype phasing to identify CNAs with high resolution. In simulated data, HiScanner consistently outperforms state-of-the-art methods across various CNA types and sizes. When applied to high-coverage scWGS data from human brain cells, HiScanner shows a superior ability to detect smaller CNAs, uncovering distinct CNA patterns between neurons and oligodendrocytes. For 179 cells we sequenced from longitudinal meningioma samples, integration of CNAs with point mutations revealed evolutionary trajectories of tumor cells. These findings show that HiScanner enables accurate characterization of frequency, clonality, and distribution of CNAs at the single-cell level in both non-neoplastic and neoplastic cells.

Cell lineage analysis with somatic mutations reveals late divergence of neuronal cell types and cortical areas in human cerebral cortex.

The mammalian cerebral cortex shows functional specialization into regions with distinct neuronal compositions, most strikingly in the human brain, but little is known in about how cellular lineages shape cortical regional variation and neuronal cell types during development. Here, we use somatic single nucleotide variants (sSNVs) to map lineages of neuronal sub-types and cortical regions. Early-occurring sSNVs rarely respect Brodmann area (BA) borders, while late-occurring sSNVs mark neuron-generating clones with modest regional restriction, though descendants often dispersed into neighboring BAs. Nevertheless, in visual cortex, BA17 contains 30-70% more sSNVs compared to the neighboring BA18, with clones across the BA17/18 border distributed asymmetrically and thus displaying different cortex-wide dispersion patterns. Moreover, we find that excitatory neuron-generating clones with modest regional restriction consistently share low-mosaic sSNVs with some inhibitory neurons, suggesting significant co-generation of excitatory and some inhibitory neurons in the dorsal cortex. Our analysis reveals human-specific cortical cell lineage patterns, with both regional inhomogeneities in progenitor proliferation and late divergence of excitatory/inhibitory lineages.

Comprehensive identification of transposable element insertions using multiple sequencing technologies.

Transposable elements (TEs) help shape the structure and function of the human genome. When inserted into some locations, TEs may disrupt gene regulation and cause diseases. Here, we present xTea (x-Transposable element analyzer), a tool for identifying TE insertions in whole-genome sequencing data. Whereas existing methods are mostly designed for short-read data, xTea can be applied to both short-read and long-read data. Our analysis shows that xTea outperforms other short read-based methods for both germline and somatic TE insertion discovery. With long-read data, we created a catalogue of polymorphic insertions with full assembly and annotation of insertional sequences for various types of retroelements, including pseudogenes and endogenous retroviruses. Notably, we find that individual genomes have an average of nine groups of full-length L1s in centromeres, suggesting that centromeres and other highly repetitive regions such as telomeres are a significant yet unexplored source of active L1s. xTea is available at https://github.com/parklab/xTea .

B1a and B2 cells are characterized by distinct CpG modification states at DNMT3A-maintained enhancers.

The B1 and B2 lineages of B cells contribute to protection from pathogens in distinct ways. The role of the DNA CpG methylome in specifying these two B-cell fates is still unclear. Here we profile the CpG modifications and transcriptomes of peritoneal B1a and follicular B2 cells, as well as their respective proB cell precursors in the fetal liver and adult bone marrow from wild-type and CD19-Cre Dnmt3a floxed mice lacking DNMT3A in the B lineage. We show that an underlying foundational CpG methylome is stably established during B lineage commitment and is overlaid with a DNMT3A-maintained dynamic methylome that is sculpted in distinct ways in B1a and B2 cells. This dynamic DNMT3A-maintained methylome is composed of novel enhancers that are closely linked to lineage-specific genes. While DNMT3A maintains the methylation state of these enhancers in both B1a and B2 cells, the dynamic methylome undergoes a prominent programmed demethylation event during B1a but not B2 cell development. We propose that the methylation pattern of DNMT3A-maintained enhancers is determined by the coincident recruitment of DNMT3A and TET enzymes, which regulate the developmental expression of B1a and B2 lineage-specific genes.