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.

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.

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.

•

Development of PharmOmics, a platform for drug repositioning and toxicity prediction

•

Contains >18000 species/tissue-specific gene signatures for 941 drugs and chemicals

•

Benchmarked and validated network-based drug repositioning and toxicity prediction

•

PharmOmics is freely accessible via an online web server to facilitate user access

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.

ScRNA-seq Reveals a Role of Mammary Luminal Epithelium in Adipocyte Adaptations

Almost four decades of research suggest a dynamic role of ductal epithelial cells in adipocyte adaptation in mammary gland white adipose tissue (mgWAT), but factors that mediate such communication are not known. Here, we identify a complex intercellular crosstalk in mgWAT revealed by single-cell RNA-seq (scRNA-seq) and comprehensive data analysis suggest that epithelial luminal cells during cold exposure undergo major transcriptomic changes that lead to the expression of an array of genes that encode for secreted factors involved in adipose metabolism such as Adropin (Enho), neuregulin 4 (Nrg4), angiopoietin-like 4 (Angptl4), lipocalin 2 (Lcn2), milk fat globule-EGF factor 8 (Mfge8), Insulin-like growth factor-binding protein 1 (Igfbp1), and haptoglobin (Hp). To define the mammary epithelial secretome, we coin the phrase “mammokines”. We validated our cluster annotations and cluster-specific transcriptomics using eight different adipose scRNA-seq datasets including Tabula Muris and Tabula Muris Senis. In situ mRNA hybridization and ex vivo isolated mgWAT luminal cells show highly localized expression of mammokines in mammary ducts. Trajectory inference demonstrates that cold-exposed luminal cells have similar transcriptional profiles to lactation post-involution (PI), a phase defined by reappearance and maintenance of adipocytes in the mammary gland. Concomitantly, we found that under cold exposure female mgWAT maintains more adipogenic and less thermogenic potential than male scWAT and ex vivo removal of luminal epithelial cells from mgWAT markedly potentiates beige adipocyte differentiation. Conditioned media from isolated mammary epithelial cells treated with isoproterenol suppressed thermogenesis in differentiated beige/brown adipocytes and treatment of beige/brown differentiated adipocyte with mammokine LCN2 suppresses thermogenesis and increases adipogenesis. Finally, we find that mice lacking LCN2 show markedly higher cold-dependent thermogenesis in mgWAT than controls, and reconstitution of LCN2 in the mgWAT of LCN2 knockout mice promotes inhibition of thermogenesis. These results show a previously unknown role of mammary epithelium in adipocyte metabolism and suggest a potentially redundant evolutionary role of mammokines in maintaining mgWAT adiposity during cold exposure. Our data highlight mammary gland epithelium as a highly active metabolic cell type and mammokines could have broader implications in mammary gland physiology and lipid metabolism.

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.

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.

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.

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.

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.

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.

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.

Changes in microRNA-mRNA Signatures Agree with Morphological, Physiological, and Behavioral Changes in Larval Mahi-Mahi Treated with Deepwater Horizon Oil

In this study, we performed a systematic evaluation of global microRNA-mRNA interactions associated with the developmental toxicity of Deepwater Horizon oil using a combination of integrated mRNA and microRNA deep sequencing, expression profiling, gene ontology enrichment, and functional predictions by a series of advanced bioinformatic tools. After exposure to water accommodated fraction (WAF) of both weathered slick oil (0.5%, 1%, and 2%) and source oil (0.125%, 0.25%, and 0.5%) from the Deep Water Horizon oil spill, four dose-dependent miRNAs were identified, including three up-regulated (miR-23b, miR-34b, and miR-181b) and one down-regulated miRNAs (miR-203a) in mahi-mahi hatchings exposed from 6 h postfertilization (hpf) to 48 hpf. Consistent with morphological, physiological, and behavioral changes, the target genes of these miRNAs were largely involved in the development of the cardiovascular, visual, nervous system and associated toxicity pathways, suggesting that miRNAs play an essential role in regulating the responses to oil exposure. The results obtained from this study improve our understanding of the role of miRNAs and their target genes in relation to dose-dependent oil toxicity and provide the potential of using miRNAs as novel biomarkers in future oil studies.

