Environmental microbiome in the home and daycare settings during the COVID-19 pandemic, and potential risk of non-communicable disease in children.

An exposure to diverse microbial population early in life is important for the development of immunity against various non-communicable diseases including asthma, childhood leukaemia and other cancers. Social mixing in daycare settings helps with exposure to a variety of microbes. However, social isolation and a high emphasis on workplace hygiene during the COVID pandemic may have affected children's exposure to diverse microbiota. The structure of microbial communities and their role in developing immunity to various diseases are not well understood. In this study, we investigated the structure of microbial communities in daycare and home settings during the pandemic. Interestingly, microbial diversity was relatively higher in dust samples collected from homes, with human-associated taxa being more prevalent compared to those from daycare settings. Environmental microbes were more abundant in dust samples from daycare providers. These results potentially suggest that cleaning practices during the pandemic may have influenced the diversity and microbial abundance of the daycare samples. Several bacterial taxa detected in both the environments are known to induce anti-inflammatory and immunomodulatory responses, conferring protection from various diseases. Therefore, exposure to diverse microbial population in early childhood may play an important role in developing immunity against various non-communicable and infectious diseases.

Single-cell multi-omics identifies chronic inflammation as a driver of TP53-mutant leukemic evolution.

Understanding the genetic and nongenetic determinants of tumor protein 53 (TP53)-mutation-driven clonal evolution and subsequent transformation is a crucial step toward the design of rational therapeutic strategies. Here we carry out allelic resolution single-cell multi-omic analysis of hematopoietic stem/progenitor cells (HSPCs) from patients with a myeloproliferative neoplasm who transform to TP53-mutant secondary acute myeloid leukemia (sAML). All patients showed dominant TP53 'multihit' HSPC clones at transformation, with a leukemia stem cell transcriptional signature strongly predictive of adverse outcomes in independent cohorts, across both TP53-mutant and wild-type (WT) AML. Through analysis of serial samples, antecedent TP53-heterozygous clones and in vivo perturbations, we demonstrate a hitherto unrecognized effect of chronic inflammation, which suppressed TP53 WT HSPCs while enhancing the fitness advantage of TP53-mutant cells and promoted genetic evolution. Our findings will facilitate the development of risk-stratification, early detection and treatment strategies for TP53-mutant leukemia, and are of broad relevance to other cancer types.

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.

Molecular characterization and clinical outcome of B-cell precursor acute lymphoblastic leukemia with IG-MYC rearrangement.

Rarely, immunophenotypically immature B-cell precursor acute lymphoblastic leukemia (BCP-ALL) carries an immunoglobulin- MYC rearrangement (IG-MYC-r). This can result in diagnostic confusion with Burkitt lymphoma/leukemia and use of individualized treatment schedules of unproven efficacy. Here we compare the molecular characteristics of these conditions and investigate historic clinical outcome data. We identified 90 cases registered in a national BCP-ALL clinical trial/registry. When present, diagnostic material underwent cytogenetic, exome, methylome and transcriptome analyses. The outcomes analyzed were 3-year event-free survival and overall survival. IG-MYC-r was identified in diverse cytogenetic backgrounds, co-existing with either established BCP-ALL-specific abnormalities (high hyperdiploidy, n=3; KMT2A-rearrangement, n=6; iAMP21, n=1; BCR-ABL1, n=1); BCL2/BCL6-rearrangements (n=15); or, most commonly, as the only defining feature (n=64). Within this final group, precursor-like V(D)J breakpoints predominated (8/9) and KRAS mutations were common (5/11). DNA methylation identified a cluster of V(D)J-rearranged cases, clearly distinct from Burkitt leukemia/lymphoma. Children with IG-MYC-r within that subgroup had a 3-year event-free survival of 47% and overall survival of 60%, representing a high-risk BCP-ALL. To develop effective management strategies this group of patients must be allowed access to contemporary, minimal residual disease-adapted, prospective clinical trial protocols.

