AI-enabled evaluation of genome-wide association relevance and polygenic risk score prediction in Alzheimer's disease.

GWAS focuses on significance loosing false positives; machine learning probes sub-significant features relying on predictivity. Yet, these are far from orthogonal. We sought to explore how these inform each other in sub-genome-wide significant situations to define relevance for predictive features. We introduce the SVM-based RubricOE that selects heavily cross-validated feature sets, and LDpred2 PRS as a strong contrast to SVM, to explore significance and predictivity. Our Alzheimer's test case notoriously lacks strong genetic signals except for few very strong phenotype-SNP associations, which suits the problem we are exploring. We found that the most significant SNPs among ML and PRS-selected SNPs captured most of the predictivity, while weaker associations tend also to contribute weakly to predictivity. SNPs with weak associations tend not to contribute to predictivity, but deletion of these features does not injure it. Significance provides a ranking that helps identify weakly predictive features.

Germline mutations and developmental mosaicism underlying EGFR-mutant lung cancer.

While the development of multiple primary tumors in smokers with lung cancer can be attributed to carcinogen-induced field cancerization, the occurrence of multiple primary tumors in individuals with EGFR-mutant lung cancer who lack known environmental exposures remains unexplained. We identified ten patients with early-stage, resectable non-small cell lung cancer who presented with multiple anatomically distinct EGFR-mutant tumors. We analyzed the phylogenetic relationships among multiple tumors from each patient using whole exome sequencing (WES) and hypermutable poly-guanine (poly-G) repeat genotyping, as orthogonal methods for lineage tracing. In two patients, we identified germline EGFR variants, which confer moderately enhanced signaling when modeled in vitro. In four other patients, developmental mosaicism is supported by the poly-G lineage tracing and WES, indicating a common non-germline cell-of-origin. Thus, developmental mosaicism and germline variants define two distinct mechanisms of genetic predisposition to multiple EGFR-mutant primary tumors, with implications for understanding their etiology and clinical management.

Genetic events associated with venetoclax resistance in CLL identified by whole-exome sequencing of patient samples.

Although BCL2 mutations are reported as later occurring events leading to venetoclax resistance, many other mechanisms of progression have been reported though remain poorly understood. Here, we analyze longitudinal tumor samples from 11 patients with disease progression while receiving venetoclax to characterize the clonal evolution of resistance. All patients tested showed increased in vitro resistance to venetoclax at the posttreatment time point. We found the previously described acquired BCL2-G101V mutation in only 4 of 11 patients, with 2 patients showing a very low variant allele fraction (0.03%-4.68%). Whole-exome sequencing revealed acquired loss(8p) in 4 of 11 patients, of which 2 patients also had gain (1q21.2-21.3) in the same cells affecting the MCL1 gene. In vitro experiments showed that CLL cells from the 4 patients with loss(8p) were more resistant to venetoclax than cells from those without it, with the cells from 2 patients also carrying gain (1q21.2-21.3) showing increased sensitivity to MCL1 inhibition. Progression samples with gain (1q21.2-21.3) were more susceptible to the combination of MCL1 inhibitor and venetoclax. Differential gene expression analysis comparing bulk RNA sequencing data from pretreatment and progression time points of all patients showed upregulation of proliferation, B-cell receptor (BCR), and NF-κB gene sets including MAPK genes. Cells from progression time points demonstrated upregulation of surface immunoglobulin M and higher pERK levels compared with those from the preprogression time point, suggesting an upregulation of BCR signaling that activates the MAPK pathway. Overall, our data suggest several mechanisms of acquired resistance to venetoclax in CLL that could pave the way for rationally designed combination treatments for patients with venetoclax-resistant CLL.

Lesion Shedding Model: unraveling site-specific contributions to ctDNA.

