Genome-wide loss of heterozygosity predicts aggressive, treatment-refractory behavior in pituitary neuroendocrine tumors.

Pituitary neuroendocrine tumors (PitNETs) exhibiting aggressive, treatment-refractory behavior are the rare subset that progress after surgery, conventional medical therapies, and an initial course of radiation and are characterized by unrelenting growth and/or metastatic dissemination. Two groups of patients with PitNETs were sequenced: a prospective group of patients (n = 66) who consented to sequencing prior to surgery and a retrospective group (n = 26) comprised of aggressive/higher risk PitNETs. A higher mutational burden and fraction of loss of heterozygosity (LOH) was found in the aggressive, treatment-refractory PitNETs compared to the benign tumors (p = 1.3 × 10-10 and p = 8.5 × 10-9, respectively). Within the corticotroph lineage, a characteristic pattern of recurrent chromosomal LOH in 12 specific chromosomes was associated with treatment-refractoriness (occurring in 11 of 14 treatment-refractory versus 1 of 14 benign corticotroph PitNETs, p = 1.7 × 10-4). Across the cohort, a higher fraction of LOH was identified in tumors with TP53 mutations (p = 3.3 × 10-8). A machine learning approach identified loss of heterozygosity as the most predictive variable for aggressive, treatment-refractory behavior, outperforming the most common gene-level alteration, TP53, with an accuracy of 0.88 (95% CI: 0.70-0.96). Aggressive, treatment-refractory PitNETs are characterized by significant aneuploidy due to widespread chromosomal LOH, most prominently in the corticotroph tumors. This LOH predicts treatment-refractoriness with high accuracy and represents a novel biomarker for this poorly defined PitNET category.

Lineage tracing reveals clonal progenitors and long-term persistence of tumor-specific T cells during immune checkpoint blockade.

Paired single-cell RNA and T cell receptor sequencing (scRNA/TCR-seq) has allowed for enhanced resolution of clonal T cell dynamics in cancer. Here, we report a scRNA/TCR-seq analysis of 187,650 T cells from 31 tissue regions, including tumor, adjacent normal tissues, and lymph nodes (LN), from three patients with non-small cell lung cancer after immune checkpoint blockade (ICB). Regions with viable cancer cells are enriched for exhausted CD8+ T cells, regulatory CD4+ T cells (Treg), and follicular helper CD4+ T cells (TFH). Tracking T cell clonotypes across tissues, combined with neoantigen specificity assays, reveals that TFH and tumor-specific exhausted CD8+ T cells are clonally linked to TCF7+SELL+ progenitors in tumor draining LNs, and progressive exhaustion trajectories of CD8+ T, Treg, and TFH cells with proximity to the tumor microenvironment. Finally, longitudinal tracking of tumor-specific CD8+ and CD4+ T cell clones reveals persistence in the peripheral blood for years after ICB therapy.

Pilot study of bempegaldesleukin in combination with nivolumab in patients with metastatic sarcoma.

PD-1 blockade (nivolumab) efficacy remains modest for metastatic sarcoma. In this paper, we present an open-label, non-randomized, non-comparative pilot study of bempegaldesleukin, a CD122-preferential interleukin-2 pathway agonist, with nivolumab in refractory sarcoma at Memorial Sloan Kettering/MD Anderson Cancer Centers (NCT03282344). We report on the primary outcome of objective response rate (ORR) and secondary endpoints of toxicity, clinical benefit, progression-free survival, overall survival, and durations of response/treatment. In 84 patients in 9 histotype cohorts, all patients experienced ≥1 adverse event and treatment-related adverse event; 1 death was possibly treatment-related. ORR was highest in angiosarcoma (3/8) and undifferentiated pleomorphic sarcoma (2/10), meeting predefined endpoints. Results of our exploratory investigation of predictive biomarkers show: CD8 + T cell infiltrates and PD-1 expression correlate with improved ORR; upregulation of immune-related pathways correlate with improved efficacy; Hedgehog pathway expression correlate with resistance. Exploration of this combination in selected sarcomas, and of Hedgehog signaling as a predictive biomarker, warrants further study in larger cohorts.

Genetic interaction effects reveal lipid-metabolic and inflammatory pathways underlying common metabolic disease risks.

