Multi-Omics Resolves a Sharp Disease-State Shift between Mild and Moderate COVID-19.

We present an integrated analysis of the clinical measurements, immune cells, and plasma multi-omics of 139 COVID-19 patients representing all levels of disease severity, from serial blood draws collected during the first week of infection following diagnosis. We identify a major shift between mild and moderate disease, at which point elevated inflammatory signaling is accompanied by the loss of specific classes of metabolites and metabolic processes. Within this stressed plasma environment at moderate disease, multiple unusual immune cell phenotypes emerge and amplify with increasing disease severity. We condensed over 120,000 immune features into a single axis to capture how different immune cell classes coordinate in response to SARS-CoV-2. This immune-response axis independently aligns with the major plasma composition changes, with clinical metrics of blood clotting, and with the sharp transition between mild and moderate disease. This study suggests that moderate disease may provide the most effective setting for therapeutic intervention.

Multiomic blood correlates of genetic risk identify presymptomatic disease alterations.

Transitions from health to disease are characterized by dysregulation of biological networks under the influence of genetic and environmental factors, often over the course of years to decades before clinical symptoms appear. Understanding these dynamics has important implications for preventive medicine. However, progress has been hindered both by the difficulty of identifying individuals who will eventually go on to develop a particular disease and by the inaccessibility of most disease-relevant tissues in living individuals. Here we developed an alternative approach using polygenic risk scores (PRSs) based on genome-wide association studies (GWAS) for 54 diseases and complex traits coupled with multiomic profiling and found that these PRSs were associated with 766 detectable alterations in proteomic, metabolomic, and standard clinical laboratory measurements (clinical labs) from blood plasma across several thousand mostly healthy individuals. We recapitulated a variety of known relationships (e.g., glutamatergic neurotransmission and inflammation with depression, IL-33 with asthma) and found associations directly suggesting therapeutic strategies (e.g., Ω-6 supplementation and IL-13 inhibition for amyotrophic lateral sclerosis) and influences on longevity (leukemia inhibitory factor, ceramides). Analytes altered in high-genetic-risk individuals showed concordant changes in disease cases, supporting the notion that PRS-associated analytes represent presymptomatic disease alterations. Our results provide insights into the molecular pathophysiology of a range of traits and suggest avenues for the prevention of health-to-disease transitions.

Longitudinal analysis reveals transition barriers between dominant ecological states in the gut microbiome.

The Pioneer 100 Wellness Project involved quantitatively profiling 108 participants' molecular physiology over time, including genomes, gut microbiomes, blood metabolomes, blood proteomes, clinical chemistries, and data from wearable devices. Here, we present a longitudinal analysis focused specifically around the Pioneer 100 gut microbiomes. We distinguished a subpopulation of individuals with reduced gut diversity, elevated relative abundance of the genus Prevotella, and reduced levels of the genus Bacteroides We found that the relative abundances of Bacteroides and Prevotella were significantly correlated with certain serum metabolites, including omega-6 fatty acids. Primary dimensions in distance-based redundancy analysis of clinical chemistries explained 18.5% of the variance in bacterial community composition, and revealed a Bacteroides/Prevotella dichotomy aligned with inflammation and dietary markers. Finally, longitudinal analysis of gut microbiome dynamics within individuals showed that direct transitions between Bacteroides-dominated and Prevotella-dominated communities were rare, suggesting the presence of a barrier between these states. One implication is that interventions seeking to transition between Bacteroides- and Prevotella-dominated communities will need to identify permissible paths through ecological state-space that circumvent this apparent barrier.

Characteristics and Factors Associated With Coronavirus Disease 2019 Infection, Hospitalization, and Mortality Across Race and Ethnicity.

