Lymphoma Driver Mutations in the Pathogenic Evolution of an Iconic Human Autoantibody.

Pathogenic autoantibodies arise in many autoimmune diseases, but it is not understood how the cells making them evade immune checkpoints. Here, single-cell multi-omics analysis demonstrates a shared mechanism with lymphoid malignancy in the formation of public rheumatoid factor autoantibodies responsible for mixed cryoglobulinemic vasculitis. By combining single-cell DNA and RNA sequencing with serum antibody peptide sequencing and antibody synthesis, rare circulating B lymphocytes making pathogenic autoantibodies were found to comprise clonal trees accumulating mutations. Lymphoma driver mutations in genes regulating B cell proliferation and V(D)J mutation (CARD11, TNFAIP3, CCND3, ID3, BTG2, and KLHL6) were present in rogue B cells producing the pathogenic autoantibody. Antibody V(D)J mutations conferred pathogenicity by causing the antigen-bound autoantibodies to undergo phase transition to insoluble aggregates at lower temperatures. These results reveal a pre-neoplastic stage in human lymphomagenesis and a cascade of somatic mutations leading to an iconic pathogenic autoantibody.

Analysis of antibodies from HCV elite neutralizers identifies genetic determinants of broad neutralization.

The high genetic diversity of hepatitis C virus (HCV) complicates effective vaccine development. We screened a cohort of 435 HCV-infected individuals and found that 2%-5% demonstrated outstanding HCV-neutralizing activity. From four of these patients, we isolated 310 HCV antibodies, including neutralizing antibodies with exceptional breadth and potency. High neutralizing activity was enabled by the use of the VH1-69 heavy-chain gene segment, somatic mutations within CDRH1, and CDRH2 hydrophobicity. Structural and mutational analyses revealed an important role for mutations replacing the serines at positions 30 and 31, as well as the presence of neutral and hydrophobic residues at the tip of the CDRH3. Based on these characteristics, we computationally created a de novo antibody with a fully synthetic VH1-69 heavy chain that efficiently neutralized multiple HCV genotypes. Our findings provide a deep understanding of the generation of broadly HCV-neutralizing antibodies that can guide the design of effective vaccine candidates.

Mosaicism in Human Health and Disease.

Mosaicism refers to the occurrence of two or more genomes in an individual derived from a single zygote. Germline mosaicism is a mutation that is limited to the gonads and can be transmitted to offspring. Somatic mosaicism is a postzygotic mutation that occurs in the soma, and it may occur at any developmental stage or in adult tissues. Mosaic variation may be classified in six ways: (a) germline or somatic origin, (b) class of DNA mutation (ranging in scale from single base pairs to multiple chromosomes), (c) developmental context, (d) body location(s), (e) functional consequence (including deleterious, neutral, or advantageous), and (f) additional sources of mosaicism, including mitochondrial heteroplasmy, exogenous DNA sources such as vectors, and epigenetic changes such as imprinting and X-chromosome inactivation. Technological advances, including single-cell and other next-generation sequencing, have facilitated improved sensitivity and specificity to detect mosaicism in a variety of biological contexts.

Somatic mutations in arteriovenous malformations in hereditary hemorrhagic telangiectasia support a bi-allelic two-hit mutation mechanism of pathogenesis.

Hereditary hemorrhagic telangiectasia (HHT) is an inherited disorder of vascular malformations characterized by mucocutaneous telangiectases and arteriovenous malformations (AVMs) in internal organs. HHT is caused by inheritance of a loss of function mutation in one of three genes. Although individuals with HHT are haploinsufficient for one of these genes throughout their entire body, rather than exhibiting a systemic vascular phenotype, vascular malformations occur as focal lesions in discrete anatomic locations. The inconsistency between genotype and phenotype has provoked debate over whether haploinsufficiency or a different mechanism gives rise to the vascular malformations. We previously showed that HHT-associated skin telangiectases develop by a two-hit mutation mechanism in an HHT gene. However, somatic mutations were identified in only half of the telangiectases, raising the question whether a second-hit somatic mutation is a necessary (required) event in HHT pathogenesis. Here, we show that another mechanism for the second hit is loss of heterozygosity across the chromosome bearing the germline mutation. Secondly, we investigate the two-hit mutation mechanism for internal organ AVMs, the source of much of the morbidity of HHT. Here, we identified somatic molecular genetic events in eight liver telangiectases, including point mutations and a loss of heterozygosity event. We also identified somatic mutations in one pulmonary AVM and two brain AVMs, confirming that mucocutaneous and internal organ vascular malformations undergo the same molecular mechanisms. Together, these data argue that bi-allelic loss of function in an HHT gene is a required event in the pathogenesis of HHT-associated vascular malformations.

