The monoclonal antibody combination REGEN-COV protects against SARS-CoV-2 mutational escape in preclinical and human studies.

Monoclonal antibodies against SARS-CoV-2 are a clinically validated therapeutic option against COVID-19. Because rapidly emerging virus mutants are becoming the next major concern in the fight against the global pandemic, it is imperative that these therapeutic treatments provide coverage against circulating variants and do not contribute to development of treatment-induced emergent resistance. To this end, we investigated the sequence diversity of the spike protein and monitored emergence of virus variants in SARS-COV-2 isolates found in COVID-19 patients treated with the two-antibody combination REGEN-COV, as well as in preclinical in vitro studies using single, dual, or triple antibody combinations, and in hamster in vivo studies using REGEN-COV or single monoclonal antibody treatments. Our study demonstrates that the combination of non-competing antibodies in REGEN-COV provides protection against all current SARS-CoV-2 variants of concern/interest and also protects against emergence of new variants and their potential seeding into the population in a clinical setting.

Common and rare variant associations with clonal haematopoiesis phenotypes.

Clonal haematopoiesis involves the expansion of certain blood cell lineages and has been associated with ageing and adverse health outcomes1-5. Here we use exome sequence data on 628,388 individuals to identify 40,208 carriers of clonal haematopoiesis of indeterminate potential (CHIP). Using genome-wide and exome-wide association analyses, we identify 24 loci (21 of which are novel) where germline genetic variation influences predisposition to CHIP, including missense variants in the lymphocytic antigen coding gene LY75, which are associated with reduced incidence of CHIP. We also identify novel rare variant associations with clonal haematopoiesis and telomere length. Analysis of 5,041 health traits from the UK Biobank (UKB) found relationships between CHIP and severe COVID-19 outcomes, cardiovascular disease, haematologic traits, malignancy, smoking, obesity, infection and all-cause mortality. Longitudinal and Mendelian randomization analyses revealed that CHIP is associated with solid cancers, including non-melanoma skin cancer and lung cancer, and that CHIP linked to DNMT3A is associated with the subsequent development of myeloid but not lymphoid leukaemias. Additionally, contrary to previous findings from the initial 50,000 UKB exomes6, our results in the full sample do not support a role for IL-6 inhibition in reducing the risk of cardiovascular disease among CHIP carriers. Our findings demonstrate that CHIP represents a complex set of heterogeneous phenotypes with shared and unique germline genetic causes and varied clinical implications.

A cellular and molecular spatial atlas of dystrophic muscle.

Asynchronous skeletal muscle degeneration/regeneration is a hallmark feature of Duchenne muscular dystrophy (DMD); however, traditional -omics technologies that lack spatial context make it difficult to study the biological mechanisms of how asynchronous regeneration contributes to disease progression. Here, using the severely dystrophic D2-mdx mouse model, we generated a high-resolution cellular and molecular spatial atlas of dystrophic muscle by integrating spatial transcriptomics and single-cell RNAseq datasets. Unbiased clustering revealed nonuniform distribution of unique cell populations throughout D2-mdx muscle that were associated with multiple regenerative timepoints, demonstrating that this model faithfully recapitulates the asynchronous regeneration observed in human DMD muscle. By probing spatiotemporal gene expression signatures, we found that propagation of inflammatory and fibrotic signals from locally damaged areas contributes to widespread pathology and that querying expression signatures within discrete microenvironments can identify targetable pathways for DMD therapy. Overall, this spatial atlas of dystrophic muscle provides a valuable resource for studying DMD disease biology and therapeutic target discovery.

Activin E-ACVR1C cross talk controls energy storage via suppression of adipose lipolysis in mice.

Body fat distribution is a heritable risk factor for cardiovascular and metabolic disease. In humans, rare Inhibin beta E (INHBE, activin E) loss-of-function variants are associated with a lower waist-to-hip ratio and protection from type 2 diabetes. Hepatic fatty acid sensing promotes INHBE expression during fasting and in obese individuals, yet it is unclear how the hepatokine activin E governs body shape and energy metabolism. Here, we uncover activin E as a regulator of adipose energy storage. By suppressing ß-agonist-induced lipolysis, activin E promotes fat accumulation and adipocyte hypertrophy and contributes to adipose dysfunction in mice. Mechanistically, we demonstrate that activin E elicits its effect on adipose tissue through ACVR1C, activating SMAD2/3 signaling and suppressing PPARG target genes. Conversely, loss of activin E or ACVR1C in mice increases fat utilization, lowers adiposity, and drives PPARG-regulated gene signatures indicative of healthy adipose function. Our studies identify activin E-ACVR1C as a metabolic rheostat promoting liver-adipose cross talk to restrain excessive fat breakdown and preserve fat mass during prolonged fasting, a mechanism that is maladaptive in obese individuals.

