Doxycycline Changes the Transcriptome Profile of mIMCD3 Renal Epithelial Cells

Abstract

Tetracycline-inducible gene expression systems have been used successfully to study gene function in vivo and in vitro renal epithelial models but the effects of the common inducing agent, doxycycline (DOX), on gene expression are not well appreciated. Here, we evaluated the DOX effects on the transcriptome of a widely used renal epithelial cell model, mIMCD3 cells, to establish a reference. Cells were grown on permeable filter supports in the absence and presence of DOX (3 or 6 days), and genome-wide transcriptome profiles were assessed using RNA-Seq. We found DOX significantly altered the transcriptome profile, changing the abundance of 1,549 transcripts at 3 days and 2,643 transcripts at 6 days. Within 3 days of treatment, DOX significantly decreased the expression of multiple signaling pathways (ERK, cAMP, and Notch) that are associated with cell proliferation and differentiation. Genes associated with cell cycle progression were subsequently downregulated in cells treated with DOX for 6 days, as were genes involved in cellular immune response processes and several cytokines and chemokines, correlating with a remarkable repression of genes encoding cell proliferation markers. The results provide new insight into responses of renal epithelial cells to DOX and a establish a resource for DOX-mediated gene expression systems.

Introduction

Gene expression systems utilizing drug-induced trans-activation provide the means to conditionally investigate gene function in a temporal manner. They can be used to activate expression of a target gene for study, or to induce Cre-recombinase for conditional knockout cell of floxed alleles. The tetracycline-inducible gene expression system has proven especially popular. Doxycycline (DOX), a synthetic derivative of tetracycline, is widely used in these systems (Kistner et al., 1996; Das et al., 2016) because of its long half live and low cell toxicity. However, effects of DOX on a global gene expression and cellular processes in renal epithelial cell models are not known.

The absence of DOX effects on global gene expression in Saccharomyces cerevisiae (Wishart et al., 2005) initially suggested DOX might be inert. However, several more recent studies revealed DOX can change gene expression in mammalian tissues. For example, in ex vivo studies with surgically removed pterygia tissue from the eye, RNA sequencing (RNA-Seq)-based transcriptomic analysis (Larrayoz et al., 2012) revealed DOX affected the expression of mitochondrial genes, the ER stress cascade, growth factors, interleukins, cell cycle regulators, integrins, and components of the extracellular matrix. In rat aortic tissue (Lu et al., 2017), DOX was found to change the expression of other sets of genes, and these were mainly enriched in pathways that control neutrophil chemotaxis, chronic inflammatory responses, and cellular responses to mechanical stimuli and negatively regulate apoptotic processes. In vivo effects of DOX have been documented in studies to characterize the trans-activator (rtTA2) system in the mouse liver (Reboredo et al., 2008), and DOX-sensitive genes were found to enriched in cellular pathways related to cell growth and death and mitochondrial electron transport (Reboredo et al., 2008). Although the sample size still remains small, these few studies question the extent to which DOX may have generalizable responses; the heterogenous nature of responses suggest DOX action may be specifically governed by cell type or experimental condition. More studies are required with other cell types to identify common pathways. It remains unknown if DOX affects gene expression in renal epithelial cells.

Here we focus on characterizing the DOX response in the renal epithelial cell line, mIMCD3. The cells have been engineered to include a tetracycline-inducible gene expression system (Schlimpert et al., 2018; Lashhab et al., 2019; Kang et al., 2019), and the model has proven popular for studying genes in a renal epithelial environment. Here we examine how DOX affects the transcriptome profiles of mIMCD3 cells with RNA-Seq. The data should provide an informative resource for future studies with the tetracycline-inducible mIMCD3 cell line and similar in vitro renal epithelial models.

Materials and Methods

Cell Culture