Changes in thyroid status of Menidia beryllina exposed to the antifouling booster irgarol: Impacts of temperature and salinity

The triazine-based herbicide irgarol is widely used in antifouling systems as an algicide and has been detected recently in multiple coastal environments. Studies evaluating sub-lethal responses of fish following exposure to irgarol are limited. Moreover, impacts of climate change on fish endocrinology may also contribute to the sublethal toxicity of irgarol. We assessed the effects of irgarol on thyroid endpoints in juveniles of Menidia beryllina under two different treatments of salinity (10 and 20 ‰) and two temperatures (10 and 20°C). Condition factor coefficients (K) of animals were significantly affected by 0.1 to 10 μg/L of irgarol at the higher temperature. Levels of T3 were changed in whole body homogenates from both temperatures at 10‰ following exposure to 1 to 10 μg/L. T4 levels were altered only at 10°C when animals were treated with 1 to 10 μg/L (10 ‰), and in 0.1 and 10 μg/L (20 ‰). Increased transcripts of deiodinase enzymes at 10 °C may be impacted by salinity and alter thyroid hormone homeostasis. Impact on gene expression of thyroid (α and β) and growth hormone receptors were also determined. Our results highlight the relevance of environmental variable that may impact the ecological risk of irgarol in estuarine systems.

Impacts of Salinity and Temperature on the Thyroidogenic Effects of the Biocide Diuron in Menidia beryllina

Diuron is a herbicide used in agricultural and urban settings and also as an antifouling agent. Recent studies have indicated sublethal responses of diuron in the endocrine system of fish and amphibians. Given the potential of climate change to also alter fish endocrinology, the combination of environmental stressors with diuron may contribute to its sublethal toxicity. In this study, the effects of temperature and salinity on thyroid targets of diuron were assessed in juveniles of the estuarine fish Menidia beryllina under different conditions of salinity (10 and 20‰) and temperature (10 and 20 °C). Environmentally relevant concentrations of diuron affected the growth, and the higher temperature reduced the condition factor of animals. Increased levels of T3 were observed in fish from all treatments, and at 10 °C, T4 levels were augmented at 10‰ but reduced at 20‰. Increased gene expression of deiodinases at 20‰ in both temperatures suggests the influence of salinity on the regulation of hormone imbalance via deiodination pathway activation. Decreased transcripts of thyroid and growth hormone receptors were also observed following diuron treatment. These results indicate that changes in environmental stressors may have significant impacts on the ecological risk of diuron in estuarine fish.

Developmental toxicity of hydroxylated chrysene metabolites in zebrafish embryos

One of the primary sources of polycyclic aromatic hydrocarbons (PAHs) in marine environments is oil. Photochemical oxidation and microbial transformation of PAH-containing oils can result in the formation of oxygenated products. Among the PAHs in crude oil, chrysene is one of the most persistent within the water column and may be transformed to 2- and 6-hydroxychrysene (OHCHR). Both of these compounds have been shown to activate (2-OHCHR) and antagonize (6-OHCHR) the estrogen receptor (ER). Previous studies in our lab have shown that estrogen can significantly alter zebrafish development. However, little is known about the developmental toxicity of hydroxylated PAHs. Zebrafish embryos were exposed to 0.5-10μM of 2- or 6-OHCHR from 2h post-fertilization (hpf) until 76hpf. A significant decrease in survival was observed following exposure to 6-OHCHR - but not 2-OHCHR. Both OHCHRs significantly increased the percentage of overall deformities after treatment. In addition to cardiac malformations, ocular and circulatory defects were also observed in embryos exposed to both compounds, while 2-OHCHR generally resulted in a higher prevalence of effect. Moreover, treatment with 2-OHCHR resulted in a significant decrease in hemoglobin levels. ER nor G-Protein coupled estrogen receptor (GPER) antagonists and agonists did not rescue the observed defects. We also analyzed the expression of cardiac-, eye- and circulation-related genes previously shown to be affected by oil. Rhodopsin mRNA expresssion was significantly decreased by both compounds equally. However, exposure to 2-OHCHR significantly increased the expression of the hematopoietic regulator, runx1 (runt related transcription factor 1). These results indicate the toxicity of oxygenated photoproducts of PAHs and suggest that other targets and signaling pathways may contribute to developmental toxicity of weathered oil. Our findings also demonstrate the regio-selective toxicity of hydroxy-PAHs in the effects on eye and circulatory development and raise the need to identify mechanisms and ecological risks of oxy-PAHs to fish populations.