Interactive analysis and assessment of single-cell copy-number variations.

We present Ginkgo (http://qb.cshl.edu/ginkgo), a user-friendly, open-source web platform for the analysis of single-cell copy-number variations (CNVs). Ginkgo automatically constructs copy-number profiles of cells from mapped reads and constructs phylogenetic trees of related cells. We validated Ginkgo by reproducing the results of five major studies. After comparing three commonly used single-cell amplification techniques, we concluded that degenerate oligonucleotide-primed PCR is the most consistent for CNV analysis.

Equitability, mutual information, and the maximal information coefficient.

How should one quantify the strength of association between two random variables without bias for relationships of a specific form? Despite its conceptual simplicity, this notion of statistical "equitability" has yet to receive a definitive mathematical formalization. Here we argue that equitability is properly formalized by a self-consistency condition closely related to Data Processing Inequality. Mutual information, a fundamental quantity in information theory, is shown to satisfy this equitability criterion. These findings are at odds with the recent work of Reshef et al. [Reshef DN, et al. (2011) Science 334(6062):1518-1524], which proposed an alternative definition of equitability and introduced a new statistic, the "maximal information coefficient" (MIC), said to satisfy equitability in contradistinction to mutual information. These conclusions, however, were supported only with limited simulation evidence, not with mathematical arguments. Upon revisiting these claims, we prove that the mathematical definition of equitability proposed by Reshef et al. cannot be satisfied by any (nontrivial) dependence measure. We also identify artifacts in the reported simulation evidence. When these artifacts are removed, estimates of mutual information are found to be more equitable than estimates of MIC. Mutual information is also observed to have consistently higher statistical power than MIC. We conclude that estimating mutual information provides a natural (and often practical) way to equitably quantify statistical associations in large datasets.

Parametric inference in the large data limit using maximally informative models.

Motivated by data-rich experiments in transcriptional regulation and sensory neuroscience, we consider the following general problem in statistical inference: when exposed to a high-dimensional signal S, a system of interest computes a representation R of that signal, which is then observed through a noisy measurement M. From a large number of signals and measurements, we wish to infer the "filter" that maps S to R. However, the standard method for solving such problems, likelihood-based inference, requires perfect a priori knowledge of the "noise function" mapping R to M. In practice such noise functions are usually known only approximately, if at all, and using an incorrect noise function will typically bias the inferred filter. Here we show that in the large data limit, this need for a precharacterized noise function can be circumvented by searching for filters that instead maximize the mutual information I[M; R] between observed measurements and predicted representations. Moreover, if the correct filter lies within the space of filters being explored, maximizing mutual information becomes equivalent to simultaneously maximizing every dependence measure that satisfies the data processing inequality. It is important to note that maximizing mutual information will typically leave a small number of directions in parameter space unconstrained. We term these directions diffeomorphic modes and present an equation that allows these modes to be derived systematically. The presence of diffeomorphic modes reflects a fundamental and nontrivial substructure within parameter space, one that is obscured by standard likelihood-based inference.

Designing meaningful continuous representations of T cell receptor sequences with deep generative models.

T Cell Receptor (TCR) antigen binding underlies a key mechanism of the adaptive immune response yet the vast diversity of TCRs and the complexity of protein interactions limits our ability to build useful low dimensional representations of TCRs. To address the current limitations in TCR analysis we develop a capacity-controlled disentangling variational autoencoder trained using a dataset of approximately 100 million TCR sequences, that we name TCR-VALID. We design TCR-VALID such that the model representations are low-dimensional, continuous, disentangled, and sufficiently informative to provide high-quality TCR sequence de novo generation. We thoroughly quantify these properties of the representations, providing a framework for future protein representation learning in low dimensions. The continuity of TCR-VALID representations allows fast and accurate TCR clustering and is benchmarked against other state-of-the-art TCR clustering tools and pre-trained language models.

Digital wearable insole-based identification of knee arthropathies and gait signatures using machine learning.

Gait is impaired in musculoskeletal conditions, such as knee arthropathy. Gait analysis is used in clinical practice to inform diagnosis and monitor disease progression or intervention response. However, clinical gait analysis relies on subjective visual observation of walking as objective gait analysis has not been possible within clinical settings due to the expensive equipment, large-scale facilities, and highly trained staff required. Relatively low-cost wearable digital insoles may offer a solution to these challenges. In this work, we demonstrate how a digital insole measuring osteoarthritis-specific gait signatures yields similar results to the clinical gait-lab standard. To achieve this, we constructed a machine learning model, trained on force plate data collected in participants with knee arthropathy and controls. This model was highly predictive of force plate data from a validation set (area under the receiver operating characteristics curve [auROC] = 0.86; area under the precision-recall curve [auPR] = 0.90) and of a separate, independent digital insole dataset containing control and knee osteoarthritis subjects (auROC = 0.83; auPR = 0.86). After showing that digital insole-derived gait characteristics are comparable to traditional gait measurements, we next showed that a single stride of raw sensor time-series data could be accurately assigned to each subject, highlighting that individuals using digital insoles can be identified by their gait characteristics. This work provides a framework for a promising alternative to traditional clinical gait analysis methods, adds to the growing body of knowledge regarding wearable technology analytical pipelines, and supports clinical development of at-home gait assessments, with the potential to improve the ease, frequency, and depth of patient monitoring.

The way we walk ­ our 'gait' ­ is a key indicator of health. Gait irregularities like limping, shuffling or a slow pace can be signs of muscle or joint problems. Assessing a patient's gait is therefore an important element in diagnosing these conditions, and in evaluating whether treatments are working. Gait is often assessed via a simple visual inspection, with patients being asked to walk back and forth in a doctor's office. While quick and easy, this approach is highly subjective and therefore imprecise. 'Objective gait analysis' is a more accurate alternative, but it relies on tests being conducted in specialised laboratories with large-scale, expensive equipment operated by highly trained staff. Unfortunately, this means that gait laboratories are not accessible for everyday clinical use. In response, Wipperman et al. aimed to develop a low-cost alternative to the complex equipment used in gait laboratories. To do this, they harnessed wearable sensor technologies ­ devices that can directly measure physiological data while embedded in clothing or attached to the user. Wearable sensors have the advantage of being cheap, easy to use, and able to provide clinically useful information without specially trained staff. Wipperman et al. analysed data from classic gait laboratory devices, as well as 'digital insoles' equipped with sensors that captured foot movements and pressure as participants walked. The analysis first 'trained' on data from gait laboratories (called force plates) and then applied the method to gait measurements obtained from digital insoles worn by either healthy participants or patients with knee problems. Analysis of the pressure data from the insoles confirmed that they could accurately predict which measurements were from healthy individuals, and which were from patients. The gait characteristics detected by the insoles were also comparable to lab-based measurements ­ in other words, the insoles provided similar type and quality of data as a gait laboratory. Further analysis revealed that information from just a single step could reveal additional information about the subject's walking. These results support the use of wearable devices as a simple and relatively inexpensive way to measure gait in everyday clinical practice, without the need for specialised laboratories and visits to the doctor's office. Although the digital insoles will require further analytical and clinical study before they can be widely used, Wipperman et al. hope they will eventually make monitoring muscle and joint conditions easier and more affordable.

Kerfuffle: a web tool for multi-species gene colocalization analysis.

BACKGROUND: The evolutionary pressures that underlie the large-scale functional organization of the genome are not well understood in eukaryotes. Recent evidence suggests that functionally similar genes may colocalize (cluster) in the eukaryotic genome, suggesting the role of chromatin-level gene regulation in shaping the physical distribution of coordinated genes. However, few of the bioinformatic tools currently available allow for a systematic study of gene colocalization across several, evolutionarily distant species. Furthermore, most tools require the user to input manually curated lists of gene position information, DNA sequence or gene homology relations between species. With the growing number of sequenced genomes, there is a need to provide new comparative genomics tools that can address the analysis of multi-species gene colocalization. RESULTS: Kerfuffle is a web tool designed to help discover, visualize, and quantify the physical organization of genomes by identifying significant gene colocalization and conservation across the assembled genomes of available species (currently up to 47, from humans to worms). Kerfuffle only requires the user to specify a list of human genes and the names of other species of interest. Without further input from the user, the software queries the e!Ensembl BioMart server to obtain positional information and discovers homology relations in all genes and species specified. Using this information, Kerfuffle performs a multi-species clustering analysis, presents downloadable lists of clustered genes, performs Monte Carlo statistical significance calculations, estimates how conserved gene clusters are across species, plots histograms and interactive graphs, allows users to save their queries, and generates a downloadable visualization of the clusters using the Circos software. These analyses may be used to further explore the functional roles of gene clusters by interrogating the enriched molecular pathways associated with each cluster. CONCLUSIONS: Kerfuffle is a new, easy-to-use and publicly available tool to aid our understanding of functional genomics and comparative genomics. This software allows for flexibility and quick investigations of a user-defined set of genes, and the results may be saved online for further analysis. Kerfuffle is freely available at http://atwallab.org/kerfuffle, is implemented in JavaScript (using jQuery and jsCharts libraries) and PHP 5.2, runs on an Apache server, and stores data in flat files and an SQLite database.

Identification of pan-cancer/testis genes and validation of therapeutic targeting in triple-negative breast cancer: Lin28a- and Siglece-based vaccination induces anti-tumor immunity and inhibits metastasis.

Background: Cancer-testis (CT) genes are targets for tumor antigen-specific immunotherapy given that their expression is normally restricted to the immune-privileged testis in healthy individuals with aberrant expression in tumor tissues. While they represent targetable germ-tissue antigens and play important functional roles in tumorigenesis, there is currently no standardized approach for identifying clinically relevant CT genes. Optimized algorithms and validated methods for accurate prediction of reliable CT antigens with high immunogenicity are also lacking. Methods: Sequencing data from the Genotype-Tissue Expression (GTEx) and The Genomic Data Commons (GDC) databases was utilized for the development of a bioinformatic pipeline to identify CT exclusive genes. A CT germness score was calculated based on the number of CT genes expressed within a tumor type and their degree of expression. The impact of tumor germness with clinical outcome was evaluated using healthy GTEx and GDC tumor samples. We then used a triple-negative breast cancer mouse model to develop and test an algorithm that predicts epitope immunogenicity based on the identification of germline sequences with strong MHCI and MHCII binding affinities. Germline sequences for CT genes were synthesized as long synthetic peptide vaccines and tested in the 4T1 triple-negative model of invasive breast cancer with Poly(I:C) adjuvant. Vaccine immunogenicity was determined by flow cytometric analysis of in vitro and in vivo T cell responses. Primary tumor growth and lung metastasis was evaluated by histopathology, flow cytometry and colony formation assay. Results: We developed a new bioinformatic pipeline to reliably identify CT exclusive genes as immunogenic targets for immunotherapy. We identified CT genes that are exclusively expressed within the testis, lack detectable thymic expression, and are significantly expressed in multiple tumor types. High tumor germness correlated with tumor progression but not with tumor mutation burden, supporting CT antigens as appealing targets in low mutation burden tumors. Importantly, tumor germness also correlated with markers of anti-tumor immunity. Vaccination of 4T1 tumor bearing mice with Siglece and Lin28a antigens resulted in increased T cell anti-tumor immunity and reduced primary tumor growth and lung metastases. Conclusion: Our results present a novel strategy for the identification of highly immunogenic CT antigens for the development of targeted vaccines that induce anti-tumor immunity and inhibit metastasis.

Identification of pan-cancer/testis genes and validation of therapeutic targeting in triple-negative breast cancer: Lin28a-based and Siglece-based vaccination induces antitumor immunity and inhibits metastasis.

BACKGROUND: Cancer-testis (CT) genes are targets for tumor antigen-specific immunotherapy given that their expression is normally restricted to the immune-privileged testis in healthy individuals with aberrant expression in tumor tissues. While they represent targetable germ tissue antigens and play important functional roles in tumorigenesis, there is currently no standardized approach for identifying clinically relevant CT genes. Optimized algorithms and validated methods for accurate prediction of reliable CT antigens (CTAs) with high immunogenicity are also lacking. METHODS: Sequencing data from the Genotype-Tissue Expression (GTEx) and The Genomic Data Commons (GDC) databases was used for the development of a bioinformatic pipeline to identify CT exclusive genes. A CT germness score was calculated based on the number of CT genes expressed within a tumor type and their degree of expression. The impact of tumor germness on clinical outcome was evaluated using healthy GTEx and GDC tumor samples. We then used a triple-negative breast cancer mouse model to develop and test an algorithm that predicts epitope immunogenicity based on the identification of germline sequences with strong major histocompatibility complex class I (MHCI) and MHCII binding affinities. Germline sequences for CT genes were synthesized as long synthetic peptide vaccines and tested in the 4T1 triple-negative model of invasive breast cancer with Poly(I:C) adjuvant. Vaccine immunogenicity was determined by flow cytometric analysis of in vitro and in vivo T-cell responses. Primary tumor growth and lung metastasis was evaluated by histopathology, flow cytometry and colony formation assay. RESULTS: We developed a new bioinformatic pipeline to reliably identify CT exclusive genes as immunogenic targets for immunotherapy. We identified CT genes that are exclusively expressed within the testis, lack detectable thymic expression, and are significantly expressed in multiple tumor types. High tumor germness correlated with tumor progression but not with tumor mutation burden, supporting CTAs as appealing targets in low mutation burden tumors. Importantly, tumor germness also correlated with markers of antitumor immunity. Vaccination of 4T1 tumor-bearing mice with Siglece and Lin28a antigens resulted in increased T-cell antitumor immunity and reduced primary tumor growth and lung metastases. CONCLUSION: Our results present a novel strategy for the identification of highly immunogenic CTAs for the development of targeted vaccines that induce antitumor immunity and inhibit metastasis.

TRAPnSeq allows high-throughput profiling of antigen-specific antibody-secreting cells.

Following activation by cognate antigen, B cells undergo fine-tuning of their antigen receptors and may ultimately differentiate into antibody-secreting cells (ASCs). While antigen-specific B cells that express surface receptors (B cell receptors [BCRs]) can be readily cloned and sequenced following flow sorting, antigen-specific ASCs that lack surface BCRs cannot be easily profiled. Here, we report an approach, TRAPnSeq (antigen specificity mapping through immunoglobulin [Ig] secretion TRAP and Sequencing), that allows capture of secreted antibodies on the surface of ASCs, which in turn enables high-throughput screening of single ASCs against large antigen panels. This approach incorporates flow cytometry, standard microfluidic platforms, and DNA-barcoding technologies to characterize antigen-specific ASCs through single-cell V(D)J, RNA, and antigen barcode sequencing. We show the utility of TRAPnSeq by profiling antigen-specific IgG and IgE ASCs from both mice and humans and highlight its capacity to accelerate therapeutic antibody discovery from ASCs.

Age-related gene expression signatures from limb skeletal muscles and the diaphragm in mice and rats reveal common and species-specific changes.

BACKGROUND: As a result of aging, skeletal muscle undergoes atrophy and a decrease in function. This age-related skeletal muscle weakness is known as "sarcopenia". Sarcopenia is part of the frailty observed in humans. In order to discover treatments for sarcopenia, it is necessary to determine appropriate preclinical models and the genes and signaling pathways that change with age in these models. METHODS AND RESULTS: To understand the changes in gene expression that occur as a result of aging in skeletal muscles, we generated a multi-time-point gene expression signature throughout the lifespan of mice and rats, as these are the most commonly used species in preclinical research and intervention testing. Gastrocnemius, tibialis anterior, soleus, and diaphragm muscles from male and female C57Bl/6J mice and male Sprague Dawley rats were analyzed at ages 6, 12, 18, 21, 24, and 27 months, plus an additional 9-month group was used for rats. More age-related genes were identified in rat skeletal muscles compared with mice; this was consistent with the finding that rat muscles undergo more robust age-related decline in mass. In both species, pathways associated with innate immunity and inflammation linearly increased with age. Pathways linked with extracellular matrix remodeling were also universally downregulated. Interestingly, late downregulated pathways were exclusively found in the rat limb muscles and these were linked to metabolism and mitochondrial respiration; this was not seen in the mouse. CONCLUSIONS: This extensive, side-by-side transcriptomic profiling shows that the skeletal muscle in rats is impacted more by aging compared with mice, and the pattern of decline in the rat may be more representative of the human. The observed changes point to potential therapeutic interventions to avoid age-related decline in skeletal muscle function.

Microenvironment-dependent growth of Sezary cells in humanized IL-15 mice.

Sezary syndrome (SS) is a rare, aggressive leukemic variant of cutaneous T-cell lymphoma (CTCL) that lacks adequate therapeutic options and representative small-animal models. Here, we demonstrate that IL-15 is a critical CTCL growth factor. Importantly, an immunodeficient knock-in mouse model genetically engineered to express human IL-15 uniquely supported the growth of SS patient samples relative to conventional immunodeficient mouse strains. SS patient-derived xenograft (PDX) models recapacitated key pathological features of the human disease, including skin infiltration and spread of leukemic cells to the periphery, and maintained the dependence on human IL-15 upon serial in vivo passaging. Detailed molecular characterization of the engrafted cells by single-cell transcriptomic analysis revealed congruent neoplastic gene expression signatures but distinct clonal engraftment patterns. Overall, we document an important dependence of Sezary cell survival and proliferation on IL-15 signaling and the utility of immunodeficient humanized IL-15 mice as hosts for SS - and potentially other T and NK cell-derived hematologic malignancies - PDX model generation. Furthermore, these studies advocate the thorough molecular understanding of the resultant PDX models to maximize their translational impact.

Structure of the Inmazeb cocktail and resistance to Ebola virus escape.

Monoclonal antibodies can provide important pre- or post-exposure protection against infectious disease for those not yet vaccinated or in individuals that fail to mount a protective immune response after vaccination. Inmazeb (REGN-EB3), a three-antibody cocktail against Ebola virus, lessened disease and improved survival in a controlled trial. Here, we present the cryo-EM structure at 3.1 Å of the Ebola virus glycoprotein, determined without symmetry averaging, in a simultaneous complex with the antibodies in the Inmazeb cocktail. This structure allows the modeling of previously disordered portions of the glycoprotein glycan cap, maps the non-overlapping epitopes of Inmazeb, and illuminates the basis for complementary activities and residues critical for resistance to escape by these and other clinically relevant antibodies. We further provide direct evidence that Inmazeb protects against the rapid emergence of escape mutants, whereas monotherapies even against conserved epitopes do not, supporting the benefit of a cocktail versus a monotherapy approach.

Intratumoral dendritic cell-CD4+ T helper cell niches enable CD8+ T cell differentiation following PD-1 blockade in hepatocellular carcinoma.

Despite no apparent defects in T cell priming and recruitment to tumors, a large subset of T cell rich tumors fail to respond to immune checkpoint blockade (ICB). We leveraged a neoadjuvant anti-PD-1 trial in patients with hepatocellular carcinoma (HCC), as well as additional samples collected from patients treated off-label, to explore correlates of response to ICB within T cell-rich tumors. We show that ICB response correlated with the clonal expansion of intratumoral CXCL13+CH25H+IL-21+PD-1+CD4+ T helper cells ("CXCL13+ TH") and Granzyme K+ PD-1+ effector-like CD8+ T cells, whereas terminally exhausted CD39hiTOXhiPD-1hiCD8+ T cells dominated in nonresponders. CD4+ and CD8+ T cell clones that expanded post-treatment were found in pretreatment biopsies. Notably, PD-1+TCF-1+ (Progenitor-exhausted) CD8+ T cells shared clones mainly with effector-like cells in responders or terminally exhausted cells in nonresponders, suggesting that local CD8+ T cell differentiation occurs upon ICB. We found that these Progenitor CD8+ T cells interact with CXCL13+ TH within cellular triads around dendritic cells enriched in maturation and regulatory molecules, or "mregDC". These results suggest that discrete intratumoral niches that include mregDC and CXCL13+ TH control the differentiation of tumor-specific Progenitor exhasuted CD8+ T cells following ICB.