BACKGROUND: Data on the characteristics of coronavirus disease 2019 (COVID-19) patients disaggregated by race/ethnicity remains limited. We evaluated the sociodemographic and clinical characteristics of patients across racial/ethnic groups and assessed their associations with COVID-19 outcomes. METHODS: This retrospective cohort study examined 629 953 patients tested for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in a large health system spanning California, Oregon, and Washington between March 1 and December 31, 2020. Sociodemographic and clinical characteristics were obtained from electronic health records. Odds of SARS-CoV-2 infection, COVID-19 hospitalization, and in-hospital death were assessed with multivariate logistic regression. RESULTS: A total of 570 298 patients with known race/ethnicity were tested for SARS-CoV-2, of whom 27.8% were non-White minorities: 54 645 individuals tested positive, with minorities representing 50.1%. Hispanics represented 34.3% of infections but only 13.4% of tests. Although generally younger than White patients, Hispanics had higher rates of diabetes but fewer other comorbidities. A total of 8536 patients were hospitalized and 1246 died, of whom 56.1% and 54.4% were non-White, respectively. Racial/ethnic distributions of outcomes across the health system tracked with state-level statistics. Increased odds of testing positive and hospitalization were associated with all minority races/ethnicities. Hispanic patients also exhibited increased morbidity, and Hispanic race/ethnicity was associated with in-hospital mortality (odds ratio [OR], 1.39; 95% confidence interval [CI], 1.14-1.70). CONCLUSION: Major healthcare disparities were evident, especially among Hispanics who tested positive at a higher rate, required excess hospitalization and mechanical ventilation, and had higher odds of in-hospital mortality despite younger age. Targeted, culturally responsive interventions and equitable vaccine development and distribution are needed to address the increased risk of poorer COVID-19 outcomes among minority populations.

Correction to: Small Molecule KRAS Agonist for Mutant KRAS Cancer Therapy.

An amendment to this paper has been published and can be accessed via the original article.

Small Molecule KRAS Agonist for Mutant KRAS Cancer Therapy.

BACKGROUND: Lung cancer patients with KRAS mutation(s) have a poor prognosis due in part to the development of resistance to currently available therapeutic interventions. Development of a new class of anticancer agents that directly targets KRAS may provide a more attractive option for the treatment of KRAS-mutant lung cancer. RESULTS: Here we identified a small molecule KRAS agonist, KRA-533, that binds the GTP/GDP-binding pocket of KRAS. In vitro GDP/GTP exchange assay reveals that KRA-533 activates KRAS by preventing the cleavage of GTP into GDP, leading to the accumulation of GTP-KRAS, an active form of KRAS. Treatment of human lung cancer cells with KRA-533 resulted in increased KRAS activity and suppression of cell growth. Lung cancer cell lines with KRAS mutation were relatively more sensitive to KRA-533 than cell lines without KRAS mutation. Mutating one of the hydrogen-bonds among the KRA-533 binding amino acids in KRAS (mutant K117A) resulted in failure of KRAS to bind KRA-533. KRA-533 had no effect on the activity of K117A mutant KRAS, suggesting that KRA-533 binding to K117 is required for KRA-533 to enhance KRAS activity. Intriguingly, KRA-533-mediated KRAS activation not only promoted apoptosis but also autophagic cell death. In mutant KRAS lung cancer xenografts and genetically engineered mutant KRAS-driven lung cancer models, KRA-533 suppressed malignant growth without significant toxicity to normal tissues. CONCLUSIONS: The development of this KRAS agonist as a new class of anticancer drug offers a potentially effective strategy for the treatment of lung cancer with KRAS mutation and/or mutant KRAS-driven lung cancer.

A bony fish immunological receptor of the NITR multigene family mediates allogeneic recognition.

Novel immune-type receptors (NITRs) comprise an exceptionally large, diversified family of activating and inhibitory receptors that has been identified in bony fish. Here, we characterized the structure of an activating NITR that is expressed by a cytotoxic natural killer (NK)-like cell line and that specifically binds an allogeneic B cell target. A single amino acid residue within the NITR immunoglobulin variable (V)-type domain accounts for specificity of the interaction. Structures solved by X-ray crystallography revealed that the V-type domains of NITRs form homodimers resembling rearranging antigen-binding receptor heterodimers. CDR1 elements of both subunits of NITR dimers form ligand-binding surfaces that determine specificity for the nonself target. In the evolution of immune function, it appears that a specific NK type of innate recognition may be mediated by a complex germline multigene family of V structures resembling those that are somatically diversified in adaptive immunological responses.

