Hyun Jun Jung
Richard Coleman1
Owen M. Woodward- 1Division of Nephrology, Department of Medicine, Johns Hopkins University School of Medicine, Baltimore, MD, United States
- 2Department of Physiology, University of Maryland School of Medicine, Baltimore, MD, United States
- 3Department of Physiology, Johns Hopkins University School of Medicine, Baltimore, MD, United States
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 QubitTM 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 platform1 (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 (GSE1715732).
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. Transcriptomic characteristics of renal epithelial cell line mIMCD3 in different culture conditions. (A) Transcriptomic correlation was assessed in mIMCD3 cells grown on the permeable filter supports for 3 and 6 days. (B,C) Transcriptome profiles of mIMCD3 grown on the permeable filter support and the solid support were compared. The transcriptome dataset of mIMCD3 cells grown on the solid support was obtained from GSE97770 (Chan et al., 2018). (D–G) Box plots of gene expression associated with epithelial cell, mesenchymal cell, tight junction formation, and extracellular matrix (ECM) in mIMCD3 cells grown on the solid support (Chan et al., 2018) and the permeable filter support.
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).
Table 1. Top 30 transcription factors (see the full transcription factor list in Supplementary Table 1).
Table 2. Top 30 ion channels (see the full ion channel list in Supplementary Table 1).
Table 3. Top 30 transporters (see the full transporter list in Supplementary Table 1).
Table 4. Top 30 GPCRs (see the full GPCR list in Supplementary Table 1).
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. DOX-responsive transcriptomic changes in mIMCD3 cells. (A,B) Identification of transcriptomic changes in mIMCD3 cells treated with DOX for 3 and 6 days using RNA-Seq-based differential expression analysis. Significantly changed genes (FDR < 0.05 over 20% expression change threshold) were labeled as green in volcano plots and used further downstream analysis. (C) Gene Ontology (GO) biological processes significantly enriched by DOX-responsive genes. Cellular processes (GO Biological Process) were identified from significantly changed genes by DOX treatment for 3 and 6 days. To avoid sample-source bias, whole transcriptome was used as a background gene set. Significant enrichment was considered as q < 0.05. (D) Bubble chart view of the GO term association using NaviGO. In the plot, the X-axis is the Resnik semantic similarity score (RSS) and the Y-axis chosen is Lin’s similarity score (LSS). GO terms are colored according to time point when DE genes are found, day 3, pink; day 6: green; and day 3 and 6: blue. The number of GO terms are listed in (C). (E) Heatmap of gene sets associated with “cell proliferation/differentiation” and “immune responses.” Full list of genes associated with GO terms are provided in Supplementary Table 4.
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. DOX-responsive genes associated with cell proliferation and differentiation. (A) DOX-responsive expressional change of transcription factors (TFs) and protein kinases in three signaling pathways (ERK, cAMP, and Notch). Full list of DOX-responsive genes involved in these signaling pathways was provided in Supplementary Table 3. (B) Expressional change of gene associated with cell cycle progression. Gene sets of each stage of cell cycle progression were obtained from Gene Ontology database (GO:0022402: cell cycle process).
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. DOX-responsive cytokines and cell proliferation markers. Significantly changed genes encoding cytokines including chemokines were identified and listed in the plot. In addition, expressional changes of genes previously known as cell proliferation markers were listed in the plot.
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.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
HJ, OW, and PW conceived and designed the studies. HJ, RC, OW, and PW performed the experiments, analyzed the data, and wrote the manuscript. All authors discussed and reviewed the manuscript.
Funding
This work was supported by grants from the National Institute of Diabetes and Digestive and Kidney Diseases (Grants DK054231, 1U54DK126114, and DK093501) and the Leducq Foundation.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
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.
Acknowledgments
Next-generation DNA sequencing for RNA-Seq transcriptomic analysis was performed by the Genomic Resource Center, Institute for Genome Sciences, University of Maryland School of Medicine (Director: Dr. Lisa DeShong Sadzewicz).
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2021.771691/full#supplementary-material
Footnotes
References
Aboudehen, K., Noureddine, L., Cobo-Stark, P., Avdulov, S., Farahani, S., Gearhart, M. D., et al. (2017). Hepatocyte nuclear factor-1beta regulates urinary concentration and response to hypertonicity. J. Am. Soc. Nephrol. 28, 2887–2900. doi: 10.1681/ASN.2016101095
Ahler, E., Sullivan, W. J., Cass, A., Braas, D., York, A. G., Bensinger, S. J., et al. (2013). Doxycycline alters metabolism and proliferation of human cell lines. PLoS One 8:e64561. doi: 10.1371/journal.pone.0064561
Amsler, K., Ghatani, S., and Hemmings, B. A. (1991). cAMP-dependent protein kinase regulates renal epithelial cell properties. Am. J. Physiol. 260(6 Pt 1), C1290–C1299. doi: 10.1152/ajpcell.1991.260.6.C1290
Aschauer, L., Gruber, L. N., Pfaller, W., Limonciel, A., Athersuch, T. J., Cavill, R., et al. (2013). Delineation of the key aspects in the regulation of epithelial monolayer formation. Mol. Cell. Biol. 33, 2535–2550. doi: 10.1128/MCB.01435-12
Barth, A. I. M., Kim, H., and Riedel-Kruse, I. H. (2018). Regulation of epithelial migration by epithelial cell adhesion molecule requires its Claudin-7 interaction domain. PLoS One 13:e0204957. doi: 10.1371/journal.pone.0204957
Bonnans, C., Chou, J., and Werb, Z. (2014). Remodelling the extracellular matrix in development and disease. Nat. Rev. Mol. Cell Biol. 15, 786–801. doi: 10.1038/nrm3904
Cai, Q., Dmitrieva, N. I., Ferraris, J. D., Brooks, H. L., van Balkom, B. W., and Burg, M. (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. 102, 503–508. doi: 10.1073/pnas.0408840102
Chambers, B. E., Gerlach, G. F., Clark, E. G., Chen, K. H., Levesque, A. E., Leshchiner, I., et al. (2019). Tfap2a is a novel gatekeeper of nephron differentiation during kidney development. Development 146:dev172387. doi: 10.1242/dev.172387
Chan, S. C., Zhang, Y., Shao, A., Avdulov, S., Herrera, J., Aboudehen, K., et al. (2018). Mechanism of fibrosis in HNF1B-related autosomal dominant tubulointerstitial kidney disease. J. Am. Soc. Nephrol. 29, 2493–2509. doi: 10.1681/ASN.2018040437
Chang, W. Y., Clements, D., and Johnson, S. R. (2010). Effect of doxycycline on proliferation, MMP production, and adhesion in LAM-related cells. Am. J. Physiol. Lung Cell. Mol. Physiol. 299, L393–L400. doi: 10.1152/ajplung.00437.2009
Chen, Y., Lun, A. T., and Smyth, G. K. (2016). From reads to genes to pathways: differential expression analysis of RNA-Seq experiments using Rsubread and the edgeR quasi-likelihood pipeline. F1000Researh 5:1438. doi: 10.12688/f1000research.8987.2
Chukwudi, C. U. (2016). rRNA binding sites and the molecular mechanism of action of the tetracyclines. Antimicrob. Agents Chemother. 60, 4433–4441. doi: 10.1128/AAC.00594-16
Chukwudi, C. U., and Good, L. (2016). Interaction of the tetracyclines with double-stranded RNAs of random base sequence: new perspectives on the target and mechanism of action. J. Antibiot. 69, 622–630. doi: 10.1038/ja.2015.145
Chung, E., Deacon, P., and Park, J. S. (2017). Notch is required for the formation of all nephron segments and primes nephron progenitors for differentiation. Development 144, 4530–4539. doi: 10.1242/dev.156661
Cohen, D. M., Chin, W. W., and Gullans, S. R. (1994). Hyperosmotic urea increases transcription and synthesis of Egr-1 in murine inner medullary collecting duct (mIMCD3) cells. J. Biol. Chem. 269, 25865–25870.
Cohen, D. M., Gullans, S. R., and Chin, W. 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. 97, 1884–1889. doi: 10.1172/JCI118619
Colombo, M., Mirandola, L., Chiriva-Internati, M., Basile, A., Locati, M., Lesma, E., et al. (2018). Cancer cells exploit Notch signaling to redefine a supportive cytokine milieu. Front. Immunol. 9:1823. doi: 10.3389/fimmu.2018.01823
Commins, S. P., Borish, L., and Steinke, J. W. (2010). Immunologic messenger molecules: cytokines, interferons, and chemokines. J. Allergy Clin. Immunol. 125(2 Suppl. 2), S53–S72. doi: 10.1016/j.jaci.2009.07.008
Coste, B., Mathur, J., Schmidt, M., Earley, T. J., Ranade, S., Petrus, M. J., et al. (2010). Piezo1 and Piezo2 are essential components of distinct mechanically activated cation channels. Science 330, 55–60. doi: 10.1126/science.1193270
Das, A. T., Tenenbaum, L., and Berkhout, B. (2016). Tet-on systems for doxycycline-inducible gene expression. Curr. Gene Ther. 16, 156–167. doi: 10.2174/1566523216666160524144041
Diaz-Coranguez, M., Liu, X., and Antonetti, D. A. (2019). Tight junctions in cell proliferation. Int. J. Mol. Sci. 20:5972. doi: 10.3390/ijms20235972
Fife, R. S., Sledge, G. W. Jr., Roth, B. J., and Proctor, C. (1998). Effects of doxycycline on human prostate cancer cells in vitro. Cancer Lett. 127, 37–41. doi: 10.1016/s0304-3835(98)00003-2
Flannery, R. J., Kleene, N. K., and Kleene, S. J. (2015). A TRPM4-dependent current in murine renal primary cilia. Am. J. Physiol. Renal Physiol. 309, F697–F707. doi: 10.1152/ajprenal.00294.2015
Garrison, J. I., Berens, M. E., Shapiro, J. R., Treasurywala, S., and Floyd-Smith, G. (1996). Interferon-beta inhibits proliferation and progression through S phase of the cell cycle in five glioma cell lines. J. Neurooncol. 30, 213–223. doi: 10.1007/BF00177272
Gomes, S. A., Hare, J. M., and Rangel, E. B. (2018). Kidney-derived c-Kit(+) cells possess regenerative potential. Stem Cells Transl. Med. 7, 317–324. doi: 10.1002/sctm.17-0232
Groenestege, W. M., Hoenderop, J. G., van den Heuvel, L., Knoers, N., and Bindels, R. J. (2006). The epithelial Mg2+ channel transient receptor potential melastatin 6 is regulated by dietary Mg2+ content and estrogens. J. Am. Soc. Nephrol. 17, 1035–1043. doi: 10.1681/ASN.2005070700
Hertzog, P. J., Hwang, S. Y., and Kola, I. (1994). Role of interferons in the regulation of cell proliferation, differentiation, and development. Mol. Reprod. Dev. 39, 226–232. doi: 10.1002/mrd.1080390216
Hilgendorf, K. I., Johnson, C. T., and Jackson, P. K. (2016). The primary cilium as a cellular receiver: organizing ciliary GPCR signaling. Curr. Opin. Cell Biol. 39, 84–92. doi: 10.1016/j.ceb.2016.02.008
Ichikawa-Tomikawa, N., Sugimoto, K., Satohisa, S., Nishiura, K., and Chiba, H. (2011). Possible involvement of tight junctions, extracellular matrix and nuclear receptors in epithelial differentiation. J. Biomed. Biotechnol. 2011:253048. doi: 10.1155/2011/253048
Jedroszka, D., Orzechowska, M., Hamouz, R., Gorniak, K., and Bednarek, A. K. (2017). Markers of epithelial-to-mesenchymal transition reflect tumor biology according to patient age and Gleason score in prostate cancer. PLoS One 12:e0188842. doi: 10.1371/journal.pone.0188842
Kang, K., Huang, L., Li, Q., Liao, X., Dang, Q., Yang, Y., et al. (2019). An improved Tet-on system in microRNA overexpression and CRISPR/Cas9-mediated gene editing. J. Anim. Sci. Biotechnol. 10:43. doi: 10.1186/s40104-019-0354-5
Kim, T. H., and Shivdasani, R. A. (2011). Notch signaling in stomach epithelial stem cell homeostasis. J. Exp. Med. 208, 677–688. doi: 10.1084/jem.20101737
Kistner, A., Gossen, M., Zimmermann, F., Jerecic, J., Ullmer, C., Lubbert, H., et al. (1996). Doxycycline-mediated quantitative and tissue-specific control of gene expression in transgenic mice. Proc. Natl. Acad. Sci. U.S.A. 93, 10933–10938. doi: 10.1073/pnas.93.20.10933
Klamp, T., Boehm, U., Schenk, D., Pfeffer, K., and Howard, J. C. (2003). A giant GTPase, very large inducible GTPase-1, is inducible by IFNs. J. Immunol. 171, 1255–1265. doi: 10.4049/jimmunol.171.3.1255
Koga, Y., Tsurumaki, H., Aoki-Saito, H., Sato, M., Yatomi, M., Takehara, K., 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. doi: 10.3390/ijms20061346
Kroon, A. M., Dontje, B. H., Holtrop, M., and Van den Bogert, C. (1984). The mitochondrial genetic system as a target for chemotherapy: tetracyclines as cytostatics. Cancer Lett. 25, 33–40. doi: 10.1016/s0304-3835(84)80023-3
Kume, T., Deng, K., and Hogan, B. L. (2000). Murine forkhead/winged helix genes Foxc1 (Mf1) and Foxc2 (Mfh1) are required for the early organogenesis of the kidney and urinary tract. Development 127, 1387–1395.
Larrayoz, I. M., de Luis, A., Rua, O., Velilla, S., Cabello, J., and Martinez, A. (2012). Molecular effects of doxycycline treatment on pterygium as revealed by massive transcriptome sequencing. PLoS One 7:e39359. doi: 10.1371/journal.pone.0039359
Lashhab, R., Rumley, A. C., Arutyunov, D., Rizvi, M., You, C., Dimke, H., et al. (2019). The kidney anion exchanger 1 affects tight junction properties via claudin-4. Sci. Rep. 9:3099. doi: 10.1038/s41598-019-39430-9
LeBleu, V. S., Taduri, G., O’Connell, J., Teng, Y., Cooke, V. G., Woda, C., et al. (2013). Origin and function of myofibroblasts in kidney fibrosis. Nat. Med. 19, 1047–1053. doi: 10.1038/nm.3218
Leiz, J., and Schmidt-Ott, K. M. (2019). Claudins in the renal collecting duct. Int. J. Mol. Sci. 21:221. doi: 10.3390/ijms21010221
Liu, X., Vien, T., Duan, J., Sheu, S. H., DeCaen, P. G., and Clapham, D. E. (2018). Polycystin-2 is an essential ion channel subunit in the primary cilium of the renal collecting duct epithelium. Elife 7:e33183. doi: 10.7554/eLife.33183
Lu, H., Jiang, W., Yang, H., Qin, Z., Guo, S. E., Hu, M., et al. (2017). Doxycycline affects gene expression profiles in aortic tissues in a rat model of vascular calcification. Microvasc. Res. 114, 12–18. doi: 10.1016/j.mvr.2017.04.007
Lubeseder-Martellato, C., Guenzi, E., Jorg, A., Topolt, K., Naschberger, E., Kremmer, E., 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. 161, 1749–1759. doi: 10.1016/S0002-9440(10)64452-5
MacMicking, J. D. (2004). IFN-inducible GTPases and immunity to intracellular pathogens. Trends Immunol. 25, 601–609. doi: 10.1016/j.it.2004.08.010
Matter, K., Aijaz, S., Tsapara, A., and Balda, M. S. (2005). Mammalian tight junctions in the regulation of epithelial differentiation and proliferation. Curr. Opin. Cell Biol. 17, 453–458. doi: 10.1016/j.ceb.2005.08.003
Mortison, J. D., Schenone, M., Myers, J. A., Zhang, Z., Chen, L., Ciarlo, C., et al. (2018). Tetracyclines modify translation by targeting key human rRNA substructures. Cell Chem. Biol. 25, 1506–1518.e13. doi: 10.1016/j.chembiol.2018.09.010
Mykytyn, K., and Askwith, C. (2017). G-protein-coupled receptor signaling in cilia. Cold Spring Harb. Perspect. Biol. 9:a028183. doi: 10.1101/cshperspect.a028183
Patro, R., Duggal, G., Love, M. I., Irizarry, R. A., and Kingsford, C. (2017). Salmon provides fast and bias-aware quantification of transcript expression. Nat. Methods 14, 417–419. doi: 10.1038/nmeth.4197
Ransick, A., Lindstrom, N. O., Liu, J., Zhu, Q., Guo, J. J., Alvarado, G. F., et al. (2019). Single-cell profiling reveals sex, lineage, and regional diversity in the mouse kidney. Dev. Cell 51, 399–413.e7. doi: 10.1016/j.devcel.2019.10.005
Rauchman, M. I., Nigam, S. K., Delpire, E., and Gullans, S. 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. doi: 10.1152/ajprenal.1993.265.3.F416
Reboredo, M., Kramer, M. G., Smerdou, C., Prieto, J., and De Las Rivas, J. (2008). Transcriptomic effects of Tet-on and mifepristone-inducible systems in mouse liver. Hum. Gene Ther. 19, 1233–1247. doi: 10.1089/hum.2008.057
Ribes, D., Fischer, E., Calmont, A., and Rossert, J. (2003). Transcriptional control of epithelial differentiation during kidney development. J. Am. Soc. Nephrol. 14(Suppl. 1), S9–S15. doi: 10.1097/01.asn.0000067647.05964.9f
Robinson, M. D., McCarthy, D. J., and Smyth, G. K. (2010). edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139–140. doi: 10.1093/bioinformatics/btp616
Schlimpert, M., Lagies, S., Budnyk, V., Muller, B., Walz, G., and Kammerer, B. (2018). Metabolic phenotyping of Anks3 depletion in mIMCD-3 cells – a putative nephronophthisis candidate. Sci. Rep. 8:9022. doi: 10.1038/s41598-018-27389-y
Schwab, K., Patterson, L. T., Aronow, B. J., Luckas, R., Liang, H. C., and Potter, S. S. (2003). A catalogue of gene expression in the developing kidney. Kidney Int. 64, 1588–1604. doi: 10.1046/j.1523-1755.2003.00276.x
Siroky, B. J., Kleene, N. K., Kleene, S. J., Varnell, C. D. Jr., Comer, R. G., Liu, J., et al. (2017). Primary cilia regulate the osmotic stress response of renal epithelial cells through TRPM3. Am. J. Physiol. Renal Physiol. 312, F791–F805. doi: 10.1152/ajprenal.00465.2015
Slaats, G. G., Wheway, G., Foletto, V., Szymanska, K., van Balkom, B. W., Logister, I., et al. (2015). Screen-based identification and validation of four new ion channels as regulators of renal ciliogenesis. J. Cell Sci. 128, 4550–4559. doi: 10.1242/jcs.176065
Sourdeval, M., Lemaire, C., Brenner, C., Boisvieux-Ulrich, E., and Marano, F. (2006). Mechanisms of doxycycline-induced cytotoxicity on human bronchial epithelial cells. Front. Biosci. 11:3036–3048. doi: 10.2741/2031
Stadnyk, A. W. (1994). Cytokine production by epithelial cells. FASEB J. 8, 1041–1047. doi: 10.1096/fasebj.8.13.7926369
Stork, P. J., and Schmitt, J. M. (2002). Crosstalk between cAMP and MAP kinase signaling in the regulation of cell proliferation. Trends Cell Biol. 12, 258–266. doi: 10.1016/s0962-8924(02)02294-8
Strait, K. A., Stricklett, P. K., Chapman, M., and Kohan, D. E. (2010). Characterization of vasopressin-responsive collecting duct adenylyl cyclases in the mouse. Am. J. Physiol. Renal Physiol. 298, F859–F867. doi: 10.1152/ajprenal.00109.2009
Sukhatme, V. P., Cao, X. M., Chang, L. C., Tsai-Morris, C. H., Stamenkovich, D., Ferreira, P. C., et al. (1988). A zinc finger-encoding gene coregulated with c-fos during growth and differentiation, and after cellular depolarization. Cell 53, 37–43. doi: 10.1016/0092-8674(88)90485-0
Torban, E., Eccles, M. R., Favor, J., and Goodyer, P. R. (2000). PAX2 suppresses apoptosis in renal collecting duct cells. Am. J. Pathol. 157, 833–842. doi: 10.1016/S0002-9440(10)64597-X
Valkova, N., and Kultz, D. (2006). Constitutive and inducible stress proteins dominate the proteome of the murine inner medullary collecting duct-3 (mIMCD3) cell line. Biochim. Biophys. Acta 1764, 1007–1020. doi: 10.1016/j.bbapap.2006.03.007
van den Bogert, C., Dontje, B. H., Holtrop, M., Melis, T. E., Romijn, J. C., van Dongen, J. W., et al. (1986). Arrest of the proliferation of renal and prostate carcinomas of human origin by inhibition of mitochondrial protein synthesis. Cancer Res. 46, 3283–3289.
van den Bogert, C., Dontje, B. H., and Kroon, A. M. (1985). The antitumour effect of doxycycline on a T-cell leukaemia in the rat. Leuk. Res. 9, 617–623. doi: 10.1016/0145-2126(85)90142-0
Vitale, G., de Herder, W. W., van Koetsveld, P. M., Waaijers, M., Schoordijk, W., Croze, E., et al. (2006). IFN-beta is a highly potent inhibitor of gastroenteropancreatic neuroendocrine tumor cell growth in vitro. Cancer Res. 66, 554–562. doi: 10.1158/0008-5472.CAN-05-3043
Wei, Q., Khan, I. K., Ding, Z., Yerneni, S., and Kihara, D. (2017). NaviGO: interactive tool for visualization and functional similarity and coherence analysis with gene ontology. BMC Bioinformatics 18:177. doi: 10.1186/s12859-017-1600-5
Wishart, J. A., Hayes, A., Wardleworth, L., Zhang, N., and Oliver, S. G. (2005). Doxycycline, the drug used to control the tet-regulatable promoter system, has no effect on global gene expression in Saccharomyces cerevisiae. Yeast 22, 565–569. doi: 10.1002/yea.1225
Yamaguchi, T., Pelling, J. C., Ramaswamy, N. T., Eppler, J. W., Wallace, D. P., Nagao, S., et al. (2000). cAMP stimulates the in vitro proliferation of renal cyst epithelial cells by activating the extracellular signal-regulated kinase pathway. Kidney Int. 57, 1460–1471. doi: 10.1046/j.1523-1755.2000.00991.x
Yu, Z., Kong, Q., and Kone, B. C. (2013). Aldosterone reprograms promoter methylation to regulate alphaENaC transcription in the collecting duct. Am. J. Physiol. Renal Physiol. 305, F1006–F1013. doi: 10.1152/ajprenal.00407.2013
Zhao, B., Ye, X., Yu, J., Li, L., Li, W., Li, S., et al. (2008). TEAD mediates YAP-dependent gene induction and growth control. Genes Dev. 22, 1962–1971. doi: 10.1101/gad.1664408
Keywords: doxycycline, RNA-seq, transcriptional response, cell proliferation, mIMCD3
Citation: Jung HJ, Coleman R, Woodward OM and Welling PA (2021) Doxycycline Changes the Transcriptome Profile of mIMCD3 Renal Epithelial Cells. Front. Physiol. 12:771691. doi: 10.3389/fphys.2021.771691
Received: 06 September 2021; Accepted: 18 October 2021;
Published: 05 November 2021.
Edited by:
Carolyn Mary Ecelbarger, Georgetown University, United StatesReviewed by:
Lihe Chen, National Institutes of Health Clinical Center (NIH), United StatesMohammad Al-bataineh, University of Pittsburgh, United States
Copyright © 2021 Jung, Coleman, Woodward and Welling. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Paul A. Welling, pwellin1@jhmi.edu