Genomic Determinants of Outcome in Acute Lymphoblastic Leukemia.

PURPOSE: Although cure rates for childhood acute lymphoblastic leukemia (ALL) exceed 90%, ALL remains a leading cause of cancer death in children. Half of relapses arise in children initially classified with standard-risk (SR) disease. MATERIALS AND METHODS: To identify genomic determinants of relapse in children with SR ALL, we performed genome and transcriptome sequencing of diagnostic and remission samples of children with SR (n = 1,381) or high-risk B-ALL with favorable cytogenetic features (n = 115) enrolled on Children's Oncology Group trials. We used a case-control study design analyzing 439 patients who relapsed and 1,057 who remained in complete remission for at least 5 years. RESULTS: Genomic subtype was associated with relapse, which occurred in approximately 50% of cases of PAX5-altered ALL (odds ratio [OR], 3.31 [95% CI, 2.17 to 5.03]; P = 3.18 × 10-8). Within high-hyperdiploid ALL, gain of chromosome 10 with disomy of chromosome 7 was associated with favorable outcome (OR, 0.27 [95% CI, 0.17 to 0.42]; P = 8.02 × 10-10; St Jude Children's Research Hospital validation cohort: OR, 0.22 [95% CI, 0.05 to 0.80]; P = .009), and disomy of chromosomes 10 and 17 with gain of chromosome 6 was associated with relapse (OR, 7.16 [95% CI, 2.63 to 21.51]; P = 2.19 × 10-5; validation cohort: OR, 21.32 [95% CI, 3.62 to 119.30]; P = .0004). Genomic alterations were associated with relapse in a subtype-dependent manner, including alterations of INO80 in ETV6::RUNX1 ALL, IKZF1, and CREBBP in high-hyperdiploid ALL and FHIT in BCR::ABL1-like ALL. Genomic alterations were also associated with the presence of minimal residual disease, including NRAS and CREBBP in high-hyperdiploid ALL. CONCLUSION: Genetic subtype, patterns of aneuploidy, and secondary genomic alterations determine risk of relapse in childhood ALL. Comprehensive genomic analysis is required for optimal risk stratification.

Blinatumomab for MRD-Negative Acute Lymphoblastic Leukemia in Adults.

BACKGROUND: Many older adults with B-cell precursor acute lymphoblastic leukemia (BCP-ALL) have a relapse despite having a measurable residual disease (MRD)-negative complete remission with combination chemotherapy. The addition of blinatumomab, a bispecific T-cell engager molecule that is approved for the treatment of relapsed, refractory, and MRD-positive BCP-ALL, may have efficacy in patients with MRD-negative remission. METHODS: In a phase 3 trial, we randomly assigned patients 30 to 70 years of age with BCR::ABL1-negative BCP-ALL (with :: indicating fusion) who had MRD-negative remission (defined as <0.01% leukemic cells in bone marrow as assessed on flow cytometry) after induction and intensification chemotherapy to receive four cycles of blinatumomab in addition to four cycles of consolidation chemotherapy or to receive four cycles of consolidation chemotherapy alone. The primary end point was overall survival, and relapse-free survival was a secondary end point. RESULTS: The data and safety monitoring committee reviewed the results from the third efficacy interim analysis and recommended that they be reported. Complete remission with or without full count recovery was observed in 395 of 488 enrolled patients (81%). Of the 224 patients with MRD-negative status, 112 were assigned to each group. The characteristics of the patients were balanced between the groups. At a median follow-up of 43 months, an advantage was observed in the blinatumomab group as compared with the chemotherapy-only group with regard to overall survival (at 3 years: 85% vs. 68%; hazard ratio for death, 0.41; 95% confidence interval [CI], 0.23 to 0.73; P = 0.002), and the 3-year relapse-free survival was 80% with blinatumomab and 64% with chemotherapy alone (hazard ratio for relapse or death, 0.53; 95% CI, 0.32 to 0.87). A higher incidence of neuropsychiatric events was reported in the blinatumomab group than in the chemotherapy-only group. CONCLUSIONS: The addition of blinatumomab to consolidation chemotherapy in adult patients in MRD-negative remission from BCP-ALL significantly improved overall survival. (Funded by the National Institutes of Health and others; E1910 ClinicalTrials.gov number, NCT02003222.).