Integrative genomic analysis of childhood acute lymphoblastic leukaemia lacking a genetic biomarker in the UKALL2003 clinical trial.

Incorporating genetics into risk-stratification for treatment of childhood B-progenitor acute lymphoblastic leukaemia (B-ALL) has contributed significantly to improved survival. In about 30% B-ALL (B-other-ALL) without well-established chromosomal changes, new genetic subtypes have recently emerged, yet their true prognostic relevance largely remains unclear. We integrated next generation sequencing (NGS): whole genome sequencing (WGS) (n = 157) and bespoke targeted NGS (t-NGS) (n = 175) (overlap n = 36), with existing genetic annotation in a representative cohort of 351 B-other-ALL patients from the childhood ALL trail, UKALL2003. PAX5alt was most frequently observed (n = 91), whereas PAX5 P80R mutations (n = 11) defined a distinct PAX5 subtype. DUX4-r subtype (n = 80) was defined by DUX4 rearrangements and/or ERG deletions. These patients had a low relapse rate and excellent survival. ETV6::RUNX1-like subtype (n = 21) was characterised by multiple abnormalities of ETV6 and IKZF1, with no reported relapses or deaths, indicating their excellent prognosis in this trial. An inferior outcome for patients with ABL-class fusions (n = 25) was confirmed. Integration of NGS into genomic profiling of B-other-ALL within a single childhood ALL trial, UKALL2003, has shown the added clinical value of NGS-based approaches, through improved accuracy in detection and classification into the range of risk stratifying genetic subtypes, while validating their prognostic significance.

Loss of ALK hotspot mutations in relapsed neuroblastoma.

ALK is the most commonly mutated oncogene in neuroblastoma with increased mutation frequency reported at relapse. Here we report the loss of an ALK mutation in two patients at relapse and a paired neuroblastoma cell line at relapse. ALK detection methods including Sanger sequencing, targeted next-generation sequencing and a new ALK Agena MassARRAY technique were used to detect common hotspot ALK variants in tumors at diagnosis and relapse from two high-risk neuroblastoma patients. Copy number analysis including single nucleotide polymorphism array and array comparative genomic hybridization confirmed adequate tumor cell content in DNA used for mutation testing. Case 1 presented with an ALK F1174L mutation at diagnosis with a variant allele frequency (VAF) ranging between 23.5% and 28.5%, but the mutation was undetectable at relapse. Case 2 presented with an ALK R1257Q mutation at diagnosis (VAF = 39%-47.4%) which decreased to <0.01% at relapse. Segmental chromosomal aberrations were maintained between diagnosis and relapse confirming sufficient tumor cell content for mutation detection. The diagnostic SKNBE1n cell line harbors an ALK F1174S mutation, which was lost in the relapsed SKNBE2c cell line. To our knowledge, these are the first reported cases of loss of ALK mutations at relapse in neuroblastoma in the absence of ALK inhibitor therapy, reflecting intra-tumoral spatial and temporal heterogeneity. As ALK inhibitors are increasingly used in the treatment of refractory/relapsed neuroblastoma, our study highlights the importance of confirming whether an ALK mutation detected at diagnosis is still present in clones leading to relapse.

Emergence and maintenance of actionable genetic drivers at medulloblastoma relapse.