Multi-study integration of brain cancer transcriptomes reveals organ-level molecular signatures.

We utilized abundant transcriptomic data for the primary classes of brain cancers to study the feasibility of separating all of these diseases simultaneously based on molecular data alone. These signatures were based on a new method reported herein--Identification of Structured Signatures and Classifiers (ISSAC)--that resulted in a brain cancer marker panel of 44 unique genes. Many of these genes have established relevance to the brain cancers examined herein, with others having known roles in cancer biology. Analyses on large-scale data from multiple sources must deal with significant challenges associated with heterogeneity between different published studies, for it was observed that the variation among individual studies often had a larger effect on the transcriptome than did phenotype differences, as is typical. For this reason, we restricted ourselves to studying only cases where we had at least two independent studies performed for each phenotype, and also reprocessed all the raw data from the studies using a unified pre-processing pipeline. We found that learning signatures across multiple datasets greatly enhanced reproducibility and accuracy in predictive performance on truly independent validation sets, even when keeping the size of the training set the same. This was most likely due to the meta-signature encompassing more of the heterogeneity across different sources and conditions, while amplifying signal from the repeated global characteristics of the phenotype. When molecular signatures of brain cancers were constructed from all currently available microarray data, 90% phenotype prediction accuracy, or the accuracy of identifying a particular brain cancer from the background of all phenotypes, was found. Looking forward, we discuss our approach in the context of the eventual development of organ-specific molecular signatures from peripheral fluids such as the blood.

Generally-healthy individuals with aberrant bowel movement frequencies show enrichment for microbially-derived blood metabolites associated with reduced kidney function.

Bowel movement frequency (BMF) has been linked to changes in the composition of the human gut microbiome and to many chronic conditions, like metabolic disorders, neurodegenerative diseases, chronic kidney disease (CKD), and other intestinal pathologies like irritable bowel syndrome and inflammatory bowel disease. Lower BMF (constipation) can lead to compromised intestinal barrier integrity and a switch from saccharolytic to proteolytic fermentation within the microbiota, giving rise to microbially-derived toxins that may make their way into circulation and cause damage to organ systems. However, the connections between BMF, gut microbial metabolism, and the early-stage development and progression of chronic disease remain underexplored. Here, we examined the phenotypic impact of BMF variation in a cohort of generally-healthy, community dwelling adults with detailed clinical, lifestyle, and multi-omic data. We showed significant differences in microbially-derived blood plasma metabolites, gut bacterial genera, clinical chemistries, and lifestyle factors across BMF groups that have been linked to inflammation, cardiometabolic health, liver function, and CKD severity and progression. We found that the higher plasma levels of 3-indoxyl sulfate (3-IS), a microbially-derived metabolite associated with constipation, was in turn negatively associated with estimated glomerular filtration rate (eGFR), a measure of kidney function. Causal mediation analysis revealed that the effect of BMF on eGFR was significantly mediated by 3-IS. Finally, we identify self-reported diet, lifestyle, and psychological factors associated with BMF variation, which indicate several common-sense strategies for mitigating constipation and diarrhea. Overall, we suggest that aberrant BMF is an underappreciated risk factor in the development of chronic diseases, even in otherwise healthy populations.

Aberrant bowel movement frequencies coincide with increased microbe-derived blood metabolites associated with reduced organ function.

Bowel movement frequency (BMF) directly impacts the gut microbiota and is linked to diseases like chronic kidney disease or dementia. In particular, prior work has shown that constipation is associated with an ecosystem-wide switch from fiber fermentation and short-chain fatty acid production to more detrimental protein fermentation and toxin production. Here, we analyze multi-omic data from generally healthy adults to see how BMF affects their molecular phenotypes, in a pre-disease context. Results show differential abundances of gut microbial genera, blood metabolites, and variation in lifestyle factors across BMF categories. These differences relate to inflammation, heart health, liver function, and kidney function. Causal mediation analysis indicates that the association between lower BMF and reduced kidney function is partially mediated by the microbially derived toxin 3-indoxyl sulfate (3-IS). This result, in a generally healthy context, suggests that the accumulation of microbiota-derived toxins associated with abnormal BMF precede organ damage and may be drivers of chronic, aging-related diseases.

