Systematic transcriptome-wide meta-analysis across endocrine disrupting chemicals reveals shared and unique liver pathways, gene networks, and disease associations.

Cardiometabolic disorders (CMD) are a growing public health problem across the world. Among the known cardiometabolic risk factors are compounds that induce endocrine and metabolic dysfunctions, such as endocrine disrupting chemicals (EDCs). To date, how EDCs influence molecular programs and cardiometabolic risks has yet to be fully elucidated, especially considering the complexity contributed by species-, chemical-, and dose-specific effects. Moreover, different experimental and analytical methodologies employed by different studies pose challenges when comparing findings across studies. To explore the molecular mechanisms of EDCs in a systematic manner, we established a data-driven computational approach to meta-analyze 30 human, mouse, and rat liver transcriptomic datasets for 4 EDCs, namely bisphenol A (BPA), bis(2-ethylhexyl) phthalate (DEHP), tributyltin (TBT), and perfluorooctanoic acid (PFOA). Our computational pipeline uniformly re-analyzed pre-processed quality-controlled microarray data and raw RNAseq data, derived differentially expressed genes (DEGs) and biological pathways, modeled gene regulatory networks and regulators, and determined CMD associations based on gene overlap analysis. Our approach revealed that DEHP and PFOA shared stable transcriptomic signatures that are enriched for genes associated with CMDs, suggesting similar mechanisms of action such as perturbations of peroxisome proliferator-activated receptor gamma (PPARγ) signaling and liver gene network regulators VNN1 and ACOT2. In contrast, TBT exhibited highly divergent gene signatures, pathways, network regulators, and disease associations from the other EDCs. In addition, we found that the rat, mouse, and human BPA studies showed highly variable transcriptomic patterns, providing molecular support for the variability in BPA responses. Our work offers insights into the commonality and differences in the molecular mechanisms of various EDCs and establishes a streamlined data-driven workflow to compare molecular mechanisms of environmental substances to elucidate the underlying connections between chemical exposure and disease risks.

Role of Matrix Gla Protein in Transforming Growth Factor-ß Signaling and Nonalcoholic Steatohepatitis in Mice.

BACKGROUND & AIMS: Nonalcoholic steatohepatitis (NASH) is a complex disease involving both genetic and environmental factors in its onset and progression. We analyzed NASH phenotypes in a genetically diverse cohort of mice, the Hybrid Mouse Diversity Panel, to identify genes contributing to disease susceptibility. METHODS: A "systems genetics" approach, involving integration of genetic, transcriptomic, and phenotypic data, was used to identify candidate genes and pathways in a mouse model of NASH. The causal role of Matrix Gla Protein (MGP) was validated using heterozygous MGP knockout (Mgp+/-) mice. The mechanistic role of MGP in transforming growth factor-beta (TGF-ß) signaling was examined in the LX-2 stellate cell line by using a loss of function approach. RESULTS: Local cis-acting regulation of MGP was correlated with fibrosis, suggesting a causal role in NASH, and this was validated using loss of function experiments in 2 models of diet-induced NASH. Using single-cell RNA sequencing, Mgp was found to be primarily expressed in hepatic stellate cells and dendritic cells in mice. Knockdown of MGP expression in stellate LX-2 cells led to a blunted response to TGF-ß stimulation. This was associated with reduced regulatory SMAD phosphorylation and TGF-ß receptor ALK1 expression as well as increased expression of inhibitory SMAD6. Hepatic MGP expression was found to be significantly correlated with the severity of fibrosis in livers of patients with NASH, suggesting relevance to human disease. CONCLUSIONS: MGP regulates liver fibrosis and TGF-ß signaling in hepatic stellate cells and contributes to NASH pathogenesis.

Mammary duct luminal epithelium controls adipocyte thermogenic programme.