BACKGROUND: Less than 5% of medulloblastoma (MB) patients survive following failure of contemporary radiation-based therapies. Understanding the molecular drivers of medulloblastoma relapse (rMB) will be essential to improve outcomes. Initial genome-wide investigations have suggested significant genetic divergence of the relapsed disease. METHODS: We undertook large-scale integrated characterization of the molecular features of rMB-molecular subgroup, novel subtypes, copy number variation (CNV), and driver gene mutation. 119 rMBs were assessed in comparison with their paired diagnostic samples (n = 107), alongside an independent reference cohort sampled at diagnosis (n = 282). rMB events were investigated for association with outcome post-relapse in clinically annotated patients (n = 54). RESULTS: Significant genetic evolution occurred over disease-course; 40% of putative rMB drivers emerged at relapse and differed significantly between molecular subgroups. Non-infant MBSHH displayed significantly more chromosomal CNVs at relapse (TP53 mutation-associated). Relapsed MBGroup4 demonstrated the greatest genetic divergence, enriched for targetable (eg, CDK amplifications) and novel (eg, USH2A mutations) events. Importantly, many hallmark features of MB were stable over time; novel subtypes (>90% of tumors) and established genetic drivers (eg, SHH/WNT/P53 mutations; 60% of rMB events) were maintained from diagnosis. Critically, acquired and maintained rMB events converged on targetable pathways which were significantly enriched at relapse (eg, DNA damage signaling) and specific events (eg, 3p loss) predicted survival post-relapse. CONCLUSIONS: rMB is characterised by the emergence of novel events and pathways, in concert with selective maintenance of established genetic drivers. Together, these define the actionable genetic landscape of rMB and provide a basis for improved clinical management and development of stratified therapeutics, across disease-course.

Persistent SARS-CoV-2 infection in patients with secondary antibody deficiency: successful clearance following combination casirivimab and imdevimab (REGN-COV2) monoclonal antibody therapy.

BACKGROUND: There is growing evidence that antibody responses play a role in the resolution of SARS-CoV-2 infection. Patients with primary or secondary antibody deficiency are at increased risk of persistent infection. This challenging clinical scenario is associated with adverse patient outcome and potentially creates an ecological niche for the evolution of novel SARS-CoV-2 variants with immune evasion capacity. Case reports and/or series have implied a therapeutic role for convalescent plasma (CP) to secure virological clearance, although concerns have been raised about the effectiveness of CP and its potential to drive viral evolution, and it has largely been withdrawn from clinical use in the UK. CASE PRESENTATION: We report two cases in which persistent SARS-CoV-2 infection was cleared following administration of the monoclonal antibody combination casirivimab and imdevimab (REGN-COV2, Ronapreve). A 55-year-old male with follicular lymphoma, treated with B cell depleting therapy, developed SARS-CoV-2 infection in September 2020 which then persisted for over 200 days. He was hospitalised on four occasions with COVID-19 and suffered debilitating fatigue and malaise throughout. There was no clinical response to antiviral therapy with remdesivir or CP, and SARS-CoV-2 was consistently detected in nasopharyngeal swabs. Intrahost evolution of several spike variants of uncertain significance was identified by viral sequence analysis. Delivery of REGN-COV2, in combination with remdesivir, was associated with clinical improvement and viral clearance within 6 days, which was sustained for over 150 days despite immunotherapy for relapsed follicular lymphoma. The second case, a 68-year-old female with chronic lymphocytic leukaemia on ibrutinib, also developed persistent SARS-CoV-2 infection. Despite a lack of response to remdesivir, infection promptly cleared following REGN-COV2 in combination with remdesivir, accompanied by resolution of inflammation and full clinical recovery that has been maintained for over 290 days. CONCLUSIONS: These cases highlight the potential benefit of REGN-COV2 as therapy for persistent SARS-CoV-2 infection in antibody deficient individuals, including after failure of CP treatment. Formal clinical studies are warranted to assess the effectiveness of REGN-COV2 in antibody-deficient patients, especially in light of the emergence of variants of concern, such as Omicron, that appear to evade REGN-COV2 neutralisation.

Concordance of copy number abnormality detection using SNP arrays and Multiplex Ligation-dependent Probe Amplification (MLPA) in acute lymphoblastic leukaemia.

