Pathway identification through transcriptome analysis

GLIS3: A novel transcriptional regulator of mitochondrial functions and metabolic reprogramming in postnatal kidney and polycystic kidney disease

OBJECTIVES

Deficiency in the transcription factor (TF) GLI-Similar 3 (GLIS3) in humans and mice leads to the development of polycystic kidney disease (PKD). In this study, we investigate the role of GLIS3 in the regulation of energy metabolism and mitochondrial functions in relation to its role in normal kidney and metabolic reprogramming in PKD pathogenesis.

METHODS

Transcriptomics, cistromics, and metabolomics were used to obtain insights into the role of GLIS3 in the regulation of energy homeostasis and mitochondrial metabolism in normal kidney and PKD pathogenesis using GLIS3-deficient mice.

RESULTS

Transcriptome analysis showed that many genes critical for mitochondrial biogenesis, oxidative phosphorylation (OXPHOS), fatty acid oxidation (FAO), and the tricarboxylic acid (TCA) cycle, including Tfam, Tfb1m, Tfb2m, Ppargc1a, Ppargc1b, Atp5j2, Hadha, and Sdha, are significantly suppressed in kidneys from both ubiquitous and tissue-specific Glis3-deficient mice. ChIP-Seq analysis demonstrated that GLIS3 is associated with the regulatory region of many of these genes, indicating that their transcription is directly regulated by GLIS3. Cistrome analyses revealed that GLIS3 binding loci frequently located near those of hepatocyte nuclear factor 1-Beta (HNF1B) and nuclear respiratory factor 1 (NRF1) suggesting GLIS3 regulates transcription of many metabolic and mitochondrial function-related genes in coordination with these TFs. Seahorse analysis and untargeted metabolomics corroborated that mitochondrial OXPHOS utilization is suppressed in GLIS3-deficient kidneys and showed that key metabolites in glycolysis, TCA cycle, and glutamine pathways were altered indicating increased reliance on aerobic glycolysis and glutamine anaplerosis. These features of metabolic reprogramming may contribute to a bioenergetic environment that supports renal cyst formation and progression in Glis3-deficient mice kidneys.

CONCLUSIONS

We identify GLIS3 as a novel positive regulator of the transition from aerobic glycolysis to OXPHOS in normal early postnatal kidney development by directly regulating the transcription of mitochondrial metabolic genes. Loss of GLIS3 induces several features of renal cell metabolic reprogramming. Our study identifies GLIS3 as a new participant in an interconnected transcription regulatory network, that includes HNF1B and NRF1, critical in the regulation of mitochondrial-related gene expression and energy metabolism in normal postnatal kidneys and PKD pathogenesis in Glis3-deficient mice.

Role of NHERF1 in MicroRNA Landscape Changes in Aging Mouse Kidneys

MicroRNAs (miRNAs) play important roles in the regulation of cellular function and fate via post-transcriptional regulation of gene expression. Although several miRNAs are associated with physiological processes and kidney diseases, not much is known about changes in miRNAs in aging kidneys. We previously demonstrated that sodium hydrogen exchanger 1 (NHERF1) expression regulates cellular responses to cisplatin, age-dependent salt-sensitive hypertension, and sodium-phosphate cotransporter trafficking. However, the mechanisms driving these regulatory effects of NHERF1 on cellular processes are unknown. Here, we hypothesize that dysregulation of miRNA-mediated gene regulatory networks that induce fibrosis and cytokines may depend on NHERF1 expression. To address this hypothesis, we compared miRNA expression in kidneys from both male and female old (12-18-month-old) and young (4-7-month-old) wild-type (WT) and NHERF1 knockout (NHERF1-/-) mice. Our results identified that miRNAs significantly decreased in NHERF1-/- mice included miR-669m, miR-590-3p, miR-153, miR-673-3p, and miR-127. Only miR-702 significantly decreased in aged WT mice, while miR-678 decreased in both WT and NHERF1-/- old versus young mice. miR-153 was shown to downregulate transcription factors NFATc2 and NFATc3 which regulate the transcription of several cytokines. Immunohistochemistry and western blotting revealed a significant increase in nuclear NFATc2 and NFATc3 in old NHERF1-/- mice compared to old WT mice. Our data further show that expression of the cytokines IL-1β, IL-6, IL-17A, MCP1, and TNF-α significantly increased in the old NHERF1-/- mice compared to the WT mice. We conclude that loss of NHERF1 expression induces cytokine expression in the kidney through interactive regulation between miR-153 and NFATc2/NFATc3 expression.