Sympathetic activation during cold exposure increases adipocyte thermogenesis via the expression of mitochondrial protein uncoupling protein 1 (UCP1)1. The propensity of adipocytes to express UCP1 is under a critical influence of the adipose microenvironment and varies between sexes and among various fat depots2-7. Here we report that mammary gland ductal epithelial cells in the adipose niche regulate cold-induced adipocyte UCP1 expression in female mouse subcutaneous white adipose tissue (scWAT). Single-cell RNA sequencing shows that glandular luminal epithelium subtypes express transcripts that encode secretory factors controlling adipocyte UCP1 expression under cold conditions. We term these luminal epithelium secretory factors 'mammokines'. Using 3D visualization of whole-tissue immunofluorescence, we reveal sympathetic nerve-ductal contact points. We show that mammary ducts activated by sympathetic nerves limit adipocyte UCP1 expression via the mammokine lipocalin 2. In vivo and ex vivo ablation of mammary duct epithelium enhance the cold-induced adipocyte thermogenic gene programme in scWAT. Since the mammary duct network extends throughout most of the scWAT in female mice, females show markedly less scWAT UCP1 expression, fat oxidation, energy expenditure and subcutaneous fat mass loss compared with male mice, implicating sex-specific roles of mammokines in adipose thermogenesis. These results reveal a role of sympathetic nerve-activated glandular epithelium in adipocyte UCP1 expression and suggest that mammary duct luminal epithelium has an important role in controlling glandular adiposity.

Conserved multi-tissue transcriptomic adaptations to exercise training in humans and mice.

Physical activity is associated with beneficial adaptations in human and rodent metabolism. We studied over 50 complex traits before and after exercise intervention in middle-aged men and a panel of 100 diverse strains of female mice. Candidate gene analyses in three brain regions, muscle, liver, heart, and adipose tissue of mice indicate genetic drivers of clinically relevant traits, including volitional exercise volume, muscle metabolism, adiposity, and hepatic lipids. Although ∼33% of genes differentially expressed in skeletal muscle following the exercise intervention are similar in mice and humans independent of BMI, responsiveness of adipose tissue to exercise-stimulated weight loss appears controlled by species and underlying genotype. We leveraged genetic diversity to generate prediction models of metabolic trait responsiveness to volitional activity offering a framework for advancing personalized exercise prescription. The human and mouse data are publicly available via a user-friendly Web-based application to enhance data mining and hypothesis development.

Prdm16 and Vcam1 regulate the postnatal disappearance of embryonic radial glia and the ending of cortical neurogenesis.

Embryonic neural stem cells (NSCs, i.e., radial glia) in the ventricular-subventricular zone (V-SVZ) generate the majority of neurons and glia in the forebrain. Postnatally, embryonic radial glia disappear and a subpopulation of radial glia transition into adult NSCs. As this transition occurs, widespread neurogenesis in brain regions such as the cerebral cortex ends. The mechanisms that regulate the postnatal disappearance of radial glia and the ending of embryonic neurogenesis remain poorly understood. Here, we show that PR domain-containing 16 (Prdm16) promotes the disappearance of radial glia and the ending of neurogenesis in the cerebral cortex. Genetic deletion of Prdm16 from NSCs leads to the persistence of radial glia in the adult V-SVZ and prolonged postnatal cortical neurogenesis. Mechanistically, Prdm16 induces the postnatal reduction in Vascular Cell Adhesion Molecule 1 (Vcam1). The postnatal disappearance of radial glia and the ending of cortical neurogenesis occur normally in Prdm16-Vcam1 double conditional knockout mice. These observations reveal novel molecular regulators of the postnatal disappearance of radial glia and the ending of embryonic neurogenesis, filling a key knowledge gap in NSC biology.

Systems spatiotemporal dynamics of traumatic brain injury at single-cell resolution reveals humanin as a therapeutic target.