AUREA: an open-source software system for accurate and user-friendly identification of relative expression molecular signatures.

BACKGROUND: Public databases such as the NCBI Gene Expression Omnibus contain extensive and exponentially increasing amounts of high-throughput data that can be applied to molecular phenotype characterization. Collectively, these data can be analyzed for such purposes as disease diagnosis or phenotype classification. One family of algorithms that has proven useful for disease classification is based on relative expression analysis and includes the Top-Scoring Pair (TSP), k-Top-Scoring Pairs (k-TSP), Top-Scoring Triplet (TST) and Differential Rank Conservation (DIRAC) algorithms. These relative expression analysis algorithms hold significant advantages for identifying interpretable molecular signatures for disease classification, and have been implemented previously on a variety of computational platforms with varying degrees of usability. To increase the user-base and maximize the utility of these methods, we developed the program AUREA (Adaptive Unified Relative Expression Analyzer)-a cross-platform tool that has a consistent application programming interface (API), an easy-to-use graphical user interface (GUI), fast running times and automated parameter discovery. RESULTS: Herein, we describe AUREA, an efficient, cohesive, and user-friendly open-source software system that comprises a suite of methods for relative expression analysis. AUREA incorporates existing methods, while extending their capabilities and bringing uniformity to their interfaces. We demonstrate that combining these algorithms and adaptively tuning parameters on the training sets makes these algorithms more consistent in their performance and demonstrate the effectiveness of our adaptive parameter tuner by comparing accuracy across diverse datasets. CONCLUSIONS: We have integrated several relative expression analysis algorithms and provided a unified interface for their implementation while making data acquisition, parameter fixing, data merging, and results analysis 'point-and-click' simple. The unified interface and the adaptive parameter tuning of AUREA provide an effective framework in which to investigate the massive amounts of publically available data by both 'in silico' and 'bench' scientists. AUREA can be found at http://price.systemsbiology.net/AUREA/.

Multiomic signatures of body mass index identify heterogeneous health phenotypes and responses to a lifestyle intervention.

Multiomic profiling can reveal population heterogeneity for both health and disease states. Obesity drives a myriad of metabolic perturbations and is a risk factor for multiple chronic diseases. Here we report an atlas of cross-sectional and longitudinal changes in 1,111 blood analytes associated with variation in body mass index (BMI), as well as multiomic associations with host polygenic risk scores and gut microbiome composition, from a cohort of 1,277 individuals enrolled in a wellness program (Arivale). Machine learning model predictions of BMI from blood multiomics captured heterogeneous phenotypic states of host metabolism and gut microbiome composition better than BMI, which was also validated in an external cohort (TwinsUK). Moreover, longitudinal analyses identified variable BMI trajectories for different omics measures in response to a healthy lifestyle intervention; metabolomics-inferred BMI decreased to a greater extent than actual BMI, whereas proteomics-inferred BMI exhibited greater resistance to change. Our analyses further identified blood analyte-analyte associations that were modified by metabolomics-inferred BMI and partially reversed in individuals with metabolic obesity during the intervention. Taken together, our findings provide a blood atlas of the molecular perturbations associated with changes in obesity status, serving as a resource to quantify metabolic health for predictive and preventive medicine.

The top-scoring 'N' algorithm: a generalized relative expression classification method from small numbers of biomolecules.

BACKGROUND: Relative expression algorithms such as the top-scoring pair (TSP) and the top-scoring triplet (TST) have several strengths that distinguish them from other classification methods, including resistance to overfitting, invariance to most data normalization methods, and biological interpretability. The top-scoring 'N' (TSN) algorithm is a generalized form of other relative expression algorithms which uses generic permutations and a dynamic classifier size to control both the permutation and combination space available for classification. RESULTS: TSN was tested on nine cancer datasets, showing statistically significant differences in classification accuracy between different classifier sizes (choices of N). TSN also performed competitively against a wide variety of different classification methods, including artificial neural networks, classification trees, discriminant analysis, k-Nearest neighbor, naïve Bayes, and support vector machines, when tested on the Microarray Quality Control II datasets. Furthermore, TSN exhibits low levels of overfitting on training data compared to other methods, giving confidence that results obtained during cross validation will be more generally applicable to external validation sets. CONCLUSIONS: TSN preserves the strengths of other relative expression algorithms while allowing a much larger permutation and combination space to be explored, potentially improving classification accuracies when fewer numbers of measured features are available.