Functions of the primary cilium in the kidney and its connection with renal diseases

The nonmotile primary cilium is a sensory structure found on most mammalian cell types that integrates multiple signaling pathways involved in tissue development and postnatal function. As such, mutations disrupting cilia activities cause a group of disorders referred to as ciliopathies. These disorders exhibit a wide spectrum of phenotypes impacting nearly every tissue. In the kidney, primary cilia dysfunction caused by mutations in polycystin 1 (Pkd1), polycystin 2 (Pkd2), or polycystic kidney and hepatic disease 1 (Pkhd1), result in polycystic kidney disease (PKD), a progressive disorder causing renal functional decline and end-stage renal disease. PKD affects nearly 1 in 1000 individuals and as there is no cure for PKD, patients frequently require dialysis or renal transplantation. Pkd1, Pkd2, and Pkhd1 encode membrane proteins that all localize in the cilium. Pkd1 and Pkd2 function as a nonselective cation channel complex while Pkhd1 protein function remains uncertain. Data indicate that the cilium may act as a mechanosensor to detect fluid movement through renal tubules. Other functions proposed for the cilium and PKD proteins in cyst development involve regulation of cell cycle and oriented division, regulation of renal inflammation and repair processes, maintenance of epithelial cell differentiation, and regulation of mitochondrial structure and metabolism. However, how loss of cilia or cilia function leads to cyst development remains elusive. Studies directed at understanding the roles of Pkd1, Pkd2, and Pkhd1 in the cilium and other locations within the cell will be important for developing therapeutic strategies to slow cyst progression.

Kidney transcriptome and cystic kidney disease genes in zebrafish

Introduction: Polycystic kidney disease (PKD) is a condition where fluid filled cysts form on the kidney which leads to overall renal failure. Zebrafish has been recently adapted to study polycystic kidney disease, because of its powerful embryology and genetics. However, there are concerns on the conservation of this lower vertebrate in modeling polycystic kidney disease. Methods: Here, we aim to assess the molecular conservation of zebrafish by searching homologues polycystic kidney disease genes and carrying transcriptome studies in this animal. Results and Discussion: We found that out of 82 human cystic kidney disease genes, 81 have corresponding zebrafish homologs. While 75 of the genes have a single homologue, only 6 of these genes have two homologs. Comparison of the expression level of the transcripts enabled us to identify one homolog over the other homolog with >70% predominance, which would be prioritized for future experimental studies. Prompted by sexual dimorphism in human and rodent kidneys, we studied transcriptome between different sexes and noted significant differences in male vs. female zebrafish, indicating that sex dimorphism also occurs in zebrafish. Comparison between zebrafish and mouse identified 10% shared genes and 38% shared signaling pathways. String analysis revealed a cluster of genes differentially expressed in male vs. female zebrafish kidneys. In summary, this report demonstrated remarkable molecular conservation, supporting zebrafish as a useful animal model for cystic kidney disease.

A compilation of bioinformatic approaches to identify novel downstream targets for the detection and prophylaxis of cancer

The paradigm of cancer genomics has been radically changed by the development in next-generation sequencing (NGS) technologies making it possible to envisage individualized treatment based on tumor and stromal cells genome in a clinical setting within a short timeframe. The abundance of data has led to new avenues for studying coordinated alterations that impair biological processes, which in turn has increased the demand for bioinformatic tools for pathway analysis. While most of this work has been concentrated on optimizing certain algorithms to obtain quicker and more accurate results. Large volumes of these existing algorithm-based data are difficult for the biologists and clinicians to access, download and reanalyze them. In the present study, we have listed the bioinformatics algorithms and user-friendly graphical user interface (GUI) tools that enable code-independent analysis of big data without compromising the quality and time. We have also described the advantages and drawbacks of each of these platforms. Additionally, we emphasize the importance of creating new, more user-friendly solutions to provide better access to open data and talk about relevant problems like data sharing and patient privacy.

Caspase-1 and the inflammasome promote polycystic kidney disease progression