Sampling circulating tumor DNA (ctDNA) using liquid biopsies offers clinically important benefits for monitoring cancer progression. A single ctDNA sample represents a mixture of shed tumor DNA from all known and unknown lesions within a patient. Although shedding levels have been suggested to hold the key to identifying targetable lesions and uncovering treatment resistance mechanisms, the amount of DNA shed by any one specific lesion is still not well characterized. We designed the Lesion Shedding Model (LSM) to order lesions from the strongest to the poorest shedding for a given patient. By characterizing the lesion-specific ctDNA shedding levels, we can better understand the mechanisms of shedding and more accurately interpret ctDNA assays to improve their clinical impact. We verified the accuracy of the LSM under controlled conditions using a simulation approach as well as testing the model on three cancer patients. The LSM obtained an accurate partial order of the lesions according to their assigned shedding levels in simulations and its accuracy in identifying the top shedding lesion was not significantly impacted by number of lesions. Applying LSM to three cancer patients, we found that indeed there were lesions that consistently shed more than others into the patients' blood. In two of the patients, the top shedding lesion was one of the only clinically progressing lesions at the time of biopsy suggesting a connection between high ctDNA shedding and clinical progression. The LSM provides a much needed framework with which to understand ctDNA shedding and to accelerate discovery of ctDNA biomarkers. The LSM source code has been available in the IBM BioMedSciAI Github (https://github.com/BiomedSciAI/Geno4SD).

Evolutionary history of transformation from chronic lymphocytic leukemia to Richter syndrome.

Richter syndrome (RS) arising from chronic lymphocytic leukemia (CLL) exemplifies an aggressive malignancy that develops from an indolent neoplasm. To decipher the genetics underlying this transformation, we computationally deconvoluted admixtures of CLL and RS cells from 52 patients with RS, evaluating paired CLL-RS whole-exome sequencing data. We discovered RS-specific somatic driver mutations (including IRF2BP2, SRSF1, B2M, DNMT3A and CCND3), recurrent copy-number alterations beyond del(9p21)(CDKN2A/B), whole-genome duplication and chromothripsis, which were confirmed in 45 independent RS cases and in an external set of RS whole genomes. Through unsupervised clustering, clonally related RS was largely distinct from diffuse large B cell lymphoma. We distinguished pathways that were dysregulated in RS versus CLL, and detected clonal evolution of transformation at single-cell resolution, identifying intermediate cell states. Our study defines distinct molecular subtypes of RS and highlights cell-free DNA analysis as a potential tool for early diagnosis and monitoring.

Pirtobrutinib targets BTK C481S in ibrutinib-resistant CLL but second-site BTK mutations lead to resistance.

Covalent inhibitors of Bruton tyrosine kinase (BTK) have transformed the therapy of chronic lymphocytic leukemia (CLL), but continuous therapy has been complicated by the development of resistance. The most common resistance mechanism in patients whose disease progresses on covalent BTK inhibitors (BTKis) is a mutation in the BTK 481 cysteine residue to which the inhibitors bind covalently. Pirtobrutinib is a highly selective, noncovalent BTKi with substantial clinical activity in patients whose disease has progressed on covalent BTKi, regardless of BTK mutation status. Using in vitro ibrutinib-resistant models and cells from patients with CLL, we show that pirtobrutinib potently inhibits BTK-mediated functions including B-cell receptor (BCR) signaling, cell viability, and CCL3/CCL4 chemokine production in both BTK wild-type and C481S mutant CLL cells. We demonstrate that primary CLL cells from responding patients on the pirtobrutinib trial show reduced BCR signaling, cell survival, and CCL3/CCL4 chemokine secretion. At time of progression, these primary CLL cells show increasing resistance to pirtobrutinib in signaling inhibition, cell viability, and cytokine production. We employed longitudinal whole-exome sequencing on 2 patients whose disease progressed on pirtobrutinib and identified selection of alternative-site BTK mutations, providing clinical evidence that secondary BTK mutations lead to resistance to noncovalent BTKis.

Distinct cellular dynamics associated with response to CAR-T therapy for refractory B cell lymphoma.