Machine learning to optimize automated RH genotyping using whole-exome sequencing data.

ABSTRACT: Rh phenotype matching reduces but does not eliminate alloimmunization in patients with sickle cell disease (SCD) due to RH genetic diversity that is not distinguishable by serological typing. RH genotype matching can potentially mitigate Rh alloimmunization but comprehensive and accessible genotyping methods are needed. We developed RHtyper as an automated algorithm to predict RH genotypes using whole-genome sequencing (WGS) data with high accuracy. Here, we adapted RHtyper for whole-exome sequencing (WES) data, which are more affordable but challenged by uneven sequencing coverage and exacerbated sequencing read misalignment, resulting in uncertain predictions for (1) RHD zygosity and hybrid alleles, (2) RHCE∗C vs. RHCE∗c alleles, (3) RHD c.1136C>T zygosity, and (4) RHCE c.48G>C zygosity. We optimized RHtyper to accurately predict RHD and RHCE genotypes using WES data by leveraging machine learning models and improved the concordance of WES with WGS predictions from 90.8% to 97.2% for RHD and 96.3% to 98.2% for RHCE among 396 patients in the Sickle Cell Clinical Research and Intervention Program. In a second validation cohort of 3030 cancer survivors (15.2% Black or African Americans) from the St. Jude Lifetime Cohort Study, the optimized RHtyper reached concordance rates between WES and WGS predications to 96.3% for RHD and 94.6% for RHCE. Machine learning improved the accuracy of RH predication using WES data. RHtyper has the potential, once implemented, to provide a precision medicine-based approach to facilitate RH genotype-matched transfusion and improve transfusion safety for patients with SCD. This study used data from clinical trials registered at ClinicalTrials.gov as #NCT02098863 and NCT00760656.

Genomic determinants of response and resistance to inotuzumab ozogamicin in B-cell ALL.

ABSTRACT: Inotuzumab ozogamicin (InO) is an antibody-drug conjugate that delivers calicheamicin to CD22-expressing cells. In a retrospective cohort of InO-treated patients with B-cell acute lymphoblastic leukemia, we sought to understand the genomic determinants of the response and resistance to InO. Pre- and post-InO-treated patient samples were analyzed by whole genome, exome, and/or transcriptome sequencing. Acquired CD22 mutations were observed in 11% (3/27) of post-InO-relapsed tumor samples, but not in refractory samples (0/16). There were multiple CD22 mutations per sample and the mechanisms of CD22 escape included epitope loss (protein truncation and destabilization) and epitope alteration. Two CD22 mutant cases were post-InO hyper-mutators resulting from error-prone DNA damage repair (nonhomologous/alternative end-joining repair, or mismatch repair deficiency), suggesting that hypermutation drove escape from CD22-directed therapy. CD22-mutant relapses occurred after InO and subsequent hematopoietic stem cell transplantation (HSCT), suggesting that InO eliminated the predominant clones, leaving subclones with acquired CD22 mutations that conferred resistance to InO and subsequently expanded. Acquired loss-of-function mutations in TP53, ATM, and CDKN2A were observed, consistent with a compromise of the G1/S DNA damage checkpoint as a mechanism for evading InO-induced apoptosis. Genome-wide CRISPR/Cas9 screening of cell lines identified DNTT (terminal deoxynucleotidyl transferase) loss as a marker of InO resistance. In conclusion, genetic alterations modulating CD22 expression and DNA damage response influence InO efficacy. Our findings highlight the importance of defining the basis of CD22 escape and eradication of residual disease before HSCT. The identified mechanisms of escape from CD22-targeted therapy extend beyond antigen loss and provide opportunities to improve therapeutic approaches and overcome resistance. These trials were registered at www.ClinicalTrials.gov as NCT01134575, NCT01371630, and NCT03441061.