BACKGROUND: Common metabolic diseases, including type 2 diabetes, coronary artery disease, and hypertension, arise from disruptions of the body's metabolic homeostasis, with relatively strong contributions from genetic risk factors and substantial comorbidity with obesity. Although genome-wide association studies have revealed many genomic loci robustly associated with these diseases, biological interpretation of such association is challenging because of the difficulty in mapping single-nucleotide polymorphisms (SNPs) onto the underlying causal genes and pathways. Furthermore, common diseases are typically highly polygenic, and conventional single variant-based association testing does not adequately capture potentially important large-scale interaction effects between multiple genetic factors. METHODS: We analyzed moderately sized case-control data sets for type 2 diabetes, coronary artery disease, and hypertension to characterize the genetic risk factors arising from non-additive, collective interaction effects, using a recently developed algorithm (discrete discriminant analysis). We tested associations of genes and pathways with the disease status while including the cumulative sum of interaction effects between all variants contained in each group. RESULTS: In contrast to non-interacting SNP mapping, which produced few genome-wide significant loci, our analysis revealed extensive arrays of pathways, many of which are involved in the pathogenesis of these metabolic diseases but have not been directly identified in genetic association studies. They comprised cell stress and apoptotic pathways for insulin-producing ß-cells in type 2 diabetes, processes covering different atherosclerotic stages in coronary artery disease, and elements of both type 2 diabetes and coronary artery disease risk factors (cell cycle, apoptosis, and hemostasis) associated with hypertension. CONCLUSIONS: Our results support the view that non-additive interaction effects significantly enhance the level of common metabolic disease associations and modify their genetic architectures and that many of the expected genetic factors behind metabolic disease risks reside in smaller genotyping samples in the form of interacting groups of SNPs.

Collective interaction effects associated with mammalian behavioral traits reveal genetic factors connecting fear and hemostasis.

BACKGROUND: Investigation of the genetic architectures that influence the behavioral traits of animals can provide important insights into human neuropsychiatric phenotypes. These traits, however, are often highly polygenic, with individual loci contributing only small effects to the overall association. The polygenicity makes it challenging to explain, for example, the widely observed comorbidity between stress and cardiac disease. METHODS: We present an algorithm for inferring the collective association of a large number of interacting gene variants with a quantitative trait. Using simulated data, we demonstrate that by taking into account the non-uniform distribution of genotypes within a cohort, we can achieve greater power than regression-based methods for high-dimensional inference. RESULTS: We analyzed genome-wide data sets of outbred mice and pet dogs, and found neurobiological pathways whose associations with behavioral traits arose primarily from interaction effects: γ-carboxylated coagulation factors and downstream neuronal signaling were highly associated with conditioned fear, consistent with our previous finding in human post-traumatic stress disorder (PTSD) data. Prepulse inhibition in mice was associated with serotonin transporter and platelet homeostasis, and noise-induced fear in dogs with hemostasis. CONCLUSIONS: Our findings suggest a novel explanation for the observed comorbidity between PTSD/anxiety and cardiovascular diseases: key coagulation factors modulating hemostasis also regulate synaptic plasticity affecting the learning and extinction of fear.

A strategy for evaluating pathway analysis methods.

BACKGROUND: Researchers have previously developed a multitude of methods designed to identify biological pathways associated with specific clinical or experimental conditions of interest, with the aim of facilitating biological interpretation of high-throughput data. Before practically applying such pathway analysis (PA) methods, we must first evaluate their performance and reliability, using datasets where the pathways perturbed by the conditions of interest have been well characterized in advance. However, such 'ground truths' (or gold standards) are often unavailable. Furthermore, previous evaluation strategies that have focused on defining 'true answers' are unable to systematically and objectively assess PA methods under a wide range of conditions. RESULTS: In this work, we propose a novel strategy for evaluating PA methods independently of any gold standard, either established or assumed. The strategy involves the use of two mutually complementary metrics, recall and discrimination. Recall measures the consistency of the perturbed pathways identified by applying a particular analysis method to an original large dataset and those identified by the same method to a sub-dataset of the original dataset. In contrast, discrimination measures specificity-the degree to which the perturbed pathways identified by a particular method to a dataset from one experiment differ from those identifying by the same method to a dataset from a different experiment. We used these metrics and 24 datasets to evaluate six widely used PA methods. The results highlighted the common challenge in reliably identifying significant pathways from small datasets. Importantly, we confirmed the effectiveness of our proposed dual-metric strategy by showing that previous comparative studies corroborate the performance evaluations of the six methods obtained by our strategy. CONCLUSIONS: Unlike any previously proposed strategy for evaluating the performance of PA methods, our dual-metric strategy does not rely on any ground truth, either established or assumed, of the pathways perturbed by a specific clinical or experimental condition. As such, our strategy allows researchers to systematically and objectively evaluate pathway analysis methods by employing any number of datasets for a variety of conditions.

Collective Genetic Interaction Effects and the Role of Antigen-Presenting Cells in Autoimmune Diseases.