GNA14 and GNAQ somatic mutations cause spinal and intracranial extra-axial cavernous hemangiomas.

Extra-axial cavernous hemangiomas (ECHs) are complex vascular lesions mainly found in the spine and cavernous sinus. Their removal poses significant risk due to their vascularity and diffuse nature, and their genetic underpinnings remain incompletely understood. Our approach involved genetic analyses on 31 tissue samples of ECHs employing whole-exome sequencing and targeted deep sequencing. We explored downstream signaling pathways, gene expression changes, and resultant phenotypic shifts induced by these mutations, both in vitro and in vivo. In our cohort, 77.4% of samples had somatic missense variants in GNA14, GNAQ, or GJA4. Transcriptomic analysis highlighted significant pathway upregulation, with the GNAQ c.626A>G (p.Gln209Arg) mutation elevating PI3K-AKT-mTOR and angiogenesis-related pathways, while GNA14 c.614A>T (p.Gln205Leu) mutation led to MAPK and angiogenesis-related pathway upregulation. Using a mouse xenograft model, we observed enlarged vessels from these mutations. Additionally, we initiated rapamycin treatment in a 14-year-old individual harboring the GNAQ c.626A>G (p.Gln209Arg) variant, resulting in gradual regression of cutaneous cavernous hemangiomas and improved motor strength, with minimal side effects. Understanding these mutations and their pathways provides a foundation for developing therapies for ECHs resistant to current therapies. Indeed, the administration of rapamycin in an individual within this study highlights the promise of targeted treatments in treating these complex lesions.

Somatic mutations in tumor-infiltrating lymphocytes impact on antitumor immunity.

Immune checkpoint inhibitors (ICIs) exert clinical efficacy against various types of cancers by reinvigorating exhausted CD8+ T cells that can expand and directly attack cancer cells (cancer-specific T cells) among tumor-infiltrating lymphocytes (TILs). Although some reports have identified somatic mutations in TILs, their effect on antitumor immunity remains unclear. In this study, we successfully established 18 cancer-specific T cell clones, which have an exhaustion phenotype, from the TILs of four patients with melanoma. We conducted whole-genome sequencing for these T cell clones and identified various somatic mutations in them with high clonality. Among the somatic mutations, an SH2D2A loss-of-function frameshift mutation and TNFAIP3 deletion could activate T cell effector functions in vitro. Furthermore, we generated CD8+ T cell-specific Tnfaip3 knockout mice and showed that Tnfaip3 function loss in CD8+ T cell increased antitumor immunity, leading to remarkable response to PD-1 blockade in vivo. In addition, we analyzed bulk CD3+ T cells from TILs in additional 12 patients and identified an SH2D2A mutation in one patient through amplicon sequencing. These findings suggest that somatic mutations in TILs can affect antitumor immunity and suggest unique biomarkers and therapeutic targets.

BRN2 suppresses apoptosis, reprograms DNA damage repair, and is associated with a high somatic mutation burden in melanoma.

Whether cell types exposed to a high level of environmental insults possess cell type-specific prosurvival mechanisms or enhanced DNA damage repair capacity is not well understood. BRN2 is a tissue-restricted POU domain transcription factor implicated in neural development and several cancers. In melanoma, BRN2 plays a key role in promoting invasion and regulating proliferation. Here we found, surprisingly, that rather than interacting with transcription cofactors, BRN2 is instead associated with DNA damage response proteins and directly binds PARP1 and Ku70/Ku80. Rapid PARP1-dependent BRN2 association with sites of DNA damage facilitates recruitment of Ku80 and reprograms DNA damage repair by promoting Ku-dependent nonhomologous end-joining (NHEJ) at the expense of homologous recombination. BRN2 also suppresses an apoptosis-associated gene expression program to protect against UVB-, chemotherapy- and vemurafenib-induced apoptosis. Remarkably, BRN2 expression also correlates with a high single-nucleotide variation prevalence in human melanomas. By promoting error-prone DNA damage repair via NHEJ and suppressing apoptosis of damaged cells, our results suggest that BRN2 contributes to the generation of melanomas with a high mutation burden. Our findings highlight a novel role for a key transcription factor in reprogramming DNA damage repair and suggest that BRN2 may impact the response to DNA-damaging agents in BRN2-expressing cancers.