Interaction between host G3BP and viral nucleocapsid protein regulates SARS-CoV-2 replication and pathogenicity.

G3BP1/2 are paralogous proteins that promote stress granule formation in response to cellular stresses, including viral infection. The nucleocapsid (N) protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) inhibits stress granule assembly and interacts with G3BP1/2 via an ITFG motif, including residue F17, in the N protein. Prior studies examining the impact of the G3PB1-N interaction on SARS-CoV-2 replication have produced inconsistent findings, and the role of this interaction in pathogenesis is unknown. Here, we use structural and biochemical analyses to define the residues required for G3BP1-N interaction and structure-guided mutagenesis to selectively disrupt this interaction. We find that N-F17A mutation causes highly specific loss of interaction with G3BP1/2. SARS-CoV-2 N-F17A fails to inhibit stress granule assembly in cells, has decreased viral replication, and causes decreased pathology in vivo. Further mechanistic studies indicate that the N-F17-mediated G3BP1-N interaction promotes infection by limiting sequestration of viral genomic RNA (gRNA) into stress granules.

Metabolic reprogramming of cancer cells by JMJD6-mediated pre-mRNA splicing associated with therapeutic response to splicing inhibitor.

Dysregulated pre-mRNA splicing and metabolism are two hallmarks of MYC-driven cancers. Pharmacological inhibition of both processes has been extensively investigated as potential therapeutic avenues in preclinical and clinical studies. However, how pre-mRNA splicing and metabolism are orchestrated in response to oncogenic stress and therapies is poorly understood. Here, we demonstrate that jumonji domain containing 6, arginine demethylase, and lysine hydroxylase, JMJD6, acts as a hub connecting splicing and metabolism in MYC-driven human neuroblastoma. JMJD6 cooperates with MYC in cellular transformation of murine neural crest cells by physically interacting with RNA binding proteins involved in pre-mRNA splicing and protein homeostasis. Notably, JMJD6 controls the alternative splicing of two isoforms of glutaminase (GLS), namely kidney-type glutaminase (KGA) and glutaminase C (GAC), which are rate-limiting enzymes of glutaminolysis in the central carbon metabolism in neuroblastoma. Further, we show that JMJD6 is correlated with the anti-cancer activity of indisulam, a 'molecular glue' that degrades splicing factor RBM39, which complexes with JMJD6. The indisulam-mediated cancer cell killing is at least partly dependent on the glutamine-related metabolic pathway mediated by JMJD6. Our findings reveal a cancer-promoting metabolic program is associated with alternative pre-mRNA splicing through JMJD6, providing a rationale to target JMJD6 as a therapeutic avenue for treating MYC-driven cancers.

A spatial cell atlas of neuroblastoma reveals developmental, epigenetic and spatial axis of tumor heterogeneity.

Neuroblastoma is a pediatric cancer arising from the developing sympathoadrenal lineage with complex inter- and intra-tumoral heterogeneity. To chart this complexity, we generated a comprehensive cell atlas of 55 neuroblastoma patient tumors, collected from two pediatric cancer institutions, spanning a range of clinical, genetic, and histologic features. Our atlas combines single-cell/nucleus RNA-seq (sc/scRNA-seq), bulk RNA-seq, whole exome sequencing, DNA methylation profiling, spatial transcriptomics, and two spatial proteomic methods. Sc/snRNA-seq revealed three malignant cell states with features of sympathoadrenal lineage development. All of the neuroblastomas had malignant cells that resembled sympathoblasts and the more differentiated adrenergic cells. A subset of tumors had malignant cells in a mesenchymal cell state with molecular features of Schwann cell precursors. DNA methylation profiles defined four groupings of patients, which differ in the degree of malignant cell heterogeneity and clinical outcomes. Using spatial proteomics, we found that neuroblastomas are spatially compartmentalized, with malignant tumor cells sequestered away from immune cells. Finally, we identify spatially restricted signaling patterns in immune cells from spatial transcriptomics. To facilitate the visualization and analysis of our atlas as a resource for further research in neuroblastoma, single cell, and spatial-omics, all data are shared through the Human Tumor Atlas Network Data Commons at www.humantumoratlas.org.