Chimeric antigen receptor (CAR)-T cell therapy has revolutionized the treatment of hematologic malignancies. Approximately half of patients with refractory large B cell lymphomas achieve durable responses from CD19-targeting CAR-T treatment; however, failure mechanisms are identified in only a fraction of cases. To gain new insights into the basis of clinical response, we performed single-cell transcriptome sequencing of 105 pretreatment and post-treatment peripheral blood mononuclear cell samples, and infusion products collected from 32 individuals with large B cell lymphoma treated with either of two CD19 CAR-T products: axicabtagene ciloleucel (axi-cel) or tisagenlecleucel (tisa-cel). Expansion of proliferative memory-like CD8 clones was a hallmark of tisa-cel response, whereas axi-cel responders displayed more heterogeneous populations. Elevations in CAR-T regulatory cells among nonresponders to axi-cel were detected, and these populations were capable of suppressing conventional CAR-T cell expansion and driving late relapses in an in vivo model. Our analyses reveal the temporal dynamics of effective responses to CAR-T therapy, the distinct molecular phenotypes of CAR-T cells with differing designs, and the capacity for even small increases in CAR-T regulatory cells to drive relapse.

A common methodological phylogenomics framework for intra-patient heteroplasmies to infer SARS-CoV-2 sublineages and tumor clones.

BACKGROUND: All diseases containing genetic material undergo genetic evolution and give rise to heterogeneity including cancer and infection. Although these illnesses are biologically very different, the ability for phylogenetic retrodiction based on the genomic reads is common between them and thus tree-based principles and assumptions are shared. Just as the different frequencies of tumor genomic variants presupposes the existence of multiple tumor clones and provides a handle to computationally infer them, we postulate that the different variant frequencies in viral reads offers the means to infer multiple co-infecting sublineages. RESULTS: We present a common methodological framework to infer the phylogenomics from genomic data, be it reads of SARS-CoV-2 of multiple COVID-19 patients or bulk DNAseq of the tumor of a cancer patient. We describe the Concerti computational framework for inferring phylogenies in each of the two scenarios.To demonstrate the accuracy of the method, we reproduce some known results in both scenarios. We also make some additional discoveries. CONCLUSIONS: Concerti successfully extracts and integrates information from multi-point samples, enabling the discovery of clinically plausible phylogenetic trees that capture the heterogeneity known to exist both spatially and temporally. These models can have direct therapeutic implications by highlighting "birth" of clones that may harbor resistance mechanisms to treatment, "death" of subclones with drug targets, and acquisition of functionally pertinent mutations in clones that may have seemed clinically irrelevant. Specifically in this paper we uncover new potential parallel mutations in the evolution of the SARS-CoV-2 virus. In the context of cancer, we identify new clones harboring resistant mutations to therapy.

K-Mer Analyses Reveal Different Evolutionary Histories of Alpha, Beta, and Gamma Papillomaviruses.

Papillomaviruses (PVs) are a heterogeneous group of DNA viruses that can infect fish, birds, reptiles, and mammals. PVs infecting humans (HPVs) phylogenetically cluster into five genera (Alpha-, Beta-, Gamma-, Mu- and Nu-PV), with differences in tissue tropism and carcinogenicity. The evolutionary features associated with the divergence of Papillomaviridae are not well understood. Using a combination of k-mer distributions, genetic metrics, and phylogenetic algorithms, we sought to evaluate the characteristics and differences of Alpha-, Beta- and Gamma-PVs constituting the majority of HPV genomes. A total of 640 PVs including 442 HPV types, 27 non-human primate PV types, and 171 non-primate animal PV types were evaluated. Our analyses revealed the highest genetic diversity amongst Gamma-PVs compared to the Alpha and Beta PVs, suggesting reduced selective pressures on Gamma-PVs. Using a sequence alignment-free trimer (k = 3) phylogeny algorithm, we reconstructed a phylogeny that grouped most HPV types into a monophyletic clade that was further split into three branches similar to alignment-based classifications. Interestingly, a subset of low-risk Alpha HPVs (the species Alpha-2, 3, 4, and 14) split from other HPVs and were clustered with non-human primate PVs. Surprisingly, the trimer-constructed phylogeny grouped the Gamma-6 species types originally isolated from the cervicovaginal region with the main Alpha-HPV clade. These data indicate that characterization of papillomavirus heterogeneity via orthogonal approaches reveals novel insights into the biological understanding of HPV genomes.

