Understanding disease-associated metabolic changes in human colon epithelial cells using i ColonEpithelium metabolic reconstruction.

The colon epithelium plays a key role in the host-microbiome interactions, allowing uptake of various nutrients and driving important metabolic processes. To unravel detailed metabolic activities in the human colon epithelium, our present study focuses on the generation of the first cell-type specific genome-scale metabolic model (GEM) of human colonic epithelial cells, named iColonEpithelium. GEMs are powerful tools for exploring reactions and metabolites at systems level and predicting the flux distributions at steady state. Our cell-type-specific iColonEpithelium metabolic reconstruction captures genes specifically expressed in the human colonic epithelial cells. The iColonEpithelium is also capable of performing metabolic tasks specific to the cell type. A unique transport reaction compartment has been included to allow simulation of metabolic interactions with the gut microbiome. We used iColonEpithelium to identify metabolic signatures associated with inflammatory bowel disease. We integrated single-cell RNA sequencing data from Crohn's Diseases (CD) and ulcerative colitis (UC) samples with the iColonEpithelium metabolic network to predict metabolic signatures of colonocytes between CD and UC compared to healthy samples. We identified reactions in nucleotide interconversion, fatty acid synthesis and tryptophan metabolism were differentially regulated in CD and UC conditions, which were in accordance with experimental results. The iColonEpithelium metabolic network can be used to identify mechanisms at the cellular level, and our network has the potential to be integrated with gut microbiome models to explore the metabolic interactions between host and gut microbiota under various conditions.

Effects of Perfluorinated Alkyl Substances (PFAS) on Amphibian Body and Liver Conditions: Is Lipid Metabolism Being Perturbed throughout Metamorphosis?

Per- and polyfluoroalkyl substances (PFAS) may interact with peroxisome proliferator activated receptors (PPARs) and alter lipid homeostasis. Using Xenopus laevis, we investigated the effect of PFAS on (a) lipid homeostasis and whether this correlated to changes in body and hepatic condition; (b) the expression of hepatic genes regulated by PPAR; and (c) the hepatic lipidome. We chronically exposed tadpoles to 0.5 µg/L of either PFOS, PFHxS, PFOA, PFHxA, a binary mixture of PFOS and PFHxS (0.5 µg/L of each), or a control, from NF stage 52 through metamorphic climax. Growth, development, and survival were not affected, but we detected a sex-specific decrease in body condition at NF 66 (6.8%) and in hepatic condition (16.6%) across metamorphic climax for male tadpoles exposed to PFOS. We observed weak evidence for the transient downregulation of apolipoprotein-V (apoa5) at NF 62 in tadpoles exposed to PFHxA. Acyl-CoA oxidase 1 (acox1) was downregulated only in males exposed to PFHxS (Ln(Fold Change) = -0.54). We detected PFAS-specific downregulation of structural glycerophospholipids, while semi-quantitative profiling detected the upregulation in numerous glycerophospholipids, sphingomyelins, and diglycerides. Overall, our findings indicate that PFAS can induce sex-specific effects that change across larval development and metamorphosis. We demonstrate that PFAS alter lipid metabolism at environmentally relevant concentrations through divergent mechanisms that may not be related to PPARs, with an absence of effects on body condition, demonstrating the need for more molecular studies to elucidate mechanisms of PFAS-induced lipid dysregulation in amphibians and in other taxa.

A genome scale transcriptional regulatory model of the human placenta.

Gene regulation is essential to placental function and fetal development. We built a genome-scale transcriptional regulatory network (TRN) of the human placenta using digital genomic footprinting and transcriptomic data. We integrated 475 transcriptomes and 12 DNase hypersensitivity datasets from placental samples to globally and quantitatively map transcription factor (TF)-target gene interactions. In an independent dataset, the TRN model predicted target gene expression with an out-of-sample R2 greater than 0.25 for 73% of target genes. We performed siRNA knockdowns of four TFs and achieved concordance between the predicted gene targets in our TRN and differences in expression of knockdowns with an accuracy of >0.7 for three of the four TFs. Our final model contained 113,158 interactions across 391 TFs and 7712 target genes and is publicly available. We identified 29 TFs which were significantly enriched as regulators for genes previously associated with preterm birth, and eight of these TFs were decreased in preterm placentas.

Timing of selective serotonin reuptake inhibitor use and risk for preterm birth and related adverse events: with a consideration of the COVID-19 pandemic period.