Acute myeloid leukemias with UBTF tandem duplications are sensitive to menin inhibitors.

ABSTRACT: UBTF tandem duplications (UBTF-TDs) have recently emerged as a recurrent alteration in pediatric and adult acute myeloid leukemia (AML). UBTF-TD leukemias are characterized by a poor response to conventional chemotherapy and a transcriptional signature that mirrors NUP98-rearranged and NPM1-mutant AMLs, including HOX-gene dysregulation. However, the mechanism by which UBTF-TD drives leukemogenesis remains unknown. In this study, we investigated the genomic occupancy of UBTF-TD in transformed cord blood CD34+ cells and patient-derived xenograft models. We found that UBTF-TD protein maintained genomic occupancy at ribosomal DNA loci while also occupying genomic targets commonly dysregulated in UBTF-TD myeloid malignancies, such as the HOXA/HOXB gene clusters and MEIS1. These data suggest that UBTF-TD is a gain-of-function alteration that results in mislocalization to genomic loci dysregulated in UBTF-TD leukemias. UBTF-TD also co-occupies key genomic loci with KMT2A and menin, which are known to be key partners involved in HOX-dysregulated leukemias. Using a protein degradation system, we showed that stemness, proliferation, and transcriptional signatures are dependent on sustained UBTF-TD localization to chromatin. Finally, we demonstrate that primary cells from UBTF-TD leukemias are sensitive to the menin inhibitor SNDX-5613, resulting in markedly reduced in vitro and in vivo tumor growth, myeloid differentiation, and abrogation of the UBTF-TD leukemic expression signature. These findings provide a viable therapeutic strategy for patients with this high-risk AML subtype.

Genomic determinants of response and resistance to inotuzumab ozogamicin in B-cell ALL.

Inotuzumab ozogamicin (InO) is an antibody-drug conjugate that delivers calicheamicin to CD22-expressing cells. In a retrospective cohort of InO treated patients with B-cell acute lymphoblastic leukemia, we sought to understand the genomic determinants of response to InO. Acquired CD22 mutations were observed in 11% (3/27) of post-InO relapsed tumor samples. There were multiple CD22 mutations per sample and the mechanisms of CD22 escape included protein truncation, protein destabilization, and epitope alteration. Hypermutation by error-prone DNA damage repair (alternative end-joining, mismatch repair deficiency) drove CD22 escape. Acquired loss-of-function mutations in TP53 , ATM and CDKN2A were observed, suggesting compromise of the G1/S DNA damage checkpoint as a mechanism of evading InO-induced apoptosis. In conclusion, genetic alterations modulating CD22 expression and DNA damage response influence InO efficacy. The escape strategies within and beyond antigen loss to CD22-targeted therapy elucidated in this study provide insights into improving therapeutic approaches and overcoming resistance. KEY POINTS: We identified multiple mechanisms of CD22 antigen escape from inotuzumab ozogamicin, including protein truncation, protein destabilization, and epitope alteration.Hypermutation caused by error-prone DNA damage repair was a driver of CD22 mutation and escape.

Identification and targeting of treatment resistant progenitor populations in T-cell Acute Lymphoblastic Leukemia.