We and others have previously shown that the presence of renal innate immune cells can promote polycystic kidney disease (PKD) progression. In this study, we examined the influence of the inflammasome, a key part of the innate immune system, on PKD. The inflammasome is a system of molecular sensors, receptors, and scaffolds that responds to stimuli like cellular damage or microbes by activating Caspase-1, and generating critical mediators of the inflammatory milieu, including IL-1β and IL-18. We provide evidence that the inflammasome is primed in PKD, as multiple inflammasome sensors were upregulated in cystic kidneys from human ADPKD patients, as well as in kidneys from both orthologous (PKD1 RC/RC or RC/RC) and non-orthologous (jck) mouse models of PKD. Further, we demonstrate that the inflammasome is activated in female RC/RC mice kidneys, and this activation occurs in renal leukocytes, primarily in CD11c+ cells. Knock-out of Casp1, the gene encoding Caspase-1, in the RC/RC mice significantly restrained cystic disease progression in female mice, implying sex-specific differences in the renal immune environment. RNAseq analysis implicated the promotion of MYC/YAP pathways as a mechanism underlying the pro-cystic effects of the Caspase-1/inflammasome in females. Finally, treatment of RC/RC mice with hydroxychloroquine, a widely used immunomodulatory drug that has been shown to inhibit the inflammasome, protected renal function specifically in females and restrained cyst enlargement in both male and female RC/RC mice. Collectively, these results provide evidence for the first time that the activated Caspase-1/inflammasome promotes cyst expansion and disease progression in PKD, particularly in females. Moreover, the data suggest that this innate immune pathway may be a relevant target for therapy in PKD.

Exposure to Veterinary Antibiotics via Food Chain Disrupts Gut Microbiota and Drives Increased Escherichia coli Virulence and Drug Resistance in Young Adults

Exposure to veterinary antibiotics (VAs) and preferred as veterinary antibiotics (PVAs) via the food chain is unavoidable for their extensive use not only for treating bacterial infections, but also for use as growth promoters in livestock and aquaculture. One of the consequences is the disturbance of gut microbiota. However, its impact on the virulence and drug resistance of opportunistic pathogens is still unclear. In this study, a total of 26 antibiotics were detected in the urine of 300 young undergraduates in Anhui Province. We found that excessive intake of milk was positively correlated to high levels of VAs and PVAs. It led to the dysbiosis of gut microbiota characterized by high abundance of Bacteroidetes and Proteobacteria. The increase in Proteobacteria was mainly due to a single operational taxonomic unit (OTU) of Escherichia coli (E. coli). We isolated several E. coli strains from participants and compared their drug resistance and virulence using PCR assay and virulence-related assays. We observed that exposure to high levels of VAs and PVAs induced more resistant genes and drove E. coli strain to become more virulent. At last, we conducted transcriptome analysis to investigate the molecular mechanism of virulent and drug-resistant regulators in the highly virulent E. coli strain. We noted that there were multiple pathways involved in the drug resistance and virulence of the highly virulent strain. Our results demonstrated that participants with high-level VAs and PVAs exposure have a disrupted gut microbiota following the appearance of highly drug-resistant and virulent E. coli and, therefore may be at elevated risk for long-term health complications.

Single-Cell RNA Sequencing with Spatial Transcriptomics of Cancer Tissues

Single-cell RNA sequencing (RNA-seq) techniques can perform analysis of transcriptome at the single-cell level and possess an unprecedented potential for exploring signatures involved in tumor development and progression. These techniques can perform sequence analysis of transcripts with a better resolution that could increase understanding of the cellular diversity found in the tumor microenvironment and how the cells interact with each other in complex heterogeneous cancerous tissues. Identifying the changes occurring in the genome and transcriptome in the spatial context is considered to increase knowledge of molecular factors fueling cancers. It may help develop better monitoring strategies and innovative approaches for cancer treatment. Recently, there has been a growing trend in the integration of RNA-seq techniques with contemporary omics technologies to study the tumor microenvironment. There has been a realization that this area of research has a huge scope of application in translational research. This review article presents an overview of various types of single-cell RNA-seq techniques used currently for analysis of cancer tissues, their pros and cons in bulk profiling of transcriptome, and recent advances in the techniques in exploring heterogeneity of various types of cancer tissues. Furthermore, we have highlighted the integration of single-cell RNA-seq techniques with other omics technologies for analysis of transcriptome in their spatial context, which is considered to revolutionize the understanding of tumor microenvironment.