Parallel Genomic Alterations of Antigen and Payload Targets Mediate Polyclonal Acquired Clinical Resistance to Sacituzumab Govitecan in Triple-Negative Breast Cancer.

Sacituzumab govitecan (SG), the first antibody-drug conjugate (ADC) approved for triple-negative breast cancer, incorporates the anti-TROP2 antibody hRS7 conjugated to a topoisomerase-1 (TOP1) inhibitor payload. We sought to identify mechanisms of SG resistance through RNA and whole-exome sequencing of pretreatment and postprogression specimens. One patient exhibiting de novo progression lacked TROP2 expression, in contrast to robust TROP2 expression and focal genomic amplification of TACSTD2/TROP2 observed in a patient with a deep, prolonged response to SG. Analysis of acquired genomic resistance in this case revealed one phylogenetic branch harboring a canonical TOP1 E418K resistance mutation and subsequent frameshift TOP1 mutation, whereas a distinct branch exhibited a novel TACSTD2/TROP2 T256R missense mutation. Reconstitution experiments demonstrated that TROP2T256R confers SG resistance via defective plasma membrane localization and reduced cell-surface binding by hRS7. These findings highlight parallel genomic alterations in both antibody and payload targets associated with resistance to SG. SIGNIFICANCE: These findings underscore TROP2 as a response determinant and reveal acquired SG resistance mechanisms involving the direct antibody and drug payload targets in distinct metastatic subclones of an individual patient. This study highlights the specificity of SG and illustrates how such mechanisms will inform therapeutic strategies to overcome ADC resistance.This article is highlighted in the In This Issue feature, p. 2355.

Functional profiling of COVID-19 respiratory tract microbiomes.

In response to the ongoing global pandemic, characterizing the molecular-level host interactions of the new coronavirus SARS-CoV-2 responsible for COVID-19 has been at the center of unprecedented scientific focus. However, when the virus enters the body it also interacts with the micro-organisms already inhabiting the host. Understanding the virus-host-microbiome interactions can yield additional insights into the biological processes perturbed by viral invasion. Alterations in the gut microbiome species and metabolites have been noted during respiratory viral infections, possibly impacting the lungs via gut-lung microbiome crosstalk. To better characterize microbial functions in the lower respiratory tract during COVID-19 infection, we carry out a functional analysis of previously published metatranscriptome sequencing data of bronchoalveolar lavage fluid from eight COVID-19 cases, twenty-five community-acquired pneumonia patients, and twenty healthy controls. The functional profiles resulting from comparing the sequences against annotated microbial protein domains clearly separate the cohorts. By examining the associated metabolic pathways, distinguishing functional signatures in COVID-19 respiratory tract microbiomes are identified, including decreased potential for lipid metabolism and glycan biosynthesis and metabolism pathways, and increased potential for carbohydrate metabolism pathways. The results include overlap between previous studies on COVID-19 microbiomes, including decrease in the glycosaminoglycan degradation pathway and increase in carbohydrate metabolism. The results also suggest novel connections to consider, possibly specific to the lower respiratory tract microbiome, calling for further research on microbial functions and host-microbiome interactions during SARS-CoV-2 infection.

Re-purposing software for functional characterization of the microbiome.

