Search

Abstract Missing covariate data is a common problem that has not been addressed in observational studies of gene expression. Here, we present a multiple imputation method that accommodates high dimensional gene expression data by incorporating principal component analysis of the transcriptome into the multiple imputation prediction models to avoid bias. Simulation studies using three datasets show that this method outperforms complete case and single imputation analyses at uncovering true positive differentially expressed genes, limiting false discovery rates, and minimizing bias. This method is easily implemented via an R Bioconductor package, RNAseqCovarImpute that integrates with the limma-voom pipeline for differential expression analysis.

Abstract The microgravity environment aboard the International Space Station (ISS) provides a unique stressor that can help understand underlying cellular and molecular drivers of pathological changes observed in astronauts with the ultimate goals of developing strategies to enable long- term spaceflight and better treatment of diseases on Earth. We used this unique environment to evaluate the effects of microgravity on kidney proximal tubule epithelial cell (PTEC) response to serum exposure and vitamin D biotransformation capacity. To test if microgravity alters the pathologic response of the proximal tubule to serum exposure, we treated PTECs cultured in a microphysiological system (PT-MPS) with human serum and measured biomarkers of toxicity and inflammation (KIM-1 and IL-6) and conducted global transcriptomics via RNAseq on cells undergoing flight (microgravity) and respective controls (ground). Given the profound bone loss observed in microgravity and PTECs produce the active form of vitamin D, we treated 3D cultured PTECs with 25(OH)D3 (vitamin D) and monitored vitamin D metabolite formation, conducted global transcriptomics via RNAseq, and evaluated transcript expression of CYP27B1, CYP24A1, or CYP3A5 in PTECs undergoing flight (microgravity) and respective ground controls. We demonstrated that microgravity neither altered PTEC metabolism of vitamin D nor did it induce a unique response of PTECs to human serum, suggesting that these fundamental biochemical pathways in the kidney proximal tubule are not significantly altered by short-term exposure to microgravity. Given the prospect of extended spaceflight, more study is needed to determine if these responses are consistent with extended (>6 months) exposure to microgravity.

Background: The first spatially resolved transcriptomics platforms, GeoMx (Nanostring) and Visium (10x Genomics) were launched in 2019 and were recognized as the method of the year by Nature Methods in 2020. The subsequent refinement and expansion of these and other technologies to increase -plex, work with formalin-fixed paraffin-embedded tissue, and analyze protein in addition to gene expression have only added to their significance and impact on the biomedical sciences. In this perspective, we focus on two platforms for spatial transcriptomics, GeoMx and Visium, and how these platforms have been used to provide novel insight into kidney disease. The choice of platform will depend largely on experimental questions and design. The application of these technologies to clinically sourced biopsies presents the opportunity to identify specific tissue biomarkers that help define disease etiology and more precisely target therapeutic interventions in the future. Summary: In this review, we provide a description of the existing and emerging technologies that can be used to capture spatially resolved gene and protein expression data from tissue. These technologies have provided new insight into the spatial heterogeneity of diseases, how reactions to disease are distributed within a tissue, which cells are affected, and molecular pathways that predict disease and response to therapy. Key Message: The upcoming years will see intense use of spatial transcriptomics technologies to better define the pathophysiology of kidney diseases and develop novel diagnostic tests to guide personalized treatments for patients.