Refractoriness to initial chemotherapy and relapse after remission are the main obstacles to cure in T-cell Acute Lymphoblastic Leukemia (T-ALL). Biomarker guided risk stratification and targeted therapy have the potential to improve outcomes in high-risk T-ALL; however, cellular and genetic factors contributing to treatment resistance remain unknown. Previous bulk genomic studies in T-ALL have implicated tumor heterogeneity as an unexplored mechanism for treatment failure. To link tumor subpopulations with clinical outcome, we created an atlas of healthy pediatric hematopoiesis and applied single-cell multiomic (CITE-seq/snATAC-seq) analysis to a cohort of 40 cases of T-ALL treated on the Children's Oncology Group AALL0434 clinical trial. The cohort was carefully selected to capture the immunophenotypic diversity of T-ALL, with early T-cell precursor (ETP) and Near/Non-ETP subtypes represented, as well as enriched with both relapsed and treatment refractory cases. Integrated analyses of T-ALL blasts and normal T-cell precursors identified a bone-marrow progenitor-like (BMP-like) leukemia sub-population associated with treatment failure and poor overall survival. The single-cell-derived molecular signature of BMP-like blasts predicted poor outcome across multiple subtypes of T-ALL within two independent patient cohorts using bulk RNA-sequencing data from over 1300 patients. We defined the mutational landscape of BMP-like T-ALL, finding that NOTCH1 mutations additively drive T-ALL blasts away from the BMP-like state. We transcriptionally matched BMP-like blasts to early thymic seeding progenitors that have low NR3C1 expression and high stem cell gene expression, corresponding to a corticosteroid and conventional cytotoxic resistant phenotype we observed in ex vivo drug screening. To identify novel targets for BMP-like blasts, we performed in silico and in vitro drug screening against the BMP-like signature and prioritized BMP-like overexpressed cell-surface (CD44, ITGA4, LGALS1) and intracellular proteins (BCL-2, MCL-1, BTK, NF-κB) as candidates for precision targeted therapy. We established patient derived xenograft models of BMP-high and BMP-low leukemias, which revealed vulnerability of BMP-like blasts to apoptosis-inducing agents, TEC-kinase inhibitors, and proteasome inhibitors. Our study establishes the first multi-omic signatures for rapid risk-stratification and targeted treatment of high-risk T-ALL.

Interaction between host G3BP and viral nucleocapsid protein regulates SARS-CoV-2 replication.

G3BP1/2 are paralogous proteins that promote stress granule formation in response to cellular stresses, including viral infection. G3BP1/2 are prominent interactors of the nucleocapsid (N) protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). However, the functional consequences of the G3BP1-N interaction in the context of viral infection remain unclear. Here we used structural and biochemical analyses to define the residues required for G3BP1-N interaction, followed by structure-guided mutagenesis of G3BP1 and N to selectively and reciprocally disrupt their interaction. We found that mutation of F17 within the N protein led to selective loss of interaction with G3BP1 and consequent failure of the N protein to disrupt stress granule assembly. Introduction of SARS-CoV-2 bearing an F17A mutation resulted in a significant decrease in viral replication and pathogenesis in vivo, indicating that the G3BP1-N interaction promotes infection by suppressing the ability of G3BP1 to form stress granules.

Metabolic reprogramming of cancer cells by JMJD6-mediated pre-mRNA splicing is associated with therapeutic response to splicing inhibitor.

Dysregulated pre-mRNA splicing and metabolism are two hallmarks of MYC-driven cancers. Pharmacological inhibition of both processes has been extensively investigated as potential therapeutic avenues in preclinical and clinical studies. However, how pre-mRNA splicing and metabolism are orchestrated in response to oncogenic stress and therapies is poorly understood. Here, we demonstrate that Jumonji Domain Containing 6, Arginine Demethylase and Lysine Hydroxylase, JMJD6, acts as a hub connecting splicing and metabolism in MYC-driven neuroblastoma. JMJD6 cooperates with MYC in cellular transformation by physically interacting with RNA binding proteins involved in pre-mRNA splicing and protein homeostasis. Notably, JMJD6 controls the alternative splicing of two isoforms of glutaminase (GLS), namely kidney-type glutaminase (KGA) and glutaminase C (GAC), which are rate-limiting enzymes of glutaminolysis in the central carbon metabolism in neuroblastoma. Further, we show that JMJD6 is correlated with the anti-cancer activity of indisulam, a "molecular glue" that degrades splicing factor RBM39, which complexes with JMJD6. The indisulam-mediated cancer cell killing is at least partly dependent on the glutamine-related metabolic pathway mediated by JMJD6. Our findings reveal a cancer-promoting metabolic program is associated with alternative pre-mRNA splicing through JMJD6, providing a rationale to target JMJD6 as a therapeutic avenue for treating MYC-driven cancers.