BACKGROUND: Widespread bioinformatic resource development generates a constantly evolving and abundant landscape of workflows and software. For analysis of the microbiome, workflows typically begin with taxonomic classification of the microorganisms that are present in a given environment. Additional investigation is then required to uncover the functionality of the microbial community, in order to characterize its currently or potentially active biological processes. Such functional analysis of metagenomic data can be computationally demanding for high-throughput sequencing experiments. Instead, we can directly compare sequencing reads to a functionally annotated database. However, since reads frequently match multiple sequences equally well, analyses benefit from a hierarchical annotation tree, e.g. for taxonomic classification where reads are assigned to the lowest taxonomic unit. RESULTS: To facilitate functional microbiome analysis, we re-purpose well-known taxonomic classification tools to allow us to perform direct functional sequencing read classification with the added benefit of a functional hierarchy. To enable this, we develop and present a tree-shaped functional hierarchy representing the molecular function subset of the Gene Ontology annotation structure. We use this functional hierarchy to replace the standard phylogenetic taxonomy used by the classification tools and assign query sequences accurately to the lowest possible molecular function in the tree. We demonstrate this with simulated and experimental datasets, where we reveal new biological insights. CONCLUSIONS: We demonstrate that improved functional classification of metagenomic sequencing reads is possible by re-purposing a range of taxonomic classification tools that are already well-established, in conjunction with either protein or nucleotide reference databases. We leverage the advances in speed, accuracy and efficiency that have been made for taxonomic classification and translate these benefits for the rapid functional classification of microbiomes. While we focus on a specific set of commonly used methods, the functional annotation approach has broad applicability across other sequence classification tools. We hope that re-purposing becomes a routine consideration during bioinformatic resource development. Video abstract.

Impact of Clinical and Genomic Factors on COVID-19 Disease Severity.

To date, there have been 180 million confirmed cases of COVID-19, with more than 3.8 million deaths, reported to WHO worldwide. In this paper we address the problem of understanding the host genome's influence, in concert with clinical variables, on the severity of COVID-19 manifestation in the patient. Leveraging positive-unlabeled machine learning algorithms coupled with RubricOE, a state-of-the-art genomic analysis framework, on UK BioBank data we extract novel insights on the complex interplay. The algorithm is also sensitive enough to detect the changing influence of the emergent B.1.1.7 SARS-CoV-2 (alpha) variant on disease severity, and, changing treatment protocols. The genomic component also implicates biological pathways that can help in understanding the disease etiology. Our work demonstrates that it is possible to build a robust and sensitive model despite significant bias, noise and incompleteness in both clinical and genomic data by a careful interleaving of clinical and genomic methodologies.

Hierarchically Labeled Database Indexing Allows Scalable Characterization of Microbiomes.

Increasingly available microbial reference data allow interpreting the composition and function of previously uncharacterized microbial communities in detail, via high-throughput sequencing analysis. However, efficient methods for read classification are required when the best database matches for short sequence reads are often shared among multiple reference sequences. Here, we take advantage of the fact that microbial sequences can be annotated relative to established tree structures, and we develop a highly scalable read classifier, PRROMenade, by enhancing the generalized Burrows-Wheeler transform with a labeling step to directly assign reads to the corresponding lowest taxonomic unit in an annotation tree. PRROMenade solves the multi-matching problem while allowing fast variable-size sequence classification for phylogenetic or functional annotation. Our simulations with 5% added differences from reference indicated only 1.5% error rate for PRROMenade functional classification. On metatranscriptomic data PRROMenade highlighted biologically relevant functional pathways related to diet-induced changes in the human gut microbiome.

Author Correction: Liquid versus tissue biopsy for detecting acquired resistance and tumor heterogeneity in gastrointestinal cancers.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Liquid versus tissue biopsy for detecting acquired resistance and tumor heterogeneity in gastrointestinal cancers.

During cancer therapy, tumor heterogeneity can drive the evolution of multiple tumor subclones harboring unique resistance mechanisms in an individual patient1-3. Previous case reports and small case series have suggested that liquid biopsy (specifically, cell-free DNA (cfDNA)) may better capture the heterogeneity of acquired resistance4-8. However, the effectiveness of cfDNA versus standard single-lesion tumor biopsies has not been directly compared in larger-scale prospective cohorts of patients following progression on targeted therapy. Here, in a prospective cohort of 42 patients with molecularly defined gastrointestinal cancers and acquired resistance to targeted therapy, direct comparison of postprogression cfDNA versus tumor biopsy revealed that cfDNA more frequently identified clinically relevant resistance alterations and multiple resistance mechanisms, detecting resistance alterations not found in the matched tumor biopsy in 78% of cases. Whole-exome sequencing of serial cfDNA, tumor biopsies and rapid autopsy specimens elucidated substantial geographic and evolutionary differences across lesions. Our data suggest that acquired resistance is frequently characterized by profound tumor heterogeneity, and that the emergence of multiple resistance alterations in an individual patient may represent the 'rule' rather than the 'exception'. These findings have profound therapeutic implications and highlight the potential advantages of cfDNA over tissue biopsy in the setting of acquired resistance.

