Novel MAPK-dependent and -independent tubulogenes identified via microarray analysis of 3D-cultured Madin-Darby canine kidney cells

A Qualitative Change in the Transcriptome Occurs after the First Cell Cycle and Coincides with Lumen Establishment during MDCKII Cystogenesis

Madin-Darby canine kidney II (MDCKII) cells are widely used to study epithelial morphogenesis. To better understand this process, we performed time course RNA-seq analysis of MDCKII 3D cystogenesis, along with polarized 2D cells for comparison. Our study reveals a biphasic change in the transcriptome that occurs after the first cell cycle and coincides with lumen establishment. This change appears to be linked to translocation of β-catenin, supported by analyses with AVL9- and DENND5A-knockdown clones, and regulation by HNF1B, supported by ATAC-seq study. These findings indicate a qualitative change model for transcriptome remodeling during epithelial morphogenesis, leading to cell proliferation decrease and cell polarity establishment. Furthermore, our study reveals that active mitochondria are retained and chromatin accessibility decreases in 3D cysts but not in 2D polarized cells. This indicates that 3D culture is a better model than 2D culture for studying epithelial morphogenesis.

•

The transcriptome switches after the first cell cycle and during MDCKII lumenogenesis

•

The transcriptome switch is linked to β-catenin translocation and HNF1B activation

•

Chromatin accessibility decreases during MDCKII cystogenesis

•

Active mitochondria are maintained in 3D, but not 2D, epithelial morphogenesis

Identifying ultrasensitive HGF dose-response functions in a 3D mammalian system for synthetic morphogenesis

Nonlinear responses to signals are widespread natural phenomena that affect various cellular processes. Nonlinearity can be a desirable characteristic for engineering living organisms because it can lead to more switch-like responses, similar to those underlying the wiring in electronics. Steeper functions are described as ultrasensitive, and can be applied in synthetic biology by using various techniques including receptor decoys, multiple co-operative binding sites, and sequential positive feedbacks. Here, we explore the inherent non-linearity of a biological signaling system to identify functions that can potentially be exploited using cell genome engineering. For this, we performed genome-wide transcription profiling to identify genes with ultrasensitive response functions to Hepatocyte Growth Factor (HGF). We identified 3,527 genes that react to increasing concentrations of HGF, in Madin-Darby canine kidney (MDCK) cells, grown as cysts in 3D collagen cell culture. By fitting a generic Hill function to the dose-responses of these genes we obtained a measure of the ultrasensitivity of HGF-responsive genes, identifying a subset with higher apparent Hill coefficients (e.g. MMP1, TIMP1, SNORD75, SNORD86 and ERRFI1). The regulatory regions of these genes are potential candidates for future engineering of synthetic mammalian gene circuits requiring nonlinear responses to HGF signalling.

Grainyhead-like 2 inhibits the coactivator p300, suppressing tubulogenesis and the epithelial-mesenchymal transition

GRHL2 suppresses EMT to give a default epithelial phenotype. GRHL2 inhibits this process through the histone acetyltransferase coactivator p300, repressing the partial EMT and preventing induction of MMPs. The results demonstrate novel roles for p300 and GRHL2 in promoting or suppressing EMT in morphogenesis and tumor progression.

Developmental morphogenesis and tumor progression require a transient or stable breakdown of epithelial junctional complexes to permit programmed migration, invasion, and anoikis resistance, characteristics endowed by the epithelial–mesenchymal transition (EMT). The epithelial master-regulatory transcription factor Grainyhead-like 2 (GRHL2) suppresses and reverses EMT, causing a mesenchymal–epithelial transition to the default epithelial phenotype. Here we investigated the role of GRHL2 in tubulogenesis of Madin–Darby canine kidney cells, a process requiring transient, partial EMT. GRHL2 was required for cystogenesis, but it suppressed tubulogenesis in response to hepatocyte growth factor. Surprisingly, GRHL2 suppressed this process by inhibiting the histone acetyltransferase coactivator p300, preventing the induction of matrix metalloproteases and other p300-dependent genes required for tubulogenesis. A 13–amino acid region of GRHL2 was necessary for inhibition of p300, suppression of tubulogenesis, and interference with EMT. The results demonstrate that p300 is required for partial or complete EMT occurring in tubulogenesis or tumor progression and that GRHL2 suppresses EMT in both contexts through inhibition of p300.

Stable Ectopic Expression of ST6GALNAC5 Induces Autocrine MET Activation and Anchorage-Independence in MDCK Cells