Multiomic identification of factors associated with progression to cystic kidney disease in mice with nephron Ift88 disruption

Ift88 gene mutations cause primary cilia loss and polycystic kidney disease (PKD) in mice. Nephron intraflagellar transport protein 88 (Ift88) knockout (KO) at 2 mo postnatal does not affect renal histology at 4 mo postnatal and causes PKD only in males by 11 mo postnatal. To identify factors associated with PKD development, kidneys from 4-mo-old male and female control and Ift88 KO mice underwent transcriptomic, proteomic, Western blot, metabolomic, and lipidomic analyses. mRNAs involved in extracellular matrix (ECM) synthesis and degradation were selectively upregulated in male KO mice. Proteomic analysis was insufficiently sensitive to detect most ECM components, while Western blot analysis paradoxically revealed reduced fibronectin and collagen type I in male KO mice. Only male KO mice had upregulated mRNAs encoding fibrinogen subunits and receptors for vascular endothelial growth factor and platelet-derived growth factor; period 2, period 3, and nuclear receptor subfamily 1 group D member 1 clock mRNAs were selectively decreased in male KO mice. Proteomic, metabolomic, and lipidomic analyses detected a relative (vs. the same-sex control) decrease in factors involved in fatty acid β-oxidation in female KO mice, while increased or unchanged levels in male KO mice, including medium-chain acyl-CoA dehydrogenase, 3-hydroxybutyrate, and acylcarnitine. Three putative mRNA biomarkers of cystogenesis in male Ift88 KO mice (similar control levels between sexes and uniquely altered by KO in males) were identified, including high levels (fibrinogen α-chain and stromal cell-derived factor 2-like 1) and low levels (BTG3-associated nuclear protein) in male KO mice. These findings suggest that relative alterations in renal ECM metabolism, fatty acid β-oxidation, and other pathways precede cystogenesis in Ift88 KO mice. In addition, potential novel biomarkers of cystogenesis in Ift88 KO mice have been identified.NEW & NOTEWORTHY Male, but not female, mice with nephron intraflagellar transport protein 88 (Ift88) gene knockout (KO) develop polycystic kidneys by ∼1 yr postnatal. We performed multiomic analysis of precystic male and female Ift88 KO and control kidneys. Precystic male Ift88 KO mice exhibited differential alterations (vs. females) in mRNA, proteins, metabolites, and/or lipids associated with renal extracellular matrix metabolism, fatty acid β-oxidation, circadian rhythm, and other pathways. These findings suggest targets for evaluation in the pathogenesis of Ift88 KO polycystic kidneys.

Identification of pathological transcription in autosomal dominant polycystic kidney disease epithelia

Autosomal dominant polycystic kidney disease (ADPKD) affects more than 12 million people worldwide. Mutations in PKD1 and PKD2 cause cyst formation through unknown mechanisms. To unravel the pathogenic mechanisms in ADPKD, multiple studies have investigated transcriptional mis-regulation in cystic kidneys from patients and mouse models, and numerous dysregulated genes and pathways have been described. Yet, the concordance between studies has been rather limited. Furthermore, the cellular and genetic diversity in cystic kidneys has hampered the identification of mis-expressed genes in kidney epithelial cells with homozygous PKD mutations, which are critical to identify polycystin-dependent pathways. Here we performed transcriptomic analyses of Pkd1- and Pkd2-deficient mIMCD3 kidney epithelial cells followed by a meta-analysis to integrate all published ADPKD transcriptomic data sets. Based on the hypothesis that Pkd1 and Pkd2 operate in a common pathway, we first determined transcripts that are differentially regulated by both genes. RNA sequencing of genome-edited ADPKD kidney epithelial cells identified 178 genes that are concordantly regulated by Pkd1 and Pkd2. Subsequent integration of existing transcriptomic studies confirmed 31 previously described genes and identified 61 novel genes regulated by Pkd1 and Pkd2. Cluster analyses then linked Pkd1 and Pkd2 to mRNA splicing, specific factors of epithelial mesenchymal transition, post-translational protein modification and epithelial cell differentiation, including CD34, CDH2, CSF2RA, DLX5, HOXC9, PIK3R1, PLCB1 and TLR6. Taken together, this model-based integrative analysis of transcriptomic alterations in ADPKD annotated a conserved core transcriptomic profile and identified novel candidate genes for further experimental studies.