OBJECTIVE: There is uncertainty around the safety of SSRIs for treating depression during pregnancy. Nevertheless, the use of SSRIs has been gradually increasing, especially during the COVID-19 pandemic period. We aimed to (1) characterize maternal depression rate and use of SSRIs in a recent 10-year period, (2) address confounding by indication, as well as socioeconomic and environmental factors, and (3) evaluate associations of the timing of SSRI exposure in pregnancy with risk for preterm birth (PTB), low birthweight (LBW), and small for gestational age (SGA) infants among women with depression before pregnancy. METHODS: We conducted propensity score-adjusted regression to calculate odds ratios (ORs) of PTB, LBW, and SGA. We accounted for maternal/pregnancy characteristics, comorbidity, depression severity, time of delivery, social vulnerability, and rural residence. RESULTS: There were 50.3% and 40.3% increases in the prevalence rate of prenatal depression and prenatal SSRI prescription rate during the pandemic. We identified women with depression ≤180 days before pregnancy (n = 8406). Women with no SSRI order during pregnancy (n = 3760) constituted the unexposed group. The late SSRI exposure group consisted of women with an SSRI order after the first trimester (n = 3759). The early-only SSRI exposure group consisted of women with SSRI orders only in the first trimester (n = 887). The late SSRI exposure group had an increased risk of PTB of OR = 1.5 ([1.2,1.8]) and LBW of OR = 1.5 ([1.2,2.0]), relative to the unexposed group. Associations between late SSRI exposure and risk of PTB/LBW were similar among a subsample of patients who delivered during the pandemic. CONCLUSIONS: These findings suggest an association between PTB/LBW and SSRI exposure is dependent on exposure timing during pregnancy. Small for gestational age is not associated with SSRI exposure.

Lifespan-extending interventions induce consistent patterns of fatty acid oxidation in mouse livers.

Aging manifests as progressive deteriorations in homeostasis, requiring systems-level perspectives to investigate the gradual molecular dysregulation of underlying biological processes. Here, we report systemic changes in the molecular regulation of biological processes under multiple lifespan-extending interventions. Differential Rank Conservation (DIRAC) analyses of mouse liver proteomics and transcriptomics data show that mechanistically distinct lifespan-extending interventions (acarbose, 17α-estradiol, rapamycin, and calorie restriction) generally tighten the regulation of biological modules. These tightening patterns are similar across the interventions, particularly in processes such as fatty acid oxidation, immune response, and stress response. Differences in DIRAC patterns between proteins and transcripts highlight specific modules which may be tightened via augmented cap-independent translation. Moreover, the systemic shifts in fatty acid metabolism are supported through integrated analysis of liver transcriptomics data with a mouse genome-scale metabolic model. Our findings highlight the power of systems-level approaches for identifying and characterizing the biological processes involved in aging and longevity.

A metabolomic signature of the APOE2 allele.

With the goal of identifying metabolites that significantly correlate with the protective e2 allele of the apolipoprotein E (APOE) gene, we established a consortium of five studies of healthy aging and extreme human longevity with 3545 participants. This consortium includes the New England Centenarian Study, the Baltimore Longitudinal Study of Aging, the Arivale study, the Longevity Genes Project/LonGenity studies, and the Long Life Family Study. We analyzed the association between APOE genotype groups E2 (e2e2 and e2e3 genotypes, N = 544), E3 (e3e3 genotypes, N = 2299), and E4 (e3e4 and e4e4 genotypes, N = 702) with metabolite profiles in the five studies and used fixed effect meta-analysis to aggregate the results. Our meta-analysis identified a signature of 19 metabolites that are significantly associated with the E2 genotype group at FDR < 10%. The group includes 10 glycerolipids and 4 glycerophospholipids that were all higher in E2 carriers compared to E3, with fold change ranging from 1.08 to 1.25. The organic acid 6-hydroxyindole sulfate, previously linked to changes in gut microbiome that were reflective of healthy aging and longevity, was also higher in E2 carriers compared to E3 carriers. Three sterol lipids and one sphingolipid species were significantly lower in carriers of the E2 genotype group. For some of these metabolites, the effect of the E2 genotype opposed the age effect. No metabolites reached a statistically significant association with the E4 group. This work confirms and expands previous results connecting the APOE gene to lipid regulation and suggests new links between the e2 allele, lipid metabolism, aging, and the gut-brain axis.

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.

Multi-Omic analyses characterize the ceramide/sphingomyelin pathway as a therapeutic target in Alzheimer's disease.