Somatic Mutagenesis in Mammals and Its Implications for Human Disease and Aging.

DNA mutations as a consequence of errors during DNA damage repair, replication, or mitosis are the substrate for evolution. In multicellular organisms, mutations can occur in the germline and also in somatic tissues, where they are associated with cancer and other chronic diseases and possibly with aging. Recent advances in high-throughput sequencing have made it relatively easy to study germline de novo mutations, but in somatic cells, the vast majority of mutations are low-abundant and can be detected only in clonal lineages, such as tumors, or single cells. Here we review recent results on somatic mutations in normal human and animal tissues with a focus on their possible functional consequences.

Molecular mechanisms insulating proliferation from genotoxic stress in B lymphocytes.

In mammals, B cells strictly segregate proliferation from somatic mutation as they develop within the bone marrow and then mature through germinal centers (GCs) in the periphery. Failure to do so risks autoimmunity and neoplastic transformation. Recent work has described how B cell progenitors transition between proliferation and mutation via cytokine signaling pathways, epigenetic chromatin regulation, and remodeling of 3D chromatin conformation. We propose a three-zone model of the GC that describes how proliferation and mutation are regulated. Using this model, we consider how recent mechanistic discoveries in B cell progenitors inform models of GC B cell function and reveal fundamental mechanisms underpinning humoral immunity, autoimmunity, and lymphomagenesis.

Transcriptome-based measurement of CD8+ T cell age and its applications.

The ability of T cells to undergo robust cell division in response to antigenic stimulation is essential for competent T cell function. However, this ability is reduced with aging and contributes to increased susceptibility to infectious diseases, cancers, and other diseases among older adults. To better understand T cell aging, improved measurements of age-related cellular changes in T cells are necessary. The recent development of machine learning (ML)-assisted transcriptome-based quantification of individual CD8+ T cell age represents a significant step forward in this regard. It reveals both prominent and subtle changes in gene expression and points to potential functional alterations of CD8+ T cells with aging. I argue that single-cell transcriptome-based age prediction in the immune system may have promising future applications.

Transforming mutations in the development of pathogenic B cell clones and autoantibodies.

Autoimmune diseases are characterized by serum autoantibodies, some of which are pathogenic, causing severe manifestations and organ injury. However, autoantibodies of the same antigenic reactivity are also present in the serum of asymptomatic people years before they develop any clinical signs of autoimmunity. Autoantibodies can arise during multiple stages of B cell development, and various genetic and environmental factors drive their production. However, what drives the development of pathogenic autoantibodies is poorly understood. Advances in single-cell technology have enabled the deep analysis of rare B cell clones producing pathogenic autoantibodies responsible for vasculitis in patients with primary Sjögren's syndrome complicated by mixed cryoglobulinaemia. These findings demonstrated a cascade of genetic events involving stereotypic immunoglobulin V(D)J recombination and transforming somatic mutations in lymphoma genes and V(D)J regions that disrupted antibody quality control mechanisms and decreased autoantibody solubility. Most studies consider V(D)J mutations that enhance autoantibody affinity to drive pathology; however, V(D)J mutations that increase autoantibody propensity to form insoluble complexes could be a major contributor to autoantibody pathogenicity. Defining the molecular characteristics of pathogenic autoantibodies and failed tolerance checkpoints driving their formation will improve prognostication, enabling early treatment to prevent escalating organ damage and B cell malignancy.

Mosaicism in Cutaneous Disorders.