In acute lymphoblastic leukaemia, MLPA has been used in research studies to identify clinically relevant copy number abnormality (CNA) profiles. However, in diagnostic settings other techniques are often employed. We assess whether equivalent CNA profiles are called using SNP arrays, ensuring platform independence. We demonstrate concordance between SNP6.0 and MLPA CNA calling on 143 leukaemia samples from two UK trials; comparing 1,287 calls within eight genes and a region. The techniques are 99% concordant using manually augmented calling, and 98% concordant using an automated pipeline. We classify these discordant calls and examine reasons for discordance. In nine cases the circular binary segmentation (CBS) algorithm failed to detect focal abnormalities or those flanking gaps in IKZF1 probe coverage. Eight cases were discordant due to probe design differences, with focal abnormalities detectable using one technique not observable by the other. Risk classification using manually augmented array calling resulted in four out of 143 patients being assigned to a different CNA risk group and eight patients using the automated pipeline. We conclude that MLPA defined CNA profiles can be accurately mirrored by SNP6.0 or similar array platforms. Automated calling using the CBS algorithm proved successful, except for IKZF1 which should be manually inspected.

IKZF1 Deletions with COBL Breakpoints Are Not Driven by RAG-Mediated Recombination Events in Acute Lymphoblastic Leukemia.

IKZF1 deletion (ΔIKZF1) is an important predictor of relapse in both childhood and adult B-cell precursor acute lymphoblastic leukemia (B-ALL). Previously, we revealed that COBL is a hotspot for breakpoints in leukemia and could promote IKZF1 deletions. Through an international collaboration, we provide a detailed genetic and clinical picture of B-ALL with COBL rearrangements (COBL-r). Patients with B-ALL and IKZF1 deletion (n = 133) were included. IKZF1 ∆1-8 were associated with large alterations within chromosome 7: monosomy 7 (18%), isochromosome 7q (10%), 7p loss (19%), and interstitial deletions (53%). The latter included COBL-r, which were found in 12% of the IKZF1 ∆1-8 cohort. Patients with COBL-r are mostly classified as intermediate cytogenetic risk and frequently harbor ETV6, PAX5, CDKN2A/B deletions. Overall, 56% of breakpoints were located within COBL intron 5. Cryptic recombination signal sequence motifs were broadly distributed within the sequence of COBL, and no enrichment for the breakpoint cluster region was found. In summary, a diverse spectrum of alterations characterizes ΔIKZF1 and they also include deletion breakpoints within COBL. We confirmed that COBL is a hotspot associated with ΔIKZF1, but these rearrangements are not driven by RAG-mediated recombination.

SH2B3 inactivation through CN-LOH 12q is uniquely associated with B-cell precursor ALL with iAMP21 or other chromosome 21 gain.

In more than 30% of B-cell precursor acute lymphoblastic leukaemia (B-ALL), chromosome 21 sequence is overrepresented through aneuploidy or structural rearrangements, exemplified by intrachromosomal amplification of chromosome 21 (iAMP21). Although frequent, the mechanisms by which these abnormalities promote B-ALL remain obscure. Intriguingly, we found copy number neutral loss of heterozygosity (CN-LOH) of 12q was recurrent in iAMP21-ALL, but never observed in B-ALL without some form of chromosome 21 gain. As a consequence of CN-LOH 12q, mutations or deletions of the adaptor protein, SH2B3, were converted to homozygosity. In patients without CN-LOH 12q, bi-allelic abnormalities of SH2B3 occurred, but only in iAMP21-ALL, giving an overall incidence of 18% in this sub-type. Review of published data confirmed a tight association between overrepresentation of chromosome 21 and both CN-LOH 12q and SH2B3 abnormalities in B-ALL. Despite relatively small patient numbers, preliminary analysis linked 12q abnormalities to poor outcome in iAMP21-ALL (p = 0.03). Homology modelling of a leukaemia-associated SH2 domain mutation and in vitro analysis of patient-derived xenograft cells implicated the JAK/STAT pathway as one likely target for SH2B3 tumour suppressor activity in iAMP21-ALL.

High prevalence of focal and multi-focal somatic genetic variants in the human brain.