Graphics processing unit implementations of relative expression analysis algorithms enable dramatic computational speedup.

SUMMARY: The top-scoring pair (TSP) and top-scoring triplet (TST) algorithms are powerful methods for classification from expression data, but analysis of all combinations across thousands of human transcriptome samples is computationally intensive, and has not yet been achieved for TST. Implementation of these algorithms for the graphics processing unit results in dramatic speedup of two orders of magnitude, greatly increasing the searchable combinations and accelerating the pace of discovery. AVAILABILITY: http://www.igb.illinois.edu/labs/price/downloads/.

A shift in the salt bridge interaction of residues D620 and E621 mediates the constitutive activation of Jak2-H538Q/K539L.

Jak2 mutations in the exon 14 and exon 12 regions that cause constitutive activation have been associated with myeloproliferative neoplasms. We have previously shown that a pi stacking interaction between F617 and F595 is important for the constitutive activation of Jak2-V617F (Gnanasambandan et al., Biochemistry 49:9972-9984, 2010). Here, using a combination of molecular dynamics (MD) simulations and in vitro mutagenesis, we studied the molecular mechanism for the constitutive activation of the Jak2 exon 12 mutation, H538Q/K539L. The activation levels of Jak2-H538Q/K539L were found to be similar to that of Jak2-V617F, and Jak2-H538Q/K539L/V617F. Data from MD simulations indicated a shift in the salt bridge interactions of D620 and E621 with K539 in Jak2-WT to R541 in Jak2-H538Q/K539L. When compared to Jak2-WT, K539A mutation resulted in increased activation, while K539D or K539E mutations diminished Jak2 activation by 50 %. In the context of Jak2-H538Q/K539L, R541A mutation reduced its activation by 50 %, while R541D and R541E mutations returned its activation levels to that of Jak2-WT. Collectively, these results indicate that a shift in the salt bridge interaction of D620 and E621 with K539 in Jak2-WT to R541 in Jak2-H538Q/K539L is critical for constitutive activation of this Jak2 exon 12 mutant.

Personal Dense Dynamic Data Clouds Connect Systems Biomedicine to Scientific Wellness.

The dramatic convergence of molecular biology, genomics, proteomics, metabolomics, bioinformatics, and artificial intelligence has provided a substrate for deep understanding of the biological basis of health and disease. Systems biology is a holistic, dynamic, integrative, cross-disciplinary approach to biological complexity that embraces experimentation, technology, computation, and clinical translation. Systems Medicine integrates genome analyses and longitudinal deep phenotyping with biological pathways and networks to understand mechanisms of disease, identify relevant blood biomarkers, define druggable molecular targets, and enhance the maintenance or restoration of wellness. Two programs initiated our understanding of data-driven population-based wellness. The Pioneer 100 Study of Scientific Wellness and the much larger Arivale commercial program that followed had two spectacular results: demonstrating the feasibility and utility of collecting longitudinal multiomic data, and then generating dense, dynamic data clouds for each individual to utilize actionable metrics for promoting health and preventing disease when combined with personalized coaching. Future developments in these domains will enable better population health and personal, preventive, predictive, participatory (P4) health care.

Genome-microbiome interplay provides insight into the determinants of the human blood metabolome.

Variation in the blood metabolome is intimately related to human health. However, few details are known about the interplay between genetics and the microbiome in explaining this variation on a metabolite-by-metabolite level. Here, we perform analyses of variance for each of 930 blood metabolites robustly detected across a cohort of 1,569 individuals with paired genomic and microbiome data while controlling for a number of relevant covariates. We find that 595 (64%) of these blood metabolites are significantly associated with either host genetics or the gut microbiome, with 69% of these associations driven solely by the microbiome, 15% driven solely by genetics and 16% under hybrid genome-microbiome control. Additionally, interaction effects, where a metabolite-microbe association is specific to a particular genetic background, are quite common, albeit with modest effect sizes. This knowledge will help to guide targeted interventions designed to alter the composition of the human blood metabolome.