Mouse kidney epithelial mIMCD3 cells derived from the inner medulla of a simian virus SV40 transgenic mouse were obtained from Maryland Polycystic Kidney Disease Research and Translation Core Center at the University of Maryland, and documented to be mycoplasma free. Cells were grown in T75 flasks with the RenaLife Epithelial Basal Medium (LM-0010, LIFELINE CELL TECHNOLOGY) supplemented with a RenaLife LifeFactors Kit (LS-1048, LIFELINE CELL TECHNOLOGY) and 5% FBS. This culture medium, which contains a nutrient blend of amino acids, vitamins, organic and inorganic supplements and salts, growth factors (0.5 μg/mL Insulin, 1 μM Epinephrine, 0.1 μg/mL Hydrocortisone, 10 nM Triiodothyronine, 10 ng/mL EGF, 5 μg/mL Transferrin, 2.4 mM L-Alanyl-L-Glutamine), and antibiotics (30 μg/mL Gentamicin, 15 ng/mL Amphotericin B), has become the choice for in vitro studies in the Maryland Cell engineering core of the NIH U54 funded Polycystic Kidney Disease Research Consortium. It was used here to provide a reference for on-going and future transcriptome profiling studies in the core. For RNA-Seq studies, cells were plated permeable filters (Polyester membrane) of a 6-well transwell plate (#3450, CORNING) and grown to confluence (10 days) before doxycycline treatment. Doxycycline (2 μg/mL) or vehicle (DMSO) added to both sides of the permeable filter and changed daily for 3 or 6 days.

Total RNA Isolation and RNA Sequencing-Based Transcriptome Profiling

Cells were lysed in Trizol reagent (15596018, Invitrogen). Total RNA was isolated from the Trizol lysate using Direct-Zol RNA Miniprep plus kit (R2070, ZYMO RESEARCH) and eluted in RNase-free water, and RNA concentration was measured using Qubit^TM RNA HS Assay Kit (Q32852, Invitrogen). To enrich mRNA, 1–3 μg of total RNA was applied to oligo dT-based mRNA isolation using NEBNext^® Poly(A) mRNA Magnetic Isolation Module (E7490, NEW ENGLAND BioLabs Inc.) according to manufacturer’s instructions. mRNA (20 ng) was used to create the cDNA libraries, using NEBNext^® Ultra II Directional RNA Library Prep Kit for Illumina (E7760, NEW ENGLAND BioLabs Inc.) and NEBNext^® Multiplex Oligos (E7335; E7500, NEW ENGLAND BioLabs Inc.). cDNA libraries were sequenced on the Illumina HiSeq 4,000 platform. Sequence reads (2 × 75 bp, paired-end) were aligned on Ensembl genome GRCm38p6 using STAR (2.6.0c).

Differential Expression Analysis and Bioinformatic Analysis

To identify differentially expressed genes between vehicle- and doxycycline-treated cells, transcript abundance was quantified using salmon (Patro et al., 2017). Differential expression (DE) analysis was carried out using edgeR (Robinson et al., 2010) on R (3.6.0). Low abundant genes with CPM (Counts Per Million) less than 1 were removed from the data set for the DE analysis. Significance of DE was determined using a modified statistical test (edgeR “glmTreat”) with a threshold of expression changes above 20% at FDR < 0.05 (Chen et al., 2016). Plots were generated using R package ggplot2. Pathway enrichment analysis of DE genes was carried out using Gene Ontology (Biological process) on Metascape platform¹ (Zhou et al., 2019). The GO term analysis was performed using NaviGO (Wei et al., 2017).

Statistics

Significance of DE was evaluated using edgeR and p-values were corrected using the Benjamini-Hochberg method. For significance of gene enrichment in the pathway analysis, q < 0.05 was considered as significant enrichment.

Data Availability

All fastq files and a raw count file from RNA-Seq were deposited in GEO (GSE171573²).

Results

Transcriptome Profiles of mIMCD3 on Different Culture Environments

RNA-Seq analysis was performed to provide a transcript database of fully polarized mIMCD3 cells on permeable filter supports in the RenalLife medium which includes growth factors and antibiotics. For these studies, the cells were grown to confluence on the permeable filter supports for 10 days, and then were treated with DMSO or DOX in the RenalLife medium for 3 or 6 days.

To establish the baseline, we first examined the control cells, treated with DMSO. We found a high correlation between the transcriptome of cells treated with DMSO at 3 and 6 days (Pearson correlation: 0.9956, Figure 1A), consistent with stable gene expression and lack of transient DMSO responses. Based comparison of normalized expression values (transcripts per million, TPM), the transcriptome exhibited a high correlation with native mIMCD3 cells grown on plastic solid supports (GSE97770) as reported by Chan et al. (2018) (Figures 1B,C). Whether grown on permeable supports or plastic solid supports, mIMCD3 equally express the conventional epithelial cell markers (Jedroszka et al., 2017), such as Tjp1, Krt18, Dsp, Muc1, and Sdc1 (Figure 1D). Other extracellular matrix components, Col4a2 and Sparc, were also more abundantly expressed in the mIMCD3 cells grown on the solid support compared to mIMCD3 cells on the filter support (Figure 1G). Interestingly, mesenchymal cell marker genes, namely Cdh2, Fn1, Sparc, and Vim, were also more abundant in mIMCD3 cells grown on the solid support compared to mIMCD3 cells grown on the filter support (Figure 1E). By contrast, in mIMCD3 cells grown on the filter supports, tight junction proteins, Cldn4, Cldn7, and Epcam, were more abundant than in cells on the solid support (Figure 1F). Thus, transcriptomic signature of mIMCD3 cells is most compatible with a differentiated epithelial state when cells are grown on permeable supports compared to cells grown on plastic supports. Yet the transcriptome profiles of mIMCD3 cells grown on permeable and non-permeable are more similar than different. Interestingly, we did not identify the enrichment of specific transcripts of medullary collecting duct in this cell line.

FIGURE 1

Nevertheless, the dataset of comprehensive transcriptome profile in mIMCD3 provides an informative resource for further utilization as the in vitro kidney epithelial cell model. To expand understanding of intracellular signaling pathways in mIMCD3 cells grown on solid supports (Valkova and Kultz, 2006; Chan et al., 2018) compared to cells grown on filters, we classified transcriptome profile into genes encoding transcription factors, ion channels, transporters, and G protein-coupled receptors (GPCRs) in Tables 1–4 (Supplementary Table 1). Many of the listed genes have been identified in previous studies. For example, Hnf1b (hepatocyte nuclear factor-1 beta, TPM: 43.0) (Aboudehen et al., 2017), Pax2 (Paired box 2, TPM: 20.8) (Torban et al., 2000; Cai et al., 2005), and Egr1 (Early Growth Response 1, TPM: 49.7) (Cohen et al., 1994, 1996) are transcription factors reliably expressed in the mIMCD3 cells (Supplementary Table 1, “Transcription factors”). In addition to known expression of transcription factors, expression of transcription factors that could be involved in the kidney nephron development, such as Emx2, Pax8, Tfap2a, Hmga2, Hmgb2, and Hoxa11, were identified (Ribes et al., 2003; Schwab et al., 2003; Chambers et al., 2019).

As shown in the previous studies (Slaats et al., 2015; Liu et al., 2018), expression of two ion transport proteins, Pkd2 (Polycystic kidney disease 2, TPM: 20.6) and Clcn4 (H⁺/Cl^– exchange transporter 4, TPM: 9.5), were found in the current transcriptome (Table 2). We also identified the mechanosensitive cation channel Piezo1 (piezo-type mechanosensitive ion channel component 1) (Coste et al., 2010; Table 2). Expression of ion channels, Trpm4 (Transient Receptor Potential Cation Channel Subfamily M Member 4, TPM: 13.9), Trpm6 (Transient Receptor Potential Cation Channel Subfamily M Member 6, TPM: 3.3), and Trpm7 (Transient Receptor Potential Cation Channel Subfamily M Member 7, TPM: 23.7), as well as a Mg⁺ transporter, Magt1 (TPM: 20.1), suggests further utilization as an in vitro model for studying transepithelial Mg⁺ transport mechanisms in kidney epithelial cells (Groenestege et al., 2006; Tables 2, 3 and Supplementary Table 1, “Ion channels” and “Transporters”).

The mIMCD3 cell line has been widely used in several studies of primary cilia conformation and function. The current transcriptome profile provides information about ciliary GPCRs expressed in mIMCD3 cells. Among the ciliary GPCRs, as previously reported (Hilgendorf et al., 2016; Mykytyn and Askwith, 2017), we identified expression of Gpr161 (G protein-coupled receptor 161, TPM: 11.5), Ptger4 (prostaglandin E receptor 4, TPM: 11.7), and Smo (Smoothened homolog, TPM: 22.86) in mIMCD3 cells (Table 4). In addition to ciliary GPCRs, expression of two adenylyl cyclases, Adcy1 (TPM: 8.9) and Adcy6 (TPM: 31.5), was found in the transcriptome profile. These adenylyl cyclases could be considered in further studies for examining cAMP responses in mIMCD3 cells (Strait et al., 2010).

Doxycycline-Responsive Transcriptome Changes in mIMCD3 Cells

To characterize the DOX response in mIMCD3 cells, we compared the transcriptomes of DOX-treated mIMCD3 cells and DMSO-treated mIMCD3 cells (Figure 2 and Supplementary Table 2). We found DOX changed the gene expression profile, reflecting a change in abundance of 1,662 genes at 3 days (Figure 2A) and 2,858 genes at 6 days of treatment (Figure 2B). A total of 1,157 genes were consistently changed at both times. Downstream analysis identified GO Biological Processes (GOBPs) that are enriched in DOX-treated cells (Figure 2C). The NaviGO based GO term association analysis, which classifies similarity within GO term hierarchies, revealed three major clusters of GOBPs: (1) cell proliferation/differentiation; (2) signal transduction; and (3) immune responses (Figure 2D). In the cells treated with DOX for 3 days, differentially expressed genes were enriched for signal transduction pathways and cell proliferation/differentiation (Figures 2C,D). In particular, this involved significant changes of genes associated with ERK, cAMP, and Notch signaling pathways, known to mainly involve cell proliferation and development processes (Figure 2C). Differentially expressed genes from the 6 days-dataset were also especially enriched in biological processes associated with cell proliferation/differentiation but were also enriched in immune response pathways (Figures 2C,D). As shown in the heatmap (Figure 2E), expression of genes involved in cell proliferation/differentiation processes and immune responses were mostly decreased at both times. The results indicate that DOX attenuates gene expression associated with cell proliferation and immune responses. Full list of genes associated with these pathways are provided in Supplementary Table 3.

FIGURE 2

Examination of genes encoding downstream factors in signal transduction processes, including transcription factors and kinases, revealed alterations in three different signaling pathways (cAMP, ERK, Notch). Remarkable downregulation of several transcription factors and protein kinases was found in DOX-treated cells (Figure 3A). Especially, Hnf1b (Chung et al., 2017), Egr1 (Sukhatme et al., 1988), Tead1 (Zhao et al., 2008), and members of the Fox family (Kume et al., 2000; Aschauer et al., 2013), which have been proposed previously as key players in epithelial cell proliferation. Additionally, several protein kinases downregulated in response to DOX, namely Tgfbr2 (LeBleu et al., 2013), Kit (Gomes et al., 2018), and protein kinase A (Amsler et al., 1991), are highly related to regulation of epithelial cell function (Figure 3A). Furthermore, the effects of DOX to decrease the expression of genes involved in each step of the cell cycle and cell cycle progression suggest that DOX may have proclivity to inhibit cell proliferation (Figure 3B).

FIGURE 3

The transcriptome dataset showing the attenuated cell cycle progression also exhibited that several genes known as cell proliferation markers were consistently downregulated by DOX treatment at 6 days (Figure 4). The result confirmed that DOX treatment suppresses cell proliferation. In addition to cell proliferation markers, DOX treatment for 3 and 6 days induced significant changes of immune response-associated cellular pathways (Figure 2), largely reflecting DOX-responsive reduction in cytokine production at day 3 and comprehensive alteration of cellular inflammatory response at day 6. In particularly, we identified several chemokines that were significantly downregulated at 3 days after DOX treatment from the literature-based gene sets of immunologic mediators (Commins et al., 2010; Figure 4), corresponding the pathways involved in the repression of epithelial cell proliferation, which are associated with epithelial responses to cytokines (Stadnyk, 1994).

FIGURE 4

Discussion

The mIMCD3 cell line, an inner medullary collecting duct (IMCD) cell line derived from a simian virus 40 (SV40) T-antigen transgenic mice, has been widely used as an in vitro model to study renal epithelial physiology and cell biology (Rauchman et al., 1993; Yu et al., 2013; Flannery et al., 2015; Siroky et al., 2017). With inclusion of a DOX-inducible gene expression system, the cell line has become an invaluable model to study gene function in renal epithelial background. Here, we applied next-generation sequencing technology to characterize how the mIMCD3 transcriptome is affected by DOX with the goal of creating a database to guide interpretation of these studies.

Appropriate differentiation of epithelial cells and development of apical-basolateral polarity in vitro is commonly believed to require growth on a permeable surface. Comparison of the mIMCD3 transcriptome of cells grown on permeable polyester (PE) filter membrane surface, as defined here, with the published the transcriptome profile of mIMCD3 cells grown on plastic (Chan et al., 2018) lends credence to this idea. Although we found transcriptome profiles of mIMCD3 cells grown on filters exhibit a surprising degree of similarity with the transcriptome profile of cells grown on a solid surface (Reboredo et al., 2008), the profile of cells grown on filters is consistent with a more differentiated state. We found genes encoding key junctional stability components (Cldn4, Cldn7, Epcam) are induced when cells are grown on filters, whereas mesenchymal markers and ECM components were less abundant in filter-grown cells. Cldn4 (Claudin-4) and Cldn7 (Claudin-7) are among claudin family proteins that are specifically expressed in the loop of Henle (LOH) and distal nephron (Leiz and Schmidt-Ott, 2019), dictating the unique pericellular solute permeation profile and the electrical resistance of these nephron segments. Epcam (Epithelial Cell Adhesion Molecule) is best known as a marker of epithelial oncogenesis, but it is highly expressed in the LOH and collecting duct, where its natural role in cell adhesion and morphogenesis may determine the differentiated state of these nephron segments by interaction with Claudin-7 to negatively regulate epithelial migration by inhibiting ERK and actomyosin contractility (Barth et al., 2018). The increased abundances of ECM components may identify a gene expression signature of solid support grown cells, reflecting a futile cellular response to establish a permeable basement membrane, and thereby influence cell proliferation and differentiation (Matter et al., 2005; Ichikawa-Tomikawa et al., 2011; Bonnans et al., 2014; Diaz-Coranguez et al., 2019). Supporting this idea, the upregulated ECM are ubiquitous components of most basement membranes, including members of type IV collagen family proteins, Col4a1 and Col4a2, and the Secreted Protein Acidic and Rich in Cysteine (SPARC, aka Osteonectin).

We found the subtle growth substrate-dependent changes occur against an otherwise, remarkably stable mIMCD3 transcriptome profile. Constancy of the gene expression profile may explain why the phenotypic properties of the cell line have been so stable over time and numerous different studies. Comparison of the mIMCD3 transcriptome to the single cell RNA profiles of mouse cell kidney (Ransick et al., 2019) indicates that stability has come at the cost of lost nephron cell definition, however. The transcriptome profile of mIMCD3 does not uniquely resemble transcriptomes of any cell type in the kidney. Instead, it is compatible with a generic renal epithelial cell model, presumably reflecting a dedifferentiated state.

Tetracycline drugs, like Doxycycline (DOX), are widely used in eukaryotic cell gene expression systems, in part, because they have been assumed to have limited toxicity profiles. They inhibit protein synthesis in a broad spectrum of bacteria by binding to rRNAs of ribosomal complexes and double-strand RNAs (Chukwudi, 2016; Chukwudi and Good, 2016), but do not usually accumulate enough in mammalian cells to affect protein synthesis. We found DOX, at least at the higher concentrations commonly used in vitro cell models to drive tet-o gene expression, has a larger than anticipated toxic profile. Differentially expressed genes in response to DOX were especially enriched in the suppression of cell proliferation processes, including a decrease in cell cycle progression genes (Figures 3, 4). DOX also has been reported to suppress cell proliferation in epithelial-derived cancers (Fife et al., 1998), providing reason to speculate that tetracycline derivatives might be developed as anti-cancer therapeutics (Kroon et al., 1984; van den Bogert et al., 1985, 1986; Mortison et al., 2018). In human epithelial cell lines, DOX reduces cell proliferation, associated with alterations in mitochondrial function (Sourdeval et al., 2006; Chang et al., 2010; Ahler et al., 2013).

Although future studies will be required to unravel the mechanisms by which DOX suppresses cell proliferation, the change in the mIMCD3 transcriptome profile provides important insights. Decreased expression of cytokines and chemokines (Colombo et al., 2018; Koga et al., 2019), together with the attenuation of intracellular signal transduction pathways (cAMP, ERK, and Notch) that are known to regulate cell proliferation (Yamaguchi et al., 2000; Stork and Schmitt, 2002; Kim and Shivdasani, 2011) provide reason to suggest that DOX suppresses the production of autocrine cell growth factors that are required for mIMCD3 proliferation. Additionally, DOX activated the interferon-beta pathway, involving guanylate-binding proteins (GBPs) and interferon activated genes (Lubeseder-Martellato et al., 2002; Klamp et al., 2003; MacMicking, 2004), which is known to inhibit cell proliferation and differentiation (Hertzog et al., 1994; Garrison et al., 1996; Vitale et al., 2006). It will be interesting to learn if either of these pathways can be manipulated to reduce the anti-proliferation effects of DOX. It will be important to corroborate that changes in protein abundance track changes in transcript levels.

The effects of DOX to rapidly induce an anti-proliferative transcriptomic signature in mIMCD3 cells should not be taken as a reason to discontinue use of the otherwise highly effective DOX-based gene induction system. Because it seems likely that all induction agents will have some off-target effects, it more practical to develop strategies to work around the off-target responses. Washing out DOX after gene induction is an obvious remedy. However, future studies will be required to determine the longevity of the DOX response, after withdrawal.

Growth and differentiation of cells in vitro can be profoundly influenced by the cell culture medium, nutrients, growth factors, hormones, and other supplements. They were not varied in the present study, and thus should not contribute to differences in the transcriptome profiles between vehicle and DOX treated cells. However, the basal transcriptome profile is likely to be influenced by the supplements in the media. Moreover, potential synergistic actions between DOX and any one of the many supplements is unknowable. Consequently, caution should be exercised in generalizing our results with other studies that use different culture media.

In summary, we have assessed the transcriptomic response to DOX in mIMCD3 at genome-wide resolution. The data base provides an atlas to guide future interpretation of studies using the DOX-inducible gene expression system in renal epithelial cells.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

• 1

  AboudehenK.NoureddineL.Cobo-StarkP.AvdulovS.FarahaniS.GearhartM. D.et al (2017). Hepatocyte nuclear factor-1beta regulates urinary concentration and response to hypertonicity.J. Am. Soc. Nephrol.282887–2900. 10.1681/ASN.2016101095

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 2

  AhlerE.SullivanW. J.CassA.BraasD.YorkA. G.BensingerS. J.et al (2013). Doxycycline alters metabolism and proliferation of human cell lines.PLoS One8:e64561. 10.1371/journal.pone.0064561

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 3

  AmslerK.GhataniS.HemmingsB. A. (1991). cAMP-dependent protein kinase regulates renal epithelial cell properties.Am. J. Physiol.260(6 Pt 1)C1290–C1299. 10.1152/ajpcell.1991.260.6.C1290

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 4

  AschauerL.GruberL. N.PfallerW.LimoncielA.AthersuchT. J.CavillR.et al (2013). Delineation of the key aspects in the regulation of epithelial monolayer formation.Mol. Cell. Biol.332535–2550. 10.1128/MCB.01435-12

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 5

  BarthA. I. M.KimH.Riedel-KruseI. H. (2018). Regulation of epithelial migration by epithelial cell adhesion molecule requires its Claudin-7 interaction domain.PLoS One13:e0204957. 10.1371/journal.pone.0204957

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 6

  BonnansC.ChouJ.WerbZ. (2014). Remodelling the extracellular matrix in development and disease.Nat. Rev. Mol. Cell Biol.15786–801. 10.1038/nrm3904

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 7

  CaiQ.DmitrievaN. I.FerrarisJ. D.BrooksH. L.van BalkomB. W.BurgM. (2005). Pax2 expression occurs in renal medullary epithelial cells in vivo and in cell culture, is osmoregulated, and promotes osmotic tolerance.Proc. Natl. Acad. Sci. U.S.A.102503–508. 10.1073/pnas.0408840102

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 8

  ChambersB. E.GerlachG. F.ClarkE. G.ChenK. H.LevesqueA. E.LeshchinerI.et al (2019). Tfap2a is a novel gatekeeper of nephron differentiation during kidney development.Development146:dev172387. 10.1242/dev.172387

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 9

  ChanS. C.ZhangY.ShaoA.AvdulovS.HerreraJ.AboudehenK.et al (2018). Mechanism of fibrosis in HNF1B-related autosomal dominant tubulointerstitial kidney disease.J. Am. Soc. Nephrol.292493–2509. 10.1681/ASN.2018040437

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 10

  ChangW. Y.ClementsD.JohnsonS. R. (2010). Effect of doxycycline on proliferation, MMP production, and adhesion in LAM-related cells.Am. J. Physiol. Lung Cell. Mol. Physiol.299L393–L400. 10.1152/ajplung.00437.2009

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 11

  ChenY.LunA. T.SmythG. K. (2016). From reads to genes to pathways: differential expression analysis of RNA-Seq experiments using Rsubread and the edgeR quasi-likelihood pipeline.F1000Researh5:1438. 10.12688/f1000research.8987.2

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 12

  ChukwudiC. U. (2016). rRNA binding sites and the molecular mechanism of action of the tetracyclines.Antimicrob. Agents Chemother.604433–4441. 10.1128/AAC.00594-16

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 13

  ChukwudiC. U.GoodL. (2016). Interaction of the tetracyclines with double-stranded RNAs of random base sequence: new perspectives on the target and mechanism of action.J. Antibiot.69622–630. 10.1038/ja.2015.145

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 14

  ChungE.DeaconP.ParkJ. S. (2017). Notch is required for the formation of all nephron segments and primes nephron progenitors for differentiation.Development1444530–4539. 10.1242/dev.156661

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 15

  CohenD. M.ChinW. W.GullansS. R. (1994). Hyperosmotic urea increases transcription and synthesis of Egr-1 in murine inner medullary collecting duct (mIMCD3) cells.J. Biol. Chem.26925865–25870.

  • Google Scholar

• 16

  CohenD. M.GullansS. R.ChinW. W. (1996). Urea signaling in cultured murine inner medullary collecting duct (mIMCD3) cells involves protein kinase C, inositol 1,4,5-trisphosphate (IP3), and a putative receptor tyrosine kinase.J. Clin. Invest.971884–1889. 10.1172/JCI118619

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 17

  ColomboM.MirandolaL.Chiriva-InternatiM.BasileA.LocatiM.LesmaE.et al (2018). Cancer cells exploit Notch signaling to redefine a supportive cytokine milieu.Front. Immunol.9:1823. 10.3389/fimmu.2018.01823

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 18

  ComminsS. P.BorishL.SteinkeJ. W. (2010). Immunologic messenger molecules: cytokines, interferons, and chemokines.J. Allergy Clin. Immunol.125(2 Suppl. 2)S53–S72. 10.1016/j.jaci.2009.07.008

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 19

  CosteB.MathurJ.SchmidtM.EarleyT. J.RanadeS.PetrusM. J.et al (2010). Piezo1 and Piezo2 are essential components of distinct mechanically activated cation channels.Science33055–60. 10.1126/science.1193270

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 20

  DasA. T.TenenbaumL.BerkhoutB. (2016). Tet-on systems for doxycycline-inducible gene expression.Curr. Gene Ther.16156–167. 10.2174/1566523216666160524144041

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 21

  Diaz-CoranguezM.LiuX.AntonettiD. A. (2019). Tight junctions in cell proliferation.Int. J. Mol. Sci.20:5972. 10.3390/ijms20235972

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 22

  FifeR. S.SledgeG. W.Jr.RothB. J.ProctorC. (1998). Effects of doxycycline on human prostate cancer cells in vitro.Cancer Lett.12737–41. 10.1016/s0304-3835(98)00003-2

  • CrossRef
  • Google Scholar

• 23

  FlanneryR. J.KleeneN. K.KleeneS. J. (2015). A TRPM4-dependent current in murine renal primary cilia.Am. J. Physiol. Renal Physiol.309F697–F707. 10.1152/ajprenal.00294.2015

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 24

  GarrisonJ. I.BerensM. E.ShapiroJ. R.TreasurywalaS.Floyd-SmithG. (1996). Interferon-beta inhibits proliferation and progression through S phase of the cell cycle in five glioma cell lines.J. Neurooncol.30213–223. 10.1007/BF00177272

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 25

  GomesS. A.HareJ. M.RangelE. B. (2018). Kidney-derived c-Kit(+) cells possess regenerative potential.Stem Cells Transl. Med.7317–324. 10.1002/sctm.17-0232

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 26

  GroenestegeW. M.HoenderopJ. G.van den HeuvelL.KnoersN.BindelsR. J. (2006). The epithelial Mg2+ channel transient receptor potential melastatin 6 is regulated by dietary Mg2+ content and estrogens.J. Am. Soc. Nephrol.171035–1043. 10.1681/ASN.2005070700

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 27

  HertzogP. J.HwangS. Y.KolaI. (1994). Role of interferons in the regulation of cell proliferation, differentiation, and development.Mol. Reprod. Dev.39226–232. 10.1002/mrd.1080390216

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 28

  HilgendorfK. I.JohnsonC. T.JacksonP. K. (2016). The primary cilium as a cellular receiver: organizing ciliary GPCR signaling.Curr. Opin. Cell Biol.3984–92. 10.1016/j.ceb.2016.02.008

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 29

  Ichikawa-TomikawaN.SugimotoK.SatohisaS.NishiuraK.ChibaH. (2011). Possible involvement of tight junctions, extracellular matrix and nuclear receptors in epithelial differentiation.J. Biomed. Biotechnol.2011:253048. 10.1155/2011/253048

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 30

  JedroszkaD.OrzechowskaM.HamouzR.GorniakK.BednarekA. K. (2017). Markers of epithelial-to-mesenchymal transition reflect tumor biology according to patient age and Gleason score in prostate cancer.PLoS One12:e0188842. 10.1371/journal.pone.0188842

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 31

  KangK.HuangL.LiQ.LiaoX.DangQ.YangY.et al (2019). An improved Tet-on system in microRNA overexpression and CRISPR/Cas9-mediated gene editing.J. Anim. Sci. Biotechnol.10:43. 10.1186/s40104-019-0354-5

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 32

  KimT. H.ShivdasaniR. A. (2011). Notch signaling in stomach epithelial stem cell homeostasis.J. Exp. Med.208677–688. 10.1084/jem.20101737

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 33

  KistnerA.GossenM.ZimmermannF.JerecicJ.UllmerC.LubbertH.et al (1996). Doxycycline-mediated quantitative and tissue-specific control of gene expression in transgenic mice.Proc. Natl. Acad. Sci. U.S.A.9310933–10938. 10.1073/pnas.93.20.10933

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 34

  KlampT.BoehmU.SchenkD.PfefferK.HowardJ. C. (2003). A giant GTPase, very large inducible GTPase-1, is inducible by IFNs.J. Immunol.1711255–1265. 10.4049/jimmunol.171.3.1255

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35

  KogaY.TsurumakiH.Aoki-SaitoH.SatoM.YatomiM.TakeharaK.et al (2019). Roles of cyclic AMP response element binding activation in the ERK1/2 and p38 MAPK signalling pathway in central nervous system, cardiovascular system, osteoclast differentiation and mucin and cytokine production.Int. J. Mol. Sci.20:1346. 10.3390/ijms20061346

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 36

  KroonA. M.DontjeB. H.HoltropM.Van den BogertC. (1984). The mitochondrial genetic system as a target for chemotherapy: tetracyclines as cytostatics.Cancer Lett.2533–40. 10.1016/s0304-3835(84)80023-3

  • CrossRef
  • Google Scholar

• 37

  KumeT.DengK.HoganB. L. (2000). Murine forkhead/winged helix genes Foxc1 (Mf1) and Foxc2 (Mfh1) are required for the early organogenesis of the kidney and urinary tract.Development1271387–1395.

  • Google Scholar

• 38

  LarrayozI. M.de LuisA.RuaO.VelillaS.CabelloJ.MartinezA. (2012). Molecular effects of doxycycline treatment on pterygium as revealed by massive transcriptome sequencing.PLoS One7:e39359. 10.1371/journal.pone.0039359

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 39

  LashhabR.RumleyA. C.ArutyunovD.RizviM.YouC.DimkeH.et al (2019). The kidney anion exchanger 1 affects tight junction properties via claudin-4.Sci. Rep.9:3099. 10.1038/s41598-019-39430-9

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 40

  LeBleuV. S.TaduriG.O’ConnellJ.TengY.CookeV. G.WodaC.et al (2013). Origin and function of myofibroblasts in kidney fibrosis.Nat. Med.191047–1053. 10.1038/nm.3218

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 41

  LeizJ.Schmidt-OttK. M. (2019). Claudins in the renal collecting duct.Int. J. Mol. Sci.21:221. 10.3390/ijms21010221

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 42

  LiuX.VienT.DuanJ.SheuS. H.DeCaenP. G.ClaphamD. E. (2018). Polycystin-2 is an essential ion channel subunit in the primary cilium of the renal collecting duct epithelium.Elife7:e33183. 10.7554/eLife.33183

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 43

  LuH.JiangW.YangH.QinZ.GuoS. E.HuM.et al (2017). Doxycycline affects gene expression profiles in aortic tissues in a rat model of vascular calcification.Microvasc. Res.11412–18. 10.1016/j.mvr.2017.04.007

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 44

  Lubeseder-MartellatoC.GuenziE.JorgA.TopoltK.NaschbergerE.KremmerE.et al (2002). Guanylate-binding protein-1 expression is selectively induced by inflammatory cytokines and is an activation marker of endothelial cells during inflammatory diseases.Am. J. Pathol.1611749–1759. 10.1016/S0002-9440(10)64452-5

  • CrossRef
  • Google Scholar

• 45

  MacMickingJ. D. (2004). IFN-inducible GTPases and immunity to intracellular pathogens.Trends Immunol.25601–609. 10.1016/j.it.2004.08.010

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 46

  MatterK.AijazS.TsaparaA.BaldaM. S. (2005). Mammalian tight junctions in the regulation of epithelial differentiation and proliferation.Curr. Opin. Cell Biol.17453–458. 10.1016/j.ceb.2005.08.003

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 47

  MortisonJ. D.SchenoneM.MyersJ. A.ZhangZ.ChenL.CiarloC.et al (2018). Tetracyclines modify translation by targeting key human rRNA substructures.Cell Chem. Biol.251506–1518.e13. 10.1016/j.chembiol.2018.09.010

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 48

  MykytynK.AskwithC. (2017). G-protein-coupled receptor signaling in cilia.Cold Spring Harb. Perspect. Biol.9:a028183. 10.1101/cshperspect.a028183

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 49

  PatroR.DuggalG.LoveM. I.IrizarryR. A.KingsfordC. (2017). Salmon provides fast and bias-aware quantification of transcript expression.Nat. Methods14417–419. 10.1038/nmeth.4197

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 50

  RansickA.LindstromN. O.LiuJ.ZhuQ.GuoJ. J.AlvaradoG. F.et al (2019). Single-cell profiling reveals sex, lineage, and regional diversity in the mouse kidney.Dev. Cell51399–413.e7. 10.1016/j.devcel.2019.10.005

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 51

  RauchmanM. I.NigamS. K.DelpireE.GullansS. R. (1993). An osmotically tolerant inner medullary collecting duct cell line from an SV40 transgenic mouse.Am. J. Physiol.265(3 Pt 2)F416–F424. 10.1152/ajprenal.1993.265.3.F416

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 52

  ReboredoM.KramerM. G.SmerdouC.PrietoJ.De Las RivasJ. (2008). Transcriptomic effects of Tet-on and mifepristone-inducible systems in mouse liver.Hum. Gene Ther.191233–1247. 10.1089/hum.2008.057

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 53

  RibesD.FischerE.CalmontA.RossertJ. (2003). Transcriptional control of epithelial differentiation during kidney development.J. Am. Soc. Nephrol.14(Suppl. 1)S9–S15. 10.1097/01.asn.0000067647.05964.9f

  • CrossRef
  • Google Scholar

• 54

  RobinsonM. D.McCarthyD. J.SmythG. K. (2010). edgeR: a bioconductor package for differential expression analysis of digital gene expression data.Bioinformatics26139–140. 10.1093/bioinformatics/btp616

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 55

  SchlimpertM.LagiesS.BudnykV.MullerB.WalzG.KammererB. (2018). Metabolic phenotyping of Anks3 depletion in mIMCD-3 cells – a putative nephronophthisis candidate.Sci. Rep.8:9022. 10.1038/s41598-018-27389-y

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 56

  SchwabK.PattersonL. T.AronowB. J.LuckasR.LiangH. C.PotterS. S. (2003). A catalogue of gene expression in the developing kidney.Kidney Int.641588–1604. 10.1046/j.1523-1755.2003.00276.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 57

  SirokyB. J.KleeneN. K.KleeneS. J.VarnellC. D.Jr.ComerR. G.LiuJ.et al (2017). Primary cilia regulate the osmotic stress response of renal epithelial cells through TRPM3.Am. J. Physiol. Renal Physiol.312F791–F805. 10.1152/ajprenal.00465.2015

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 58

  SlaatsG. G.WhewayG.FolettoV.SzymanskaK.van BalkomB. W.LogisterI.et al (2015). Screen-based identification and validation of four new ion channels as regulators of renal ciliogenesis.J. Cell Sci.1284550–4559. 10.1242/jcs.176065

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 59

  SourdevalM.LemaireC.BrennerC.Boisvieux-UlrichE.MaranoF. (2006). Mechanisms of doxycycline-induced cytotoxicity on human bronchial epithelial cells.Front. Biosci.11:3036–3048. 10.2741/2031

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 60

  StadnykA. W. (1994). Cytokine production by epithelial cells.FASEB J.81041–1047. 10.1096/fasebj.8.13.7926369

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 61

  StorkP. J.SchmittJ. M. (2002). Crosstalk between cAMP and MAP kinase signaling in the regulation of cell proliferation.Trends Cell Biol.12258–266. 10.1016/s0962-8924(02)02294-8

  • CrossRef
  • Google Scholar

• 62

  StraitK. A.StricklettP. K.ChapmanM.KohanD. E. (2010). Characterization of vasopressin-responsive collecting duct adenylyl cyclases in the mouse.Am. J. Physiol. Renal Physiol.298F859–F867. 10.1152/ajprenal.00109.2009

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 63

  SukhatmeV. P.CaoX. M.ChangL. C.Tsai-MorrisC. H.StamenkovichD.FerreiraP. C.et al (1988). A zinc finger-encoding gene coregulated with c-fos during growth and differentiation, and after cellular depolarization.Cell5337–43. 10.1016/0092-8674(88)90485-0

  • CrossRef
  • Google Scholar

• 64

  TorbanE.EcclesM. R.FavorJ.GoodyerP. R. (2000). PAX2 suppresses apoptosis in renal collecting duct cells.Am. J. Pathol.157833–842. 10.1016/S0002-9440(10)64597-X

  • CrossRef
  • Google Scholar

• 65

  ValkovaN.KultzD. (2006). Constitutive and inducible stress proteins dominate the proteome of the murine inner medullary collecting duct-3 (mIMCD3) cell line.Biochim. Biophys. Acta17641007–1020. 10.1016/j.bbapap.2006.03.007

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 66

  van den BogertC.DontjeB. H.HoltropM.MelisT. E.RomijnJ. C.van DongenJ. W.et al (1986). Arrest of the proliferation of renal and prostate carcinomas of human origin by inhibition of mitochondrial protein synthesis.Cancer Res.463283–3289.

  • Google Scholar

• 67

  van den BogertC.DontjeB. H.KroonA. M. (1985). The antitumour effect of doxycycline on a T-cell leukaemia in the rat.Leuk. Res.9617–623. 10.1016/0145-2126(85)90142-0

  • CrossRef
  • Google Scholar

• 68

  VitaleG.de HerderW. W.van KoetsveldP. M.WaaijersM.SchoordijkW.CrozeE.et al (2006). IFN-beta is a highly potent inhibitor of gastroenteropancreatic neuroendocrine tumor cell growth in vitro.Cancer Res.66554–562. 10.1158/0008-5472.CAN-05-3043

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 69

  WeiQ.KhanI. K.DingZ.YerneniS.KiharaD. (2017). NaviGO: interactive tool for visualization and functional similarity and coherence analysis with gene ontology.BMC Bioinformatics18:177. 10.1186/s12859-017-1600-5

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 70

  WishartJ. A.HayesA.WardleworthL.ZhangN.OliverS. G. (2005). Doxycycline, the drug used to control the tet-regulatable promoter system, has no effect on global gene expression in Saccharomyces cerevisiae.Yeast22565–569. 10.1002/yea.1225

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 71

  YamaguchiT.PellingJ. C.RamaswamyN. T.EpplerJ. W.WallaceD. P.NagaoS.et al (2000). cAMP stimulates the in vitro proliferation of renal cyst epithelial cells by activating the extracellular signal-regulated kinase pathway.Kidney Int.571460–1471. 10.1046/j.1523-1755.2000.00991.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 72

  YuZ.KongQ.KoneB. C. (2013). Aldosterone reprograms promoter methylation to regulate alphaENaC transcription in the collecting duct.Am. J. Physiol. Renal Physiol.305F1006–F1013. 10.1152/ajprenal.00407.2013

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 73

  ZhaoB.YeX.YuJ.LiL.LiW.LiS.et al (2008). TEAD mediates YAP-dependent gene induction and growth control.Genes Dev.221962–1971. 10.1101/gad.1664408

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 74

  ZhouY.ZhouB.PacheL.ChangM.KhodabakhshiA. H.TanaseichukO.et al (2019). Metascape provides a biologist-oriented resource for the analysis of systems-level datasets.Nat. Commun.10:1523. 10.1038/s41467-019-09234-6

  • Pubmed Abstract
  • CrossRef
  • Google Scholar