The genomic landscape of acute lymphoblastic leukemia with intrachromosomal amplification of chromosome 21.

Intrachromosomal amplification of chromosome 21 defines a subtype of high-risk childhood acute lymphoblastic leukemia (iAMP21-ALL) characterized by copy number changes and complex rearrangements of chromosome 21. The genomic basis of iAMP21-ALL and the pathogenic role of the region of amplification of chromosome 21 to leukemogenesis remains incompletely understood. In this study, using integrated whole genome and transcriptome sequencing of 124 patients with iAMP21-ALL, including rare cases arising in the context of constitutional chromosomal aberrations, we identified subgroups of iAMP21-ALL based on the patterns of copy number alteration and structural variation. This large data set enabled formal delineation of a 7.8 Mb common region of amplification harboring 71 genes, 43 of which were differentially expressed compared with non-iAMP21-ALL ones, including multiple genes implicated in the pathogenesis of acute leukemia (CHAF1B, DYRK1A, ERG, HMGN1, and RUNX1). Using multimodal single-cell genomic profiling, including single-cell whole genome sequencing of 2 cases, we documented clonal heterogeneity and genomic evolution, demonstrating that the acquisition of the iAMP21 chromosome is an early event that may undergo progressive amplification during disease ontogeny. We show that UV-mutational signatures and high mutation load are characteristic secondary genetic features. Although the genomic alterations of chromosome 21 are variable, these integrated genomic analyses and demonstration of an extended common minimal region of amplification broaden the definition of iAMP21-ALL for more precise diagnosis using cytogenetic or genomic methods to inform clinical management.

Genomic landscape of Down syndrome-associated acute lymphoblastic leukemia.

Trisomy 21, the genetic cause of Down syndrome (DS), is the most common congenital chromosomal anomaly. It is associated with a 20-fold increased risk of acute lymphoblastic leukemia (ALL) during childhood and results in distinctive leukemia biology. To comprehensively define the genomic landscape of DS-ALL, we performed whole-genome sequencing and whole-transcriptome sequencing (RNA-Seq) on 295 cases. Our integrated genomic analyses identified 15 molecular subtypes of DS-ALL, with marked enrichment of CRLF2-r, IGH::IGF2BP1, and C/EBP altered (C/EBPalt) subtypes compared with 2257 non-DS-ALL cases. We observed abnormal activation of the CEBPD, CEBPA, and CEBPE genes in 10.5% of DS-ALL cases via a variety of genomic mechanisms, including chromosomal rearrangements and noncoding mutations leading to enhancer hijacking. A total of 42.3% of C/EBP-activated DS-ALL also have concomitant FLT3 point mutations or insertions/deletions, compared with 4.1% in other subtypes. CEBPD overexpression enhanced the differentiation of mouse hematopoietic progenitor cells into pro-B cells in vitro, particularly in a DS genetic background. Notably, recombination-activating gene-mediated somatic genomic abnormalities were common in DS-ALL, accounting for a median of 27.5% of structural alterations, compared with 7.7% in non-DS-ALL. Unsupervised hierarchical clustering analyses of CRLF2-rearranged DS-ALL identified substantial heterogeneity within this group, with the BCR::ABL1-like subset linked to an inferior event-free survival, even after adjusting for known clinical risk factors. These results provide important insights into the biology of DS-ALL and point to opportunities for targeted therapy and treatment individualization.

Discovery of novel predisposing coding and noncoding variants in familial Hodgkin lymphoma.