Contribution of G protein-coupled estrogen receptor 1 (GPER) to 17β-estradiol-induced developmental toxicity in zebrafish

Exposure to 17β-estradiol (E2) influences the regulation of multiple signaling pathways, and E2-mediated disruption of signaling events during early development can lead to malformations such as cardiac defects. In this study, we investigated the potential role of the G-protein estrogen receptor 1 (GPER) in E2-induced developmental toxicity. Zebrafish embryos were exposed to E2 from 2h post fertilization (hpf) to 76 hpf with subsequent transcriptional measurements of heart and neural crest derivatives expressed 2 (hand2), leucine rich repeat containing 10 (lrrc10), and gper at 12, 28 and 76 hpf. Alteration in the expression of lrrc10, hand2 and gper was observed at 12 hpf and 76 hpf, but not at 28 hpf. Expression of these genes was also altered after exposure to G1 (a GPER agonist) at 76 hpf. Expression of lrrc10, hand2 and gper all coincided with the formation of cardiac edema at 76 hpf as well as other developmental abnormalities. While co-exposure of G1 with G36 (a GPER antagonist) rescued G1-induced abnormalities and altered gene expression, co-exposure of E2 with G36, or ICI 182,780 (an estrogen receptor antagonist) did not rescue E2-induced cardiac deformities or gene expression. In addition, no effects on the concentrations of downstream ER and GPER signaling molecules (cAMP or calcium) were observed in embryo homogenates after E2 treatment. These data suggest that the impacts of E2 on embryonic development at this stage are complex and may involve multiple receptor and/or signaling pathways.

Measuring the mass, volume, and density of microgram-sized objects in fluid

Measurements of an object’s fundamental physical properties like mass, volume, and density can offer valuable insights into the composition and state of the object. However, many important biological samples reside in a liquid environment where it is difficult to accurately measure their physical properties. We show that by using a simple piece of glass tubing and some inexpensive off-the-shelf electronics, we can create a sensor that can measure the mass, volume, and density of microgram-sized biological samples in their native liquid environment. As a proof-of-concept, we use this sensor to measure mass changes in zebrafish embryos reacting to toxicant exposure, density changes in seeds undergoing rehydration and germination, and degradation rates of biomaterials used in medical implants. Since all objects have these physical properties, this sensor has immediate applications in a wide variety of different fields including developmental biology, toxicology, materials science, plant science, and many others.

Influence of Temperature on the Thyroidogenic Effects of Diuron and Its Metabolite 3,4-DCA in Tadpoles of the American Bullfrog (Lithobates catesbeianus)

Temperature is a key variable affecting the timing of amphibian metamorphosis from tadpoles to tetrapods, through the production and subsequent function of thyroid hormones (TH). Thyroid function can be impaired by environmental contaminants as well as temperature. Tadpoles can experience large temperature fluctuations in their habitats and many species are distributed in areas that may be impacted by agriculture. Diuron is a widely used herbicide detected in freshwater ecosystems and may impact endocrine function in aquatic organisms. We evaluated the influence of temperature (28 and 34 °C) on the action of diuron and its metabolite 3,4-dichloroaniline (3,4-DCA) on thyroid function and metamorphosis in tadpoles of Lithobates catesbeianus. Exposure to both compounds induced more pronounced changes in gene expression and plasma 3,3',5-triiodothyronine (T₃) concentrations in tadpoles treated at higher temperature. T₃ concentrations were increased in tadpoles exposed to 200 ng/L of diuron at 34 °C and an acceleration of metamorphosis was observed for the same group. Transcriptomic responses included alteration of thyroid hormone induced bZip protein (thibz), deiodinases (dio2, dio3), thyroid receptors (trα, trβ) and Krüppel-like factor 9 (klf9), suggesting regulation by temperature on TH-gene expression. These results suggest that environmental temperature should be considered in risk assessments of environmental contaminants for amphibian species.