Genetic mosaicism arises when a zygote harbors two or more distinct genotypes, typically due to de novo, somatic mutation during embryogenesis. The clinical manifestations largely depend on the differentiation status of the mutated cell; earlier mutations target pluripotent cells and generate more widespread disease affecting multiple organ systems. If gonadal tissue is spared-as in somatic genomic mosaicism-the mutation and its effects are limited to the proband, whereas mosaicism also affecting the gametes, such as germline or gonosomal mosaicism, is transmissible. Mosaicism is easily appreciated in cutaneous disorders, as phenotypically distinct mutant cells often give rise to lesions in patterns determined by the affected cell type. Genetic investigation of cutaneous mosaic disorders has identified pathways central to disease pathogenesis, revealing novel therapeutic targets. In this review, we discuss examples of cutaneous mosaicism, approaches to gene discovery in these disorders, and insights into molecular pathobiology that have potential for clinical translation.

Integrating somatic mutation profiles with structural deep clustering network for metabolic stratification in pancreatic cancer: a comprehensive analysis of prognostic and genomic landscapes.

Pancreatic cancer is a globally recognized highly aggressive malignancy, posing a significant threat to human health and characterized by pronounced heterogeneity. In recent years, researchers have uncovered that the development and progression of cancer are often attributed to the accumulation of somatic mutations within cells. However, cancer somatic mutation data exhibit characteristics such as high dimensionality and sparsity, which pose new challenges in utilizing these data effectively. In this study, we propagated the discrete somatic mutation data of pancreatic cancer through a network propagation model based on protein-protein interaction networks. This resulted in smoothed somatic mutation profile data that incorporate protein network information. Based on this smoothed mutation profile data, we obtained the activity levels of different metabolic pathways in pancreatic cancer patients. Subsequently, using the activity levels of various metabolic pathways in cancer patients, we employed a deep clustering algorithm to establish biologically and clinically relevant metabolic subtypes of pancreatic cancer. Our study holds scientific significance in classifying pancreatic cancer based on somatic mutation data and may provide a crucial theoretical basis for the diagnosis and immunotherapy of pancreatic cancer patients.

Threshold of somatic mosaicism leading to brain dysfunction with focal epilepsy.

Somatic mosaicism in a fraction of brain cells causes neurodevelopmental disorders, including childhood intractable epilepsy. However, the threshold for somatic mosaicism leading to brain dysfunction is unknown. In this study, we induced various mosaic burdens in focal cortical dysplasia type II (FCD II) mice, featuring mTOR somatic mosaicism and spontaneous behavioural seizures. The mosaic burdens ranged from approximately 1000 to 40 000 neurons expressing the mTOR mutant in the somatosensory or medial prefrontal cortex. Surprisingly, approximately 8000-9000 neurons expressing the MTOR mutant, extrapolated to constitute 0.08%-0.09% of total cells or roughly 0.04% of variant allele frequency in the mouse hemicortex, were sufficient to trigger epileptic seizures. The mutational burden was correlated with seizure frequency and onset, with a higher tendency for electrographic inter-ictal spikes and beta- and gamma-frequency oscillations in FCD II mice exceeding the threshold. Moreover, mutation-negative FCD II patients in deep sequencing of their bulky brain tissues revealed somatic mosaicism of the mTOR pathway genes as low as 0.07% in resected brain tissues through ultra-deep targeted sequencing (up to 20 million reads). Thus, our study suggests that extremely low levels of somatic mosaicism can contribute to brain dysfunction.

Focal cortical dysplasia II caused by brain somatic mutation of IRS-1 is associated with ERK signaling pathway activation.

Somatic mutations have been identified in 10% to 63% of focal cortical dysplasia type II samples, primarily linked to the mTOR pathway. When the causative genetic mutations are not identified, this opens the possibility of discovering new pathogenic genes or pathways that could be contributing to the condition. In our previous study, we identified a novel candidate pathogenic somatic variant of IRS-1 c.1791dupG in the brain tissue of a child with focal cortical dysplasia type II. This study further explored the variant's role in causing type II focal cortical dysplasia through in vitro overexpression in 293T and SH-SY5Y cells and in vivo evaluation via in utero electroporation in fetal brains, assessing effects on neuronal migration, morphology, and network integrity. It was found that the mutant IRS-1 variant led to hyperactivity of p-ERK, increased cell volume, and was predominantly associated with the MAPK signaling pathway. In vivo, the IRS-1 c.1791dupG variant induced abnormal neuron migration, cytomegaly, and network hyperexcitability. Notably, the ERK inhibitor GDC-0994, rather than the mTOR inhibitor rapamycin, effectively rescued the neuronal defects. This study directly highlighted the ERK signaling pathway's role in the pathogenesis of focal cortical dysplasia II and provided a new therapeutic target for cases of focal cortical dysplasia II that are not treatable by rapamycin analogs.