Dark-matter matters: Discriminating subtle blood cancers using the darkest DNA.

The confluence of deep sequencing and powerful machine learning is providing an unprecedented peek at the darkest of the dark genomic matter, the non-coding genomic regions lacking any functional annotation. While deep sequencing uncovers rare tumor variants, the heterogeneity of the disease confounds the best of machine learning (ML) algorithms. Here we set out to answer if the dark-matter of the genome encompass signals that can distinguish the fine subtypes of disease that are otherwise genomically indistinguishable. We introduce a novel stochastic regularization, ReVeaL, that empowers ML to discriminate subtle cancer subtypes even from the same 'cell of origin'. Analogous to heritability, implicitly defined on whole genome, we use predictability (F1 score) definable on portions of the genome. In an effort to distinguish cancer subtypes using dark-matter DNA, we applied ReVeaL to a new WGS dataset from 727 patient samples with seven forms of hematological cancers and assessed the predictivity over several genomic regions including genic, non-dark, non-coding, non-genic, and dark. ReVeaL enabled improved discrimination of cancer subtypes for all segments of the genome. The non-genic, non-coding and dark-matter had the highest F1 scores, with dark-matter having the highest level of predictability. Based on ReVeaL's predictability of different genomic regions, dark-matter contains enough signal to significantly discriminate fine subtypes of disease. Hence, the agglomeration of rare variants, even in the hitherto unannotated and ill-understood regions of the genome, may play a substantial role in the disease etiology and deserve much more attention.

Correction to: Sequencing and curation strategies for identifying candidate glioblastoma treatments.

Following publication of the original article [1], it was reported that the given name of the fourteenth author was incorrectly published. The incorrect and the correct names are given below.

Transcriptome characterization and differentially expressed genes under flooding and drought stress in the biomass grasses Phalaris arundinacea and Dactylis glomerata.

BACKGROUND AND AIMS: Perennial grasses are a global resource as forage, and for alternative uses in bioenergy and as raw materials for the processing industry. Marginal lands can be valuable for perennial biomass grass production, if perennial biomass grasses can cope with adverse abiotic environmental stresses such as drought and waterlogging. METHODS: In this study, two perennial grass species, reed canary grass (Phalaris arundinacea) and cocksfoot (Dactylis glomerata) were subjected to drought and waterlogging stress to study their responses for insights to improving environmental stress tolerance. Physiological responses were recorded, reference transcriptomes established and differential gene expression investigated between control and stress conditions. We applied a robust non-parametric method, RoDEO, based on rank ordering of transcripts to investigate differential gene expression. Furthermore, we extended and validated vRoDEO for comparing samples with varying sequencing depths. KEY RESULTS: This allowed us to identify expressed genes under drought and waterlogging whilst using only a limited number of RNA sequencing experiments. Validating the methodology, several differentially expressed candidate genes involved in the stage 3 step-wise scheme in detoxification and degradation of xenobiotics were recovered, while several novel stress-related genes classified as of unknown function were discovered. CONCLUSIONS: Reed canary grass is a species coping particularly well with flooding conditions, but this study adds novel information on how its transcriptome reacts under drought stress. We built extensive transcriptomes for the two investigated C3 species cocksfoot and reed canary grass under both extremes of water stress to provide a clear comparison amongst the two species to broaden our horizon for comparative studies, but further confirmation of the data would be ideal to obtain a more detailed picture.