Somatic mutations during stem cell division are responsible for several cancers. In principle, a similar process could occur during the intense cell proliferation accompanying human brain development, leading to the accumulation of regionally distributed foci of mutations. Using dual platform >5000-fold depth sequencing of 102 genes in 173 adult human brain samples, we detect and validate somatic mutations in 27 of 54 brains. Using a mathematical model of neurodevelopment and approximate Bayesian inference, we predict that macroscopic islands of pathologically mutated neurons are likely to be common in the general population. The detected mutation spectrum also includes DNMT3A and TET2 which are likely to have originated from blood cell lineages. Together, these findings establish developmental mutagenesis as a potential mechanism for neurodegenerative disorders, and provide a novel mechanism for the regional onset and focal pathology in sporadic cases.

Dynamic clonal progression in xenografts of acute lymphoblastic leukemia with intrachromosomal amplification of chromosome 21.

Intrachromosomal amplification of chromosome 21 is a heterogeneous chromosomal rearrangement occurring in 2% of cases of childhood precursor B-cell acute lymphoblastic leukemia. These abnormalities are too complex to engineer faithfully in animal models and are unrepresented in leukemia cell lines. As a resource for future functional and preclinical studies, we have created xenografts from the leukemic blasts of patients with intrachromosomal amplification of chromosome 21 and characterized them by in-vivo and ex-vivo luminescent imaging, flow immunophenotyping, and histological and ultrastructural analyses of bone marrow and the central nervous system. Investigation of up to three generations of xenografts revealed phenotypic evolution, branching genomic architecture and, compared with other B-cell acute lymphoblastic leukemia genetic subtypes, greater clonal diversity of leukemia-initiating cells. In support of intrachromosomal amplification of chromosome 21 as a primary genetic abnormality, it was always retained through generations of xenografts, although we also observed the first example of structural evolution of this rearrangement. Clonal segregation in xenografts revealed convergent evolution of different secondary genomic abnormalities implicating several known tumor suppressor genes and a region, containing the B-cell adaptor, PIK3AP1, and nuclear receptor co-repressor, LCOR, in the progression of B-cell acute lymphoblastic leukemia. Tracking of mutations in patients and derived xenografts provided evidence for co-operation between abnormalities activating the RAS pathway in B-cell acute lymphoblastic leukemia and for their aggressive clonal expansion in the xeno-environment. Bi-allelic loss of the CDKN2A/B locus was recurrently maintained or emergent in xenografts and also strongly selected as RNA sequencing demonstrated a complete absence of reads for genes associated with the deletions.

Cognate ligand domain mapping for enzymes.

Here, we present an automatic assignment of potential cognate ligands to domains of enzymes in the CATH and SCOP protein domain classifications on the basis of structural data available in the wwPDB. This procedure involves two steps; firstly, we assign the binding of particular ligands to particular domains; secondly, we compare the chemical similarity of the PDB ligands to ligands in KEGG in order to assign cognate ligands. We find that use of the Enzyme Commission (EC) numbers is necessary to enable efficient and accurate cognate ligand assignment. The PROCOGNATE database currently has cognate ligand mapping for 3277 (4118) protein structures and 351 (302) superfamilies, as described by the CATH and (SCOP) databases, respectively. We find that just under half of all ligands are only and always bound by a single domain, with 16% bound by more than one domain and the remainder of the ligands showing a variety of binding modes. This finding has implications for domain recombination and the evolution of new protein functions. Domain architecture or context is also found to affect substrate specificity of particular domains, and we discuss example cases. The most popular PDB ligands are all found to be generic components of crystallisation buffers, highlighting the non-cognate ligand problem inherent in the PDB. In contrast, the most popular cognate ligands are all found to be universal cellular currencies of reducing power and energy such as NADH, FADH2 and ATP, respectively, reflecting the fact that the vast majority of enzymatic reactions utilise one of these popular co-factors. These ligands all share a common adenine ribonucleotide moiety, suggesting that many different domain superfamilies have converged to bind this chemical framework.