Time- and Oil-Dependent Transcriptomic and Physiological Responses to Deepwater Horizon Oil in Mahi-Mahi (Coryphaena hippurus) Embryos and Larvae

The Deepwater Horizon (DWH) oil spill contaminated the spawning habitats for numerous commercially and ecologically important fishes. Exposure to the water accommodated fraction (WAF) of oil from the spill has been shown to cause cardiac toxicity during early developmental stages across fishes. To better understand the molecular events and explore new pathways responsible for toxicity, RNA sequencing was performed in conjunction with physiological and morphological assessments to analyze the time-course (24, 48, and 96 h post fertilization (hpf)) of transcriptional and developmental responses in embryos/larvae of mahi-mahi exposed to WAF of weathered (slick) and source DWH oils. Slick oil exposure induced more pronounced changes in gene expression over time than source oil exposure. Predominant transcriptomic responses included alteration of EIF2 signaling, steroid biosynthesis, ribosome biogenesis and activation of the cytochrome P450 pathway. At 96 hpf, slick oil exposure resulted in significant perturbations in eye development and peripheral nervous system, suggesting novel targets in addition to the heart may be involved in the developmental toxicity of DHW oil. Comparisons of changes of cardiac genes with phenotypic responses were consistent with reduced heart rate and increased pericardial edema in larvae exposed to slick oil but not source oil.

Insights into the mechanism of cell death induced by saporin delivered into cancer cells by an antibody fusion protein targeting the transferrin receptor 1

We previously developed an antibody-avidin fusion protein (ch128.1Av) that targets the human transferrin receptor 1 (TfR1) and exhibits direct cytotoxicity against malignant B cells in an iron-dependent manner. ch128.1Av is also a delivery system and its conjugation with biotinylated saporin (b-SO6), a plant ribosome-inactivating toxin, results in a dramatic iron-independent cytotoxicity, both in malignant cells that are sensitive or resistant to ch128.1Av alone, in which the toxin effectively inhibits protein synthesis and triggers caspase activation. We have now found that the ch128.1Av/b-SO6 complex induces a transcriptional response consistent with oxidative stress and DNA damage, a response that is not observed with ch128.1Av alone. Furthermore, we show that the antioxidant N-acetylcysteine partially blocks saporin-induced apoptosis suggesting that oxidative stress contributes to DNA damage and ultimately saporin-induced cell death. Interestingly, the toxin was detected in nuclear extracts by immunoblotting, suggesting the possibility that saporin might induce direct DNA damage. However, confocal microscopy did not show a clear and consistent pattern of intranuclear localization. Finally, using the long-term culture-initiating cell assay we found that ch128.1Av/b-SO6 is not toxic to normal human hematopoietic stem cells suggesting that this critical cell population would be preserved in therapeutic interventions using this immunotoxin.

Effect of NRG1, GDNF, EGF and NGF in the migration of a Schwann cell precursor line

The Schwann cells are the myelinating glia of the peripheral nervous system that originated during development from the highly motile neural crest. However, we do not know what the guidance signals for the Schwann cell precursors are. Therefore, we set to test some of the known neurotrophins that are expressed early in developing embryos and have been shown to be critical for the survival and patterning of developing glia and neurons. The goal of this study was to determine more specifically if GDNF, NRG1 and NGF are chemoattractants and/or chemokinetic molecules for a Schwann cell precursor line, the Spl201. We performed live chemoattraction assays, with imaging and also presented these molecules as part of their growing substrate. Our results show for the first time that GDNF and NRG1 are potent chemoattractive and chemokinetic molecules for these cells while NGF is a chemokinetic molecule stimulating their motility.