The centrosomal protein 83 (CEP83) regulates human pluripotent stem cell differentiation toward the kidney lineage
Abstract
During embryonic development, the mesoderm undergoes patterning into diverse lineages including axial, paraxial, and lateral plate mesoderm (LPM). Within the LPM, the so-called intermediate mesoderm (IM) forms kidney and urogenital tract progenitor cells, while the remaining LPM forms cardiovascular, hematopoietic, mesothelial, and additional progenitor cells. The signals that regulate these early lineage decisions are incompletely understood. Here, we found that the centrosomal protein 83 (CEP83), a centriolar component necessary for primary cilia formation and mutated in pediatric kidney disease, influences the differentiation of human-induced pluripotent stem cells (hiPSCs) toward IM. We induced inactivating deletions of CEP83 in hiPSCs and applied a 7-day in vitro protocol of IM kidney progenitor differentiation, based on timed application of WNT and FGF agonists. We characterized induced mesodermal cell populations using single-cell and bulk transcriptomics and tested their ability to form kidney structures in subsequent organoid culture. While hiPSCs with homozygous CEP83 inactivation were normal regarding morphology and transcriptome, their induced differentiation into IM progenitor cells was perturbed. Mesodermal cells induced after 7 days of monolayer culture of CEP83-deficient hiPCS exhibited absent or elongated primary cilia, displayed decreased expression of critical IM genes (PAX8, EYA1, HOXB7), and an aberrant induction of LPM markers (e.g. FOXF1, FOXF2, FENDRR, HAND1, HAND2). Upon subsequent organoid culture, wildtype cells differentiated to form kidney tubules and glomerular-like structures, whereas CEP83-deficient cells failed to generate kidney cell types, instead upregulating cardiomyocyte, vascular, and more general LPM progenitor markers. Our data suggest that CEP83 regulates the balance of IM and LPM formation from human pluripotent stem cells, identifying a potential link between centriolar or ciliary function and mesodermal lineage induction.
Editor's evaluation
This paper describes a novel role of the centrosomal protein CEP83 in mesoderm patterning. Specifically, CEP83 is shown to be required for the formation of the intermediate mesoderm and kidney progenitor tissue derived from it. In a CEP83 knockout situation in human pluripotent stem cells, a shift to lateral plate mesoderm at the expense of intermediate mesoderm occurs, which is convincingly demonstrated. This work, therefore, presents important new insights into the processes involved in mesodermal lineage induction and the fine-tuning of kidney differentiation.
https://doi.org/10.7554/eLife.80165.sa0Introduction
During mammalian embryonic development, the mesoderm forms axial, paraxial, and lateral plate domains that harbor precursor cells for distinct organ systems. Forming as a major part of the lateral plate mesoderm (LPM), the intermediate mesoderm (IM) harbors progenitor cells of all kidney epithelial cells (Davidson et al., 2019), whereas the remaining LPM contributes progenitors of various cell types, including cells of the cardiovascular system (Prummel et al., 2020). The molecular and cellular mechanisms that drive induction of the IM and distinct LPM domains during embryonic development are not fully understood.
The centrosomal protein 83 (CEP83) is a component of distal appendages (DAPs) of centrioles. DAPs are involved in the anchoring of the mother centriole to the cell membrane, an early and critical step in ciliogenesis (Lo et al., 2019; Tanos et al., 2013; Yang et al., 2018; Kurtulmus et al., 2018; Wheway et al., 2015; Bowler et al., 2019; Failler et al., 2014; Shao et al., 2020; Mansour et al., 2021). CEP83 recruits other DAP components to the ciliary base, and loss of CEP83 disrupts ciliogenesis (Tanos et al., 2013). In radial glial progenitors, removal of CEP83 disrupts DAP assembly and impairs the anchoring of the centrosome to the apical membrane as well as primary ciliogenesis (Yang et al., 2018; Shao et al., 2020). Mutations of CEP83 in humans have been associated with infantile nephronophthisis (Failler et al., 2014), an early onset kidney disease that results in end-stage renal disease before the age of 3 years (Hildebrandt, 2004; Luo and Tao, 2018) and additional organ anomalies (Failler et al., 2014). To date, how the loss of CEP83 function contributes to aberrant kidney development remains unclear.
Human-induced pluripotent stem cells (hiPSCs) provide useful tools to study molecular mechanisms of cellular differentiation. Protocols for the induction of kidney organoids from iPSC have been successfully developed (Takahashi et al., 2007; Morizane et al., 2015; Taguchi et al., 2014; Takasato et al., 2015; Freedman et al., 2015; Kumar et al., 2019). The protocol by Takasato et al. uses stepwise exposure of iPSC to WNT and FGF agonists in a monolayer culture system for a 7-day period, which results in the induction of cells with a transcriptional phenotype resembling kidney progenitors in the IM (Takasato et al., 2015). Transfer of these cells to an organoid culture system followed by another series of WNT and FGF signals results in the differentiation of three-dimensional (3D) kidney organoids composed of different kidney cell types, including glomerular and tubular cells. Genome editing studies have previously been used to study the effects of genetic defects associated with kidney diseases on kidney differentiation in human iPSC systems (Freedman et al., 2015; Tan et al., 2018; Boyle et al., 2008; Kobayashi et al., 2008; Howden et al., 2019; Kuraoka et al., 2020). Here, we studied the effect of an induced knockout of CEP83 in human iPSCs on kidney organoid differentiation. We uncovered a novel role of CEP83 in determining the balance of IM versus LPM differentiation, implicating a centrosomal protein in early mesodermal lineage decisions.
Results
CEP83 is essential for the differentiation of hiPSCs into kidney cells
To investigate the effect of CEP83 loss on the differentiation of hiPSCs into IM kidney progenitors, we applied CRISPR-Cas9 technology to induce a null mutation in the CEP83 gene in hiPSCs (Figure 1A). Three hiPSCs clones designated CEP83−/− (KO1, KO2, and KO3) carried deletions within CEP83 exon 7, each of which led to an induction of a premature stop codon resulting in a predicted truncated protein (Figure 1B–D and Figure 1—figure supplement 1A). These clones exhibited a complete loss of CEP83 protein by immunoblotting (Figure 1E). Three wildtype clones were derived as controls (WT1, WT2, and WT3). All six clones were morphologically indistinguishable (by brightfield microscopy) and had similar overall gene expression profiles (by bulk RNA-seq and qRT-PCR), including pluripotency and lineage marker expression (Figure 1—figure supplement 1B, C, and Figure 1—figure supplement 2A, B). In KO clones, the anticipated altered transcripts of CEP83 were detectable based on bulk RNA-seq (data not shown). Single nucleotide polymorphism - analysis confirmed identical karyotypes of all six clones (Figure 1—figure supplement 2C).
-
Figure 1—source data 1
- https://cdn.elifesciences.org/articles/80165/elife-80165-fig1-data1-v2.pdf
-
Figure 1—source data 2
- https://cdn.elifesciences.org/articles/80165/elife-80165-fig1-data2-v2.xlsx
-
Figure 1—source data 3
- https://cdn.elifesciences.org/articles/80165/elife-80165-fig1-data3-v2.zip
Together, these findings confirmed the successful deletion of CEP83 in iPSCs without any overt direct cellular phenotypic consequences. We applied a 7-day monolayer protocol using timed application of WNT and FGF agonists as reported by Takasato et al., 2015 to differentiate WT and KO hiPSCs into IM kidney progenitors (Takahashi et al., 2007; Morizane et al., 2015; Taguchi et al., 2014; Figure 2A).
After 7 days of culture (D7), WT and KO cells exhibited an indistinguishable morphology by bright field microscopy (Figure 2B and C). Immunostaining for acetylated tubulin, however, indicated abnormal primary cilia formation in CEP83-deficient cells (Figure 2D and E). The number of ciliated cells was reduced from approximately 30% (in WT clones) to less than 10% (in KO clones) (Figure 2F). Among ciliated cells, the length of cilia was increased from 2 to 5 µm (in WT clones) to 5–13 µm (in KO clones) (Figure 2G). This indicated that CEP83−/− hiPSCs differentiated toward IM progenitors exhibited ciliary abnormalities. To analyze the induced IM kidney progenitor cells functionally, we collected D7 WT and CEP83−/− cells and placed them into an organoid culture system again applying timed WNT and FGF agonists to foster differentiation of mature kidney cell types, as previously reported (Takasato et al., 2015; Figure 2A). Organoids harvested from WT clones after a total of 25 days of culture (D25) had formed patterned kidney epithelial-like structures, including NPHS1-positive glomerulus-like structures, Lotus tetragonolobus lectin-positive proximal tubule-like, and E-cadherin (CDH1)-positive distal tubule-like structures (Figure 3A, C and E, and Figure 3—figure supplement 1). In contrast, CEP83−/− organoids at day 25 were composed of monomorphic cells with a mesenchyme-like appearance, which stained negative for an array of kidney cell markers (Figure 3B, D and F).
Kidney epithelial-like structures formed only in WT but not in CEP83−/− organoids (Figure 3G). Similar to the findings in day 7 cells reported above, primary cilia were found in fewer cells of CEP83−/− organoids (<5% of cells) and were abnormally elongated (Figure 2—figure supplement 1).
Next, bulk RNA sequencing of WT (WT1, WT2, WT3) and CEP83−/− (KO1, KO2, KO3) organoids was carried out to evaluate differential gene expression on a genome-wide level, and RT-PCR was used to validate selected genes. Hierarchical clustering of the samples indicated strong gene expression differences between WT and CEP83−/− samples (Figure 3—figure supplement 2). Several genes associated with kidney development and kidney epithelial differentiation were differentially expressed with high expression in WT organoids but showed comparatively low or absent expression in CEP83−/− organoids: these genes included kidney-specific lineage genes (PAX2, PAX8), and lineage/differentiation markers of glomerular cells (NPHS1, PODXL, WT1, PTPRO), proximal (HNF1B, LRP2, CUBN), and distal (EMX2, MAL2, EPCAM, GATA3) kidney epithelial cells (Figure 3H–L, Figure 3—figure supplement 2). This indicated that CEP83−/− IM progenitors failed to differentiate into kidney cells, suggesting that CEP83 function is necessary to complete essential steps in the process of differentiation from pluripotent stem cells to kidney cells.
CEP83 deficiency associates with molecular defects of nephron progenitor cells
We next aimed to gain molecular insights into the lineage impact of CEP83 deficiency during the course of kidney epithelial differentiation. Since no global transcriptomic differences were detectable between WT and CEP83−/− hiPSCs prior to differentiation (see above), we focused on mesodermal cell stages induced at D7, which displayed mild overall gene expression differences between WT- and CEP83-deficient cells as detected by bulk RNA sequencing (Figure 4—figure supplement 1).
A marked upregulation of nephron progenitor marker genes (GATA3, HOXB7, HOXD11, EYA1) (Bilous et al., 1992; Grote et al., 2006; Kress et al., 1990; Srinivas et al., 1999; Wellik et al., 2002; Mugford et al., 2008a; Ruf et al., 2004; Sajithlal et al., 2005) was observed in both WT and CEP83−/− cells at day 7 (Figure 4—figure supplement 2), suggesting that the differentiation path of pluripotent CEP83−/− cells to IM nephron progenitors was largely intact. To understand the potential molecular defects at the IM stage in more detail, we performed single-cell RNA (scRNA) sequencing on D7 WT and CEP83−/− cells (representing two different hiPSC clones for each condition differentiated in two separate experiments). We obtained transcriptomes from 27,328 cells, representing clones WT1 (experiment 1: 3768 cells), WT2 (experiment 2: 5793 cells), KO1 (experiment 1: 8503 cells), and KO2 (experiment 2: 9264 cells). Principal component analysis (PCA) using pseudo-bulk expression data of the top 1000 highly variable genes (HVGs) indicated that the first principal component (dimension 1, explaining 54% of expression variation) was driven by the genotype (WT vs. KO), while the second principal component (dimension 2, explaining 51% of expression variation) was driven by a batch effect of the two experiments (Figure 4A). We combined all cells and generated a uniform manifold approximation and projection plot uncovering 10 different cell states/clusters (0–9; Figure 4B). We identified marker genes for each cluster (Figure 4C), indicating that clusters 1, 3, and 4 represented kidney progenitors/nascent nephrons (expressing, e.g., PAX8, EYA1, HOXB7) in different phases of the cell cycle. Other clusters represented as-of-yet uncharacterized cell types, which were consistent with previous single-cell transcriptome analyses from iPSC-derived cells induced by the same induction protocol (Subramanian et al., 2019; Low et al., 2019). Each of the four samples (WT1, WT2, KO1, and KO2) contributed to each cluster (Figure 4D). However, one cluster representing damaged cells (cluster 5) was observed at numerically higher percentages in KO cells compared to WT cells. Cluster 5 cells expressed high levels of mitochondrial RNAs, and staining for active caspase 3 demonstrated an increased percentage of apoptotic cells in KO samples compared to controls (Figure 4—figure supplement 3). We focused on kidney progenitors (clusters 1, 3, and 4) and found that a numerically lower percentage of KO cells (11.9 and 12.5% in KO clones) contributed to cluster 1 when compared with WT cells (25.9 and 36.3% in WT clones) (Figure 5A). In contrast, similar percentages of WT and KO cells were represented in kidney progenitor clusters 3 and 4 (Figure 5B and C). Differential gene expression analysis in these three clusters indicated significantly lower expression of kidney progenitor markers PAX8, EYA1, CITED1, and HOXB7 in KO cells from clusters 1, 3, and 4 when compared to WT cells (Figure 5D, E and F; Figure 5—figure supplements 1–2). Interestingly, scRNA sequencing data also showed downregulated expression of genes encoding ciliary proteins, including OFD1, PCM1, and RAB11A (Ferrante et al., 2006; Dammermann and Merdes, 2002; Knödler et al., 2010; Figure 5—figure supplement 3), consistent with the ciliogenesis defects in CEP83 knockout cells. These results indicate that CEP83 deficiency remained permissive with initial kidney progenitor induction, but that these cells exhibited mild molecular defects detectable by differential expression of kidney progenitor genes, which potentially contributed to the failure of CEP83-deficient cells to further differentiate toward mature kidney cell types.
CEP83 deficiency promotes ectopic induction of lateral plate mesoderm-like cells followed by an expansion of cardiac and vascular progenitors
We next inspected single-cell transcriptomes and bulk RNA sequencing data from D7 cells for genes that were upregulated in CEP83−/− cells compared to WT cells. From this analysis, we observed a consistent upregulation of genes that are normally expressed in early LPM, including OSR1, FOXF1, FOXF2, FENDRR, HAND1, HAND2, CXCL12, GATA5, and GATA6 (Prummel et al., 2020; Mugford et al., 2008b; Mae et al., 2013; Mahlapuu et al., 2001; Ormestad et al., 2004; Wilm et al., 2004; Wotton et al., 2008; Grote et al., 2013; Schindler et al., 2014; Tsuchihashi et al., 2011; McFadden et al., 2005; Firulli et al., 1998; Risebro et al., 2006; Angelo et al., 2000; Perens et al., 2016; Salcedo and Oppenheim, 2003; Liekens et al., 2010; Loh et al., 2016; Koutsourakis et al., 1999; Holtzinger and Evans, 2007; Reiter et al., 1999; Pikkarainen et al., 2004; Laverriere et al., 1994; Zhao et al., 2008; Prummel et al., 2021; Figure 6A–I). This suggested that CEP83−/− cells entered an aberrant differentiation path assuming a phenotype indicative of broader LPM instead of more specific IM. To further substantiate this idea, we restricted the analysis to progenitor cells of clusters 1, 3, and 4 and to cells from cluster 0, which exhibited a mesenchymal transcriptome fingerprint (see Figure 4C). Within each cell, we analyzed the expression of LPM markers (FOXF1, HAND1, HAND2, and CXCL12) and of more restricted IM markers (PAX8, EYA1, and HOXB7) (Figure 6—figure supplement 1). This analysis indicated that WT cells of these clusters exhibited an IM-like phenotype, while KO cells were shifted toward an LPM-like phenotype. The common IM/LPM marker OSR1 was expressed at higher level in KO cells comparing to the WT cells.
We then inspected RNA-seq data from WT and KO organoids at day 25 for the expression of LPM genes and markers of LPM derivatives. The expression of several LPM genes (OSR1, FOXF1, FOXF2, FENDRR, HAND1, HAND2, and CXCL12) was strongly upregulated in KO cells compared to WT cells suggesting that an LPM-like cell pool persisted in D25 KO organoids (Figure 6A–I). To further substantiate the potential differentiation of the CEP83-mutant cells into broadly LPM-like cells, we compared genes that were upregulated genes in D25 organoids (in total, 397 genes) with LPM genes that were previously identified by single-cell transcriptomics of sorted post-gastrulation LPM cells from developing zebrafish (Prummel et al., 2021; Prummel et al., 2020; Prummel et al., 2022). Our targeted comparison documented that CEP83−/− organoids showed significant enrichment for expression of orthologs of early LPM genes (p=0.006) (Figure 6—figure supplement 2), including OSR1, CXCL12, HAND1/2, KCTD12, PIK3R3, and ZBTB2. A subset of LPM genes enriched for expression in CEP83-mutant cells at D25 of differentiation was indicative of cardiac or cardiopharyngeal (ISL1, TBX1) as well as vascular progenitor (SOX7, SOX11, NAP1L3, LMO2, GATA2) differentiation (Morikawa and Cserjesi, 2004; Cai et al., 2003; Kwon et al., 2009; Laugwitz et al., 2005; Moretti et al., 2006; Gao et al., 2019; Stennard and Harvey, 2005; Baldini, 2005; Chen et al., 2009; Figure 7A–G, Figure 6—figure supplement 2). Interestingly, three genes of the upregulated LPM genes, namely OSR1, FOXF1, and FOXF2, are downstream genes of hedgehog signaling (Han et al., 2017; Kugler et al., 2015). In addition, HAND2 is upstream of hedgehog signaling (Anderson et al., 2012), essential for early embryonic development and regulated by primary cilium (Kugler et al., 2015). Furthermore, bulk RNA sequencing data showed significant upregulation of hedgehog signaling components GLI1 and PTCH1 in CEP83−/− cells (Lee et al., 1997; Villavicencio et al., 2000; Figure 6—figure supplement 3).
Taken together, these observations document that hiPSCs deficient in CEP83 respond to an in vitro differentiation program toward kidney progenitors, yet diverge toward a broader LPM progenitor composition without significant IM instead.
Discussion
This study indicates a novel contribution of CEP83 in regulating the differentiation path from human pluripotent stem cells to kidney progenitors. We pinpoint a stage at day 7 of IM induction where CEP83 loss of function results in a decreased nephron progenitor pool with downregulation of critical kidney progenitor genes (PAX8, EYA1, HOXB7). At the same stage, genes typical of LPM specification (including FOXF1, FOXF2, FENDRR, HAND1, and HAND2) are upregulated (Figure 8). Functionally, these alterations are associated with an inability of CEP83-deficient cells to form kidney epithelia. Organoids derived from CEP83-deficient cells fail to induce any detectable nephron structures, suggesting a novel role for CEP83 during the specification of functional kidney progenitors in the mesoderm.
Our findings are relevant to understanding the cellular and molecular functions of CEP83 and might be relevant to the pathophysiology of human genetic diseases. To date, 11 patients with biallelic mutations of CEP83 have been reported, 8 of which displayed kidney phenotypes (Failler et al., 2014; Veldman et al., 2021; Haer-Wigman et al., 2017). Available kidney histologies identified microcystic tubular dilatations, tubular atrophy, thickened basement membranes, and interstitial fibrosis.
Extrarenal phenotypes included speech delay, intellectual disability, hydrocephalus, strabismus, retinal degeneration, retinitis pigmentosa, hepatic cytolysis, cholestasis, and portal septal fibrosis with mild thickening of arterial walls and an increase in the number of the biliary canalicules on liver biopsy. Among individuals with CEP83 mutations, all but one carried at least one missense mutation or short in-frame deletion, suggesting that CEP83 function may have been partially preserved. One individual with presumed full loss of CEP83 displayed a more severe phenotype with multiple organ dysfunction. It will be interesting to await future reports of additional CEP83 mutations in humans and whether complete loss of function alleles will result in broader mesoderm defects or renal agenesis. In this regard, it is interesting that mice with a targeted homozygous loss-of-function mutation of their CEP83 ortholog (Cep83tm1.1(KOMP)Vlcg) display midembryonic lethality (at E12.5) with evidence of severe developmental delay as early as E9.5 (https://www.mousephenotype.org/data/genes/MGI:1924298). These phenotypes are potentially consistent with the role of CEP83 in germ layer patterning and mesoderm development, but a more detailed phenotypical characterization of Cep83 knockout embryos would be required to substantiate this possibility.
The precise molecular and cellular mechanisms underlying our observations remain to be determined. CEP83 is a protein that is necessary for the assembly of DAPs and primary cilia formation in several cell types (Tanos et al., 2013; Yang et al., 2018; Shao et al., 2020; Kumar et al., 2021; Stinchcombe et al., 2015; Joo et al., 2013). A potential involvement of CEP83-mediated primary cilia formation in the findings reported here is suggested by obvious ciliary defects in CEP83-deficient cells at the D7 and the organoid stage (Figure 2D–G, Figure 2—figure supplement 1). These defects include reduced percentages of ciliated cells and elongated primary cilia in those cells that continue to form a primary cilium. In addition, CEP83-deficient cells displayed downregulated expression of several transcripts encoding ciliary components (Figure 5—figure supplement 3) and evidence of an activation of several hedgehog signaling associated genes. The primary cilium is critically involved in hedgehog signaling (Ho and Stearns, 2021). Moreover, Hedgehog signaling is important for mesodermal lineage decisions during gastrulation (Guzzetta et al., 2020). This raises the possibility that CEP83 controls mesodermal cell fate decisions by modulating hedgehog signaling in the mesoderm. Nevertheless, additional studies will be necessary to address this possibility. In addition, it remains unknown whether abnormal cilia formation in CEP83-deficient cells causally contributes to the cell fate phenotype.
We observed downregulated expression of the key nephron progenitor genes PAX8, EYA1, and HOXB7 in CEP83−/− cells at day 7, which might explain their failure to differentiate into kidney cells since each of these genes is essential for normal kidney development (Vincent et al., 1997; Pfeffer et al., 1998; Kumar et al., 1998; Bouchard et al., 2002; Xu et al., 1999; Patterson and Potter, 2004; Rojek et al., 2006). Furthermore, inductive signals from HOXB7-positive ureteric bud are known to maintain viability of nephron progenitor cells in the IM (Barasch et al., 1997), which may contribute to increased numbers of apoptotic cells we observed in CEP83−/− cells. Defects during nephron progenitor differentiation in the IM would be expected to result in severe kidney phenotypes such as renal agenesis or renal hypodysplasia. Defects of centriolar components or cilia have previously been linked to such phenotypes: in mice, centrosome amplification, i.e., the formation of excess centrosomes per cell severely disrupts kidney development, resulting in depletion of renal progenitors and renal hypoplasia (Dionne et al., 2018). In humans, loss of KIF14, a protein necessary for proper DAP assembly and cilium formation, has been associated with kidney malformations, including renal agenesis and renal dysplasia (Filges et al., 2014; Reilly et al., 2019; Pejskova et al., 2020). Furthermore, Kif3a, a ciliary protein involved in intraflagellar transport, is necessary for normal mesoderm formation, and kidney progenitor-specific defects of Kif3a have been associated with reduced nephron numbers (Takeda et al., 1999; Chi et al., 2013). Similarly, mouse genes encoding the ciliary intraflagellar transport proteins IFT25 and IFT27 have been associated with renal agenesis or renal hypoplasia (Desai et al., 2018; Quélin et al., 2018). Together, these studies highlight the importance of molecules involved in ciliogenesis for mesoderm and kidney progenitor development and suggest that CEP83 contributes to such processes by facilitating an early step of ciliogenesis. Nevertheless, the detailed molecular processes that link CEP83 function, cilia formation, and kidney progenitor specification remain to be determined.
The finding of various upregulated LPM markers in CEP83−/− cells starting from day 7 suggests that CEP83 function may be involved in fine-tuning the balance of LPM and IM, thereby contributing to lineage decisions during mesoderm formation. Crosstalk of LPM and IM has been reported previously in zebrafish, overexpression of LPM transcription factors Scl/Tal1 and Lmo2 induces ectopic vessel and blood specification while inhibiting IM formation (Gering et al., 2003). Furthermore, the LPM transcription factor Hand2 is critical in determining the size of the IM, while natively expressed in the IM-adjacent LPM progenitors that form mesothelia (Barbosa et al., 2007Perens et al., 2016; Prummel et al., 2022). Loss of Hand2 in zebrafish results in an expanded IM, whereas Hand2 overexpression reduces or abolishes the IM. Interestingly, HAND2 was among the most strongly induced transcripts in our CEP83−/− cells at day 7; connecting with the developmental role of Hand2 in IM formation, these observations suggest that HAND2 expression in CEP83-deficient cells may have contributed to the reduced numbers of nephron progenitor cells at this stage. Of note, CEP83-deficient cells at D25 expressed increased levels of LPM genes expressed in mesothelial (including OSR1, CXCL12, HAND1/2), cardiopharyngeal (including ISL1, TBX1), and endothelial/hematopoietic (including TAL1, LMO2, GATA2) progenitors (Prummel et al., 2020; Prummel et al., 2022). In sum, we propose a novel role for CEP83 in regulating the development of IM nephron progenitors, which may involve direct effects of CEP83 in the nephron progenitor differentiation program and indirect LPM-mediated effects on the IM. Future studies are warranted to delineate the molecular and cellular mechanisms underlying CEP83 function in LPM and specifically IM patterning.
Materials and methods
hiPSCs cell line
Request a detailed protocolWe used the human iPSC cell line BIHi005-A, which was generated by the Berlin Institute of Health (BIH). The hiPSCs were maintained in six-well plates (Corning, 353046) coated with Matrigel (Corning, 354277) and cultured in Essential 8 medium (E8, A1517001, Gibco-Thermo Fisher Scientific) supplemented with 10 µM Y-27632 (Rocki, Wako, 253–00513). Cells were authenticated and tested for the mycoplasma infection.
CRISPR CAS9 technology to generate CEP83−/− hiPSCs clones
Request a detailed protocolClustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 technology was used to generate CEP83−/− hiPSCs clones. We designed two CRISPR RNAs (crRNAs) (5'-GGCTGAAGTAGCGGAATTAA-AGG-3' and 5'-AAGAATACAGGTGCGGCAGT-TGG-3') using CRISPOR software (Concordet and Haeussler, 2018). The two crRNAs were annealed with trans-activating CRISPR RNA (tracrRNA) to form two guide RNAs (gRNA1 and gRNA2) and then formed a ribonucleoprotein (RNP) complex by incubating gRNA1 and gRNA2 separately with Alt-R S.p. Cas9 Nuclease V3 (1 μM concentration, IDT, 1081058). The hiPSCs were transfected with RNP complexes using Neon transfection system (Thermo Fisher Scientific, MPK5000; Yumlu et al., 2017) and Neon transfection 10 μl kit (Thermo Fisher Scientific, MPK10025) according to the manufacturer’s instructions. After 48 hr of transfection, we analyzed the editing efficiency in the pool by PCR genotyping.
For PCR genotyping, we isolated genomic DNA from the pool of transfected cells followed by PCR using Phire Tissue Direct PCR Master Mix (Thermo Scientific, F170S) according to the manufacturer’s instructions (Figure 1B). After confirming the editing efficiency in the pool, we generated single-cell clones by the clonal dilution method. We plated 500 single cells per well of a 6 well plate and picked 24 clones using a picking hood S1 (Max Delbrück Centre Stem Cell Core Facility). Then, clones were screened for homozygous deletions of CEP83 by PCR using Phire Tissue Direct PCR Master Mix. Selected knockout clones were further characterized for CEP83 loss of function on the DNA, RNA, and protein level. CEP83−/− clones (KO1, KO2, and KO3) were registered as (BIHi005-A-71, BIHi005-A-72, and BIHi005-A-73) in the European Human Pluripotent Stem Cell Registry (https://hpscreg.eu).
Single nucleotide polymorphism - karyotype
Request a detailed protocolTo assess karyotype integrity, copy number variation (CNV) analysis on the human Illumina OMNI-EXPRESS-8v1.6 BeadChip was used. In brief, genomic DNA was isolated from three WT (WT1, WT2, and WT3) and three KO (KO1, KO2, and KO3) clones using the DNeasy blood and tissue kit (Qiagen, Valencia, CA, United States), hybridized to the human Illumina OMNI-EXPRESS-8v1.6 BeadChip (Illumina), stained, and scanned using the Illumina iScan system according to a standard protocol (LaFramboise, 2009; Arsham et al., 2017; Haraksingh et al., 2017). The genotyping was initially investigated using the GenomeStudio 1 genotyping module (Illumina). Following that, KaryoStudio 1.3 (Illumina) was used to perform automatic normalization and identify genomic aberrations in detected regions by generating B-allele frequency and smoothed Log R ratio plots. To detect CNVs, the stringency parameters were set to 75 kb (loss), 100 kb (gain), and CN-LOH (loss of heterozygosity). KaryoStudio generates reports and displays chromosome, length, list of cytobands, and genes in CNV-affected regions.
Differentiation protocol
Request a detailed protocolWe used the protocol of Takasato to differentiate the hiPSCs into nephron organoids (Takasato et al., 2015). Briefly, hiPSCs were cultured first in APEL2 medium (Stem Cell Technologies, 05270) supplemented with 5% Protein Free Hybridoma Medium II (PFHMII, GIBCO, 12040077), and 8 μM CHIR99021 (R&D, 4423/10) for 5 days, with medium changes every 2 days. Then, the cells were cultured in APEL2 medium supplemented with 200 ng/ml FGF9 (R&D, 273-F9-025) and 1 μg/ml heparin (Sigma Aldrich, H4784-250MG) for 2 days. On day 7, the cells were washed with 1× Dulbecco’s PBS (DPBS, Thermo Fisher Scientific,14190–250), then trypsinized using trypsin EDTA-0.05% (Thermo Fisher Scientific, 25300–062) at 37°C for 3 min. The cells were counted and divided to achieve 1×106 cells per organoid and cultured into 3D organoid culture on 0.4-μm-pore polyester membrane of Corning 6-well Transwell cell culture plate (Corning-Sigma Aldrich, CLS3450-24EA). Four to five organoids were seeded on one membrane using a P100 wide-bore tip and cultured in APEL2 with 5 μM CHIR99021 at 37°C for 1 hr (CHIR99021 pulse). After the CHIR pulse, we changed the medium to APEL2 medium supplemented with 200 ng/ml FGF9 + 1 μg/ml heparin for 5 days with medium refreshing every 2 days. The organoids were then cultured only in APEL2 medium with 1 μg/ml heparin for additional 13 days. The total differentiation time is 25 days (7+18).
DNA isolation and PCR
Request a detailed protocolDNA was isolated from cells using DNeasy Blood & Tissue Kits (Qiagen, 69504). CEP83 primers were designed using Primer3 webtool (Supplementary file 1). PCR was done using Phusion high-fidelity DNA polymerase (Biolabs, New England, M0530) according to the manufacturer’s instructions. PCR results were visualized on 1.5% agarose gel using a BioDoc Analyze dark hood and software system (Biometra).
RNA isolation, RNA sequencing, and qPCR
Request a detailed protocolTotal RNA was isolated from the cells using RNAasy Mini Kit (QIAGEN, Hilden, Germany, 74104) following the manufacture’s instructions. The concentration, quality, and integrity of the extracted RNA were evaluated using Nanodrop (Thermo Scientific, Waltham, MA; USA), an Agilent 2100 Bioanalyzer, and the Agilent RNA 6000 Nano kit (Agilent Technologies, 5067–1511). 0.4 μg total RNA was used to obtain a poly A–enriched RNA library by Novogene (Cambridge, United Kingdom). Library concentration was performed using a Qubit fluorometer (HS RNA assay kit, Agilent Technologies). Library size was measured by Agilent 2100 bioanalyzer. The libraries were then subjected to 150 bp paired-end next-generation sequencing (Illumina NovaSeq 6000 S4 flow cells). Mutation visualization was performed using the Integrative Genomic Viewer tool (Robinson et al., 2011). Read counts of the sequenced RNA were normalized to transcripts per million (TPM). The TPM values of the variables were used to plot heatmaps and for principle component analysis (PCA) based on Pearson correlation, using self-written scripts in R (R Development Core Team (2011)) (version 4.0.4).
RNA was reverse transcribed using the RevertAid First Strand cDNA Synthesis Kit (Thermo Scientific). qPCR was performed using the FastStart Universal SYBR Green Master (Rox) mix (Hoffmann-La Roche) according to the manufacturer’s instructions. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) mRNA expression was calculated according to the ΔΔCt method. All primer pairs were designed using Primer3, purchased at BioTeZ (Berlin, Germany), and sequences are shown in Supplementary file 1.
Single-cell RNA sequencing (scRNA-seq)
Cells isolation and preparation
Request a detailed protocolDifferentiated cells at day 7 were washed twice with 1× DPBS, dissociated with Accumax solution, and resuspended in 1× DPBS. Then, cells were filtered, counted, and checked for viability.
Library preparation and single-cell sequencing
Request a detailed protocolSingle-cell 3’ RNA sequencing was performed using the 10× Genomics toolkit version v3.1 (Alles et al., 2017), according to the manufacturer’s instructions aiming for 10,000 cells. Obtained libraries were sequenced on Illumina NextSeq 500 sequencers.
Single-cell sequencing data analysis and clustering
Request a detailed protocolAfter sequencing and demultiplexing, fastq files were analyzed using Cellranger version 3.0.2. Gene expression matrices were then imported into R, and Seurat objects were created using the Seurat R package (version 4.0.5) (Stuart et al., 2019). The gene expression matrices were initially filtered by applying lower and upper cut-offs for the number of detected genes (500 and 6000, respectively). The filtered data were then log normalized and scaled according to the number of unique molecular identifiers. The normalized and scaled data derived from the four samples were then merged into one Seurat object. Clustering was performed using the first 20 principal components. We used the Seurat FindAllMarkers function to extract marker gene lists that differentiate between clusters with log fold change threshold ± 0.25 using only positive markers expressed in a minimum of 25% of cells. PCA was done using the first 20 principle components in R using the following libraries factoextra, FactoMineR, and ggplot2.
Protein extraction and immunoblotting
Request a detailed protocolProteins were extracted from hiPSCs using radioimmunoprecipitation assay (RIPA) buffer (Sigma-Aldrich, R0278) as described in details in supplementary data. 30 µg protein in RIPA buffer were mixed with 1× reducing (10% b-mercaptoethanol) NuPAGE loading buffer (Life Technologies, Carlsbad, CA), loaded on a precast polyacrylamide NuPage 4–12% Bis-Tris protein gel (Invitrogen, Carlsbad, CA, USA) and blotted on 0.45 µm pore size Immobilon-P Polyvinylidene difluoride membrane (EMD Millipore, Billerica, MA; USA). The membrane was blocked in 5% bovine serum albumin for 1 hr at room temperature and incubated overnight at 4°C with primary antibodies: anti-CEP83 produced in rabbit (1:500, Sigma-Aldrich) and anti-α-Tubulin produced in mouse (1:500, Sigma-Aldrich, T9026). Then, the membrane was incubated for 1 hr at room temperature with horseradish peroxidase-conjugated secondary antibodies (1:2000, Sigma-Aldrich, Saint Louis, MO, USA). Chemiluminescent reagent (Super Signal–West Pico; Thermo Scientific, Waltham, MA; USA) was used to detect the proteins. The spectra Multicolor Broad Range Protein Ladder (Thermo Fisher Scientific, USA) was used to evaluate the molecular weight of corresponding protein bands.
Histology and immunofluorescence staining
Request a detailed protocolCells at different time points were checked regularly under a confocal microscope (Leica DMI 6000 CEL) for differentiation progress. Quantitative analysis of nephron-like structure formation within each organoid (D25) was done on tile scanning images of each organoid by estimating the percentage of the organoid area composed of nephron-like structures using 13 WT and 9 KO organoids. Organoids were fixed in BD Cytofix buffer (554655, BD Biosciences) for 1 hr on ice. Then organoids were gradually dehydrated in increasing ethanol concentrations, cleared in xylene, and embedded in paraffin. Organoids were cut into 3.5-µm thick sections. The sections were deparaffinized, dehydrated, and stained in hematoxylin (Sigma-Aldrich, Saint Louis, MO) for 3 min and in 1% eosin (Sigma-Aldrich) for 2 min. For immunostaining, cells (day 7) and organoids (day 25) were fixed with BD Cytofix, permeabilized with BD Perm/Wash (554723, BD Biosciences), and blocked with blocking solution (1% BSA + 0.3% triton-X-100 in 1× DPBS) for 2 hr. Cells and organoids were incubated overnight at 4°C with primary antibodies (Supplementary file 2), then incubated with fluorescence-labeled secondary antibodies with 1:500 dilution including Cy3, Cy5, Alexa488, and Alexa647 (Jackson ImmunoResearch, Newmarket, UK) and Cy3 Streptavidin (Vector lab, Burlingame, USA) overnight at 4°C. DAPI was then used for nuclear staining (Cell signaling Technology, Danvers, MA, USA) with 1:300,000 dilution for 1 hr at RT. Finally, cells were mounted with Dako fluorescent mounting medium (Agilent Technologies). Images were taken using a SP8 confocal microscope (Carl Zeiss GmbH, Oberkochen, Germany). Quantitative analyses of acquired images were performed using ImageJ software (1.48 v; National Institutes of Health, Bethesda, MD).
Comparison to zebrafish LPM
Request a detailed protocolThe upregulated genes in CEP83−/− cells on day 7 and day 25 were compared with the top 20 orthologous genes identified in subclusters of zebrafish LPM identified by scRNA-seq, as deposited on ArrayExpress (E-MTAB-9727; Prummel et al., 2022).
Statistical analysis
Request a detailed protocolscRNA-seq was done on two biological replicates representing two different clones of CEP83−/− and control cells, respectively. All other experiments were performed using three biological replicates representing three independent clones of CEP83−/− and control cells at different time points. A common excel sheet for the genes present in both bulk RNA and scRNA sequencing was generated in R. The sheet includes a total of 20,894 genes and represents the TPM values of both groups (WT and KO) on day 0, day 7, and day 25 of differentiation. The maximum TPM (TPMmax) and the minimum TPM (TPMmin) were calculated for each gene across all samples. HVGs were calculated based on the ratio of TPMmax and TPMmin. For heatmaps and PCA analysis, the top 1000 HVGs were plotted with a selection of TPMmax >2 for each gene. Deregulated (upregulated and downregulated) genes between WT and KO groups were selected with expression criteria of TPM >2, fold change >1.5, and p-value calculated on log10 TPM <0.05. The unpaired two-tailed t-test was used to compare the two groups. All graphs were generated using GraphPad Prism 7.04 (GraphPad Software, San Diego, CA). Data are presented as mean ± SD.
Appendix 1
Supplemental methods
hiPSCs culture
We used the human iPSC cell line BIHi005-A, which was generated from a healthy donor by the Berlin Institute of Health (BIH) and supplied by the stem cell core facility at Max Delbrück Center for Molecular Medicine (Berlin). The hiPSCs were maintained in six-well plates (Corning, 353046), coated with Matrigel (Corning, 354277), and cultured in Essential 8 medium (E8, Gibco-Thermo Fisher Scientific, A1517001) supplemented with 10 µM Y-27632 (Rocki, Wako, 253–00513). The cells were split twice per week using EDTA/PBS or Accumax (Stem cell technology, 07921).
CRISPR CAS9 technology to generate CEP83−/− hiPSCs clones
Clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 technology was used to generate CEP83−/− hiPSCs clones. Two 20 bp-long CRISPR RNAs (crRNAs) were designed using CRISPOR software (Concordet and Haeussler, 2018) to selectively target exon7: (5'- GGCTGAAGTAGCGGAATTAA-AGG-3'); (5'-AAGAATACAGGTGCGGCAGT-TGG-3'). The crRNAs were ordered from Integrated DNA Technologies (IDT). The crRNAs (IDT) were annealed in equimolar concentrations with trans-activating CRISPR RNA (tracrRNA) to form two guide RNAs (gRNA1 and gRNA2), which were then conjugated separately with Alt-R S.p. Cas9 Nuclease V3 (1 μM concentration, IDT, 1081058) at room temperature for 1 hr to form ribonucleoprotein (RNP) complexes (RNP1 and RNP2).
One day prior to transfection, hiPSCs were split using Accumax solution and cultured in an equal proportion of E8 medium and StemFlex medium (Thermo Fisher Scientific, A3349401). The hiPSCs were transfected using a Neon transfection system (Thermo Fisher Scientific, MPK5000). Immediately before the transfection, the cells were dissociated, collected, and resuspended in Resuspension Buffer (Buffer R) that included in Neon transfection 100 μl kit (Thermo Fisher Scientific, MPK10025; Yumlu et al., 2017). Cells were transfected in 3 ml Electrolytic buffer (Buffer E) that included in the neon transfection kit and by using the Neon transfection system 10 µl tip. The used Neon Transfection parameters were voltage (1200 V), width (30 ms), and pulse (1). The transfected cells were cultured in StemFlex Medium with Rocki for 48 hr.
Then, quick DNA extraction and PCR were then done to test transfection efficiency according to the manufacturer’s instructions using Phire Tissue Direct PCR Master Mix (Thermo Scientific, F170S). The size of the PCR products was visualized on 1.5% agarose gel. After confirming the transfection’s success in the knockout cells, as shown in Figure 1B, the cells were dissociated and seeded at low densities for 24 hr. Then, 24 single cells were picked under a picking hood S1 (stem cell core facility-MDC, Buch) and cultured in StemFlex medium for 2 weeks. The clones were then tested for CEP83 mutation on the DNA level by PCR using Phire Tissue Direct PCR Master Mix. Finally, the selected clones were expanded and frozen in the Bambanker medium (Nippon Genetics, BB01) for further characterization. The selected clones were characterized for the mutation induction on the DNA, protein, and RNA level.
Differentiation protocol
We used the protocol of Takasato to differentiate the hiPSCs into nephron organoids (Takasato et al., 2015), the experiment was performed using three replicates per each group wildtype hiPSCs (WT1, WT2, and WT3) and CEP83−/− hiPSCs (KO1, KO2, and KO3). Two days prior to the differentiation, cultured hiPSCs on matrigel with 70–80% density were prepared for the differentiation. The cells were washed twice with 1× Dulbecco’s PBS (Thermo Fisher Scientific,14190–250), then cells were trypsinized using 1× Trypl E Select (Thermo Fisher Scientific, 12563011). Cells were incubated at 37°C for 3 min. Then, DMEM/F-12 medium (Thermo Fisher Scientific, 11320–033) was added to the cells to neutralize Trypl E. The cell suspension was mixed by pipetting (pipetting is maximum twice), then centrifuged at 300 g for 5 min. The cell pellet was washed and resuspended in 1 ml of E8 medium. Then the cells were counted using Countess chamber slide and Countess II Automated Cell Counter (Thermo Fisher Scientific). Then cells were centrifuged and resuspended in E8 media supplemented with 10 µM Rocki. Last, cells were cultured on a prepared coated matrigel six-well culture plates to obtain a density of 15×103 cells per cm2 and incubated overnight at 37°C CO2 incubator for 48 hr with daily medium change.
Immediately before the differentiation, the cells were checked under the microscope. Cells with 40–50% confluency were used for the differentiation. The E8 medium was changed into APEL2 medium (Stem Cell Technologies, 05270) with 5% Protein Free Hybridoma Medium II (PFHMII, GIBCO, 12040077) and 8 μM CHIR99021 (2 ml medium per a well of 6- well plate). Cells were incubated in a 37°C CO2 incubator for 5 days, with medium refreshing every 2 dys. Following the CHIR99021 phase, the medium was changed into a double volume of APEL2 medium (4 ml medium per a well of six-well plate) supplemented with 200 ng/ml FGF9 (R&D, 273-F9-025) and 1 μg/ml heparin (Sigma Aldrich, H4784-250MG) and were incubated in a 37°C CO2 incubator.
On day 7 of differentiation, the cells were washed, trypsinized with trypsin EDTA (0.05%), and incubated at 37°C for 3 min. Then, the cell suspension was transferred to a 50 ml tube containing 9 ml of MEF conditioned medium (R&D, AR005) to neutralize the trypsin. The cells were centrifuged and resuspended in APEL2 medium. Using a hemocytometer, the cells were counted, and the cell suspension was divided to achieve 1×106 cells (organoid) per 1.5 ml tube. All the tubes were centrifuged at 400× g for 3 min at RT. During centrifugation, six-well transwell cell culture plates (Corning-Sigma Aldrich, CLS3450-24EA) were prepared by adding 1.2 ml of APEL2 supplemented with 5 μM CHIR99021 to each well. Cell pellets were picked up by using a P1000 or P200 wide-bore tip. Pellets were carefully seeded onto the six-well transwell membrane with minimal APEL2 medium carryover and incubated at 37°C for 1 hr. Then the medium was changed into APEL2 medium supplemented with 200 ng/ml FGF9 plus 1 μg/ml heparin for further 5 days, with medium refreshing every 2 days. Finally, the medium to APEL2 medium with only heparin for further 13 days, with medium refreshing every 2 days.
DNA isolation and PCR
Cultured WT and CEP83−/− hiPSCs were washed, scrapped gently using a cell scraper (VWR, part of Avantor, 734–2602), and collected with a maximum 5×106 cell/ml for proper DNA extraction. The DNA was extracted using DNeasy Blood & Tissue Kits (Qiagen, 69504) following the manufacturer’s instructions. The concentrations and quality of the DNA were evaluated using Nanodrop (Thermo Scientific, Waltham, MA; USA). To detect CEP83 expression, 200 µg DNA was amplified by a standard PCR using Phusion high-fidelity DNA polymerase (Biolabs, New England, M0530). The master mix was calculated according to the manufacturer’s instructions. Primers are designed using Primer3 webtool, Supplementary file 2. PCR was carried out in a thermocycler as follow: initial denaturation at 98°C for 30 s, 35–40 cycles of 30 s at 98°C, 30 s at 63.5°C, and 30 s at 72°C; final elongation step at 72°C for 10 min. The PCR results were checked on 1.5% agarose gel and analyzed using a BioDoc Analyze dark hood and software system (Biometra).
RNA isolation, RNA sequencing, and qPCR
Total RNA was isolated from the cells at three time points: day 0 (hiPSCs), day 7 (IM), and day 25 (organoids) with a maximum of 1×107 cells using RNAasy Mini Kit (QIAGEN, Hilden, Germany, 74104,) following the manufacture instructions. The RNA was treated with RNase-free DNase I (QIAGEN, 79254) for 15 min at room temperature during the extraction. The concentration, quality, and integrity of the extracted RNA were evaluated using Nanodrop (Thermo Scientific, Waltham, MA; USA), an Agilent 2100 Bioanalyzer, and the Agilent RNA 6000 Nano kit (5067–1511, Agilent Technologies). More than 0.4 μg total RNA with high integrity (more than 6.8) and high purity (OD260/280=1.8–2.2 and OD260/230≥1.8) were collected and sent to Illumina NovaSeq 6000RNA sequencing by Novogene. RNA-Seq library preparation and next-generation sequencing: cDNA libraries with paired-end 150 bp enriched were prepared by Novogene. First, the mRNA was randomly fragmented and supplemented with oligo (dT) beads. Then cDNA synthesis was done using the random hexamers and reverse transcriptase. Second, second-strand synthesis was done using: a custom second-strand synthesis buffer from Illumina, deoxyribose nucleoside triphosphates (dNTPs), RNase H, and Escherichia coli polymerase I. The final obtained cDNA library was purified, terminally repaired, A-tailed ligated to sequencing adapters, size-selected, and PCR-enriched. Quantification of library concentration was performed using a Qubit 2.0 fluorometer. Library size was measured by Agilent 2100 bioanalyzer and was quantified by qPCR (library activity >2 nM). Libraries were sequenced on Illumina NovaSeq 6000 S4 flow cells (paired end, 150 bp).
Raw data were transformed into sequenced reads and recorded in a FASTQ file. FASTQ files were aligned to build 19 of the human genome provided by the Genome Reference Consortium (GRCh19), performed by Christian Hinze using TOPHAT2 aligner tool (Kim et al., 2013). Up to four mismatches with the reference genome were accepted. Raw counts were obtained using featureCounts (Liao et al., 2014). Mutation visualization in the knockout samples was performed using the Integrative Genomic Viewer tool (Robinson et al., 2011). For gene expression analysis reads were normalized to the sequence length and transcripts per million (TPM) values were calculated (Wagner et al., 2012). TPM values of the samples were used to plot heatmaps and for PCA analysis based on Pearson correlation, using R (R Development Core Team, 4.0.4)
500 ng of RNA was reverse transcribed using the RevertAid First Strand cDNA synthesis kit (Thermo Scientific, K1622) according to the manufacturer’s instructions. The qPCR was carried out using the Fast Universal SYBR Green Master Mix (ROX, Roche Diagnostics, 04 913 850 001,) according to the manufacturer’s instructions. For expression analysis, relative mRNA expression levels were normalized for GAPDH mRNA expression and calculated according to the ΔΔCt method. All primer pairs were designed using the free-online primer design tool Primer3, purchased at BioTeZ (Berlin, Germany), and sequences are shown in (Supplementary file 1). The statistical significance of differences between two groups (WT and KO) was analyzed using a two-sided Student’s t-test.
Single-cell RNA (scRNA) experiment
Cells isolation and preparation
The differentiated hiPSCs to intermediate mesodermal cells were collected at day 7 of the differentiation from two different experiments. The cells of the first experiment were derived from WT1 and KO1 differentiated hiPSCs, while the second experiment comprises the differentiated cells of WT2 and KO2 cells. The cells were washed twice with 1× DPBS and dissociated with Accumax for 7 min at 37 °C. Cells were centrifuged at 350× g for 5 min and resuspended in 1× DPBS. Then, cells were filtered with a 40 µm filter (Corning, 352340), counted (10,000 cells per sample), and checked for viability using Trypan blue staining.
Protein extraction and Immunoblotting
Protein extraction
Up to 1×106 hiPSCs per sample (WT1, WT2, WT3, KO1, KO2, and KO3) were washed with cold 1× PBS, then centrifuged at 3500 g for 5 min. Next, the cell pellet was resuspended in pre-ice cold 100 µl of radioimmunoprecipitation assay (RIPA) buffer (Sigma-Aldrich, R0278) supplemented with protease inhibitor (Roche, 11697498001) and maintained with constant agitation for 30 min at 4°C. Then the suspension was centrifuged at 4°C for 20 min at 12,000 rpm. The supernatant was collected as protein extract and quantified using BCA Protein Assay (Thermo Scientific, 23228).
Immunoblotting
30 µg of the extracted protein in RIPA buffer were mixed with 1RIPA reducing (10% b-mercaptoethanol) NuPAGE loading buffer (Life Technologies, Carlsbad, CA). After denaturation at 70°C for 10 min. The protein was loaded on a precast polyacrylamide NuPage 4–12% Bis-Tris protein gel (Invitrogen, Carlsbad, CA, USA) and 1× MOPS (1 M MOPS, 1 M TrisBase, 69.3 mM SDS, 20.5 mM EDTA) to be separated according to the length using SDS -PAGE (100 V, 200mA, 2 hr). Proteins were blotted on 0.45 µm pore size Immobilon-P polyvinylidene difluoride membrane (EMD Millipore, Billerica, MA; USA). The membrane was pre-activated for 20 s in methanol and equilibrated in 1× NuPage Transfer buffer (1.25 mM bicine, 1.25 mM BisTris, 0.05 mM EDTA, and 10% ethanol) for 30 min at RT. The membrane was blocked in 5% bovine serum albumin for 1 hr at RT and incubated overnight at 4°C with primary antibodies: Anti-CEP83 produced in rabbit (1:500, Sigma-Aldrich) and Anti-α-Tubulin produced in mouse (1:500, T9026, Sigma-Aldrich). The membrane was incubated for 1 hr at RT with horseradish peroxidase-conjugated secondary antibodies (Sigma-Aldrich, Saint Louis, MO, USA) with 1:2000 dilution. Chemiluminescent reagent (Super Signal–West Pico; Thermo Scientific, Waltham, MA; USA) was used to detect the proteins. The spectra Multicolor Broad Range Protein Ladder (Thermo Fisher Scientific, USA) was used to evaluate the molecular weight of corresponding protein bands.
Histology and immunofluorescence staining
After organoid fixation in BD Cytofix buffer (554655, BD Biosciences) for 1 hr on ice, the organoids were gradually dehydrated in increasing ethanol concentrations for 15 min each. Then organoids were cleared in xylene for three times 20 min each. After infiltration with melted paraffin at 65°C three times for 30 min each, the organoids were embedded in paraffin and processed in 3.5-µm thick sections using a HM 355 S microtome. The sections were deparaffinized, dehydrated, and stained with hematoxylin (Sigma-Aldrich, Saint Louis, MO) for 3 min and in 1% eosin (Sigma-Aldrich) for 2 min. The sections were mounted using Kaiser’s glycerol gelatin-based mounting medium. Images were captured with a Leica CTR 6000 microscope (Leica Biosystems, Wetzlar, Germany).
For immunostaining, cultured cells (D7) and organoids (D25) were fixed with BD Cytofix for 10 min on ice. Then cells were permeabilized with BD Perm/Wash (554723, BD Biosciences), twice, 15 min per each. Cells were blocked with a blocking solution (1% BSA + 0.3% triton-X-100 in 1× DPBS) for 2 hr at RT or overnight at 4°C. Cells were incubated overnight at 4°C with primary antibodies (Supplementary file 2). Cells were then washed twice (10 min each) and incubated with fluorescence-labeled secondary antibodies with 1:500 dilution including Cy3, Cy5, Alexa488, and Alexa647 (Jackson ImmunoResearch, Newmarket, UK) and Cy3 Streptavidin (Vector lab, Burlingame, USA) overnight at 4°C. DAPI was then used for nuclear staining (Cell signaling Technology, Danvers, MA, USA) with 1:300,000 dilution for 1 hour at RT. Finally, cells were mounted with Dako fluorescent mounting medium (Agilent Technologies). The images were taken using a SP8 confocal microscope (Carl Zeiss GmbH, Oberkochen, Germany). All the quantitative analyses of the taken images were performed using ImageJ (1.48 v; National Institutes of Health, Bethesda, MD) software.
Data availability
All data supporting the findings of this study are available within the article and its supplementary files. Source data files have been provided for Figures 1–6. Sequencing data have been deposited in GEO under accession code GSE205978.
-
NCBI Gene Expression OmnibusID GSE205978. The centrosomal protein 83 (CEP83) regulates human pluripotent stem cell differentiation toward the kidney lineage.
-
NCBI Gene Expression OmnibusID GSE75711. Parallel generation of easily selectable multiple nephronal cell types from human pluripotent stem cells.
References
-
Human limb abnormalities caused by disruption of Hedgehog signalingTrends in Genetics 28:364–373.https://doi.org/10.1016/j.tig.2012.03.012
-
BookThe AGT Cytogenetics Laboratory ManualHoboken, NJ, USA: John Wiley & Sons.https://doi.org/10.1002/9781119061199
-
Dissecting contiguous gene defects: Tbx1Current Opinion in Genetics & Development 15:279–284.https://doi.org/10.1016/j.gde.2005.03.001
-
Ureteric bud cells secrete multiple factors, including bfgf, which rescue renal progenitors from apoptosisThe American Journal of Physiology 273:F757–F767.https://doi.org/10.1152/ajprenal.1997.273.5.F757
-
Brief report: autosomal dominant familial hypoparathyroidism, sensorineural deafness, and renal dysplasiaThe New England Journal of Medicine 327:1069–1074.https://doi.org/10.1056/NEJM199210083271506
-
Nephric lineage specification by Pax2 and Pax8Genes & Development 16:2958–2970.https://doi.org/10.1101/gad.240102
-
Tbx1 regulates proliferation and differentiation of multipotent heart progenitorsCirculation Research 105:842–851.https://doi.org/10.1161/CIRCRESAHA.109.200295
-
CRISPOR: intuitive guide selection for CRISPR/Cas9 genome editing experiments and screensNucleic Acids Research 46:W242–W245.https://doi.org/10.1093/nar/gky354
-
Assembly of centrosomal proteins and microtubule organization depends on PCM-1The Journal of Cell Biology 159:255–266.https://doi.org/10.1083/jcb.200204023
-
Turning mesoderm into kidneySeminars in Cell & Developmental Biology 91:86–93.https://doi.org/10.1016/j.semcdb.2018.08.016
-
Ift25 is not a cystic kidney disease gene but is required for early steps of kidney developmentMechanisms of Development 151:10–17.https://doi.org/10.1016/j.mod.2018.04.001
-
Centrosome amplification disrupts renal development and causes cystogenesisThe Journal of Cell Biology 217:2485–2501.https://doi.org/10.1083/jcb.201710019
-
Mutations of CEP83 cause infantile nephronophthisis and intellectual disabilityAmerican Journal of Human Genetics 94:905–914.https://doi.org/10.1016/j.ajhg.2014.05.002
-
Diagnostic exome sequencing in 266 Dutch patients with visual impairmentEuropean Journal of Human Genetics 25:591–599.https://doi.org/10.1038/ejhg.2017.9
-
Osr1 functions downstream of hedgehog pathway to regulate foregut developmentDevelopmental Biology 427:72–83.https://doi.org/10.1016/j.ydbio.2017.05.005
-
Parallel generation of easily selectable multiple nephronal cell types from human pluripotent stem cellsCellular and Molecular Life Sciences 76:179–192.https://doi.org/10.1007/s00018-018-2929-2
-
Sonic hedgehog signaling in the lung from development to diseaseAmerican Journal of Respiratory Cell and Molecular Biology 52:1–13.https://doi.org/10.1165/rcmb.2014-0132TR
-
A ciliopathy complex builds distal appendages to initiate ciliogenesisThe Journal of Cell Biology 220:e202011133.https://doi.org/10.1083/jcb.202011133
-
PKD1-dependent renal cystogenesis in human induced pluripotent stem cell-derived ureteric bud/collecting duct organoidsJournal of the American Society of Nephrology 31:2355–2371.https://doi.org/10.1681/ASN.2020030378
-
LRRC45 contributes to early steps of axoneme extensionJournal of Cell Science 131:jcs223594.https://doi.org/10.1242/jcs.223594
-
Single nucleotide polymorphism arrays: a decade of biological, computational and technological advancesNucleic Acids Research 37:4181–4193.https://doi.org/10.1093/nar/gkp552
-
Gata-4/5/6, a subfamily of three transcription factors transcribed in developing heart and gutThe Journal of Biological Chemistry 269:23177–23184.https://doi.org/10.1016/S0021-9258(17)31636-8
-
CXCL12-CXCR4 axis in angiogenesis, metastasis and stem cell mobilizationCurrent Pharmaceutical Design 16:3903–3920.https://doi.org/10.2174/138161210794455003
-
Phosphorylation of CEP83 by TTBK2 is necessary for cilia initiationThe Journal of Cell Biology 218:3489–3505.https://doi.org/10.1083/jcb.201811142
-
The role of centrosome distal appendage proteins (daps) in nephronophthisis and ciliogenesisInternational Journal of Molecular Sciences 22:12253.https://doi.org/10.3390/ijms222212253
-
Nephron organoids derived from human pluripotent stem cells model kidney development and injuryNature Biotechnology 33:1193–1200.https://doi.org/10.1038/nbt.3392
-
Atlas of Hox gene expression in the developing kidneyDevelopmental Dynamics 229:771–779.https://doi.org/10.1002/dvdy.10474
-
Kif14 controls ciliogenesis via regulation of Aurora A and is important for hedgehog signalingThe Journal of Cell Biology 219:e201904107.https://doi.org/10.1083/jcb.201904107
-
Gata transcription factors in the developing and adult heartCardiovascular Research 63:196–207.https://doi.org/10.1016/j.cardiores.2004.03.025
-
Loss of function IFT27 variants associated with an unclassified lethal fetal ciliopathy with renal agenesisAmerican Journal of Medical Genetics. Part A 176:1610–1613.https://doi.org/10.1002/ajmg.a.38685
-
Gata5 is required for the development of the heart and endoderm in zebrafishGenes & Development 13:2983–2995.https://doi.org/10.1101/gad.13.22.2983
-
Eya 1 acts as a critical regulator for specifying the metanephric mesenchymeDevelopmental Biology 284:323–336.https://doi.org/10.1016/j.ydbio.2005.05.029
-
Centriole distal appendages promote membrane docking, leading to cilia initiationGenes & Development 27:163–168.https://doi.org/10.1101/gad.207043.112
-
Hand2 function in second heart field progenitors is essential for cardiogenesisDevelopmental Biology 351:62–69.https://doi.org/10.1016/j.ydbio.2010.12.023
-
Beyond nephronophthisis: retinal dystrophy in the absence of kidney dysfunction in childhood expands the clinical spectrum of CEP83 deficiencyAmerican Journal of Medical Genetics. Part A 185:2204–2210.https://doi.org/10.1002/ajmg.a.62225
-
The sonic hedgehog-patched-gli pathway in human development and diseaseAmerican Journal of Human Genetics 67:1047–1054.https://doi.org/10.1016/S0002-9297(07)62934-6
-
Bor and bo syndromes are allelic defects of EYA1European Journal of Human Genetics 5:242–246.
-
Measurement of mRNA abundance using RNA-Seq data: RPKM measure is inconsistent among samplesTheory in Biosciences = Theorie in Den Biowissenschaften 131:281–285.https://doi.org/10.1007/s12064-012-0162-3
-
Hox11 paralogous genes are essential for metanephric kidney inductionGenes & Development 16:1423–1432.https://doi.org/10.1101/gad.993302
Article and author information
Author details
Funding
Deutsche Forschungsgemeinschaft (GRK 2318)
- Kai M Schmidt-Ott
Deutsche Forschungsgemeinschaft (FOR 2841)
- Kai M Schmidt-Ott
Deutsche Forschungsgemeinschaft (SFB 1365)
- Kai M Schmidt-Ott
National Science Foundation
- Jelena Kresoja
University of Colorado
- Christian Mosimann
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Tatjana Luganskaja for her excellent technical support. This work was supported by grants to KMS-O from the Deutsche Forschungsgemeinschaft (DFG; SFB 1365, GRK 2318, and FOR 2841), by stipends to FM by the Egyptian government, by the Urological Research Foundation (Berlin), a Swiss National Science Foundation postdoctoral fellowship to JK-R, and the University of Colorado School of Medicine, Anschutz Medical Campus, and the Children’s Hospital Colorado Foundation to CM.
Copyright
© 2022, Mansour et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 924
- views
-
- 226
- downloads
-
- 3
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cell Biology
Anionic lipid molecules, including phosphatidylinositol-4,5-bisphosphate (PI(4,5)P2), are implicated in the regulation of epidermal growth factor receptor (EGFR). However, the role of the spatiotemporal dynamics of PI(4,5)P2 in the regulation of EGFR activity in living cells is not fully understood, as it is difficult to visualize the local lipid domains around EGFR. Here, we visualized both EGFR and PI(4,5)P2 nanodomains in the plasma membrane of HeLa cells using super-resolution single-molecule microscopy. The EGFR and PI(4,5)P2 nanodomains aggregated before stimulation with epidermal growth factor (EGF) through transient visits of EGFR to the PI(4,5)P2 nanodomains. The degree of coaggregation decreased after EGF stimulation and depended on phospholipase Cγ, the EGFR effector hydrolyzing PI(4,5)P2. Artificial reduction in the PI(4,5)P2 content of the plasma membrane reduced both the dimerization and autophosphorylation of EGFR after stimulation with EGF. Inhibition of PI(4,5)P2 hydrolysis after EGF stimulation decreased phosphorylation of EGFR-Thr654. Thus, EGFR kinase activity and the density of PI(4,5)P2 around EGFR molecules were found to be mutually regulated.
-
- Cell Biology
Classical G-protein-coupled receptor (GPCR) signaling takes place in response to extracellular stimuli and involves receptors and heterotrimeric G proteins located at the plasma membrane. It has recently been established that GPCR signaling can also take place from intracellular membrane compartments, including endosomes that contain internalized receptors and ligands. While the mechanisms of GPCR endocytosis are well understood, it is not clear how well internalized receptors are supplied with G proteins. To address this gap, we use gene editing, confocal microscopy, and bioluminescence resonance energy transfer to study the distribution and trafficking of endogenous G proteins. We show here that constitutive endocytosis is sufficient to supply newly internalized endocytic vesicles with 20–30% of the G protein density found at the plasma membrane. We find that G proteins are present on early, late, and recycling endosomes, are abundant on lysosomes, but are virtually undetectable on the endoplasmic reticulum, mitochondria, and the medial-trans Golgi apparatus. Receptor activation does not change heterotrimer abundance on endosomes. Our findings provide a subcellular map of endogenous G protein distribution, suggest that G proteins may be partially excluded from nascent endocytic vesicles, and are likely to have implications for GPCR signaling from endosomes and other intracellular compartments.