A quantification method of somatic mutations in normal tissues and their accumulation in pediatric patients with chemotherapy.

Somatic mutations are accumulated in normal human tissues with aging and exposure to carcinogens. If we can accurately count any passenger mutations in any single DNA molecule, since their quantity is much larger than driver mutations, we can sensitively detect mutation accumulation in polyclonal normal tissues. Duplex sequencing, which tags both DNA strands in one DNA molecule, enables accurate count of such mutations, but requires a very large number of sequencing reads for each single sample of human-genome size. Here, we reduced the genome size to 1/90 using the BamHI restriction enzyme and established a cost-effective pipeline. The enzymatically cleaved and optimal sequencing (EcoSeq) method was able to count somatic mutations in a single DNA molecule with a sensitivity of as low as 3 × 10-8 per base pair (bp), as assessed by measuring artificially prepared mutations. Taking advantages of EcoSeq, we analyzed normal peripheral blood cells of pediatric sarcoma patients who received chemotherapy (n = 10) and those who did not (n = 10). The former had a mutation frequency of 31.2 ± 13.4 × 10-8 per base pair while the latter had 9.0 ± 4.5 × 10-8 per base pair (P < 0.001). The increase in mutation frequency was confirmed by analysis of the same patients before and after chemotherapy, and increased mutation frequencies persisted 46 to 64 mo after chemotherapy, indicating that the mutation accumulation constitutes a risk of secondary leukemia. EcoSeq has the potential to reveal accumulation of somatic mutations and exposure to environmental factors in any DNA samples and will contribute to cancer risk estimation.

Rethinking cancer initiation: The role of large-scale mutational events.

Cancer initiation is revisited in light of recent discoveries in cancer pathogenesis. Of note is the detection of mutated cancer genes in benign conditions. More significantly, somatic clones, which harbor mutations in cancer genes, arise in normal tissues from early development through adulthood, but seldom do they transform into cancer. Further, clustered mutational events-kataegis, chromothripsis and chromoplexy-are widespread in cancer, generating point mutations and chromosomal rearrangements in a single cellular catastrophe. These observations are contrary to the prevailing somatic mutation theory, which states that a cancer is caused by the gradual accumulation of mutations over time. A different perspective is proposed within the framework of Waddington's epigenetic landscape wherein tumorigenesis is viewed primarily as a disruption of cell development. Cell types are defined by their specific gene-expression profiles, determined by the gene regulatory network, and can be regarded as attractor states of the network dynamics: they represent specific, self-stabilizing patterns of gene activities across the genome. However, large-scale mutational events reshape the landscape topology, creating abnormal "unphysiological" attractors. This is the crux of the process of initiation. Initiation primes the cell for conversion into a tumor phenotype by oncogenes and tumor suppressor genes, which drive cell proliferation and clonal diversification. This view of tumorigenesis calls for a different approach to therapy.

Identification of hypoxia-related gene signatures based on multi-omics analysis in lung adenocarcinoma.

Lung adenocarcinoma (LUAD) is the most common type of lung cancer and one of the malignancies with the highest incidence rate and mortality worldwide. Hypoxia is a typical feature of tumour microenvironment (TME), which affects the progression of LUAD from multiple molecular levels. However, the underlying molecular mechanisms behind LUAD hypoxia are not fully understood. In this study, we estimated the level of hypoxia by calculating a score based on 15 hypoxia genes. The hypoxia scores were relatively high in LUAD patients with poor prognosis and were bound up with tumour node metastasis (TNM) stage, tumour size, lymph node, age and gender. By comparison of high hypoxia score group and low hypoxia score group, 1820 differentially expressed genes were identified, among which up-regulated genes were mainly about cell division and proliferation while down-regulated genes were primarily involved in cilium-related biological processes. Besides, LUAD patients with high hypoxia scores had higher frequencies of gene mutations, among which TP53, TTN and MUC16 had the highest mutation rates. As for DNA methylation, 1015 differentially methylated probes-related genes were found and may play potential roles in tumour-related neurobiological processes and cell signal transduction. Finally, a prognostic model with 25 multi-omics features was constructed and showed good predictive performance. The area under curve (AUC) values of 1-, 3- and 5-year survival reached 0.863, 0.826 and 0.846, respectively. Above all, our findings are helpful in understanding the impact and molecular mechanisms of hypoxia in LUAD.