Dysregulation of sphingomyelin and ceramide metabolism have been implicated in Alzheimer's disease. Genome-wide and transcriptome-wide association studies have identified various genes and genetic variants in lipid metabolism that are associated with Alzheimer's disease. However, the molecular mechanisms of sphingomyelin and ceramide disruption remain to be determined. We focus on the sphingolipid pathway and carry out multi-omics analyses to identify central and peripheral metabolic changes in Alzheimer's patients, correlating them to imaging features. Our multi-omics approach is based on (a) 2114 human post-mortem brain transcriptomics to identify differentially expressed genes; (b) in silico metabolic flux analysis on context-specific metabolic networks identified differential reaction fluxes; (c) multimodal neuroimaging analysis on 1576 participants to associate genetic variants in sphingomyelin pathway with Alzheimer's disease pathogenesis; (d) plasma metabolomic and lipidomic analysis to identify associations of lipid species with dysregulation in Alzheimer's; and (e) metabolite genome-wide association studies to define receptors within the pathway as a potential drug target. We validate our hypothesis in amyloidogenic APP/PS1 mice and show prolonged exposure to fingolimod alleviated synaptic plasticity and cognitive impairment in mice. Our integrative multi-omics approach identifies potential targets in the sphingomyelin pathway and suggests modulators of S1P metabolism as possible candidates for Alzheimer's disease treatment.

Constraint-Based Reconstruction and Analyses of Metabolic Models: Open-Source Python Tools and Applications to Cancer.

The influence of metabolism on signaling, epigenetic markers, and transcription is highly complex yet important for understanding cancer physiology. Despite the development of high-resolution multi-omics technologies, it is difficult to infer metabolic activity from these indirect measurements. Fortunately, genome-scale metabolic models and constraint-based modeling provide a systems biology framework to investigate the metabolic states and define the genotype-phenotype associations by integrations of multi-omics data. Constraint-Based Reconstruction and Analysis (COBRA) methods are used to build and simulate metabolic networks using mathematical representations of biochemical reactions, gene-protein reaction associations, and physiological and biochemical constraints. These methods have led to advancements in metabolic reconstruction, network analysis, perturbation studies as well as prediction of metabolic state. Most computational tools for performing these analyses are written for MATLAB, a proprietary software. In order to increase accessibility and handle more complex datasets and models, community efforts have started to develop similar open-source tools in Python. To date there is a comprehensive set of tools in Python to perform various flux analyses and visualizations; however, there are still missing algorithms in some key areas. This review summarizes the availability of Python software for several components of COBRA methods and their applications in cancer metabolism. These tools are evolving rapidly and should offer a readily accessible, versatile way to model the intricacies of cancer metabolism for identifying cancer-specific metabolic features that constitute potential drug targets.

Multiple early factors anticipate post-acute COVID-19 sequelae.

Post-acute sequelae of COVID-19 (PASC) represent an emerging global crisis. However, quantifiable risk factors for PASC and their biological associations are poorly resolved. We executed a deep multi-omic, longitudinal investigation of 309 COVID-19 patients from initial diagnosis to convalescence (2-3 months later), integrated with clinical data and patient-reported symptoms. We resolved four PASC-anticipating risk factors at the time of initial COVID-19 diagnosis: type 2 diabetes, SARS-CoV-2 RNAemia, Epstein-Barr virus viremia, and specific auto-antibodies. In patients with gastrointestinal PASC, SARS-CoV-2-specific and CMV-specific CD8+ T cells exhibited unique dynamics during recovery from COVID-19. Analysis of symptom-associated immunological signatures revealed coordinated immunity polarization into four endotypes, exhibiting divergent acute severity and PASC. We find that immunological associations between PASC factors diminish over time, leading to distinct convalescent immune states. Detectability of most PASC factors at COVID-19 diagnosis emphasizes the importance of early disease measurements for understanding emergent chronic conditions and suggests PASC treatment strategies.

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.

Systems modeling of metabolic dysregulation in neurodegenerative diseases.

Neurodegenerative diseases (NDDs) encompass a wide range of conditions that arise owing to progressive degeneration and the ultimate loss of nerve cells in the brain and peripheral nervous system. NDDs such as Alzheimer's, Parkinson's, and Huntington's diseases negatively impact both length and quality of life, due to lack of effective disease-modifying treatments. Herein, we review the use of genome-scale metabolic models, network-based approaches, and integration with multiomics data to identify key biological processes that characterize NDDs. We describe powerful systems biology approaches for modeling NDD pathophysiology by leveraging in silico models that are informed by patient-derived multiomics data. These approaches can enable mechanistic insights into NDD-specific metabolic dysregulations that can be leveraged to identify potential metabolic markers of disease and predisease states.

Microglia show differential transcriptomic response to Aß peptide aggregates ex vivo and in vivo.

