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.

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.

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.

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.

Integrated analysis of plasma and single immune cells uncovers metabolic changes in individuals with COVID-19.

A better understanding of the metabolic alterations in immune cells during severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection may elucidate the wide diversity of clinical symptoms experienced by individuals with coronavirus disease 2019 (COVID-19). Here, we report the metabolic changes associated with the peripheral immune response of 198 individuals with COVID-19 through an integrated analysis of plasma metabolite and protein levels as well as single-cell multiomics analyses from serial blood draws collected during the first week after clinical diagnosis. We document the emergence of rare but metabolically dominant T cell subpopulations and find that increasing disease severity correlates with a bifurcation of monocytes into two metabolically distinct subsets. This integrated analysis reveals a robust interplay between plasma metabolites and cell-type-specific metabolic reprogramming networks that is associated with disease severity and could predict survival.

Baseline Gut Metagenomic Functional Gene Signature Associated with Variable Weight Loss Responses following a Healthy Lifestyle Intervention in Humans.

Recent human feeding studies have shown how the baseline taxonomic composition of the gut microbiome can determine responses to weight loss interventions. However, the functional determinants underlying this phenomenon remain unclear. We report a weight loss response analysis on a cohort of 105 individuals selected from a larger population enrolled in a commercial wellness program, which included healthy lifestyle coaching. Each individual in the cohort had baseline blood metabolomics, blood proteomics, clinical labs, dietary questionnaires, stool 16S rRNA gene sequencing data, and follow-up data on weight change. We generated additional targeted proteomics data on obesity-associated proteins in blood before and after intervention, along with baseline stool metagenomic data, for a subset of 25 individuals who showed the most extreme weight change phenotypes. We built regression models to identify baseline blood, stool, and dietary features associated with weight loss, independent of age, sex, and baseline body mass index (BMI). Many features were independently associated with baseline BMI, but few were independently associated with weight loss. Baseline diet was not associated with weight loss, and only one blood analyte was associated with changes in weight. However, 31 baseline stool metagenomic functional features, including complex polysaccharide and protein degradation genes, stress-response genes, respiration-related genes, and cell wall synthesis genes, along with gut bacterial replication rates, were associated with weight loss responses after controlling for age, sex, and baseline BMI. Together, these results provide a set of compelling hypotheses for how commensal gut microbiota influence weight loss outcomes in humans. IMPORTANCE Recent human feeding studies have shown how the baseline taxonomic composition of the gut microbiome can determine responses to dietary interventions, but the exact functional determinants underlying this phenomenon remain unclear. In this study, we set out to better understand interactions between baseline BMI, metabolic health, diet, gut microbiome functional profiles, and subsequent weight changes in a human cohort that underwent a healthy lifestyle intervention. Overall, our results suggest that the microbiota may influence host weight loss responses through variable bacterial growth rates, dietary energy harvest efficiency, and immunomodulation.

Discovery of Small Molecule Bak Activator for Lung Cancer Therapy.

Rationale: Bak is a major proapoptotic Bcl2 family member and a required molecule for apoptotic cell death. High levels of endogenous Bak were observed in both small cell lung cancer (SCLC) and non-small cell lung cancer (NSCLC) cell lines. Increased Bak expression was correlated with poor prognosis of NSCLC patients, suggesting that Bak protein is an attractive target for lung cancer therapy. The BH3 domain functions as death domain and is required for Bak to initiate apoptotic cell death. Thus, the BH3 domain is attractive target for discovery of Bak agonist. Methods: The BH3 death domain binding pocket (aa75-88) of Bak was chosen as a docking site for screening of small molecule Bak activators using the UCSF DOCK 6.1 program suite and the NCI chemical library (300,000 small molecules) database. The top 500 compounds determined to have the highest affinity for the BH3 domain were obtained from the NCI and tested for cytotoxicity for further screening. We identified a small molecule Bak activator BKA-073 as the lead compound. The binding affinity of BKA-073 with Bak protein was analyzed by isothermal titration calorimetry (ITC) assay. BKA-073-mediated Bak activation via oligomerization was analyzed by a cross-linking with Bis (maleimido) hexane (BMH). Sensitivity of BKA-073 to lung cancer cells in vitro was evaluated by dynamic BH3 profiling (DBP) and apoptotic cell death assay. The potency of BKA-073 alone or in combination with radiotherapy or Bcl2 inhibitor was evaluated in animal models. Results: We found that BKA-073 binds Bak at BH3 domain with high affinity and selectivity. BKA-073/Bak binding promotes Bak oligomerization and mitochondrial priming that activates its proapoptotic function. BKA-073 potently suppresses tumor growth without significant normal tissue toxicity in small cell lung cancer (SCLC) and NSCLC xenografts, patient-derived xenografts, and genetically engineered mouse models of mutant KRAS-driven cancer. Bak accumulates in radioresistant lung cancer cells and BKA-073 reverses radioresistance. Combination of BKA-073 with Bcl-2 inhibitor venetoclax exhibits strong synergy against lung cancer in vivo. Conclusions: Development of small molecule Bak activator may provide a new class of anticancer agents to treat lung cancer.

The geometry of clinical labs and wellness states from deeply phenotyped humans.