Abstract Preterm birth remains the leading identifiable risk factor for cerebral palsy (CP), a devastating form of motor impairment due to developmental brain injury occurring around the time of birth. We performed genome wide methylation and whole transcriptome analyses to elucidate the early pathogenesis of CP in extremely low gestational age neonates (ELGANs). We evaluated peripheral blood cell specimens collected during a randomized trial of erythropoietin for neuroprotection in the ELGAN (PENUT Trial, NCT# 01378273). DNA methylation data were generated from 94 PENUT subjects (n = 47 CP vs. n = 47 Control) on day 1 and 14 of life. Gene expression data were generated from a subset of 56 subjects. Only one differentially methylated region was identified for the day 1 to 14 change between CP versus no CP, without evidence for differential gene expression of the associated gene RNA Pseudouridine Synthase Domain Containing 2. iPathwayGuide meta-analyses identified a relevant upregulation of JAK1 expression in the setting of decreased methylation that was observed in control subjects but not CP subjects. Evaluation of whole transcriptome data identified several top pathways of potential clinical relevance including thermogenesis, ferroptossis, ribosomal activity and other neurodegenerative conditions that differentiated CP from controls.

Air pollution might affect atherosclerosis through DNA methylation changes in cells crucial to atherosclerosis, such as monocytes. We conducted an epigenome-wide study of DNA methylation in CD14+ monocytes and long-term ambient air pollution exposure in adults participating in the Multi-Ethnic Study of Atherosclerosis (MESA). We also assessed the association between differentially methylated signals and cis-gene expression. Using spatiotemporal models, one-year average concentrations of outdoor fine particulate matter (PM2.5) and oxides of nitrogen (NOX) were estimated at participants’ homes. We assessed DNA methylation and gene expression using Illumina 450k and HumanHT-12 v4 Expression BeadChips, respectively (n = 1,207). We used bump hunting and site-specific approaches to identify differentially methylated signals (false discovery rate of 0.05) and used linear models to assess associations between differentially methylated signals and cis-gene expression. Four differentially methylated regions (DMRs) located on chromosomes 5, 6, 7, and 16 (within or near SDHAP3, ZFP57, HOXA5, and PRM1, respectively) were associated with PM2.5. The DMRs on chromosomes 5 and 6 also associated with NOX. The DMR on chromosome 5 had the smallest p-value for both PM2.5 (p = 1.4×10−6) and NOX (p = 7.7×10−6). Three differentially methylated CpGs were identified for PM2.5, and cg05926640 (near TOMM20) had the smallest p-value (p = 5.6×10−8). NOX significantly associated with cg11756214 within ZNF347 (p = 5.6×10−8). Several differentially methylated signals were also associated with cis-gene expression. The DMR located on chromosome 7 was associated with the expression of HOXA5, HOXA9, and HOXA10. The DMRs located on chromosomes 5 and 16 were associated with expression of MRPL36 and DEXI, respectively. The CpG cg05926640 was associated with expression of ARID4B, IRF2BP2, and TOMM20. We identified differential DNA methylation in monocytes associated with long-term air pollution exposure. Methylation signals associated with gene expression might help explain how air pollution contributes to cardiovascular disease.

Summary: Fibrosis-driven solid organ failure is an enormous burden on global health. Spiny mice (Acomys) are terrestrial mammals that can regenerate severe skin wounds without scars to avoid predation. Whether spiny mice also regenerate internal organ injuries is unknown. Here, we show that despite equivalent acute obstructive or ischemic kidney injury, spiny mice fully regenerate nephron structure and organ function without fibrosis, whereas C57Bl/6 or CD1 mice progress to complete organ failure with extensive renal fibrosis. Two mechanisms for vertebrate regeneration have been proposed that emphasize either extrinsic (pro-regenerative macrophages) or intrinsic (surviving cells of the organ itself) controls. Comparative transcriptome analysis revealed that the Acomys genome appears poised at the time of injury to initiate regeneration by surviving kidney cells, whereas macrophage accumulation was not detected until about day 7. Thus, we provide evidence for rapid activation of a gene expression signature for regenerative wound healing in the spiny mouse kidney.

Abstract Background Developmental exposure to particulate matter air pollution is harmful to cardiovascular health, but the mechanisms by which this exposure mediates susceptibility to heart disease is poorly understood. We have previously shown, in a mouse model, that gestational exposure to diesel exhaust (DE) results in increased cardiac hypertrophy, fibrosis and susceptibility to heart failure in the adult offspring following transverse aortic constriction. Results In this study, we have analyzed gene expression in neonatal cardiomyocytes after gestational exposure by RNA-sequencing and have identified 300 genes that are dysregulated, including many involved in cardiac metabolism. We subsequently determined that these cardiomyocytes exhibit reduced metabolic activity as measured by Seahorse extracellular flux analysis. We also surveyed for modifications in DNA methylation at global regulatory regions using reduced representation bisulfite sequencing and found hypomethylation of DNA in neonatal cardiomyocytes isolated from in utero DE exposed neonates. Conclusion We have demonstrated that in utero exposure to diesel exhaust alters the neonatal cardiomyocyte transcriptional and epigenetic landscapes, as well as the metabolic capability of these cells. Understanding how exposure alters the developing heart through dysregulation of gene expression, metabolism and DNA methylation is vital for identifying therapeutic interventions for air pollution-related heart failure.

Summary: The endothelium first forms in the blood islands in the extra-embryonic yolk sac and then throughout the embryo to establish circulatory networks that further acquire organ-specific properties during development to support diverse organ functions. Here, we investigated the properties of endothelial cells (ECs), isolated from four human major organs—the heart, lung, liver, and kidneys—in individual fetal tissues at three months' gestation, at gene expression, and at cellular function levels. We showed that organ-specific ECs have distinct expression patterns of gene clusters, which support their specific organ development and functions. These ECs displayed distinct barrier properties, angiogenic potential, and metabolic rate and support specific organ functions. Our findings showed the link between human EC heterogeneity and organ development and can be exploited therapeutically to contribute in organ regeneration, disease modeling, as well as guiding differentiation of tissue-specific ECs from human pluripotent stem cells. : Stem Cells Research; Developmental Biology; Biology of Human Development Subject Areas: Stem Cells Research, Developmental Biology, Biology of Human Development

Neuroinflammation in utero may result in lifelong neurological disabilities. Astrocytes play a pivotal role in this process, but the mechanisms are poorly understood. No early postnatal treatment strategies exist to enhance neuroprotective potential of astrocytes. We hypothesized that agonism on α7 nicotinic acetylcholine receptor (α7nAChR) in fetal astrocytes will augment their neuroprotective transcriptome profile, while the inhibition of α7nAChR will achieve the opposite. Using an in vivo–in vitro model of developmental programming of neuroinflammation induced by lipopolysaccharide (LPS), we validated this hypothesis in primary fetal sheep astrocytes cultures re-exposed to LPS in the presence of a selective α7nAChR agonist or antagonist. Our RNAseq findings show that a pro-inflammatory astrocyte transcriptome phenotype acquired in vitro by LPS stimulation is reversed with α7nAChR agonistic stimulation. Conversely, α7nAChR inhibition potentiates the pro-inflammatory astrocytic transcriptome phenotype. Furthermore, we conducted a secondary transcriptome analysis against the identical α7nAChR experiments in fetal sheep primary microglia cultures. Similar to findings in fetal microglia, in fetal astrocytes we observed a memory effect of in vivo exposure to inflammation, expressed in a perturbation of the iron homeostasis signaling pathway (hemoxygenase 1, HMOX1), which persisted under pre-treatment with α7nAChR antagonist but was reversed with α7nAChR agonist. For both glia cell types, common pathways activated due to LPS included neuroinflammation signaling and NF-κB signaling in some, but not all comparisons. However, overall, the overlap on the level of signaling pathways was rather minimal. Astrocytes, not microglia—the primary immune cells of the brain, were characterized by unique inhibition patterns of STAT3 pathway due to agonistic stimulation of α7nAChR prior to LPS exposure. Lastly, we discuss the implications of our findings for fetal and postnatal brain development.

Epigenetic factors (EFs) regulate multiple aspects of cerebral cortex development, including proliferation, differentiation, laminar fate, and regional identity. The same neurodevelopmental processes are also regulated by transcription factors (TFs), notably the Pax6→ Tbr2→ Tbr1 cascade expressed sequentially in radial glial progenitors (RGPs), intermediate progenitors, and postmitotic projection neurons, respectively. Here, we studied the EF landscape and its regulation in embryonic mouse neocortex. Microarray and in situ hybridization assays revealed that many EF genes are expressed in specific cortical cell types, such as intermediate progenitors, or in rostrocaudal gradients. Furthermore, many EF genes are directly bound and transcriptionally regulated by Pax6, Tbr2, or Tbr1, as determined by chromatin immunoprecipitation-sequencing and gene expression analysis of TF mutant cortices. Our analysis demonstrated that Pax6, Tbr2, and Tbr1 form a direct feedforward genetic cascade, with direct feedback repression. Results also revealed that each TF regulates multiple EF genes that control DNA methylation, histone marks, chromatin remodeling, and non-coding RNA. For example, Tbr1 activates Rybp and Auts2 to promote the formation of non-canonical Polycomb repressive complex 1 (PRC1). Also, Pax6, Tbr2, and Tbr1 collectively drive massive changes in the subunit isoform composition of BAF chromatin remodeling complexes during differentiation: for example, a novel switch from Bcl7c (Baf40c) to Bcl7a (Baf40a), the latter directly activated by Tbr2. Of 11 subunits predominantly in neuronal BAF, 7 were transcriptionally activated by Pax6, Tbr2, or Tbr1. Using EFs, Pax6→ Tbr2→ Tbr1 effect persistent changes of gene expression in cell lineages, to propagate features such as regional and laminar identity from progenitors to neurons.

Mechanisms driving the initiation of brain folding are incompletely understood. We have previously characterized mouse models recapitulating human PIK3CA-related brain overgrowth, epilepsy, dysplastic gyrification and hydrocephalus (Roy et al., 2015). Using the same, highly regulatable brain-specific model, here we report PI3K-dependent mechanisms underlying gyrification of the normally smooth mouse cortex, and hydrocephalus. We demonstrate that a brief embryonic Pik3ca activation was sufficient to drive subtle changes in apical cell adhesion and subcellular Yap translocation, causing focal proliferation and subsequent initiation of the stereotypic ‘gyrification sequence’, seen in naturally gyrencephalic mammals. Treatment with verteporfin, a nuclear Yap inhibitor, restored apical surface integrity, normalized proliferation, attenuated gyrification and rescued the associated hydrocephalus, highlighting the interrelated role of regulated PI3K-Yap signaling in normal neural-ependymal development. Our data defines apical cell-adhesion as the earliest known substrate for cortical gyrification. In addition, our preclinical results support the testing of Yap-related small-molecule therapeutics for developmental hydrocephalus.

Embryonal rhabdomyosarcoma (ERMS) is the most common soft tissue cancer in children. The prognosis of patients with relapsed or metastatic disease remains poor. ERMS genomes show few recurrent mutations, suggesting that other molecular mechanisms such as epigenetic regulation might play a major role in driving ERMS tumor biology. In this study, we have demonstrated the diverse roles of histone deacetylases (HDACs) in the pathogenesis of ERMS by characterizing effects of HDAC inhibitors, trichostatin A (TSA) and suberoylanilide hydroxamic acid (SAHA; also known as vorinostat) in vitro and in vivo. TSA and SAHA suppress ERMS tumor growth and progression by inducing myogenic differentiation as well as reducing the self-renewal and migratory capacity of ERMS cells. Differential expression profiling and pathway analysis revealed downregulation of key oncogenic pathways upon HDAC inhibitor treatment. By gain-of-function, loss-of-function, and chromatin immunoprecipitation (ChIP) studies, we show that Notch1- and EphrinB1-mediated pathways are regulated by HDACs to inhibit differentiation and enhance migratory capacity of ERMS cells, respectively. Our study demonstrates that aberrant HDAC activity plays a major role in ERMS pathogenesis. Druggable targets in the molecular pathways affected by HDAC inhibitors represent novel therapeutic options for ERMS patients.

The goal of this study was to compare the effects of two inflammatory modulators, erythropoietin (EPO) and anakinra, on functional recovery and brain gene expression following a cortical contusion impact (CCI) injury. Dosage regimens were designed to provide serum concentrations in the range obtained with clinically approved doses. Functional recovery was assessed using both motor and spatial learning tasks and neuropathological measurements conducted in the cortex and hippocampus. Microarray-based transcriptional profiling was used to determine the effect on gene expression at 24 h, 72 h and 7 days post-CCI. Ingenuity Pathway Analysis was used to evaluate the effect on relevant functional categories. EPO and anakinra treatment resulted in significant changes in brain gene expression in the CCI model demonstrating acceptable brain penetration. At all three time points, EPO treatment resulted in significantly more differentially expressed genes than anakinra. For anakinra at 24 h and EPO at 24 h, 72 h and 7 days, the genes in the top 3 functional categories were involved in cellular movement, inflammatory response and cell-to-cell signaling. For EPO, the majority of the genes in the top 10 canonical pathways identified were associated with inflammatory and immune signaling processes. This was true for anakinra only at 24 h post-traumatic brain injury (TBI). The immunomodulation effects of EPO and anakinra did not translate into positive effects on functional behavioral and lesion studies. Treatment with either EPO or anakinra failed to induce significant beneficial effects on recovery of function or produce any significant effects on the prevention of injury induced tissue loss at 30 days post-injury. In conclusion, treatment with EPO or anakinra resulted in significant effects on gene expression in the brain without affecting functional outcome. This suggests that targeting these inflammatory processes alone may not be sufficient for preventing secondary injuries