Familial aggregation of Hodgkin lymphoma (HL) has been demonstrated in large population studies, pointing to genetic predisposition to this hematological malignancy. To understand the genetic variants associated with the development of HL, we performed whole genome sequencing on 234 individuals with and without HL from 36 pedigrees that had 2 or more first-degree relatives with HL. Our pedigree selection criteria also required at least 1 affected individual aged <21 years, with the median age at diagnosis of 21.98 years (3-55 years). Family-based segregation analysis was performed for the identification of coding and noncoding variants using linkage and filtering approaches. Using our tiered variant prioritization algorithm, we identified 44 HL-risk variants in 28 pedigrees, of which 33 are coding and 11 are noncoding. The top 4 recurrent risk variants are a coding variant in KDR (rs56302315), a 5' untranslated region variant in KLHDC8B (rs387906223), a noncoding variant in an intron of PAX5 (rs147081110), and another noncoding variant in an intron of GATA3 (rs3824666). A newly identified splice variant in KDR (c.3849-2A>C) was observed for 1 pedigree, and high-confidence stop-gain variants affecting IRF7 (p.W238∗) and EEF2KMT (p.K116∗) were also observed. Multiple truncating variants in POLR1E were found in 3 independent pedigrees as well. Whereas KDR and KLHDC8B have previously been reported, PAX5, GATA3, IRF7, EEF2KMT, and POLR1E represent novel observations. Although there may be environmental factors influencing lymphomagenesis, we observed segregation of candidate germline variants likely to predispose HL in most of the pedigrees studied.

Innate frequency-discrimination hyperacuity in Williams-Beuren syndrome mice.

Williams-Beuren syndrome (WBS) is a rare disorder caused by hemizygous microdeletion of ∼27 contiguous genes. Despite neurodevelopmental and cognitive deficits, individuals with WBS have spared or enhanced musical and auditory abilities, potentially offering an insight into the genetic basis of auditory perception. Here, we report that the mouse models of WBS have innately enhanced frequency-discrimination acuity and improved frequency coding in the auditory cortex (ACx). Chemogenetic rescue showed frequency-discrimination hyperacuity is caused by hyperexcitable interneurons in the ACx. Haploinsufficiency of one WBS gene, Gtf2ird1, replicated WBS phenotypes by downregulating the neuropeptide receptor VIPR1. VIPR1 is reduced in the ACx of individuals with WBS and in the cerebral organoids derived from human induced pluripotent stem cells with the WBS microdeletion. Vipr1 deletion or overexpression in ACx interneurons mimicked or reversed, respectively, the cellular and behavioral phenotypes of WBS mice. Thus, the Gtf2ird1-Vipr1 mechanism in ACx interneurons may underlie the superior auditory acuity in WBS.

Targeted Long-Read Bisulfite Sequencing Identifies Differences in the TERT Promoter Methylation Profiles between TERT Wild-Type and TERT Mutant Cancer Cells.

Background: TERT promoter methylation, located several hundred base pairs upstream of the transcriptional start site, is cancer specific and correlates with increased TERT mRNA expression and poorer patient outcome. Promoter methylation, however, is not mutually exclusive to TERT activating genetic alterations, as predicted for functionally redundant mechanisms. To annotate the altered patterns of TERT promoter methylation and their relationship with gene expression, we applied a Pacific Biosciences-based, long-read, bisulfite-sequencing technology and compared the differences in the methylation marks between wild-type and mutant cancers in an allele-specific manner. Results: We cataloged TERT genetic alterations (i.e., promoter point mutations or structural variations), allele-specific promoter methylation patterns, and allele-specific expression levels in a cohort of 54 cancer cell lines. In heterozygous mutant cell lines, the mutant alleles were significantly less methylated than their silent, mutation-free alleles (p < 0.05). In wild-type cell lines, by contrast, both epialleles were equally methylated to high levels at the TERT distal promoter, but differentially methylated in the proximal regions. ChIP analysis showed that epialleles with the hypomethylated proximal and core promoter were enriched in the active histone mark H3K4me2/3, whereas epialleles that were methylated in those regions were enriched in the repressive histone mark H3K27me3. Decitabine therapy induced biallelic expression in the wild-type cancer cells, whereas the mutant cell lines were unaffected. Conclusions: Long-read bisulfite sequencing analysis revealed differences in the methylation profiles and responses to demethylating agents between TERT wild-type and genetically altered cancer cell lines. The causal relation between TERT promoter methylation and gene expression remains to be established.