Longitudinal multi-omics measurements are highly valuable in studying heterogeneity in health and disease phenotypes. For thousands of people, we have collected longitudinal multi-omics data. To analyze, interpret and visualize this extremely high-dimensional data, we use the Pareto Task Inference (ParTI) method. We find that the clinical labs data fall within a tetrahedron. We then use all other data types to characterize the four archetypes. We find that the tetrahedron comprises three wellness states, defining a wellness triangular plane, and one aberrant health state that captures aspects of commonality in movement away from wellness. We reveal the tradeoffs that shape the data and their hierarchy, and use longitudinal data to observe individual trajectories. We then demonstrate how the movement on the tetrahedron can be used for detecting unexpected trajectories, which might indicate transitions from health to disease and reveal abnormal conditions, even when all individual blood measurements are in the norm.

DUOX2 variants associate with preclinical disturbances in microbiota-immune homeostasis and increased inflammatory bowel disease risk.

A primordial gut-epithelial innate defense response is the release of hydrogen peroxide by dual NADPH oxidase (DUOX). In inflammatory bowel disease (IBD), a condition characterized by an imbalanced gut microbiota-immune homeostasis, DUOX2 isoenzyme is the highest induced gene. Performing multiomic analyses using 2872 human participants of a wellness program, we detected a substantial burden of rare protein-altering DUOX2 gene variants of unknown physiologic significance. We identified a significant association between these rare loss-of-function variants and increased plasma levels of interleukin-17C, which is induced also in mucosal biopsies of patients with IBD. DUOX2-deficient mice replicated increased IL-17C induction in the intestine, with outlier high Il17c expression linked to the mucosal expansion of specific Proteobacteria pathobionts. Integrated microbiota/host gene expression analyses in patients with IBD corroborated IL-17C as a marker for epithelial activation by gram-negative bacteria. Finally, the impact of DUOX2 variants on IL-17C induction provided a rationale for variant stratification in case control studies that substantiated DUOX2 as an IBD risk gene. Thus, our study identifies an association of deleterious DUOX2 variants with a preclinical hallmark of disturbed microbiota-immune homeostasis that appears to precede the manifestation of IBD.

Gut microbiome pattern reflects healthy ageing and predicts survival in humans.

The gut microbiome has important effects on human health, yet its importance in human ageing remains unclear. In the present study, we demonstrate that, starting in mid-to-late adulthood, gut microbiomes become increasingly unique to individuals with age. We leverage three independent cohorts comprising over 9,000 individuals and find that compositional uniqueness is strongly associated with microbially produced amino acid derivatives circulating in the bloodstream. In older age (over ~80 years), healthy individuals show continued microbial drift towards a unique compositional state, whereas this drift is absent in less healthy individuals. The identified microbiome pattern of healthy ageing is characterized by a depletion of core genera found across most humans, primarily Bacteroides. Retaining a high Bacteroides dominance into older age, or having a low gut microbiome uniqueness measure, predicts decreased survival in a 4-year follow-up. Our analysis identifies increasing compositional uniqueness of the gut microbiome as a component of healthy ageing, which is characterized by distinct microbial metabolic outputs in the blood.

Characteristics and Factors Associated with COVID-19 Infection, Hospitalization, and Mortality Across Race and Ethnicity.

Background: Data on the characteristics of COVID-19 patients disaggregated by race/ethnicity remain 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 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: 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. While generally younger than White patients, Hispanics had higher rates of diabetes but fewer other comorbidities. 8,536 patients were hospitalized and 1,246 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 (OR: 1.39 [95% 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.

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.

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.

Untargeted longitudinal analysis of a wellness cohort identifies markers of metastatic cancer years prior to diagnosis.

We analyzed 1196 proteins in longitudinal plasma samples from participants in a commercial wellness program, including samples collected pre-diagnosis from ten cancer patients and 69 controls. For three individuals ultimately diagnosed with metastatic breast, lung, or pancreatic cancer, CEACAM5 was a persistent longitudinal outlier as early as 26.5 months pre-diagnosis. CALCA, a biomarker for medullary thyroid cancer, was hypersecreted in metastatic pancreatic cancer at least 16.5 months pre-diagnosis. ERBB2 levels spiked in metastatic breast cancer between 10.0 and 4.0 months pre-diagnosis. Our results support the value of deep phenotyping seemingly healthy individuals in prospectively inferring disease transitions.

Health and disease markers correlate with gut microbiome composition across thousands of people.

Variation in the human gut microbiome can reflect host lifestyle and behaviors and influence disease biomarker levels in the blood. Understanding the relationships between gut microbes and host phenotypes are critical for understanding wellness and disease. Here, we examine associations between the gut microbiota and ~150 host phenotypic features across ~3,400 individuals. We identify major axes of taxonomic variance in the gut and a putative diversity maximum along the Firmicutes-to-Bacteroidetes axis. Our analyses reveal both known and unknown associations between microbiome composition and host clinical markers and lifestyle factors, including host-microbe associations that are composition-specific. These results suggest potential opportunities for targeted interventions that alter the composition of the microbiome to improve host health. By uncovering the interrelationships between host diet and lifestyle factors, clinical blood markers, and the human gut microbiome at the population-scale, our results serve as a roadmap for future studies on host-microbe interactions and interventions.

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.