Heterogeneity in statin responses explained by variation in the human gut microbiome.

BACKGROUND: Statins remain one of the most prescribed medications worldwide. While effective in decreasing atherosclerotic cardiovascular disease risk, statin use is associated with adverse effects for a subset of patients, including disrupted metabolic control and increased risk of type 2 diabetes. METHODS: We investigated the potential role of the gut microbiome in modifying patient responses to statin therapy across two independent cohorts (discovery n = 1,848, validation n = 991). Microbiome composition was assessed in these cohorts using stool 16S rRNA amplicon and shotgun metagenomic sequencing, respectively. Microbiome associations with markers of statin on-target and adverse effects were tested via a covariate-adjusted interaction analysis framework, utilizing blood metabolomics, clinical laboratory tests, genomics, and demographics data. FINDINGS: The hydrolyzed substrate for 3-hydroxy-3-methylglutarate-coenzyme-A (HMG-CoA) reductase, HMG, emerged as a promising marker for statin on-target effects in cross-sectional cohorts. Plasma HMG levels reflected both statin therapy intensity and known genetic markers for variable statin responses. Through exploring gut microbiome associations between blood-derived measures of statin effectiveness and adverse metabolic effects of statins, we find that heterogeneity in statin responses was consistently associated with variation in the gut microbiome across two independent cohorts. A Bacteroides-enriched and diversity-depleted gut microbiome was associated with more intense statin responses, both in terms of on-target and adverse effects. CONCLUSIONS: With further study and refinement, gut microbiome monitoring may help inform precision statin treatment. FUNDING: This research was supported by the M.J. Murdock Charitable Trust, WRF, NAM Catalyst Award, and NIH grant U19AG023122 awarded by the NIA.

Manifestations of Alzheimer's disease genetic risk in the blood are evident in a multiomic analysis in healthy adults aged 18 to 90.

Genetics play an important role in late-onset Alzheimer's Disease (AD) etiology and dozens of genetic variants have been implicated in AD risk through large-scale GWAS meta-analyses. However, the precise mechanistic effects of most of these variants have yet to be determined. Deeply phenotyped cohort data can reveal physiological changes associated with genetic risk for AD across an age spectrum that may provide clues to the biology of the disease. We utilized over 2000 high-quality quantitative measurements obtained from blood of 2831 cognitively normal adult clients of a consumer-based scientific wellness company, each with CLIA-certified whole-genome sequencing data. Measurements included: clinical laboratory blood tests, targeted chip-based proteomics, and metabolomics. We performed a phenome-wide association study utilizing this diverse blood marker data and 25 known AD genetic variants and an AD-specific polygenic risk score (PGRS), adjusting for sex, age, vendor (for clinical labs), and the first four genetic principal components; sex-SNP interactions were also assessed. We observed statistically significant SNP-analyte associations for five genetic variants after correction for multiple testing (for SNPs in or near NYAP1, ABCA7, INPP5D, and APOE), with effects detectable from early adulthood. The ABCA7 SNP and the APOE2 and APOE4 encoding alleles were associated with lipid variability, as seen in previous studies; in addition, six novel proteins were associated with the e2 allele. The most statistically significant finding was between the NYAP1 variant and PILRA and PILRB protein levels, supporting previous functional genomic studies in the identification of a putative causal variant within the PILRA gene. We did not observe associations between the PGRS and any analyte. Sex modified the effects of four genetic variants, with multiple interrelated immune-modulating effects associated with the PICALM variant. In post-hoc analysis, sex-stratified GWAS results from an independent AD case-control meta-analysis supported sex-specific disease effects of the PICALM variant, highlighting the importance of sex as a biological variable. Known AD genetic variation influenced lipid metabolism and immune response systems in a population of non-AD individuals, with associations observed from early adulthood onward. Further research is needed to determine whether and how these effects are implicated in early-stage biological pathways to AD. These analyses aim to complement ongoing work on the functional interpretation of AD-associated genetic variants.