BACKGROUND: The etiology of mild traumatic brain injury (mTBI) remains elusive due to the tissue and cellular heterogeneity of the affected brain regions that underlie cognitive impairments and subsequent neurological disorders. This complexity is further exacerbated by disrupted circuits within and between cell populations across brain regions and the periphery, which occur at different timescales and in spatial domains. METHODS: We profiled three tissues (hippocampus, frontal cortex, and blood leukocytes) at the acute (24-h) and subacute (7-day) phases of mTBI at single-cell resolution. RESULTS: We demonstrated that the coordinated gene expression patterns across cell types were disrupted and re-organized by TBI at different timescales with distinct regional and cellular patterns. Gene expression-based network modeling implied astrocytes as a key regulator of the cell-cell coordination following mTBI in both hippocampus and frontal cortex across timepoints, and mt-Rnr2, which encodes the mitochondrial peptide humanin, as a potential target for intervention based on its broad regional and dynamic dysregulation following mTBI. Treatment of a murine mTBI model with humanin reversed cognitive impairment caused by mTBI through the restoration of metabolic pathways within astrocytes. CONCLUSIONS: Our results offer a systems-level understanding of the dynamic and spatial regulation of gene programs by mTBI and pinpoint key target genes, pathways, and cell circuits that are amenable to therapeutics.

PharmOmics: A species- and tissue-specific drug signature database and gene-network-based drug repositioning tool.

Drug development has been hampered by a high failure rate in clinical trials due to our incomplete understanding of drug functions across organs and species. Therefore, elucidating species- and tissue-specific drug functions can provide insights into therapeutic efficacy, potential adverse effects, and interspecies differences necessary for effective translational medicine. Here, we present PharmOmics, a drug knowledgebase and analytical tool that is hosted on an interactive web server. Using tissue- and species-specific transcriptome data from human, mouse, and rat curated from different databases, we implemented a gene-network-based approach for drug repositioning. We demonstrate the potential of PharmOmics to retrieve known therapeutic drugs and identify drugs with tissue toxicity using in silico performance assessment. We further validated predicted drugs for nonalcoholic fatty liver disease in mice. By combining tissue- and species-specific in vivo drug signatures with gene networks, PharmOmics serves as a complementary tool to support drug characterization and network-based medicine.

Determination of the nanoparticle- and cell-specific toxicological mechanisms in 3D liver spheroids using scRNAseq analysis.

Engineered nanomaterials (ENMs) are commonly used in consumer products, allowing exposure to target organs such as the lung, liver, and skin that could lead to adverse health effects in humans. To better reflect on toxicological effects in liver cells, it is important to consider the contribution of hepatocyte morphology, function, and intercellular interactions in a dynamic 3D microenvironment. Herein, we used a 3D liver spheroid model containing hepatocyte and Kupffer cells (KCs) to study the effects of three different material compositions, namely vanadium pentoxide (V2O5), titanium dioxide (TiO2), or graphene oxide (GO). Additionally, we used single-cell RNA sequencing (scRNAseq) to determine the nanoparticle (NP) and cell-specific toxicological responses. A general finding was that hepatocytes exhibit more variation in gene expression and adaptation of signaling pathways than KCs. TNF-α production tied to the NF-κB pathway was a commonly affected pathway by all NPs while impacts on the metabolic function of hepatocytes were unique to V2O5. V2O5 NPs also showed the largest number of differentially expressed genes in both cell types, many of which are related to pro-inflammatory and apoptotic response pathways. There was also evidence of mitochondrial ROS generation and caspase-1 activation after GO and V2O5 treatment, in association with cytokine production. All considered, this study provides insight into the impact of nanoparticles on gene responses in key liver cell types, providing us with a scRNAseq platform that can be used for high-content screening of nanomaterial impact on the liver, for use in biosafety and biomedical applications.

Disparate Metabolomic Responses to Fructose Consumption between Different Mouse Strains and the Role of Gut Microbiota.

