Discovery of novel determinants of endothelial lineage using chimeric heterokaryons
Abstract
We wish to identify determinants of endothelial lineage. Murine embryonic stem cells (mESC) were fused with human endothelial cells in stable, non-dividing, heterokaryons. Using RNA-seq, it is possible to discriminate between human and mouse transcripts in these chimeric heterokaryons. We observed a temporal pattern of gene expression in the ESCs of the heterokaryons that recapitulated ontogeny, with early mesodermal factors being expressed before mature endothelial genes. A set of transcriptional factors not known to be involved in endothelial development was upregulated, one of which was POU class 3 homeobox 2 (Pou3f2). We confirmed its importance in differentiation to endothelial lineage via loss- and gain-of-function (LOF and GOF). Its role in vascular development was validated in zebrafish embryos using morpholino oligonucleotides. These studies provide a systematic and mechanistic approach for identifying key regulators in directed differentiation of pluripotent stem cells to somatic cell lineages.
https://doi.org/10.7554/eLife.23588.001eLife digest
Endothelial cells form the inner surface of blood vessels, acting like a non-stick coating. In addition to making substances that keep blood from sticking to the vessel wall, endothelial cells generate compounds that relax the vessel, and prevent it from thickening. Endothelial cells also form capillaries, the smallest vessels that provide oxygen and nutrients for all tissues.
A regenerating organ, or a bioengineered tissue, requires a system of capillaries and other microvessels. Thus, regenerative medicine could benefit from a knowledge of how to generate endothelial cells from pluripotent stem cells – cells that can “differentiate” to form almost any type of cell in the body.
Wong, Matrone et al. have now used a cell fusion model (named heterokaryon) to track the changes in gene expression that occur as a pluripotent stem cell differentiates to ultimately become an endothelial cell. In this model, mouse embryonic stem cells (ESCs) are fused to human endothelial cells. Over time the human endothelial cells drive gene expression in the ESCs toward that of endothelial cells.
Wong, Matrone et al. discovered changes in gene expression in many genes that have not previously been described as involved in the differentiation of endothelial cells. When one of these genes – named Pou3f2 – was inactivated in ESCs, they could not be differentiated into endothelial cells. The absence of Pou3f2 also drastically impaired how blood vessels developed in zebrafish embryos.
Thus the heterokaryon model can generate important information regarding the dynamic changes in gene expression that occur as a pluripotent cell differentiates to become an endothelial cell. This model may also be useful for discovering other genes that control the differentiation of other cell types.
https://doi.org/10.7554/eLife.23588.002Introduction
Our understanding of the genetic and epigenetic processes governing endothelial development and differentiation is limited (Yan et al., 2010; De Val and Black, 2009). Accordingly, our methodologies for obtaining endothelial cells from pluripotent stem cells are empirically driven and suboptimal (Choi et al., 2009; James et al., 2010; Huang et al., 2010a, 2010b; Wong et al., 2012). There is unexplained inconsistency in the yield of iPSC-ECs; in the stability of their phenotype; and in the fidelity of differentiation (in terms of replicating the epigenetic and genetic profile of a mature endothelial cell). Furthermore, our ability to efficiently generate specific endothelial subtypes (e.g. arterial, venous, lymphatic) is poor. Thus, a systematic approach is needed to more completely define the genetic and epigenetic programs required for differentiating pluripotent stem cells to the endothelial phenotype. Here, we propose an unbiased systematic approach to discover determinants of differentiation. We use interspecies heterokaryons, RNA sequencing and third-generation bioinformatics to discover novel candidate genes critical for proper endothelial differentiation and specification.
Results
Interspecies heterokaryons as a discovery tool
To discover new genes involved in endothelial specification, we made heterokaryons consisting of human endothelial cells (hEC) and murine embryonic stem cells (mESC) (Figure 1a–c), which expressed cell surface markers and characteristics of both cell types. We hypothesized that the factors that are actively maintaining endothelial phenotype (transcription factors, epigenetic modifiers and non-coding RNA etc) would act on the pluripotent stem cell nuclei to induce expression of key determinants of endothelial lineage. We reasoned that we could use RNA seq to monitor global changes in the transcriptome of the pluripotent nucleus as it is reprogrammed in the heterokaryon toward an endothelial fate. In 95% of cases, the species-specific nucleotide differences between the mouse and human transcripts would permit us to differentiate between reads of murine versus human transcripts when the sequences were aligned to their respective genomes.
Optimization and testing of the heterokaryon system
Reprogramming of the cell population is synchronized upon the addition of the fusagen. Since there is no nuclear fusion, chromosome rearrangement, or chromosome loss in the heterokaryons (Bhutani et al., 2010), we reasoned that this synchronization would permit us to study the temporal sequence of reprogramming to endothelial lineage using RNA seq. We optimized the cell fusion strategy using the fusagen HVJ (Sendai virus) envelope protein. By skewing the ratio of the input cells so that endothelial cells outnumbered pluripotent stem cells in the multinucleate heterokaryon, we forced reprogramming of the pluripotent stem cell nuclei toward an endothelial phenotype. To confirm that the system was working as anticipated, we assessed the expression of the mESC mRNA transcripts using murine-specific primers. In mESCs which were fused with human endothelial cells, we observed upregulation of murine endothelial genes including Kdr, Tie2, Cdh5 and Vwf (Figure 1d–g), and transcription factors involved in endothelial development such as Etv2, Ets1 and Tal1 (Figure 1l–n). Intriguingly, the expression of these genes seemed to mirror ontogeny, in that genes involved early in mesoderm specification (e.g. Kdr) were expressed earlier in the heterokaryon, whereas genes that are more specific to hemato-endothelial lineage (e.g. Von Willebrand factor) were expressed later. In parallel to the upregulation of genes involved in endothelial development, we observed down-regulation of genes encoding pluripotent factors (Oct4, Sox2 and Nanog) (Figure 1o).
As a complementary approach, we generated heterokaryons consisting of murine endothelial cells (mEC) and human-induced pluripotent stem cells (hiPSC). We fused hiPSCs with mECs, and the expression of human mRNA transcripts was assessed using human specific primers. We observed in the pluripotent cells an upregulation of human genes including Kdr, Tie2, Cdh5 and Vwf (Figure 1h–k), and transcription factors involved in endothelial development such as Etv2, Ets1 and Tal1 (Figure 1p–r). Again, we observed a temporal sequence that mirrored ontogeny, with a parallel downregulation of genes encoding pluripotent factors (Oct4, Sox2 and Nanog) (Figure 1s). These results demonstrate rapid induction of endothelial genes in the heterokaryon which appears to recapitulate endothelial development.
RNA seq for novel EC determinants
Having observed that the heterokaryon system seemed to mirror endothelial ontogeny, we applied RNA-seq to identify novel determinants of endothelial lineage (Figure 2). Heterokaryons were double-positive cells (GFP+ for mESC and CellTracker Red+ for hECs) and were harvested at 6–24 hr post-fusion. These time points were chosen based on our qPCR data, which showed rapid sequential activation of endothelial genes as early as 6 hr post-fusion. Total RNA (~2 ug) was isolated from heterokaryons or co-culture controls, and prepared for RNA-seq studies. We also did RNA-Seq for parental mESC; and for mESC that were exposed to an endothelial differentiation protocol for 4 or 8 days. We defined genes that are up- or down-regulated in each sample relative to the parental mESC. We were most interested in novel genes that encode transcriptional factors that are not known to be involved in mesodermal lineages or endothelial specification. We postulated that the heterokaryon system would be particularly useful to detect transcription factors that are transiently expressed during endothelial differentiation (Figure 2a). The RNAseq analysis of the heterokaryon samples (mESC/hEC) revealed a set of 1297 genes in the mESC that were upregulated as well as a set of 795 genes that were downregulated after fusion with hEC (Figure 2b). The upregulated genes included known endothelial related genes, for example, Kdr and Vwf, as well as hemangioblast markers, that is, angiotensin-converting enzyme (ACE, CD143). Conversely, murine pluripotency markers were down-regulated. We identified some novel transcriptional factors (such as Tbx1, Cebpd, Nrarp, Pou3f2, Maf, Tbx20, Tigd5, Pdn and Batf2), and targeted one of these (Pou3f2) for further study.
Unbiased hierarchical clustering analysis indicated that the transcriptome of the mESC in the heterokaryon is closer to that of mESCs differentiated toward EC lineage, than to the parental mESCs (Figure 2c). By contrast, the transcriptome of mESC co-cultured (but not fused) with hECs bears a closer relationship to the transcriptome of the parental mESCs. Global gene function enrichment analysis indicates that differentially up- or down-regulated genes in the mESC in heterokaryons are similar to those in mESC exposed to the endothelial differentiation protocol (Figure 2d). In particular, the upregulated genes in the heterokaryons tend to be involved in transcription regulation, RNA alternative splicing, DNA binding, embryonic organ development, differentiation and regulation of cell proliferation (Figure 2e), whereas downregulated genes are implicated in ATP and ribonucleotide binding (Figure 2f). These results suggest that the heterokaryon system may serve as an effective model for cell differentiation.
Validation of Pou3f2 in differentiation to EC lineage
To validate our approach in discovery of novel determinants of cell lineage, we focused on Pou3f2, which was identified as a candidate EC transcription factor in the RNA seq studies. Accordingly, we examined the temporal sequence of gene expression of Pou3f2 in the heterokaryon system (Figure 3a). The expression of Pou3f2 was up-regulated at 6-hr post-fusion and its expression then declined over the ensuing 48 hr. In addition, we examined Pou3f2 expression during the differentiation of murine ESC to endothelial cells using our standard differentiation protocol (Figure 3b). Intriguingly, the pattern of expression was similar, albeit with a slower time course. Specifically, the expression of Pou3f2 increased from Day 3 to Day 7 and then declined in the later phase (11 days) of the endothelial differentiation protocol. Notably, the expression of other mesodermal and EC genes followed a similar pattern, in that their regulation was accelerated when ESCs were placed into heterokaryons versus the standard endothelial differentiation protocol. This observation is consistent with the notion that current differentiation protocols could be considerably improved.
To further characterize the role of Pou3f2 in the differentiation of mESC into EC, we performed lentiviral shRNA knockdown (KD) of Pou3f2 in mESCs, and subjected the KD mESCs to the endothelial cell differentiation protocol. We found that KD of Pou3f2 reduced the expression of endothelial genes including Kdr, Tie2, Nos3, Cdh5 and Vwf at Day 8 of the differentiation protocol (Figure 3c). Similarly, the expression of endothelial transcription factors such as Ets1, Erg, Etv2 and Fli1 were also reduced in the Pou3f2 KD mESC (Figure 3d). The KD of Pou3f2 in mESC did not affect the expression of mesoderm- (Bmp4, T), endoderm- (Cxcr4, Gata4) or ectoderm- (Pax6 and Nestin) related genes (Figure 3e). Notably, the generation of mESC-derived CD31+ and CD144+ cells were reduced by over 50% in Pou3f2 KD group (Figure 3g). Furthermore, Pou3f2 KD mESC-derived endothelial cells manifested poor network formation on matrigel (Figure 3h). To summarize, Pou3f2 seems to be necessary for the full expression of genes known to be involved in endothelial development, and for the efficient generation of fully functioning endothelial cells. Amongst the factors released from endothelial cells in the heterokaryons that could control Pou3f2 expression there is Wnt and β-catenin. The Wnt/B-catenin signalling pathway is highly conserved and regulates vascular cell fate and development through Dll4/Notch signalling (Corada et al., 2010). The promoter for the Pou3f2 gene is a direct target for β-catenin/Lef1 (Goodall et al., 2004a). Endothelial cells in the heterokaryon might also contribute phosphatidylinositol 3-kinase to activate Pou3f2. The PI3K pathway mediates angiogenesis and the expression of growth factors in endothelial cells (Jiang et al., 2000) and also regulates Pou3f2 in melanoma cells (Bonvin et al., 2012).
Importance of Pou3f2 in EC development
Pou3f2 promotes neurogenesis (Jaegle et al., 2003; Castro et al., 2006; Dominguez et al., 2013; Sugitani et al., 2002). Specifically, it activates the Notch ligand Delta1, synergistically with Mash1, to maintain a subset of neural progenitors in an undifferentiated state (Castro et al., 2006), whereas it suppresses the Notch effector Hes5 (Dominguez et al., 2013) that negatively regulate transcription of neurogenesis-promoting genes Neuregulins (Imayoshi et al., 2008). In human melanoma spheres and tumor xenograft, Pou3f2 is proposed to induce the Notch pathway (Thurber et al., 2011). In our studies, the importance of Pou3f2 in endothelial development in vivo was assessed using the zebrafish model. The availability of transgenics expressing endothelial-specific fluorescent reporters, for example Tg(fli1:EGFP)y1, combined with the transparency of the embryo, facilitate visualization of vascular development and blood flow in real time (Baldessari and Mione, 2008; Ellertsdóttir et al., 2010; Holden et al., 2011; Kamei et al., 2010).
In situ hybridization for Pou3f2 showed that this transcription factor is expressed in most embryo structures, although the staining is particularly evident in the brain region (Figure 4—figure supplement 1). Because endothelial and hematopoietic cells share a common progenitor, we assessed the expression of Pou3f2 in both endothelium and hematopoietic cells. To do so, we purified GFP+ cells from tg(fli1:EGFP)y1 or tg(cmyb:GFP) embryos. We found expression of Pou3f2 in both type of cells, although more pronounced in the GFP+ cells from tg(fli1:EGFP)y1 embryos (Figure 4—figure supplement 2). Injection of morpholino targeting Pou3f2 (the relevant zebrafish analogue) resulted in 100% embryo death, as embryo could not reach 50% epiboly stage, at 5.3 hr post-fertilization. Both Pou3f2-targeted morpholinos caused the same effects. Therefore, we used light-activated cyclic caged morpholinos against Pou3f2 to study the role of this gene in vascular development (Figure 4a–d). When the Pou3f2-targeted morpholino was activated by UV light at six hpf the embryo survival was around 80%. We observed a general phenotype characterized by curved body and curly tail, and specific dysmorphogenesis of the vascular system in most of the tg(fli1:EGFP)y1 embryos (Figure 4a). A reduced number of intersegmental vessels (ISVs) was observed (Figure 4b–c), which could be rescued with Pou3f2 mRNA. In situ hybridization showed that Fli1 and Kdr were downregulated following Pou3f2 KD, particularly evident in the intersegmental vessels, confirming the role for Pou3f2 in endothelial development (Figure 4d). These results were confirmed by real-time PCR (Figure 4—figure supplement 3), where Fli1 and Kdr were significantly reduced in Pou3f2 knockdown embryos compared to control. Effective knockdown by the morpholino in zebrafish embryos was documented by Western blotting using a Pou3f2-specific antibody (Figure 4e). In addition, there was a global reduction of endothelial cells in Pou3f2 KD animals. Specifically, the percentage of GFP+ cells isolated from the tg(fli1:EGFP)y1 zebrafish embryos at 24 hr post-fertilization was 6.7% of total cells in control compared to 3.7% following Pou3f2 KD (Figure 4f–g). Injection of Pou3f2-targeted caged morpholino in tg(cmyb:GFP) did not significantly reduce the percentage of GFP+ cells, suggesting that Pou3f2 is more critical for endothelial rather than hematopoietic development. Notably, caged MO activation at 24 hpf or later did not produce a significant phenotype, indicating that Pou3f2 is necessary early in endothelial development, but is not required for maintenance of the endothelial phenotype. This observation is consistent with the lack of expression of Pou3f2 in mature endothelial cells. Our studies suggest that Pou3f2 could act upstream of Notch during endothelial, as well as neuronal development. In this regard, Pou3f2 could have a role in the human vascular diseases caused by mutations in Notch family members, such as Alagille syndrome and cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL) (Loomes et al., 1999; Joutel et al., 1997). As Notch signaling is also required for arterio-venous differentiation during vascular development (Lawson et al., 2001), this suggests that Pou3f2 could be also upstream of the eph-ephrin system and play a role in arterio-venous specification. Pou3f2 gene expression seems to be controlled by MAPK (Goodall et al., 2004b) and Wnt/β-catenin (Goodall et al., 2004a) signaling, two key pathways linked to cell proliferation. Pou3f2 function can be also modulated by posttranslational modifications including sumoylation, ubiquitinylation, glycosylation and in particular phosphorylation (Kasibhatla et al., 1999; Augustijn et al., 2002; Diamond et al., 1999; Nieto et al., 2007).
-
Figure 4—source data 1
- https://doi.org/10.7554/eLife.23588.007
Role of Pou3f2 in human EC differentiation
To determine the role of Pou3f2 during endothelial differentiation from human iPSC, we examined the expression of Pou3f2 in the heterokaryon system consisting of hiPSCs and murine endothelial cells (Figure 5a). In addition, we examined the expression of Pou3f2 in the hiPSC during our endothelial differentiation protocol (Figure 5b). The expression patterns for Pou3f2 were similar to that which we observed in mESC (Figure 3a and b). Notably, we saw the same accelerated regulation of Pou3f2 in the hiPSCs in the heterokaryon by comparison to hiPSCs exposed to the endothelial differentiation protocol.
Next, we assessed the effects of shRNA-mediated KD of Pou3f2 on differentiation of endothelial cells from human iPSC. Pou3f2 shRNA significantly reduced the expression of Pou3f2 during the endothelial differentiation protocol (Figure 5c,d). Furthermore, Pou3f2 KD significantly inhibited the expression of endothelial-related genes during the course of the differentiation protocol, including Kdr, Tie2, CDdh5 Pecam1, Nos3 and Vwf (Figure 5e–j). Notably, Pou3f2 KD also reduced iPSC-EC generation by ~50% (Figure 6a–c). In addition, protein expression of CD31, CD144 and Vwf (Figure 6d) and formation of tubular networks on matrigel (Figure 6e) were reduced in ECs derived from Pou3f2 KD iPSCs. Furthermore, Pou3f2 KD iPSC-ECs exhibited impairment of other endothelial functions including nitric oxide generation (Figure 6f), and uptake of acetylated LDL (Figure 6g). These functional impairments were associated with reduced expression of EC markers including Kdr, Tie2, Nos3, CD31, Cdh5 and Vwf (Figure 6h). We have previously shown that endothelial cells derived from iPSCs (hiPSC-ECs) expressed various markers associated with arterial, venous and lymphatic ECs and thereby represent a heterogeneous population of ECs (Rufaihah et al., 2011). We found that ECs generated from Pou3f2 KD iPSCs resemble those generated from wild-type iPSCs with respect to venous (Ephb4 and Coup-TFII) and lymphatic markers (Pdpn and Lyve1) but have a significant reduction of arterial markers (Notch4, Efnb2 and Hey2) compared to scrambled control (Figure 6i).
Pou3f2 interacts with known endothelial promoters
To determine if Pou3f2 binds to the promoters of endothelial-related transcription factors, we performed ChIP-PCR. We observed that during differentiation of iPSCs to ECs, the binding of Pou3f2 to the promoters of endothelial related transcription factors including Ets1, Lmo2, Hey1 and Hey2 (Figure 7a–d) was inhibited in Pou3f2 KD cells, in association with reduced gene expression of these factors (Figure 7f–i). Finally, we observed that the generation of CD31+CD144+ cells was reduced in the Pou3f2 KD hiPSCs and could be rescued with modified mRNA encoding Pou3f2 (Figure 7j).
Discussion
Our current understanding of the genetic and epigenetic processes governing endothelial development and differentiation is limited. We lack comprehensive knowledge regarding all endothelial lineage factors and have sparse information regarding the magnitude and temporal sequence of their expression. In this paper, we find that the bi-species heterokaryons combined with RNAseq can provide new insights into determinants of endothelial lineage. Our work suggests that transcription factors and epigenetic machinery which actively maintain endothelial phenotype can also act on the pluripotent cell nucleus to recapitulate ontogeny. This system is likely to generate useful insights to improve the yield and fidelity of reprogramming to endothelial phenotype. A tangible and immediate outcome of this line of inquiry will be a more complete knowledge of the hierarchy of genes regulating differentiation to the EC lineage. Insights into these processes will be of general interest to investigators of vascular differentiation and development and may lead to new therapeutic targets for endothelial regeneration and the treatment of vascular diseases.
Finally, and perhaps most importantly, this model system should be amenable to discovery of novel determinants of other cell lineages. We believe that our studies provide proof-of-concept for using bi-species heterokaryon technology as a tool to elucidate novel genes regulating differentiation to any somatic cell. Our work opens a new vista of exploration for the broader community of scientists working in tissue regeneration, development, differentiation and the therapeutic applications of these insights.
Materials and methods
Cell culture
Request a detailed protocolHuman-induced progenitor stem cells (Takara Bio USA, Mountain View, CA) line authentication was achieved by genetic profiling using polymorphic short tandem repeat (STR profiling) loci. Our cell cultures were tested weekly for mycoplasma by real-time PCR approach and were mycoplasma-free. HiPSC lines were generated using retroviral factors encoding Oct4, Sox2, Cmyc and Klf4 in adult dermal fibroblasts. The hiPSCs were characterized for their pluripotency using PCR and IHC for known pluripotency markers, and were maintained in mTeSR1 (Stem Cell Technology, Vancouver, Canada). Murine ESCs (D3, ATCC, Manassus, VA) of the SV129 strain were cultured on gelatin-coated dish and maintained in ESGRO media plus GSK3β inhibitors. Human microvascular endothelial cells (HMVECs) (obtained from Lonza, Walkersville, MD) and murine endothelial cells (obtained from Applied Stemcell, Menlo Park, CA) were cultured in EC growth medium EGM-2 MV (cc-3162). Cells were used for all experiments at passage 6–8.
Heterokaryon formation and isolation
Request a detailed protocolOne day before cell fusion, 4 × 105 endothelial cells were seeded on one well of a 6-well dish. The endothelial cells were confluent on the day of cell fusion. On the same day, 1 hr prior to cell fusion, the endothelial cell medium was replaced with fresh EGM-2 MV medium supplemented with 1 μM Cell Tracker Red, then cells were incubated at 37°C in darkness for 30 min. Human iPSCs or murine ESCs labeled by transduction with retroviruses encoding GFP were rinsed with PBS followed by accutase treatment at 37°C for 5 min to dissociate the pluripotent stem cells into single-cell suspension. The cells were then collected in conical tubes after neutralization by MEF media containing 10% FBS. The cells were then counted with hemocytometer and 2 × 105 pluripotent stem cells were taken. The cells were then centrifuged at 200 x g for 5 min at 4°C, the supernatant was removed, and the cell pellet was resuspended in 25 μL ice-cold cell fusion buffer with 2.5 μL ice-cold HVJ-envelope fusagen. The reaction mixture was placed on ice for 5 min with regular agitation in 2.5 min apart. After 5 min, the cells were centrifuged again and the supernatant was discarded and 2 ml cell fusion buffer was added. The pluripotent cells were then plated onto the endothelial cells. The six-well dish was then centrifuged at 200 g for 5 min at 4°C. After centrifugation the dish was placed into a 37°C incubator to induce cell fusion. Twenty minutes later, the medium was removed and EGM-2 MV medium was added. For the Co-culture Control, the described procedure was the same except HVJ-enveloped fusagen was not added. The heterokaryons (double-positive cells) can be efficiently sorted by FACS. Heterokaryons (GFP+ and CellTracker Red+) were harvested by FACS at 6, 12, 24, 48 and 72 hr post-fusion.
Preparation of RNA-seq libraries and sequencing
Request a detailed protocolThe species-specific nucleotide differences between the mouse and human transcripts enable us to differentiate between reads of transcripts from the murine ESC versus those from the human EC when the sequences are aligned to their respective genomes. Heterokaryons (GFP+ and CellTracker Red+) were harvested by FACS at 6, 12 and 24 hr post-fusion and prepared for analysis by RNA-seq. Total RNA from heterokaryons were isolated. Human and mouse mRNA transcripts were isolated from the total RNA samples using polyA-based enrichment using oligo-dT magnetic beads. The majority of contaminating ribosomal RNA was eliminated by this approach. The resulting mRNA was fragmented, reverse transcribed to cDNA, ligated to adapters, and subject to brief PCR amplification in preparation of the Illumina library. The integrity and quality of RNA and complementary DNA were monitored using an Agilent Bioanalyzer 2100. The samples were sequenced using pair-end 100 base-pair reads. For the estimation of gene expression and data analysis, any remaining ribosomal reads were discarded, and the resulting murine and human transcripts were mapped to their respective genomes. Reads that map to both transcriptomes would be discarded and the RPKMs adjusted accordingly (discarded reads represent only 5% of the total reads; furthermore, virtually all genes have at least one unique read that is different between species, so that no gene is completely discarded).
RNA-seq read mapping and gene expression
Request a detailed protocolRNA-Seq reads were aligned to the mouse genome version mm9 using TopHat version 2.1.0. We use the full set of knownGene downloaded from the UCSC Genome browser (http://genome.ucsc.edu/cgi-bin/hgTables) as reference genes. RNA-Seq read counts for each gene in each sample was calculated using Cuffdiff function in Cufflinks version 2.2.1. The Cuffdiff also calculates fragment per kilobase per million reads (FPKM) for each gene. We further subject the reads counts to EdgeR version 3.12.0 for differential expression analysis, and define differential genes based on false discovery rate (FDR) cutoff 1e-5. We subject interesting gene groups to the DAVID website (https://david.ncifcrf.gov) for functional enrichment analysis. Enriched functional terms were defined based on Benjamini adjusted p value cut-off 0.05. Hierarchical clustering of gene expression heatmap was conducted using MEV based on Pearson correlation distance metric and the average linkage method.
RNA extraction and quantitative PCR
Request a detailed protocolUsing RNeasy Mini Kit (Qiagen, Chatsworth, CA), total RNA was extracted. The Quantitect reverse transcription kit (Qiagen) was used to generate cDNA and SYBR Green PCR kit (Invitrogen, Carlsbad, CA) was used for real-time qPCR with the QuantStudio 12 k Flex system (Applied Biosystems, Foster City, CA) following the manufacturer’s instructions. Genes were analyzed with the data normalized to Gapdh and expressed as relative fold changes using the ΔCt method of analysis.
Endothelial differentiation protocol
Request a detailed protocolMurine ESC-EC Differentiation: Endothelial differentiation of ESCs was carried out using the suspension culture approach with modifications. To initiate differentiation, ESCs were cultured in ultralow nonadhesive dishes to form embryoid body aggregates in a differentiation medium that consisted of α-Minimum Eagle’s Medium, 10% FBS, 1% penicillin/streptomycin, and 0.05 mmol/L β-mercaptoethanol (Sigma, St Louis, MO). After 4 days of suspension culture, the embryoid bodies were reattached onto 0.2% gelatin-coated dishes and cultured in differentiation medium. After 3 weeks of differentiation, the cells were purified by fluorescence-activated cell sorting (FACS) using anti-mouse vascular endothelial cadherin (VE-cadherin) antibody (Ab) (BD Biosciences, Bedford, CA).
Human iPSC-EC differentiation: Confluent cultures of hiPSCs were incubated with 1 mg/ml type IV collagenase for 10 min and transferred to ultra low attachment dishes containing differentiation media for 4 days to form embryoid bodies (EBs). The differentiation media used consisted of α-Minimum Eagle’s Medium, 20% fetal bovine serum, L-glutamine, β-mercaptoethanol (0.05 mmol/L) and 1% non-essential amino acids supplemented with bone morphogenetic protein-4 (BMP-4, 50 ng/ml, Peprotech) and vascular endothelial growth factor (VEGF-A, 50 ng/ml, Peprotech). The four-day EBs were reattached to gelatin-coated dishes in the presence of VEGF-A for another 10 days before purification.
Fluorescence-activated cell sorting (FACS)
Request a detailed protocolECs derived from pluripotent stem cells were purified using FACS. Cells were dissociated into single cells with Accutase (Invitrogen) for 5 min at 37°C, washed with 1x PBS containing 5% BSA and passed through a 70-μm cell strainer. Cells were then incubated with either Alexa Fluor 488-conjugated CD31 antibody (BD Bioscience, San jose, CA) or PE-conjugated CD144 antibody (BD Bioscience) for 30 min. Isotype-matched antibody served as negative control. The purified ESC- or iPSC-ECs were expanded in EGM-2 media.
Immunofluorescent imaging
Request a detailed protocolHuman iPSC-ECs were fixed with 4% paraformaldehyde, permeabilized with 0.1% Triton X-100, blocked with 1% normal goat serum and stained for anti-human CD31 (R and D Systems), anti-human CD144 (R &D Systems, Minneapolis, MN), anti-human von Willebrand factor (vWF, Abcam, Cambridge, UK) overnight at 4°C. After washes with PBS, the cells were treated with Alexa Fluor-488 or -594 secondary antibodies. Cell nuclei were stained with Hoechst 33342 (Sigma). Images were acquired on a confocal microscope (FV1000-IX81, Olympus, Tokyo, Japan).
Western blotting
Request a detailed protocolCells were homogenized with ice-cold RIPA lysis buffer containing 1 µg/mL leupeptin, 5 µg/mL aprotonin, 100 µg/mL PMSF, 1 mmol/L sodium orthovanadate, 1 mmol/L EDTA, 1 mmol/L EGTA, 1 mmol/L sodium fluoride and 2 µg/mL β-glycerolphosphate. The protein concentration was determined by Bradford method and aliquots of 20 µg of the total proteins were separated on 10% SDS-poly-acrylamide gel. Proteins were then transferred to immobilon-P polyvinylidene difluoride (PVDF) membrane (Millipore, Billerica, MA). Membranes were blocked with 5% non-fat milk in TBS-T and subsequently exposed to Pou3f2 primary antibody (Genetex, Irvine, CA) followed by HRP-conjugated secondary antibody and developed by chemiluminescence.
Functional assays
Request a detailed protocolUptake of Ac-LDL: was evaluated by incubating cells with ac-LDL-594 at 1:200 dilution for 5 hr before washing the cells with PBS and then measuring the mean fluorescence of the cells. Endothelial network formation,the ability of cells to form tube-like structures, was assessed in vitro by seeding 1.2 × 105 cells in wells coated with matrigel in the presence of EGM-2 media containing 50 ng/ml VEGF and incubated for 24 hr.
Nitric oxide production
Request a detailed protocolThe ability of the cells to produce NO was assessed by measuring the concentration of NO in the culture medium using the NO detection kit (Molecular Probe, Carlsbad, CA) according to the manufacturer’s instructions. The amount of nitrate was determined by converting it to nitrite, followed by the colorimetric determination of the total concentration of nitrite as a colored azo dye product of the Griess reaction that absorbed visible light at 540 nm using a microplate reader.
Chromatin immunoprecipitation and ChIP-qPCR
Request a detailed protocolhiPSCs were differentiated towards EC lineage and collected at Day 8. Samples were prepared by SimpleChIP enzymatic chromatin IP kit (Cell Signaling Technology). Chromatin immunoprecipitation was performed using human Pou3f2 antibody (Genetex), rabbit IgG (CST), histone H3 antibody (CST). DNA was purified using Nucleospin PCR clean-up kit (Macherey-Nagel, Bethlehem, PA) and used for quantitative PCR with primers against regions predicted within the promoter of Ets1, Lmo2, Hey1 and Hey2. Recovery of genomic DNA as the percentage input was calculated as the ratio of copy numbers in the immunoprecipitate to the input control.
Zebrafish aquaculture and husbandry
Request a detailed protocolAdult zebrafish (wild-type Wik and tg(fli1:EGFP)y1 strains) were acquired from the Zebrafish International Resource Center and raised according to standard procedures and kept at 28°C under a 14/10 hr light/dark cycle and fed with dry meal (Gemma Micro, Westbrook, ME) twice per day. Embryos used in these studies were obtained by natural matings and cultured in E3 embryo medium at 28.5°C. Animals were housed and all experiments were carried out in accordance with the recommendations of the Institutional Animal Care and Use Committee. All surgery procedures were performed under anesthesia with Tricaine 0.02 mg/ml.
Morpholinos and caged morpholinos injections
Request a detailed protocolPou3f2 KD in zebrafish was achieved using two different antisense morpholinos (Gene Tools, Oregon) targeting the Pou3f2 mRNA AUG translational start site with sequence: (Mo1) 5’-ATGATTGGATGCTGTAGTCGCCATG-3’, and (Mo2) 5’-CGGACTGATCGCTCCTATTAAAGGA-3’. As one control we used a 5-base pair mismatch MO: sequence 5’-ATcATTcGATcCTGTAcTCcCCATG-3’. To decipher the roles of Pou3f2 transcription factor in specific stages of endothelial development, we used Pou3f2-targeted caged morpholino (cMOs) (Shestopalov et al., 2007). This chemically modified morpholino allowed temporal gene silencing by using targeted UV illumination. An optimized dose of 0.5 ng/eggs (0.5 nL bolus) of Pou3f2 targeted morpholino was injected in each embryo at 1–2 cell stage, just below the cell mass.
cMOs photoactivation
Request a detailed protocolTo photoactivate the cMOs, injected zebrafish embryos were arrayed in an agarose microinjection template (560 μm x 960 μm wells), with the animal pole facing the light source. Then, the mercury lamp light was focused onto individual embryos for 10 s, using a Leica DM4500B epifluorescence microscope equipped with an A4 filtercube (Ex: 360 nm, 40 nm bandpass) and a 20 x/0.5 NA water-immersion objective. Individual embryos were irradiated at 6 or 24 hr post-fertilization.
Morpholino phenotype rescue by modified mRNA
Request a detailed protocolAs control, we also performed rescue experiments by co-injecting Pou3f2-targeted MO together with Pou3f2 modified mRNA into one-cell-stage embryos. The Pou3f2 modified mRNA version used, produced from the RNA Core available in our Institute, was modified at the 5′ untranslated region so that it was not recognized by the morpholino. An optimized dose of 300 pg was co-injected with the morpholino in rescue experiments.
Zebrafish imaging
Request a detailed protocolThe embryos were manually dechorionated at 24 or 48 hpf. Brightfield images were acquired using a Leica M205FA fluorescence stereoscope equipped with a Leica DFC500 digital camera. For immunofluorescence imaging, bright-field images of embryos were obtained with a Leica DM4500B compound microscope equipped with a 20 x/0.12 NA water-immersion objectives and a QImaging Retiga-SRV digital camera. Fluorescence images were obtained with the DM4500B/Retiga-SRV system equipped with a mercury lamp and GFP (Ex: 470 nm, 40 nm bandpass; Em: 525 nm, 50 nm bandpass) filter sets.
Western blot analysis of Pou3f2 in zebrafish
Request a detailed protocolThe embryos were de-yolked in TM1 buffer (100 mM NaCl, 5 mM KCl, 5 mM HEPES pH 7.0, 1% (w/v) PEG-200,000). Twenty de-yolked embryos from each experimental condition were homogenized in SDS-PAGE loading buffer (50 μH 7.0, 1% (w/mM 2-mercaptoethanol, 4% (w/v) glycerol, 100 mM DTT, 100 mM Tris-HCl, pH 6.8), vortexed, and heated to 95°C for 5 min. The resulting lysates were used for gel electrophoresis followed by blotting with Pou3f2 antibody (rabbit polyclonal, Abcam 137469). β-Tubulin (rabbit polyclonal, Abcam 6046) was used as loading control.
In situ hybridization in zebrafish
Request a detailed protocolThe preparation of sense (used as control) and antisense RNA probes for Kdr and Fli1 and in situ hybridization procedure were performed according to Thisse and Thisse (Thisse and Thisse, 2008).
FACS-based analysis of GFP+ cell from zebrafish
Request a detailed protocolCells were isolated according to Shestopalov et al. (Shestopalov et al., 2012) with modifications. Briefly, Tg(fli1:EGFP)y1 embryos at the appropriate developmental stage were dechorionated, transferred in an eppendorf tube with calcium-free Ringer’s solution (200 μl for 25–30 embryos; 116 mM NaCl, 2.6 mM KCl, 5 mM HEPES, pH 7.0) and dissociated with a 200 μl pipette tip. Then 1 ml solution of 1X PBS containing trypsin (0.25%, Gibco), 50 μg collagenase P (Roche, Indianapolis, IN) and 1 mM EDTA was added and samples were incubated for 30 min at 28.5°C with further pipetting every 5 min. Enzymatic processing was quenched with stop solution (200 μl; 1X PBS containing 30% (v/v) calf serum and 10 mM CaCl2), and cells were collected by centrifugation (400 g, 5 min, 4°C). After aspirating the supernatant, cells were resuspended in a chilled solution of DMEM containing 1% (v/v) calf serum, 0.8 mM CaCl2, 50 V ml−1 penicillin/streptomycin, centrifuged and resuspended in the same medium. The cell suspension was filtered through a 40 μm cell strainer (BD Biosciences) into FACS sample tubes. Cell suspensions were analyzed using a BD FACSAria. Wild-type zebrafish (Wik) was used as GFP negative control. Viable fluorescent single cells were identified as DAPI-negative. Cells viability was confirmed under fluorescence stereomicroscope (Leica M205) by using a Neubauer chamber.
Data access
Request a detailed protocolAll RNA-Seq data have been deposited to the GEO database by the accession number GSE84558.
Statistical analyses
Request a detailed protocolStatistical analysis was performed with SPSS software (SPSS Inc., Chicago, IL, USA). Results were expressed as mean ± SEM. The Shapiro-Wilk test was used to confirm the null hypothesis that the data follow a normal distribution. Statistical comparisons were performed via Student t-test for two groups and via one-way ANOVA test for multiple groups. Bonferroni corrections test was applied for multiple comparisons. p<0.05 was considered significant.
Data availability
-
Discovery of Novel Determinants of Endothelial Lineage: Insights from Chimeric HeterokaryonsPublicly available at the NCBI Gene Expression Omnibus (accession no: GSE84558).
References
-
How to create the vascular tree? (Latest) help from the zebrafishPharmacology & Therapeutics 118:206–230.https://doi.org/10.1016/j.pharmthera.2008.02.010
-
A phosphatidylinositol 3-kinase-Pax3 axis regulates Brn-2 expression in melanomaMolecular and Cellular Biology 32:4674–4683.https://doi.org/10.1128/MCB.01067-12
-
Transcriptional control of endothelial cell developmentDevelopmental Cell 16:180–195.https://doi.org/10.1016/j.devcel.2009.01.014
-
Vascular morphogenesis in the zebrafish embryoDevelopmental Biology 341:56–65.https://doi.org/10.1016/j.ydbio.2009.10.035
-
Brn-2 expression controls melanoma proliferation and is directly regulated by beta-cateninMolecular and Cellular Biology 24:2915–2922.https://doi.org/10.1128/MCB.24.7.2915-2922.2004
-
The Brn-2 transcription factor links activated BRAF to melanoma proliferationMolecular and Cellular Biology 24:2923–2931.https://doi.org/10.1128/MCB.24.7.2923-2931.2004
-
Molecular control of vascular development in the zebrafishBirth Defects Research Part C: Embryo Today: Reviews 93:134–140.https://doi.org/10.1002/bdrc.20204
-
Embryonic stem cell-derived endothelial cells engraft into the ischemic hindlimb and restore perfusionArteriosclerosis, Thrombosis, and Vascular Biology 30:984–991.https://doi.org/10.1161/ATVBAHA.110.202796
-
Role of nitric oxide signaling in endothelial differentiation of embryonic stem cellsStem Cells and Development 19:1617–1626.https://doi.org/10.1089/scd.2009.0417
-
The POU proteins Brn-2 and Oct-6 share important functions in schwann cell developmentGenes & Development 17:1380–1391.https://doi.org/10.1101/gad.258203
-
Imaging blood vessels in the zebrafishMethods in Cell Biology 100:27–54.https://doi.org/10.1016/B978-0-12-384892-5.00002-5
-
Jun kinase phosphorylates and regulates the DNA binding activity of an octamer binding protein, T-cell factor beta1Molecular and Cellular Biology 19:2021–2031.https://doi.org/10.1128/MCB.19.3.2021
-
Notch signaling is required for arterial-venous differentiation during embryonic vascular developmentDevelopment 128:3675–3683.
-
Endothelial cells derived from human iPSCS increase capillary density and improve perfusion in a mouse model of peripheral arterial diseaseArteriosclerosis, Thrombosis, and Vascular Biology 31:e72–e79.https://doi.org/10.1161/ATVBAHA.111.230938
-
Light-controlled gene silencing in zebrafish embryosNature Chemical Biology 3:650–651.https://doi.org/10.1038/nchembio.2007.30
-
Spatiotemporal resolution of the ntla transcriptome in axial mesoderm developmentNature Chemical Biology 8:270–276.https://doi.org/10.1038/nchembio.772
-
Endothelial cells derived from nuclear reprogrammingCirculation Research 111:1363–1375.https://doi.org/10.1161/CIRCRESAHA.111.247213
-
Epigenetics of the vascular endotheliumJournal of Applied Physiology 109:916–926.https://doi.org/10.1152/japplphysiol.00131.2010
Article and author information
Author details
Funding
American Heart Association (13SDG15800004)
- Wing Tak Wong
National Institutes of Health (PCBC_JS_2012/1_02)
- Wing Tak Wong
National Institutes of Health (U01HL100397)
- John P Cooke
Cancer Prevention and Research Institute of Texas
- John P Cooke
National Institutes of Health (RC2HL103400)
- John P Cooke
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This work was supported by grants to Dr. Cooke from National Institutes of Health (U01HL100397, RC2HL103400) and Cancer Prevention and Research Institute of Texas (CPRIT). Dr. Wong was supported by a NIH PCBC Jump Start Award (PCBC_JS_2012/1_02) and American Heart Association Scientist Development Grant (AHASDG) (13SDG15800004).
Ethics
Animal experimentation: Zebrafish are kept according to the laboratory protocols described in "Zebrafish: A Practical Approach" (Oxford University Press, 2002). These protocols comply with the Guide for the Care and Use of Laboratory Animals, the American Association for the Accreditation of Laboratory Animal Care (AAALAC) standards, and the regulations set forth in the Animals Welfare Act (P.L. 89-544, as amended by P.L. 91-579 and P.L. 94-279). Veterinary care is provided on a 24 hours basis, including weekends and holidays, by a staff of veterinarians with specialties in laboratory animal medicine and anesthesiology, and licensed animal health technicians. Training classes are offered. All veterinary care is provided by Houston Methodist Research Institute, which is fully accredited by AAALAC (ID A4555-01) and holds an approved NIH Assurance and USDA License (start date 03/08/2013). Support includes quarantine rooms, sterile operating rooms, post-surgical recovery rooms, radiology and diagnostic laboratory services. All surgery procedures were performed under anesthesia with Tricaine 0.02 mg/ml.
Copyright
© 2017, Wong 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
-
- 1,486
- views
-
- 262
- downloads
-
- 7
- 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
- Chromosomes and Gene Expression
During oncogene-induced senescence there are striking changes in the organisation of heterochromatin in the nucleus. This is accompanied by activation of a pro-inflammatory gene expression programme – the senescence-associated secretory phenotype (SASP) – driven by transcription factors such as NF-κB. The relationship between heterochromatin re-organisation and the SASP has been unclear. Here, we show that TPR, a protein of the nuclear pore complex basket required for heterochromatin re-organisation during senescence, is also required for the very early activation of NF-κB signalling during the stress-response phase of oncogene-induced senescence. This is prior to activation of the SASP and occurs without affecting NF-κB nuclear import. We show that TPR is required for the activation of innate immune signalling at these early stages of senescence and we link this to the formation of heterochromatin-enriched cytoplasmic chromatin fragments thought to bleb off from the nuclear periphery. We show that HMGA1 is also required for cytoplasmic chromatin fragment formation. Together these data suggest that re-organisation of heterochromatin is involved in altered structural integrity of the nuclear periphery during senescence, and that this can lead to activation of cytoplasmic nucleic acid sensing, NF-κB signalling, and activation of the SASP.
-
- Cell Biology
Defining the cellular factors that drive growth rate and proteome composition is essential for understanding and manipulating cellular systems. In bacteria, ribosome concentration is known to be a constraining factor of cell growth rate, while gene concentration is usually assumed not to be limiting. Here, using single-molecule tracking, quantitative single-cell microscopy, and modeling, we show that genome dilution in Escherichia coli cells arrested for DNA replication limits total RNA polymerase activity within physiological cell sizes across tested nutrient conditions. This rapid-onset limitation on bulk transcription results in sub-linear scaling of total active ribosomes with cell size and sub-exponential growth. Such downstream effects on bulk translation and cell growth are near-immediately detectable in a nutrient-rich medium, but delayed in nutrient-poor conditions, presumably due to cellular buffering activities. RNA sequencing and tandem-mass-tag mass spectrometry experiments further reveal that genome dilution remodels the relative abundance of mRNAs and proteins with cell size at a global level. Altogether, our findings indicate that chromosome concentration is a limiting factor of transcription and a global modulator of the transcriptome and proteome composition in E. coli. Experiments in Caulobacter crescentus and comparison with eukaryotic cell studies identify broadly conserved DNA concentration-dependent scaling principles of gene expression.