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

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 prediction for 1) RHD zygosity and hybrid alleles, 2) RHCE*C versus 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 (SCCRIP). In a second validation cohort with 3030 cancer survivors (15.2% Black or African Americans) from the St. Jude Lifetime Cohort Study (SJLIFE), the optimized RHtyper reached concordance rates between WES and WGS predications to 96.3% for RHD, and 94.6% for RHCE. In conclusion, machine learning improved the accuracy of RH predication from 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.

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.

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 and resistance to InO. Pre- and post-InO 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, protein destabilization) and epitope alteration. Two CD22 mutant cases were post-InO hypermutators resulting from error-prone DNA damage repair (non-homologous/alternative end joining, mismatch repair deficiency), suggesting hypermutation drove escape from CD22-directed therapy. CD22-mutant relapses occurred after InO and subsequent hematopoietic stem cell transplantation (HSCT), suggesting InO eliminated 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 compromise of the G1/S DNA damage checkpoint as a mechanism of evading InO-induced apoptosis. Genome wide CRISPR/Cas9 screening in cell lines identified DNTT (TdT) 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 prior to HSCT. The identified mechanisms of escape from CD22-targeted therapy extend beyond antigen loss, and provide opportunities to improve therapeutic approaches and overcome resistance.

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.

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.

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.

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.

Methylation profiling reveals novel molecular classes of rhabdomyosarcoma.

Rhabdomyosarcomas (RMS) represent a family of aggressive soft tissue sarcomas that present in both children and adults. Pathologic risk stratification for RMS has been based on histologic subtype, with poor outcomes observed in alveolar rhabdomyosarcoma (ARMS) and the adult-type pleomorphic rhabdomyosarcoma (PRMS) compared to embryonal rhabdomyosarcoma (ERMS). Genomic sequencing studies have expanded the spectrum of RMS, with several new molecularly defined entities, including fusion-driven spindle cell/sclerosing rhabdomyosarcoma (SC/SRMS) and MYOD1-mutant SC/SRMS. Comprehensive genomic analysis has previously defined the mutational and copy number spectrum for the more common ERMS and ARMS and revealed corresponding methylation signatures. Comparatively, less is known about epigenetic correlates for the rare SC/SRMS or PRMS histologic subtypes. Herein, we present exome and RNA sequencing, copy number analysis, and methylation profiling of the largest cohort of molecularly characterized RMS samples to date. In addition to ARMS and ERMS, we identify two novel methylation subtypes, one having SC/SRMS histology and defined by MYOD1 p. L122R mutations and the other matching adult-type PRMS. Selected tumors from adolescent patients grouped with the PRMS methylation class, expanding the age range of these rare tumors. Limited follow-up data suggest that pediatric tumors with MYOD1-mutations are associated with an aggressive clinical course.

Antibiotic prophylaxis and the gastrointestinal resistome in paediatric patients with acute lymphoblastic leukaemia: a cohort study with metagenomic sequencing analysis.

BACKGROUND: Although antibiotic prophylaxis with levofloxacin can reduce the risk of serious infection in immunocompromised patients, the potential contribution of prophylaxis to antibiotic resistance is a major drawback. We aimed to identify the effects of levofloxacin prophylaxis, given to paediatric patients with acute lymphoblastic leukaemia to prevent infections during induction chemotherapy, on antibiotic resistance in gastrointestinal microbiota after completion of induction and consolidation therapy. METHODS: This prospective, single-centre (St Jude Children's Research Hospital, Memphis, TN, USA) cohort study included children (≤18 years) receiving therapy for newly diagnosed acute lymphoblastic leukaemia and who received either primary levofloxacin prophylaxis or no antibacterial prophylaxis (aside from Pneumocystis jirovecii prophylaxis with trimethoprim-sulfamethoxazole) and provided at least two stool samples, including one after completion of induction therapy. We used metagenomic sequencing to identify bacterial genes that confer resistance to fluoroquinolones, trimethoprim-sulfamethoxazole, or other antibiotics, and to identify point mutations in bacterial topoisomerases (gyrA, parC) that confer resistance to fluoroquinolones. We then used generalised linear mixed models to compare the prevalence and relative abundance of antibiotic resistance gene groups after completion of induction and consolidation therapy between participants who had received levofloxacin and those who received no prophylaxis. FINDINGS: Between Feb 1, 2012, and April 30, 2016, 118 stool samples (32 baseline, 49 after induction, and 37 after consolidation) were collected from 49 evaluable participants; of these participants, 31 (63%) received levofloxacin prophylaxis during induction therapy and 18 (37%) received no antibacterial prophylaxis. Over the course of induction therapy, there was an overall increase in the relative abundance of trimethoprim-sulfamethoxazole resistance genes (estimated mean fold change 5·9, 95% CI 3·6-9·6; p<0·0001), which was not modified by levofloxacin prophylaxis (p=0·46). By contrast, the prevalence of topoisomerase point mutations increased over the course of induction therapy in levofloxacin recipients (mean prevalence 10·4% [95% CI 3·2-25·4] after induction therapy vs 3·7% [0·2-22·5] at baseline) but not other participants (0% vs 0%; p<0·0001). There was no significant difference between prophylaxis groups with respect to changes in aminoglycoside, ß-lactam, vancomycin, or multidrug resistance genes after completion of induction or consolidation therapy. INTERPRETATION: Analysing the gastrointestinal resistome can provide insights into the effects of antibiotics on the risk of antibiotic-resistant infections. In this study, antibiotic prophylaxis with trimethoprim-sulfamethoxazole or levofloxacin during induction therapy for acute lymphoblastic leukaemia appeared to increase the short-term and medium-term risk of colonisation with bacteria resistant to these antibiotics, but not to other drugs. More research is needed to determine the longer-term effects of antibacterial prophylaxis on colonisation with antibiotic-resistant bacteria. FUNDING: Children's Infection Defense Center at St Jude Children's Research Hospital, American Lebanese Syrian Associated Charities, and National Institutes of Health.