Aggregation and accumulation of amyloid-ß (Aß) is a defining feature of Alzheimer's disease pathology. To study microglial responses to Aß, we applied exogenous Aß peptide, in either oligomeric or fibrillar conformation, to primary mouse microglial cultures and evaluated system-level transcriptional changes and then compared these with transcriptomic changes in the brains of CRND8 APP mice. We find that primary microglial cultures have rapid and massive transcriptional change in response to Aß. Transcriptomic responses to oligomeric or fibrillar Aß in primary microglia, although partially overlapping, are distinct and are not recapitulated in vivo where Aß progressively accumulates. Furthermore, although classic immune mediators show massive transcriptional changes in the primary microglial cultures, these changes are not observed in the mouse model. Together, these data extend previous studies which demonstrate that microglia responses ex vivo are poor proxies for in vivo responses. Finally, these data demonstrate the potential utility of using microglia as biosensors of different aggregate conformation, as the transcriptional responses to oligomeric and fibrillar Aß can be distinguished.

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.

Bile acid synthesis, modulation, and dementia: A metabolomic, transcriptomic, and pharmacoepidemiologic study.

BACKGROUND: While Alzheimer disease (AD) and vascular dementia (VaD) may be accelerated by hypercholesterolemia, the mechanisms underlying this association are unclear. We tested whether dysregulation of cholesterol catabolism, through its conversion to primary bile acids (BAs), was associated with dementia pathogenesis. METHODS AND FINDINGS: We used a 3-step study design to examine the role of the primary BAs, cholic acid (CA), and chenodeoxycholic acid (CDCA) as well as their principal biosynthetic precursor, 7α-hydroxycholesterol (7α-OHC), in dementia. In Step 1, we tested whether serum markers of cholesterol catabolism were associated with brain amyloid accumulation, white matter lesions (WMLs), and brain atrophy. In Step 2, we tested whether exposure to bile acid sequestrants (BAS) was associated with risk of dementia. In Step 3, we examined plausible mechanisms underlying these findings by testing whether brain levels of primary BAs and gene expression of their principal receptors are altered in AD. Step 1: We assayed serum concentrations CA, CDCA, and 7α-OHC and used linear regression and mixed effects models to test their associations with brain amyloid accumulation (N = 141), WMLs, and brain atrophy (N = 134) in the Baltimore Longitudinal Study of Aging (BLSA). The BLSA is an ongoing, community-based cohort study that began in 1958. Participants in the BLSA neuroimaging sample were approximately 46% male with a mean age of 76 years; longitudinal analyses included an average of 2.5 follow-up magnetic resonance imaging (MRI) visits. We used the Alzheimer's Disease Neuroimaging Initiative (ADNI) (N = 1,666) to validate longitudinal neuroimaging results in BLSA. ADNI is an ongoing, community-based cohort study that began in 2003. Participants were approximately 55% male with a mean age of 74 years; longitudinal analyses included an average of 5.2 follow-up MRI visits. Lower serum concentrations of 7α-OHC, CA, and CDCA were associated with higher brain amyloid deposition (p = 0.041), faster WML accumulation (p = 0.050), and faster brain atrophy mainly (false discovery rate [FDR] p = <0.001-0.013) in males in BLSA. In ADNI, we found a modest sex-specific effect indicating that lower serum concentrations of CA and CDCA were associated with faster brain atrophy (FDR p = 0.049) in males.Step 2: In the Clinical Practice Research Datalink (CPRD) dataset, covering >4 million registrants from general practice clinics in the United Kingdom, we tested whether patients using BAS (BAS users; 3,208 with ≥2 prescriptions), which reduce circulating BAs and increase cholesterol catabolism, had altered dementia risk compared to those on non-statin lipid-modifying therapies (LMT users; 23,483 with ≥2 prescriptions). Patients in the study (BAS/LMT) were approximately 34%/38% male and with a mean age of 65/68 years; follow-up time was 4.7/5.7 years. We found that BAS use was not significantly associated with risk of all-cause dementia (hazard ratio (HR) = 1.03, 95% confidence interval (CI) = 0.72-1.46, p = 0.88) or its subtypes. We found a significant difference between the risk of VaD in males compared to females (p = 0.040) and a significant dose-response relationship between BAS use and risk of VaD (p-trend = 0.045) in males.Step 3: We assayed brain tissue concentrations of CA and CDCA comparing AD and control (CON) samples in the BLSA autopsy cohort (N = 29). Participants in the BLSA autopsy cohort (AD/CON) were approximately 50%/77% male with a mean age of 87/82 years. We analyzed single-cell RNA sequencing (scRNA-Seq) data to compare brain BA receptor gene expression between AD and CON samples from the Religious Orders Study and Memory and Aging Project (ROSMAP) cohort (N = 46). ROSMAP is an ongoing, community-based cohort study that began in 1994. Participants (AD/CON) were approximately 56%/36% male with a mean age of 85/85 years. In BLSA, we found that CA and CDCA were detectable in postmortem brain tissue samples and were marginally higher in AD samples compared to CON. In ROSMAP, we found sex-specific differences in altered neuronal gene expression of BA receptors in AD. Study limitations include the small sample sizes in the BLSA cohort and likely inaccuracies in the clinical diagnosis of dementia subtypes in primary care settings. CONCLUSIONS: We combined targeted metabolomics in serum and amyloid positron emission tomography (PET) and MRI of the brain with pharmacoepidemiologic analysis to implicate dysregulation of cholesterol catabolism in dementia pathogenesis. We observed that lower serum BA concentration mainly in males is associated with neuroimaging markers of dementia, and pharmacological lowering of BA levels may be associated with higher risk of VaD in males. We hypothesize that dysregulation of BA signaling pathways in the brain may represent a plausible biologic mechanism underlying these results. Together, our observations suggest a novel mechanism relating abnormalities in cholesterol catabolism to risk of dementia.