High fructose consumption has been linked to metabolic syndrome, yet the fructose-induced phenotypes, gene expression, and gut microbiota alterations are distinct between mouse strains. In this study, we aim to investigate how fructose consumption shapes the metabolomic profiles of mice with different genetic background and microbiome. We used fructose-sensitive DBA/2J (DBA) and fructose-resistant C57BL/6J (B6) mice given 8% fructose or regular water for 12 weeks. Plasma and fecal metabolites were profiled using a liquid chromatography-tandem mass spectrometry based global metabolomic approach. We found that the baseline metabolomic profiles were different between DBA and B6 mice, particularly plasma metabolites involved in lipid metabolism and fecal metabolites related to dipeptide/amino acid metabolism. In response to fructose, DBA mice showed a distinct decrease of plasma branched chain fatty acids with concordantly increased branched chain amino acids, which were correlated with adiposity; B6 mice had significantly increased plasma cholesterol and total bile acids, accompanied by decreased fecal levels of farnesoid X receptor antagonist tauro-ß-muricholate, which were correlated with fructose-responsive bacteria Dehalobacterium, Magibacteriaceae, and/or Akkermansia. Our results demonstrate that baseline metabolomic profiles differ and respond differentially to fructose between mice with different genetic background and gut microbiota, which may play a role in individualized risks to fructose-induced metabolic syndrome.

Estrogen receptor alpha in the brain mediates tamoxifen-induced changes in physiology in mice.

Adjuvant tamoxifen therapy improves survival in breast cancer patients. Unfortunately, long-term treatment comes with side effects that impact health and quality of life, including hot flashes, changes in bone density, and fatigue. Partly due to a lack of proven animal models, the tissues and cells that mediate these negative side effects are unclear. Here, we show that mice undergoing tamoxifen treatment experience changes in temperature, bone, and movement. Single-cell RNA sequencing reveals that tamoxifen treatment induces widespread gene expression changes in the hypothalamus and preoptic area (hypothalamus-POA). These expression changes are dependent on estrogen receptor alpha (ERα), as conditional knockout of ERα in the hypothalamus-POA ablates or reverses tamoxifen-induced gene expression. Accordingly, ERα-deficient mice do not exhibit tamoxifen-induced changes in temperature, bone, or movement. These findings provide mechanistic insight into the effects of tamoxifen on the hypothalamus-POA and indicate that ERα mediates several physiological effects of tamoxifen treatment in mice.

Estrogen is a hormone often known for its role in female development and reproduction. Yet, it also has an impact on many biological processes such as immunity and the health of bones, the heart, or the brain. It usually works by attaching to receptor proteins in specific cells. For instance, estrogen-responsive cells are present in the hypothalamus, the brain area that controls energy levels as well as the body's temperature and internal clock. Breast cancer cells are also often sensitive to estrogen, with the hormone fuelling the growth of tumors. The drug tamoxifen blocks estrogen receptors, stopping cells from responding to the hormone. As such, it is often used to reduce the likelihood that estrogen-dependent breast cancer will come back after treatment. However, its use can induce hot flashes, changes in bone density, fatigue and other life-altering side effects. Here, Zhang et al. investigated how estrogen receptors in the hypothalamus and a related region known as the preoptic area could be responsible for these side effects in mice. When the rodents were given tamoxifen for 28 days, they experienced changes in temperature, bone density and movement similar to those found in humans. In fact, genetic analyses revealed that the drug altered the way genes were turned on and off in certain cells types in the hypothalamus. Crucially, mice whose hypothalamus and preoptic area lacked estrogen receptors did not experience these behavioral and biological alterations. The findings by Zhang et al. help to understand how the side effects of tamoxifen emerge, singling out estrogen receptors in particular brain regions. This result could help to develop new therapies so that breast cancer can be treated with a better quality of life.

Single-Cell Study of Two Rat Models of Pulmonary Arterial Hypertension Reveals Connections to Human Pathobiology and Drug Repositioning.