The epithelial-mesenchymal transition (EMT) is a complex cancer progression that can boost the metastatic potential of transformed cells by inducing migration, loss of cell adhesion, and promoting proliferation under anchorage-independent conditions. A DNA microarray analysis was performed comparing parental anchorage-dependent MDCK cells and anchorage-independent MDCK cells that were engineered to express human siat7e (ST6GALNAC5). The comparison identified several genes involved in the EMT process that were differentially expressed between the anchorage-dependent and the anchorage-independent MDCK cell lines. The hepatocyte growth factor gene (hgf) was found to be over-expressed in the engineered MDCK-siat7e cells at both transcription and protein expression levels. Phosphorylation analysis of the MET receptor tyrosine kinase confirmed the activation of an autocrine loop of the HGF/ MET signaling pathway in the MDCK-siat7e cells. When MET activities were suppressed by using the small-molecular inhibitor drug PF-02341066 (Crizotinib), the anchorage-independent MDCK-siat7e cells reverted to the cellular morphology of the parental anchorage-dependent MDCK cells. These observations indicate that the MET receptor plays a central role in the growth properties of the MDCK cells and its phosphorylation status is likely dependent on sialylation. Further investigation of the downstream signaling targets in the MET network showed that the degree of MDCK cell adhesion correlated with secretion levels of a matrix metalloproteinase, MMP1, suggesting a role of metalloproteinases in the EMT process. These results demonstrate that in addition to its application in biotechnology processes, MDCK-siat7e may serve as a model cell for metastasis studies to decipher the sequence of events leading up to the activation of EMT.

High-throughput mRNA and miRNA profiling of epithelial-mesenchymal transition in MDCK cells

BACKGROUND

Epithelial-mesenchymal transition (EMT) is an important process in embryonic development, especially during gastrulation and organ formation. Furthermore EMT is widely observed in pathological conditions, e.g., fibrosis, tumor progression and metastasis. Madin-Darby Canine Kidney (MDCK) cells are widely used for studies of EMT and epithelial plasticity. MDCK cells show an epithelial phenotype, while oncogenic Ras-transformed MDCK (MDCK-Ras) cells undergo EMT and show a mesenchymal phenotype.

METHODS

RNA-Seq and miRNA-Seq analyses were performed on MDCK and MDCK-Ras cells. Data were validated by qRT-PCR. Gene signature analyses were carried out to identify pathways and gene ontology terms. For selected miRNAs target prediction was performed.

RESULTS

With RNA-Seq, mRNAs of approximately half of the genes known for dog were detected. These were screened for differential regulation during Ras-induced EMT. We went further and performed gene signature analyses and found Gene Ontology (GO) terms and pathways important for epithelial polarity and implicated in EMT. Among the identified pathways, TGFβ1 emerged as a central signaling factor in many EMT related pathways and biological processes. With miRNA-Seq, approximately half of the known canine miRNAs were found expressed in MDCK and MDCK-Ras cells. Furthermore, among differentially expressed miRNAs, miRNAs that are known to be important regulators of EMT were detected and new candidates were predicted. New dog miRNAs were discovered after aligning our reads to that of other species in miRBase. Importantly, we could identify 25 completely novel miRNAs with a stable hairpin structure. Two of these novel miRNAs were differentially expressed. We validated the two novel miRNAs with the highest read counts by RT-qPCR. Target prediction of a particular novel miRNA highly expressed in mesenchymal MDCK-Ras cells revealed that it targets components of epithelial cell junctional complexes. Combining target prediction for the most upregulated miRNAs and validation of the targets in MDCK-Ras cells with pathway analysis allowed us to identify two novel pathways, e.g., JAK/STAT signaling and pancreatic cancer pathways. These pathways could not be detected solely by gene set enrichment analyses of RNA-Seq data.

CONCLUSION

With deep sequencing data of mRNAs and miRNAs of MDCK cells and of Ras-induced EMT in MDCK cells, differentially regulated mRNAs and miRNAs are identified. Many of the identified genes are within pathways known to be involved in EMT. Novel differentially upregulated genes in MDCK cells are interferon stimulated genes and genes involved in Slit and Netrin signaling. New pathways not yet linked to these processes were identified. A central pathway in Ras induced EMT is TGFβ signaling, which leads to differential regulation of many target genes, including miRNAs. With miRNA-Seq we identified miRNAs involved in either epithelial cell biology or EMT. Finally, we describe completely novel miRNAs and their target genes.