What can we learn about acid-base transporters in cancer from studying somatic mutations in their genes?

Acidosis is a chemical signature of the tumour microenvironment that challenges intracellular pH homeostasis. The orchestrated activity of acid-base transporters of the solute-linked carrier (SLC) family is critical for removing the end-products of fermentative metabolism (lactate/H+) and maintaining a favourably alkaline cytoplasm. Given the critical role of pH homeostasis in enabling cellular activities, mutations in relevant SLC genes may impact the oncogenic process, emerging as negatively or positively selected, or as driver or passenger mutations. To address this, we performed a pan-cancer analysis of The Cancer Genome Atlas simple nucleotide variation data for acid/base-transporting SLCs (ABT-SLCs). Somatic mutation patterns of monocarboxylate transporters (MCTs) were consistent with their proposed essentiality in facilitating lactate/H+ efflux. Among all cancers, tumours of uterine corpus endometrial cancer carried more ABT-SLC somatic mutations than expected from median tumour mutation burden. Among these, somatic mutations in SLC4A3 had features consistent with meaningful consequences on cellular fitness. Definitive evidence for ABT-SLCs as 'cancer essential' or 'driver genes' will have to consider microenvironmental context in genomic sequencing because bulk approaches are insensitive to pH heterogeneity within tumours. Moreover, genomic analyses must be validated with phenotypic outcomes (i.e. SLC-carried flux) to appreciate the opportunities for targeting acid-base transport in cancers.

Diagnostic Utility of Next-Generation Sequencing in Circulating Free DNA and a Comparison With Matched Tissue in Gallbladder Carcinoma.

Mutation detection for therapy monitoring in cell-free DNA (cfDNA) is used clinically for some malignancies. Gallbladder carcinoma (GBC) presents a diagnostic challenge and has limited late-stage treatment options. To our knowledge, this novel study examines, for the first time, genomic alterations in cfDNA from GBC to assess diagnostic accuracy and therapeutic options. The concordance of somatic genomic changes in cfDNA and DNA from paired tumor tissue was analyzed. Paired serum and tissue samples from 40 histologically proven GBC, 20 cholecystitis, and 4 normal (noninflamed gallbladder) controls were included. Targeted next-generation sequencing with a 22-gene panel (Colon and Lung Cancer Research Panel v2, Thermo Scientific) in cfDNA and tumor tissue with high depth and uniform coverage on ION Personal Genome Machine (ION, PGM) was performed. A spectrum of 223 mutations in cfDNA and 225 mutations in formalin-fixed paraffin-embedded tissue DNA were identified in 22 genes. Mutations ranged from 1 to 17 per case. In cfDNA frequent alterations were in TP53 (85.0%), EGFR (52.5%), MET (35%) CTNNB1, SMAD4, BRAF (32.5%), PTEN (30%), FGFR3 and PIK3CA (27.5%), NOTCH1 (25.0%), and FBXW7 and ERBB4 (22.5%). At least one clinically actionable mutation was identified in all cfDNA samples. Paired samples shared 149 of 225 genetic abnormalities (66.2%). Individual gene mutation concordance ranged from 44.44% to 82.0% and was highest for EGFR (82.0%), BRAF and NOTCH1 (80.0%), TP53 (73.08%), MET (72.22%), and ERBB4 (71.42%) with a significant level of correlation (Spearman r = 0.91, P ≤ .0001). The sensitivity and specificity of the TP53 gene at the gene level was the highest (94.44% and 100.0%, respectively). Overall survival was higher for ERBB4 and ERBB2 mutant tumors. The adenocarcinoma subtype revealed specific genetic changes in ERBB4, SMAD4, ERBB2, PTEN, KRAS, and NRAS. NGS-based cfDNA mutation profiling can be used to diagnose GBC before surgery to guide treatment decisions. Targeted therapy identified in GBC included SMAD4, ERBB2, ERBB4, EGFR, KRAS, BRAF, PIK3CA, MET, and NRAS.