Somatic and Germline Variant Calling from Next-Generation Sequencing Data.

Re-sequencing of the human genome by next-generation sequencing (NGS) has been widely applied to discover pathogenic genetic variants and/or causative genes accounting for various types of diseases including cancers. The advances in NGS have allowed the sequencing of the entire genome of patients and identification of disease-associated variants in a reasonable timeframe and cost. The core of the variant identification relies on accurate variant calling and annotation. Numerous algorithms have been developed to elucidate the repertoire of somatic and germline variants. Each algorithm has its own distinct strengths, weaknesses, and limitations due to the difference in the statistical modeling approach adopted and read information utilized. Accurate variant calling remains challenging due to the presence of sequencing artifacts and read misalignments. All of these can lead to the discordance of the variant calling results and even misinterpretation of the discovery. For somatic variant detection, multiple factors including chromosomal abnormalities, tumor heterogeneity, tumor-normal cross contaminations, unbalanced tumor/normal sample coverage, and variants with low allele frequencies add even more layers of complexity to accurate variant identification. Given the discordances and difficulties, ensemble approaches have emerged by harmonizing information from different algorithms to improve variant calling performance. In this chapter, we first introduce the general scheme of variant calling algorithms and potential challenges at distinct stages. We next review the existing workflows of variant calling and annotation, and finally explore the strategies deployed by different callers as well as their strengths and caveats. Overall, NGS-based variant identification with careful consideration allows reliable detection of pathogenic variant and candidate variant selection for precision medicine.

Induction of broadly reactive influenza antibodies increases susceptibility to autoimmunity.

Infection and vaccination repeatedly expose individuals to antigens that are conserved between influenza virus subtypes. Nevertheless, antibodies recognizing variable influenza epitopes greatly outnumber antibodies reactive against conserved epitopes. Elucidating factors contributing to the paucity of broadly reactive influenza antibodies remains a major obstacle for developing a universal influenza vaccine. Here, we report that inducing broadly reactive influenza antibodies increases autoreactive antibodies in humans and mice and exacerbates disease in four distinct models of autoimmune disease. Importantly, transferring broadly reactive influenza antibodies augments disease in the presence of inflammation or autoimmune susceptibility. Further, broadly reactive influenza antibodies spontaneously arise in mice with defects in B cell tolerance. Together, these data suggest that self-tolerance mechanisms limit the prevalence of broadly reactive influenza antibodies, which can exacerbate disease in the context of additional risk factors.

Gain-of-function mutations in RPA1 cause a syndrome with short telomeres and somatic genetic rescue.

Human telomere biology disorders (TBD)/short telomere syndromes (STS) are heterogeneous disorders caused by inherited loss-of-function mutations in telomere-associated genes. Here, we identify 3 germline heterozygous missense variants in the RPA1 gene in 4 unrelated probands presenting with short telomeres and varying clinical features of TBD/STS, including bone marrow failure, myelodysplastic syndrome, T- and B-cell lymphopenia, pulmonary fibrosis, or skin manifestations. All variants cluster to DNA-binding domain A of RPA1 protein. RPA1 is a single-strand DNA-binding protein required for DNA replication and repair and involved in telomere maintenance. We showed that RPA1E240K and RPA1V227A proteins exhibit increased binding to single-strand and telomeric DNA, implying a gain in DNA-binding function, whereas RPA1T270A has binding properties similar to wild-type protein. To study the mutational effect in a cellular system, CRISPR/Cas9 was used to knock-in the RPA1E240K mutation into healthy inducible pluripotent stem cells. This resulted in severe telomere shortening and impaired hematopoietic differentiation. Furthermore, in patients with RPA1E240K, we discovered somatic genetic rescue in hematopoietic cells due to an acquired truncating cis RPA1 mutation or a uniparental isodisomy 17p with loss of mutant allele, coinciding with stabilized blood counts. Using single-cell sequencing, the 2 somatic genetic rescue events were proven to be independently acquired in hematopoietic stem cells. In summary, we describe the first human disease caused by germline RPA1 variants in individuals with TBD/STS.