Immune Subtyping in Latent Tuberculosis.

Latent tuberculosis infection (LTBI) poses a major roadblock in the global effort to eradicate tuberculosis (TB). A deep understanding of the host responses involved in establishment and maintenance of TB latency is required to propel the development of sensitive methods to detect and treat LTBI. Given that LTBI individuals are typically asymptomatic, it is challenging to differentiate latently infected from uninfected individuals. A major contributor to this problem is that no clear pattern of host response is linked with LTBI, as molecular correlates of latent infection have been hard to identify. In this study, we have analyzed the global perturbations in host response in LTBI individuals as compared to uninfected individuals and particularly the heterogeneity in such response, across LTBI cohorts. For this, we constructed individualized genome-wide host response networks informed by blood transcriptomes for 136 LTBI cases and have used a sensitive network mining algorithm to identify top-ranked host response subnetworks in each case. Our analysis indicates that despite the high heterogeneity in the gene expression profiles among LTBI samples, clear patterns of perturbation are found in the immune response pathways, leading to grouping LTBI samples into 4 different immune-subtypes. Our results suggest that different subnetworks of molecular perturbations are associated with latent tuberculosis.

Modulation of Immune Checkpoints by Chemotherapy in Human Colorectal Liver Metastases.

Metastatic colorectal cancer (CRC) is a major cause of cancer-related death, and incidence is rising in younger populations (younger than 50 years). Current chemotherapies can achieve response rates above 50%, but immunotherapies have limited value for patients with microsatellite-stable (MSS) cancers. The present study investigates the impact of chemotherapy on the tumor immune microenvironment. We treat human liver metastases slices with 5-fluorouracil (5-FU) plus either irinotecan or oxaliplatin, then perform single-cell transcriptome analyses. Results from eight cases reveal two cellular subtypes with divergent responses to chemotherapy. Susceptible tumors are characterized by a stemness signature, an activated interferon pathway, and suppression of PD-1 ligands in response to 5-FU+irinotecan. Conversely, immune checkpoint TIM-3 ligands are maintained or upregulated by chemotherapy in CRC with an enterocyte-like signature, and combining chemotherapy with TIM-3 blockade leads to synergistic tumor killing. Our analyses highlight chemomodulation of the immune microenvironment and provide a framework for combined chemo-immunotherapies.

Metabolic Network Analysis Reveals Altered Bile Acid Synthesis and Metabolism in Alzheimer's Disease.

Increasing evidence suggests Alzheimer's disease (AD) pathophysiology is influenced by primary and secondary bile acids, the end product of cholesterol metabolism. We analyze 2,114 post-mortem brain transcriptomes and identify genes in the alternative bile acid synthesis pathway to be expressed in the brain. A targeted metabolomic analysis of primary and secondary bile acids measured from post-mortem brain samples of 111 individuals supports these results. Our metabolic network analysis suggests that taurine transport, bile acid synthesis, and cholesterol metabolism differ in AD and cognitively normal individuals. We also identify putative transcription factors regulating metabolic genes and influencing altered metabolism in AD. Intriguingly, some bile acids measured in brain tissue cannot be explained by the presence of enzymes responsible for their synthesis, suggesting that they may originate from the gut microbiome and are transported to the brain. These findings motivate further research into bile acid metabolism in AD to elucidate their possible connection to cognitive decline.