Autoimmune diseases occur when immune cells fail to develop or lose their tolerance toward self and destroy body's own tissues. Both insufficient negative selection of self-reactive T cells and impaired development of regulatory T cells preventing effector cell activation are believed to contribute to autoimmunity. Genetic predispositions center around the major histocompatibility complex (MHC) class II loci involved in antigen presentation, the key determinant of CD4+ T cell activation. Recent studies suggested that variants in the MHC region also exhibit significant non-additive interaction effects. However, collective interactions involving large numbers of single nucleotide polymorphisms (SNPs) contributing to such effects are yet to be characterized. In addition, relatively little is known about the cell-type-specificity of such interactions in the context of cellular pathways. Here, we analyzed type 1 diabetes (T1D) and rheumatoid arthritis (RA) genome-wide association data sets via large-scale, high-performance computations and inferred collective interaction effects involving MHC SNPs using the discrete discriminant analysis. Despite considerable differences in the details of SNP interactions in T1D and RA data, the enrichment pattern of interacting pairs in reference epigenomes was remarkably similar: statistically significant interactions were epigenetically active in cell-type combinations connecting B cells to T cells and intestinal epithelial cells, with both helper and regulatory T cells showing strong disease-associated interactions with B cells. Our results provide direct genetic evidence pointing to the important roles B cells play as antigen-presenting cells toward CD4+ T cells in the context of central and peripheral tolerance. In addition, they are consistent with recent experimental studies suggesting that the repertoire of B cell-specific self-antigens in the thymus are critical to the effective control of corresponding autoimmune activation in peripheral tissues.

Genotype distribution-based inference of collective effects in genome-wide association studies: insights to age-related macular degeneration disease mechanism.

BACKGROUND: Genome-wide association studies provide important insights to the genetic component of disease risks. However, an existing challenge is how to incorporate collective effects of interactions beyond the level of independent single nucleotide polymorphism (SNP) tests. While methods considering each SNP pair separately have provided insights, a large portion of expected heritability may reside in higher-order interaction effects. RESULTS: We describe an inference approach (discrete discriminant analysis; DDA) designed to probe collective interactions while treating both genotypes and phenotypes as random variables. The genotype distributions in case and control groups are modeled separately based on empirical allele frequency and covariance data, whose differences yield disease risk parameters. We compared pairwise tests and collective inference methods, the latter based both on DDA and logistic regression. Analyses using simulated data demonstrated that significantly higher sensitivity and specificity can be achieved with collective inference in comparison to pairwise tests, and with DDA in comparison to logistic regression. Using age-related macular degeneration (AMD) data, we demonstrated two possible applications of DDA. In the first application, a genome-wide SNP set is reduced into a small number (∼100) of variants via filtering and SNP pairs with significant interactions are identified. We found that interactions between SNPs with highest AMD association were epigenetically active in the liver, adipocytes, and mesenchymal stem cells. In the other application, multiple groups of SNPs were formed from the genome-wide data and their relative strengths of association were compared using cross-validation. This analysis allowed us to discover novel collections of loci for which interactions between SNPs play significant roles in their disease association. In particular, we considered pathway-based groups of SNPs containing up to ∼10, 000 variants in each group. In addition to pathways related to complement activation, our collective inference pointed to pathway groups involved in phospholipid synthesis, oxidative stress, and apoptosis, consistent with the AMD pathogenesis mechanism where the dysfunction of retinal pigment epithelium cells plays central roles. CONCLUSIONS: The simultaneous inference of collective interaction effects within a set of SNPs has the potential to reveal novel aspects of disease association.

Quantitative modeling of virus evolutionary dynamics and adaptation in serial passages using empirically inferred fitness landscapes.

We describe a stochastic virus evolution model representing genomic diversification and within-host selection during experimental serial passages under cell culture or live-host conditions. The model incorporates realistic descriptions of the virus genotypes in nucleotide and amino acid sequence spaces, as well as their diversification from error-prone replications. It quantitatively considers factors such as target cell number, bottleneck size, passage period, infection and cell death rates, and the replication rate of different genotypes, allowing for systematic examinations of how their changes affect the evolutionary dynamics of viruses during passages. The relative probability for a viral population to achieve adaptation under a new host environment, quantified by the rate with which a target sequence frequency rises above 50%, was found to be most sensitive to factors related to sequence structure (distance from the wild type to the target) and selection strength (host cell number and bottleneck size). For parameter values representative of RNA viruses, the likelihood of observing adaptations during passages became negligible as the required number of mutations rose above two amino acid sites. We modeled the specific adaptation process of influenza A H5N1 viruses in mammalian hosts by simulating the evolutionary dynamics of H5 strains under the fitness landscape inferred from multiple sequence alignments of H3 proteins. In light of comparisons with experimental findings, we observed that the evolutionary dynamics of adaptation is strongly affected not only by the tendency toward increasing fitness values but also by the accessibility of pathways between genotypes constrained by the genetic code.

Stationary distribution of self-organized states and biological information generation.

Self-organization, where spontaneous orderings occur under driven conditions, is one of the hallmarks of biological systems. We consider a statistical mechanical treatment of the biased distribution of such organized states, which become favored as a result of their catalytic activity under chemical driving forces. A generalization of the equilibrium canonical distribution describes the stationary state, which can be used to model shifts in conformational ensembles sampled by an enzyme in working conditions. The basic idea is applied to the process of biological information generation from random sequences of heteropolymers, where unfavorable Shannon entropy is overcome by the catalytic activities of selected genes. The ordering process is demonstrated with the genetic distance to a genotype with high catalytic activity as an order parameter. The resulting free energy can have multiple minima, corresponding to disordered and organized phases with first-order transitions between them.