Rationale: The cellular and molecular landscape and translational value of commonly used models of pulmonary arterial hypertension (PAH) are poorly understood. Single-cell transcriptomics can enhance molecular understanding of preclinical models and facilitate their rational use and interpretation.Objectives: To determine and prioritize dysregulated genes, pathways, and cell types in lungs of PAH rat models to assess relevance to human PAH and identify drug repositioning candidates.Methods: Single-cell RNA sequencing was performed on the lungs of monocrotaline (MCT), Sugen-hypoxia (SuHx), and control rats to identify altered genes and cell types, followed by validation using flow-sorted cells, RNA in situ hybridization, and immunofluorescence. Relevance to human PAH was assessed by histology of lungs from patients and via integration with human PAH genetic loci and known disease genes. Candidate drugs were predicted using Connectivity Map.Measurements and Main Results: Distinct changes in genes and pathways in numerous cell types were identified in SuHx and MCT lungs. Widespread upregulation of NF-κB signaling and downregulation of IFN signaling was observed across cell types. SuHx nonclassical monocytes and MCT conventional dendritic cells showed particularly strong NF-κB pathway activation. Genes altered in SuHx nonclassical monocytes were significantly enriched for PAH-associated genes and genetic variants, and candidate drugs predicted to reverse the changes were identified. An open-access online platform was developed to share single-cell data and drug candidates (http://mergeomics.research.idre.ucla.edu/PVDSingleCell/).Conclusions: Our study revealed the distinct and shared dysregulation of genes and pathways in two commonly used PAH models for the first time at single-cell resolution and demonstrated their relevance to human PAH and utility for drug repositioning.

Systems toxicogenomics of prenatal low-dose BPA exposure on liver metabolic pathways, gut microbiota, and metabolic health in mice.

Bisphenol A (BPA) is an industrial plasticizer widely found in consumer products, and exposure to BPA during early development has been associated with the prevalence of various cardiometabolic diseases including obesity, metabolic syndrome, type 2 diabetes, and cardiovascular diseases. To elucidate the molecular perturbations underlying the connection of low-dose prenatal BPA exposure to cardiometabolic diseases, we conducted a multi-dimensional systems biology study assessing the liver transcriptome, gut microbial community, and diverse metabolic phenotypes in both male and female mouse offspring exposed to 5 µg/kg/day BPA during gestation. Prenatal exposure to low-dose BPA not only significantly affected liver genes involved in oxidative phosphorylation, PPAR signaling and fatty acid metabolism, but also affected the gut microbial composition in an age- and sex-dependent manner. Bacteria such as those belonging to the S24-7 and Lachnospiraceae families were correlated with offspring phenotypes, differentially expressed liver metabolic genes such as Acadl and Dgat1, and key drivers identified in our gene network modeling such as Malat1 and Apoa2. This multiomics study provides insight into the relationship between gut bacteria and host liver genes that could contribute to cardiometabolic disease risks upon low-dose BPA exposure.

Host Genetic Background and Gut Microbiota Contribute to Differential Metabolic Responses to Fructose Consumption in Mice.

BACKGROUND: It is unclear how high fructose consumption induces disparate metabolic responses in genetically diverse mouse strains. OBJECTIVE: We aimed to investigate whether the gut microbiota contributes to differential metabolic responses to fructose. METHODS: Eight-week-old male C57BL/6J (B6), DBA/2J (DBA), and FVB/NJ (FVB) mice were given 8% fructose solution or regular water (control) for 12 wk. The gut microbiota composition in cecum and feces was analyzed using 16S ribosomal DNA sequencing, and permutational multivariate ANOVA (PERMANOVA) was used to compare community across mouse strains, treatments, and time points. Microbiota abundance was correlated with metabolic phenotypes and host gene expression in hypothalamus, liver, and adipose tissues using Biweight midcorrelation. To test the causal role of the gut microbiota in determining fructose response, we conducted fecal transplants from B6 to DBA mice and vice versa for 4 wk, as well as gavaged antibiotic-treated DBA mice with Akkermansia for 9 wk, accompanied with or without fructose treatment. RESULTS: Compared with B6 and FVB, DBA mice had significantly higher Firmicutes to Bacteroidetes ratio and lower baseline abundance of Akkermansia and S24-7 (P < 0.05), accompanied by metabolic dysregulation after fructose consumption. Fructose altered specific microbial taxa in individual mouse strains, such as a 7.27-fold increase in Akkermansia in B6 and 0.374-fold change in Rikenellaceae in DBA (false discovery rate <5%), which demonstrated strain-specific correlations with host metabolic and transcriptomic phenotypes. Fecal transplant experiments indicated that B6 microbes conferred resistance to fructose-induced weight gain in DBA mice (F = 43.1, P < 0.001), and Akkermansia colonization abrogated the fructose-induced weight gain (F = 17.8, P < 0.001) and glycemic dysfunctions (F = 11.8, P = 0.004) in DBA mice. CONCLUSIONS: Our findings support that differential microbiota composition between mouse strains is partially responsible for host metabolic sensitivity to fructose, and that Akkermansia is a key bacterium that confers resistance to fructose-induced metabolic dysregulation.

Omega-3 fatty acids increase OXPHOS energy for immune therapy of Alzheimer disease patients.

Sporadic late-onset Alzheimer disease (LOAD) preceded by mild cognitive impairment (MCI) is the most common type of dementia. Long-term studies of immunity to pathogenic amyloid-ß (Aß) in LOAD are lacking. Innate immunity of LOAD patients is malfunctioning in phagocytosis and degradation of Aß and LOAD patients' macrophage transcriptome and metabolome are deregulated. We previously showed omega-3 fatty acid (ω-3)-mediated repair of unfolded protein response and here we show much broader transcriptomic effects. ω-3 treatment in vitro and ω-3 supplementation by the drink Smartfish (SMF) in vivo increased the transcripts of the genes and pathways of immunity, glycolysis, tricarboxylic acid cycle, OX-PHOS, nicotinamide dinucleotide (NAD+ ) synthesis, and reversed the defects in Aß phagocytosis. In both peripheral blood mononuclear cells (PBMC) and macrophages, ω-3 increased ATP-linked oxygen consumption rate (OCR) and ω-3 with carnitine was superior to ω-3. ω-3 treatment in vitro and supplementation by the ω-3 drink SMF in vivo rescued macrophage phagocytosis when glycolysis or glycosylation were blocked. ω-3 provide flexible energy for immune clearance of the brain throughout the diurnal cycle, even in hypo- or hyper-glycemia. In certain LOAD patients, ω-3 may delay progression to dementia.

Single-cell RNA-seq analysis of the brainstem of mutant SOD1 mice reveals perturbed cell types and pathways of amyotrophic lateral sclerosis.

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease in which motor neurons throughout the brain and spinal cord progressively degenerate resulting in muscle atrophy, paralysis and death. Recent studies using animal models of ALS implicate multiple cell-types (e.g., astrocytes and microglia) in ALS pathogenesis in the spinal motor systems. To ascertain cellular vulnerability and cell-type specific mechanisms of ALS in the brainstem that orchestrates oral-motor functions, we conducted parallel single cell RNA sequencing (scRNA-seq) analysis using the high-throughput Drop-seq method. We isolated 1894 and 3199 cells from the brainstem of wildtype and mutant SOD1 symptomatic mice respectively, at postnatal day 100. We recovered major known cell types and neuronal subpopulations, such as interneurons and motor neurons, and trigeminal ganglion (TG) peripheral sensory neurons, as well as, previously uncharacterized interneuron subtypes. We found that the majority of the cell types displayed transcriptomic alterations in ALS mice. Differentially expressed genes (DEGs) of individual cell populations revealed cell-type specific alterations in numerous pathways, including previously known ALS pathways such as inflammation (in microglia), stress response (ependymal and an uncharacterized cell population), neurogenesis (astrocytes, oligodendrocytes, neurons), synapse organization and transmission (microglia, oligodendrocyte precursor cells, and neuronal subtypes), and mitochondrial function (uncharacterized cell populations). Other cell-type specific processes altered in SOD1 mutant brainstem include those from motor neurons (axon regeneration, voltage-gated sodium and potassium channels underlying excitability, potassium ion transport), trigeminal sensory neurons (detection of temperature stimulus involved in sensory perception), and cellular response to toxic substances (uncharacterized cell populations). DEGs consistently altered across cell types (e.g., Malat1), as well as cell-type specific DEGs, were identified. Importantly, DEGs from various cell types overlapped with known ALS genes from the literature and with top hits from an existing human ALS genome-wide association study (GWAS), implicating the potential cell types in which the ALS genes function with ALS pathogenesis. Our molecular investigation at single cell resolution provides comprehensive insights into the cell types, genes and pathways altered in the brainstem in a widely used ALS mouse model.

Therapeutic IDOL Reduction Ameliorates Amyloidosis and Improves Cognitive Function in APP/PS1 Mice.

Brain lipoprotein receptors have been shown to regulate the metabolism of ApoE and ß-amyloid (Aß) and are potential therapeutic targets for Alzheimer's disease (AD). Previously, we identified E3 ubiquitin ligase IDOL as a negative regulator of brain lipoprotein receptors. Genetic ablation of Idol increases low-density lipoprotein receptor protein levels, which facilitates Aß uptake and clearance by microglia. In this study, we utilized an antisense oligonucleotide (ASO) to reduce IDOL expression therapeutically in the brains of APP/PS1 male mice. ASO treatment led to decreased Aß pathology and improved spatial learning and memory. Single-cell transcriptomic analysis of hippocampus revealed that IDOL inhibition upregulated lysosomal/phagocytic genes in microglia. Furthermore, clustering of microglia revealed that IDOL-ASO treatment shifted the composition of the microglia population by increasing the prevalence of disease-associated microglia. Our results suggest that reducing IDOL expression in the adult brain promotes the phagocytic clearance of Aß and ameliorates Aß-dependent pathology. Pharmacological inhibition of IDOL activity in the brain may represent a therapeutic strategy for the treatment of AD.

Diesel exhaust particles dysregulate multiple immunological pathways in murine macrophages: Lessons from microarray and scRNA-seq technologies.

Exposure to ambient particulate matter has been shown to promote a variety of disorders, including cardiovascular diseases predominantly of ischemic etiology. However, the mechanisms linking inhaled particulates with systemic vascular effects, resulting in worsened atherosclerosis, are not well defined. We assessed the potential role of macrophages in translating these effects by analyzing gene expression patterns in response to diesel exhaust particles (DEP) at the average cell level, using Affymetrix microarrays in peritoneal macrophages in culture (in vitro), and at the individual cell level, using single-cell RNA sequencing (scRNA-seq) in alveolar macrophages collected from exposed mice (in vivo). Peritoneal macrophages were harvested from C57BL/6J mice and treated with 25 µg/mL of a DEP methanol extract (DEPe). These cells exhibited significant (FDR < 0.05) differential expression of a large number of genes and enrichment in pathways, especially engaged in immune responses and antioxidant defense. DEPe led to marked upregulation of heme oxygenase 1 (Hmox1), the most significantly upregulated gene (FDR = 1.75E-06), and several other antioxidant genes. For the in vivo work, C57BL/6J mice were subjected to oropharyngeal aspiration of 200 µg of whole DEP. The gene expression profiles of the alveolar macrophages harvested from these mice were analyzed at the single-cell level using scRNA-seq, which showed significant dysregulation of a broad number of genes enriched in immune system pathways as well, but with a large heterogeneity in how individual alveolar macrophages responded to DEP exposures. Altogether, DEP pollutants dysregulated immunological pathways in macrophages that may mediate the development of pulmonary and systemic vascular effects.

Single cell analysis reveals immune cell-adipocyte crosstalk regulating the transcription of thermogenic adipocytes.

Immune cells are vital constituents of the adipose microenvironment that influence both local and systemic lipid metabolism. Mice lacking IL10 have enhanced thermogenesis, but the roles of specific cell types in the metabolic response to IL10 remain to be defined. We demonstrate here that selective loss of IL10 receptor α in adipocytes recapitulates the beneficial effects of global IL10 deletion, and that local crosstalk between IL10-producing immune cells and adipocytes is a determinant of thermogenesis and systemic energy balance. Single Nuclei Adipocyte RNA-sequencing (SNAP-seq) of subcutaneous adipose tissue defined a metabolically-active mature adipocyte subtype characterized by robust expression of genes involved in thermogenesis whose transcriptome was selectively responsive to IL10Rα deletion. Furthermore, single-cell transcriptomic analysis of adipose stromal populations identified lymphocytes as a key source of IL10 production in response to thermogenic stimuli. These findings implicate adaptive immune cell-adipocyte communication in the maintenance of adipose subtype identity and function.

Prenatal Bisphenol A Exposure in Mice Induces Multitissue Multiomics Disruptions Linking to Cardiometabolic Disorders.

The health impacts of endocrine-disrupting chemicals (EDCs) remain debated, and their tissue and molecular targets are poorly understood. In this study, we leveraged systems biology approaches to assess the target tissues, molecular pathways, and gene regulatory networks associated with prenatal exposure to the model EDC bisphenol A (BPA). Prenatal BPA exposure at 5 mg/kg/d, a dose below most reported no-observed-adverse-effect levels, led to tens to thousands of transcriptomic and methylomic alterations in the adipose, hypothalamus, and liver tissues in male offspring in mice, with cross-tissue perturbations in lipid metabolism as well as tissue-specific alterations in histone subunits, glucose metabolism, and extracellular matrix. Network modeling prioritized main molecular targets of BPA, including Pparg, Hnf4a, Esr1, Srebf1, and Fasn as well as numerous less studied targets such as Cyp51 and long noncoding RNAs across tissues, Fa2h in hypothalamus, and Nfya in adipose tissue. Lastly, integrative analyses identified the association of BPA molecular signatures with cardiometabolic phenotypes in mouse and human. Our multitissue, multiomics investigation provides strong evidence that BPA perturbs diverse molecular networks in central and peripheral tissues and offers insights into the molecular targets that link BPA to human cardiometabolic disorders.

IDOL regulates systemic energy balance through control of neuronal VLDLR expression.

Liver X receptors limit cellular lipid uptake by stimulating the transcription of Inducible Degrader of the LDL Receptor (IDOL), an E3 ubiquitin ligase that targets lipoprotein receptors for degradation. The function of IDOL in systemic metabolism is incompletely understood. Here we show that loss of IDOL in mice protects against the development of diet-induced obesity and metabolic dysfunction by altering food intake and thermogenesis. Unexpectedly, analysis of tissue-specific knockout mice revealed that IDOL affects energy balance, not through its actions in peripheral metabolic tissues (liver, adipose, endothelium, intestine, skeletal muscle), but by controlling lipoprotein receptor abundance in neurons. Single-cell RNA sequencing of the hypothalamus demonstrated that IDOL deletion altered gene expression linked to control of metabolism. Finally, we identify VLDLR rather than LDLR as the primary mediator of IDOL effects on energy balance. These studies identify a role for the neuronal IDOL-VLDLR pathway in metabolic homeostasis and diet-induced obesity.