Genome-Wide Transcriptional Profiling of Human Embryonic Stem Cells Differentiat



免疫细胞治疗专区	欢迎关注干细胞微信公众号

12 3 4 5 6 7 8 9 10 ... 24 下一页

返回列表

查看: 394227\|回复: 234	go Genome-Wide Transcriptional Profiling of Human Embryonic Stem Cells Differentiat

江边孤钓

版主

Rank: 7 Rank: 7 Rank: 7

积分: 0
威望: 0
包包: 483

发表于 2009-3-5 00:03 |显示全部帖子

作者：Abdelaziz Beqqalia, Jantine Klootsa, Dorien Ward-van Oostwaarda,b, Christine Mummerya,b, Robert Passiera,b
      o1 `' e+ b" F  R% i% y' K, y
   , t/ Q6 [1 I4 g  l; A( k$ w
      ) l/ F0 ?6 f$ |# l7 u, V
                     $ b6 o  H; B+ O* m0 M4 B5 O
         8 q: e8 |4 C2 M( [0 b
         4 X7 q, s; f. _5 U* `

         % \2 w- c: K6 e, u4 M

      - Q/ t% L8 i( K" d


      【摘要】
   Mammals are unable to regenerate their heart after major cardiomyocyte loss caused by myocardial infarction. Human embryonic stem cells (hESCs) can give rise to functional cardiomyocytes and therefore have exciting potential as a source of cells for replacement therapy. Understanding the molecular regulation of cardiomyocyte differentiation from stem cells is crucial for the stepwise enhancement and scaling of cardiomyocyte production that will be necessary for transplantation therapy. Our novel hESC differentiation protocol is now efficient enough for meaningful genome-wide transcriptional profiling by microarray technology of hESCs, differentiating toward cardiomyocytes. Here, we have identified and validated time-dependent gene expression patterns and shown a reflection of early embryonic events; induction of genes of the primary mesoderm and endodermal lineages is followed by those of cardiac progenitor cells and fetal cardiomyocytes in consecutive waves of known and novel genes. Collectively, these results permit enhancement of stepwise differentiation and facilitate isolation and expansion of cardiac progenitor cells. Furthermore, these genes may provide new clinically relevant clues for identifying causes of congenital heart defects.
      【关键词】 Human embryonic stem cells Cardiomyocytes Microarrays Transcriptional profiling Gene expression Differentiation2 A7 S+ Q- ?( a" w; h; \
   INTRODUCTION* E3 ~1 I  Y* c2 E# O0 w/ R& [& i
: @6 T  {* r' M% Q
Embryonic stem cells (ESCs) are pluripotent cells derived from the inner cell mass of embryos at the blastocyst stage of development. They are characterized by their capacity to self-renew and potential to differentiate into all cells of the adult individual. Since derivation of the first human embryonic stem cell (hESC) lines in 1998 .. e  V- m9 A" z

To identify stem cell-specific features and create a molecular signature for "stemness," transcriptional profiling has been carried out on different stem cell populations. Two independent studies first compared the transcriptional profiles of mouse embryonic stem cells (mESCs), mouse neural cells, and hematopoietic stem cells  to compare up to 75 cell lines under standard growth conditions using standard assays.

Only one study to date has described the transcriptome of hESCs undergoing directed differentiation specifically to hepatic and neuronal lineages  in a culture system in which hESCs are cocultured with a visceral endoderm-like cell line, END-2. The number of beating areas in the hESC-END-2 cocultures was 24-fold higher in the absence of FCS than in the presence of 20% FCS and up to 20% of the cells within the aggregates were hESC-derived cardiomyocytes (hESC-CMs). This improved efficiency made it then feasible to study the transcriptome of hESCs differentiating toward the cardiac lineage.

Here, we describe the temporal gene expression profiles of hESCs during their directed differentiation toward cardiomyocytes by whole-genome microarray analysis. mRNA from human fetal heart (HFH) of 16 weeks of gestation (after voluntary abortus provocatis) was included as a reference source, allowing us to compare gene expression profiles of hESC-CMs with that of HFH. By microarray cluster analysis, genes were identified that were enriched in UH, transiently expressed during early differentiation, and associated with mesoderm or endoderm differentiation, or genes were identified that were enriched in hESC-CMs and HFH. The specificity of the candidate genes was confirmed by reverse transcription-polymerase chain reaction (RT-PCR), real-time RT-PCR, and whole-mount in situ hybridization (ISH) on hESC colonies or on hESC-CM-containing aggregates, mouse embryos, and HFHs. This study demonstrates that directed differentiation of hESCs toward cardiomyocytes can be made sufficiently efficient to provide an excellent model for understanding commitment to the mesoderm or endoderm lineage and the early steps in human cardiac cell differentiation. Furthermore, it is an approach to data mining to identify new genes that may be associated with (abnormal) human heart development.7 ?+ E6 G* w: W" F: o. g) a4 W% t" i
9 g4 c8 K* `0 F
MATERIALS AND METHODS9 [: d: S0 l1 U
/ L- W" k$ P; C7 B  G" y
Cell Culture
: D7 _* e, H7 y% m  O
END-2 cells and HES-2 cells were cultured as described previously . At 1, 3, 6, 9, and 12 days (d) after the start of hESC-END-2 coculture (passages 79 and 80), differentiated hESCs were carefully dissected from the END-2 cell layer to avoid contamination of END-2 cells. Beating areas normally appear from days 6¨C7 onwards and were included in the 9- and 12-day time points. Samples were collected for RNA isolation. For the isolation of undifferentiated hESCs (UH), only undifferentiated HES-2 colonies based on morphology (passages 46, 79, and 82) were removed from MEFs by dispase (10 mg/ml) treatment. After 2 minutes of incubation, HES-2 colonies were lifted and washed twice in phosphate-buffered saline (PBS) and collected for RNA isolation. For the verification of the microarray data by RT-PCR or real-time quantitative PCR (qPCR), different independent samples were used.
# `( w5 g9 z" r3 g" W1 w) \5 \
Primary Human Fetal Cardiomyocytes: g; \" o/ U9 d& \6 Y* G3 }

Primary tissue was obtained after abortion after individual permission had been obtained by use of standard informed consent procedures and approval of the ethics committee of the University Medical Center, Utrecht. Human fetal hearts (HFHs) (11¨C16 weeks of development) were directly frozen in liquid nitrogen and fixed overnight in 2% paraformaldehyde or perfused by Langendorff's method, followed by culture of cardiac cells. After 4 days in culture, fetal cardiac cells were washed with PBS and collected for RNA isolation.) x+ d8 |2 j5 _) M
0 `9 \! q9 v/ G6 C" ?5 j; P! }$ X  S
RNA Isolation and Sample Preparation* @% j: i4 T; l8 a. @: O
# W( ^" y8 x. n* R/ X- P5 X
Total RNA was isolated using TRIzol (Invitrogen, Carlsbad, CA, http://www.invitrogen.com), following the manufacturers' protocol. The quality of all RNA preparations was confirmed prior to microarray hybridization with an Agilent Bioanalyzer 2100 (Agilent Technologies, Palo Alto, CA, http://www.agilent.com) using the RNA 6000 Nano Labchip kit (Agilent Technologies). Equal amounts of total RNA from UH, all time points from hESC-END-2 differentiation (1, 3, 6, 9, and 12 days), and HFH were combined to create a common reference pool (CRP). All microarray hybridizations were performed against this CRP. The use of a CRP has several advantages. First, by comparing each sample directly to the CRP, it is possible to compare relative gene expression in all other samples. Second, it prevents loss of data for genes that are expressed in one sample but not in another. Third, it is possible to extend the data set at a later date by hybridizing a new sample to the CRP, allowing comparisons for the expression of the whole genome to all other samples. Furthermore, we hybridized UH with 12-day differentiated hESCs and performed dye swap for this experiment.9 Q  [+ u, w# U* c% a, U

Probe Preparation and Hybridization of Agilent DNA Microarrays
. @1 q9 A- N! I
Agilent 44K Whole Human genome arrays (G4112A; Agilent Technologies) were prepared and hybridized with linearly amplified and labeled total RNA following the manufacturer's protocol. Briefly, probes were prepared from 500 ng of total RNA by linear amplification and cyanine-3 (cy-3) or cyanine-5 (cy-5) labeling of cRNA using an Agilent Low Input Linear Amplification and Labeling kit. cy-3-CTP (10 mM) and cy-5-CTP (10 mM) were obtained from PerkinElmer/NEN (PerkinElmer Life and Analytical Sciences, Boston, http://www.perkinelmer.com). Fluorescent labeled probes were purified using Qiagen RNeasy Mini kit (Qiagen, Hilden, Germany, http://www1.qiagen.com) as described by the manufacturer. Dye incorporation was confirmed using UV spectrophotometer analysis. According to Agilent's ISH Kit Plus protocol, we used 0.75 µg of labeled cRNA per sample for 17 hours of hybridization at 60¡ãC. After the coverslips were removed in 6x sodium chloride/sodium phosphate/EDTA (SSPE) 0.005% N-laurylsarcosine, the slides were washed in 6x SSPE 0.005% N-laurylsarcosine for 1 minute, washed in 0.06x SSPE 0.005% N-laurylsarcosine for 1 minute, and washed in Agilent stabilization and drying solution for 30 seconds. This wash is to delay degradation of cy-5 dye, which is susceptible to ozone and conveniently dries the slides. Arrays were scanned using the Agilent DNA microarray scanner.
/ d  _  ?0 a% [# r0 u9 e9 S
Data Analysis/Statistics8 B* A6 n- e( N$ V; z. W* _$ e0 i
% Y# O( J8 [8 F; z5 L+ a' x, W
DNA microarray images are processed using Feature Extraction software. The processing that Feature Extraction does involves finding the spots, analyzing the spots (pixel population statistics), analyzing the spot population (feature population statistics), background subtraction, dye normalization (combination of linear and locally weighted linear regression), and then calculating ratios from the processed signals with error values and p values. The signal intensity from the red channel (cy-5) is divided by the signal of the green channel (cy-3), and the scale is compressed by applying log to base 10 (log10) for easier visualization. Luminator from Rosetta Biosoftware (Seattle, http://www.rosettabio.com) was used to analyze the microarray data by clustering using an agglomerative algorithm, trending, and comparing profiles at a p value of 0.01. Additional microarray software analysis was performed using Genespring version 7.2 software (Agilent Technologies).: a) n) Y  t$ r  U0 y- A

RT-PCR$ h7 U0 E! T, n% }
  E* O% O7 K+ A- J& J6 \6 u' E
Total RNA (1 µg) was transcribed by reverse transcriptase (SuperScript II; Invitrogen) and used for PCR using Silverstar DNA polymerase (Eurogentec, Seraing, Belgium, http://www.eurogentec.be). For each PCR, 1 µl of cDNA was used in a reaction volume of 50 µl. The cycling parameters were 94¡ãC for 15 seconds; 58¡ãC for 30 seconds; and 72¡ãC for 45 seconds, for 35 cycles. The PCR cycles were preceded by an initial denaturation of 2 minutes at 94¡ãC and followed by a final extension of 7 minutes at 72¡ãC. Primer sequences for PCR are given in supplementary online Table 1A and 1B. Products were analyzed on ethidium bromide-stained 1.5% agarose gel. Human acidic ribosomal phosphoprotein (hARP) was used as RNA input control. HARP levels did not significantly change during differentiation (data not shown). As a negative control, total RNA was used directly for PCR.8 g$ J; a1 }4 I: v* Z' P
  L5 p0 z6 G- s
Real-Time qPCR

Real-time PCR was performed according to standard protocols on a MyIQ single-color real-time PCR detection system (Bio-Rad, Hercules, CA, http://www.bio-rad.com). Briefly, 1 µg of total RNA was DNase-treated and transcribed to cDNA. Ten µl of a 1/10 dilution of cDNA was then added to 12.5 µl of the 2x SYBR green PCR master mix (Applied BioSystems, Foster City, CA, http://www.appliedbiosystems.com) and 500 µM of each primer. PCR cycles were 3 minutes at 95¡ãC, followed by 40 cycles of 15 seconds at 95¡ãC, 30 seconds at a specific annealing temperature, and 45 seconds at 72¡ãC. The thermal denaturation protocol was run at the end of PCR to determine the number of products. Samples were run on a 2% agarose gel to confirm the correct size of the PCR products. All reactions were run in triplicate. As negative controls, PCR was performed on water and on RNA without reverse transcription. The cycle number at which the reaction crossed an arbitrarily placed threshold (C) was determined for each gene. The relative amount of mRNA levels was determined by the formula 2¨CC. Relative gene expression was normalized to ARP expression. Primer sequences and annealing temperatures used for real-time PCR are shown in supplementary online Table 2.
1 a! a& G# B6 _3 K; s8 I
In Situ Hybridization

Whole-mount in situ hybridization was performed as described before . Probe template primers were designed with the T3-promoter sequence at the 5' end of the forward primer and T7-promoter at the 5' end of the reverse primer. Digoxigenin-labeled probes were generated using the purified PCR product as a template for probe transcription with either T3 RNA polymerase or T7 RNA polymerase (Promega, Madison, WI, http://www.promega.com). Human and mouse probe template primers are shown in supplementary online Table 3. The T3 promoter sequence was 5'-ATACAATTAACCCTCACTAAAGGG-3'. The T7 promoter sequence was 5'-ATAGGTAATACGACTCACTATAGGGC-3'.
7 e% W" |: l& @8 e- x$ X, V4 v) H( w
RESULTS

Temporal Gene Expression During Differentiation of hESCs

Whole-genome microarray analysis was carried out on total RNA from UH (HES-2 cells), hESCs isolated after coculture with END-2 cells for various times (HES-2-END-2 cocultures: 1¨C12 days), and whole HFH (16 weeks of gestation). From all samples, equal amounts of RNA were pooled to create a CRP. Each separate sample was then hybridized against the CRP on Agilent whole-genome microarray chips (Fig. 1A). The use of a CRP allows comparison of relative gene expression levels of each sample with all other samples. Dye-swap experiments confirmed the reproducibility of the microarray data (data not shown). Since we previously have shown that hESC-CMs resemble human fetal cardiomyocytes , RNA from HFH was used as a reference source for comparison to differentiating hESCs. Additional microarray experiments were performed by hybridizing UH directly against hESCs cocultured for 12 days with END-2 (12 days). Again, dye-swap control experiments were performed (Tables 1, 2)." X5 h, L  l1 P, \. @

Figure 1. Whole-genome temporal gene expression profiling of hESCs differentiating toward cardiomyocytes. (A): Experimental approach in which RNA of UH, differentiating human embryonic stem cells (hESCs) from hESC-END-2 coculture (1¨C12 days), and HFH were first isolated, followed by creating a CRP. Each sample was hybridized with the CRP on a 44K Agilent Microarray, and a direct comparison of UH and 12-day differentiated hESCs was done. The first appearance of beating areas (average starting point) is indicated by a red circle on the time line. (B): Number of significantly up- and downregulated genes of UH and all timepoints during hESC differentiation toward cardiomyocytes (p

Table 1. Human embryonic stem cell (hESC)-enriched genes, including seven genes of unknown function that are rapidly downregulated during hESC-END2 differentiation

Table 2. Cardiomyocyte-enriched genes that are upregulated at 12d and human fetal heart, including two genes of unknown function (ADPRHL1 and SYNPO2L)

Differentially expressed genes were identified using a statistical cutoff value of p  .01 (Rosetta Luminator). Figure 1B shows the number of significantly up- and downregulated genes of all samples during hESC differentiation toward cardiomyocytes and HFH. A steady increase in differentially expressed genes is seen as differentiation proceeds. The relatively high number of downregulated genes in the UH and 1-day samples compared with the CRP is likely due to differentiation being initiated but not yet to distinct lineages. The results of all microarray experiments are provided in supplementary online Data 1A¨C1J.  [- v  k3 h3 [, l8 a  r5 k6 R$ R
. H; I  F, ^+ n! j9 p- o7 {
hESCs Preferentially Differentiate Toward Mesoderm and Endoderm# f- M. {9 _) z; Q6 k, A: E
# a. `. G9 u2 a
Previously, we demonstrated that after 12 days of hESC-END-2 coculture, various cell types formed, in addition to cardiomyocytes , expressed in endoderm and trophectoderm. The formation of cardiomyocytes, which are derived from mesoderm, was preceded by transient upregulation of Brachyury T, a marker of nascent mesoderm (Fig. 1C). This suggested that hESCs grown on END-2 cells differentiated preferentially toward mesoderm and endoderm derivatives. This was confirmed by calculating the percentage of tissue-enriched genes that were 3-fold upregulated at 12 days compared directly to UH. Genes were indicated as tissue-enriched according Gene Ontology Tree Machine classifications, representing genes with a p value 10¨C10. Since most genes are significantly expressed in multiple tissues and therefore classified in different tissues, the sum of the percentages of all tissue-enriched genes exceeds 100%. As expected, the highest percentage (62%) was observed for heart-enriched genes (Fig. 1D). In addition, liver (58%) (endoderm-derived) and kidney-enriched (47%) (mesoderm-derived) genes were also highly represented, whereas brain (12%) and skin-enriched (11%) (both ectoderm-derived) genes showed the lowest percentage. Thus hESC-END-2 cocultures differentiate predominantly toward mesoderm and endoderm, with little ectoderm differentiation taking place.

hESC-Enriched Genes5 A0 [: p/ n$ `4 f
/ o2 f0 I+ W- l/ x3 A2 L% U
It is to be expected that genes important for maintenance or self-renewal of hESCs would be rapidly downregulated upon differentiation. By cluster analysis, we identified a group of 15 genes (cluster 1a) that displayed rapid downregulation (Fig. 2A, 2B) and included well-established hESC-markers such as OCT4 and NANOG. Other genes in this cluster were TDGF1, GAL, LEFTB, DNMT3B, and SOX2, also described as being enriched in UH , indicating that the approach used could identify genes that may be involved in hESC self-renewal. Seven genes in this group were novel with unknown function but with expression trend-lines during differentiation similar to those of OCT4 and NANOG (Fig. 2B). A less restrictive cluster (cluster 1b) with similar downregulating gene patterns identified a group of 145 gene entries with potential roles in self-renewal of hESCs (supplementary online data 2). To verify the microarray results, expression levels of a subset of known and unknown genes from both clusters were examined by RT-PCR in two independent UH and 12-day differentiated samples (Fig. 2C).

Figure 2. Human embryonic stem cell (hESC)-enriched genes. (A): Cluster of 15 unique genes (cluster 1a) that are rapidly downregulated upon differentiation, of which seven are unknown. Red indicates upregulation versus common reference pool (CRP); blue indicates downregulation versus CRP. *, since multiple probes are present on the microarray chip, some genes are represented multiple times in one cluster. (B): Genes of cluster 1a shown in a time-course graph with a reference line at y = 0 where the expression ratio is equal to 1. (C): Reverse transcription-polymerase chain reaction confirmation of various hESC-enriched genes, including three unknown genes. (D): In situ hybridizations on hESC colonies for four unknown hESC-enriched genes. OCT4 was used as a positive control, and OCT4 sense probe was used as a negative control. To visualize nonstained hESC colonies, different transmission light settings were used. Arrowheads indicate differentiated areas within hESC colonies. Scale bar = 1 mm for OCT4, OCT4 sense, ENST00000297751, FLJ10884, and C14ORF115. Scale bar = 500 µm for THC2173136. Abbreviations: d, day(s); RT, reverse transcription; UH, undifferentiated human embryonic stem cells.- t% h/ \% }2 }% S: F+ R# ]1 d

For four unknown genes (ENST00000297751, THC1818772, FLJ10884, and C14ORF115), we carried out whole-mount ISH on UH colonies, cultured on MEFs, and compared expression with that of OCT4 (Fig. 2D). High expression of the selected genes was observed in the parts of the colonies with characteristic undifferentiated morphology, whereas the differentiated parts and the MEFs remained negative. This confirmed association of these genes with the pluripotent state.
  J( u# `* O8 }( ?% |( S2 ]" o* Y
Activated and Repressed Signaling Pathways During hESC Differentiation
& m5 D8 A% U+ V# Y1 R2 B$ l0 u0 Q3 [
The earliest appearance of beating areas in hESC-END-2 cocultures was consistently at 6¨C7 days. Before the appearance of beating areas, cells are committed to form cardiomyocytes, and regulatory as well as structural genes are activated. Accordingly, expression of the earliest mesoderm marker Brachyury T is transiently upregulated in 3-day hESC-END-2 cocultures, 3 or 4 days before the onset of beating, as described above. Clustering analysis revealed a group of 25 genes (Table 3, cluster 2a), transiently co-upregulated with Brachyury T on 3 days of hESC-END-2 coculture (Fig. 3A, 3B). Interestingly, in addition to eight unknown genes, the clustered group contained genes shown previously to be important for early embryonic development and for mesoderm formation in particular. The Notch signaling pathway, with ligands -like 1 (DLL) and DLL3, upregulated at 3 days, has pleiotropic functions during development, including a role in differentiation of paraxial mesoderm, which later gives rise to skeletal muscle  on cardiomyocyte differentiation have been described. Other genes in this cluster were the transcription factors FOXC1 and MESP1 and the transcriptional repressor SNAI1, which are all expressed in mesoderm and important for cardiac development. The majority of these genes also play roles in epithelial-mesenchymal-transformation (EMT), important for many processes during development, including that of the heart. The transient expression of these genes during differentiation was confirmed by quantitative real-time PCR (Fig. 3C). Interestingly, for DKK1, WNT3A, FOXC1, and SNAI1 a biphasic expression pattern was observed with peaks at 3 and 12 days, suggesting a role in both early and late phases of differentiation. A larger cluster with a similar pattern (cluster 2b) is available as supplementary online Data 2. Besides activation of genes involved in mesoderm formation, genes associated with endoderm formation, such as FOXA2 and SOX17, displayed a similar transient upregulation with a peak at 3 days (data not shown).3 r' x# ]0 c" [& {

Table 3. Genes transiently upregulated at 3d of differentiation include established mesoderm markers

Figure 3. Transiently upregulated genes during human embryonic stem cell (hESC) differentiation. (A): Cluster of 25 unique genes (cluster 2a) that are transiently upregulated at day 3 of differentiation. Among these genes, we found mesoderm markers and genes that are important in early cardiac development. *, since multiple probes are present on the microarray chip, some genes are represented multiple times in one cluster. (B): Genes of cluster 2a are shown in a time-course graph with a reference line at y = 0 where the expression ratio is equal to 1. (C): Quantitative real-time PCR also demonstrates a transient upregulation of genes implicated in mesoderm and heart formation. (D): Schematic representation of the regulation of genes of the WNT family during hESC differentiation. Transiently upregulated expression is indicated by blue-red-blue coloring, downregulation by blue, and upregulation by red. (BRA-T) is downstream of LEF1 and is important for mesoderm formation, whereas SOX17 may interact with LEF1 for the formation of endoderm. Abbreviations: BRA-T, Brachyury T; ß-CAT, ß-catenin; d, day(s); UH, undifferentiated human embryonic stem cells.2 L4 z, o* ^  Q3 ^% G7 [

We next identified a cluster of genes (cluster 3) with expression transiently increased at 6 days and enhanced expression in HFH. This cluster consisted of 204 gene entries, including the transcription factors LBX1, IRX3, IRX4, FOXC2, HES7, and SRF, all of which have been implicated in cardiac development (supplementary online Data 2). In addition, MSX1 and ISL1, homeodomain transcription factors, also increased from 3¨C12 days and from 6¨C12 days, respectively (supplementary online Fig. 1). MSX1 plays a role in EMT and cardiac development, whereas ISL1 has recently been described as a marker for a cardiac progenitor cell (CPC) population .1 t8 N$ x5 _! W9 s/ O3 k  i

The WNT, fibroblast growth factor (FGF), transforming growth factor ß (TGFß), and bone morphogenetic protein (BMP) signaling pathways have been shown to play important roles during embryonic development and early cardiogenesis. The early transient increase of members of the WNT family, DKK1, DKK3, and WNT3a have been discussed above, but other WNT family members, such as FRZ10, LEF1, WNT5A, WNT5B, and WNT11, showed a steady increase in expression during differentiation, whereas secreted FRZ1, FRZ2, Cerberus, WNT8B, FRZ5, and SRY17 were downregulated during differentiation (Fig. 3D). For the TGFß and FGF families, FGF9, FGF receptor 2, TGFß1, TGFß2, TGFß3, SMAD3, TGFß receptor III, and Activin A receptor type I and II were upregulated, whereas nodal, FGF receptor 1, and FGF2 (basic FGF) were downregulated. Interestingly, BMP10, enriched in heart tissue and essential for cardiogenesis in mouse , were increased during differentiation (supplementary online Data 3; supplementary online Fig. 2).  E1 y9 L: w! G
$ A6 d' `: k8 p6 t! K: N& z4 K
Cardiac-Enriched Genes

Since the number of beating areas are increasing from day 9 to day 12 of differentiation, genes responsible for cardiomyocyte differentiation and function are expected to be upregulated during this time period. As reference source, HFH was used to increase the possibility that upregulated genes may be associated with human cardiac development. A cluster of 18 genes (cluster 4a), upregulated at 9 or 12 days of differentiation and in HFH was identified (Fig. 4A, 4B). This group consisted of known cardiac genes, such as troponin T 2 (TNNT2), myosin light chain 4 (MYL4), -actinin-2 (ACTN2), and myosin light chain 2B (HUMMLC2B), but also four novel genes and several genes that have been annotated but not previously associated with cardiac development. A larger cluster of 145 gene entries with similar expression pattern (cluster 4b) is represented in supplementary online Data 2. Verification of a subset of the cardiac-enriched genes was performed by RT-PCR (Fig. 4C). Many transcription factors known to play a role in cardiac development were upregulated during differentiation. Real-time RT-PCR for cardiac transcription factors MEF2C, TBX2, and TBX5 demonstrated increased expression at 9 and 12 days of differentiation. By comparison, -actinin, a sarcomeric protein, was only increased after 12 days of differentiation (Fig. 4D). These data demonstrate distinct molecular regulation during each step of differentiation (Fig. 4E). Furthermore, ISH for a set of the cardiac-enriched genes, including phospholamban and two novel genes annotated as SRD5A2L2 and SYNPO2L (from clusters 4a and 4b), on microdissected beating areas from 12-day hESC-END-2 cocultures showed strong expression (Fig. 5A¨C5C). Expression of these genes was not detected in the nonbeating areas derived from the same hESC-END-2 cocultures. Mouse orthologs were found for both SRD5A2L2 and SYNPO2L and whole-mount ISH on mouse embryos at different development stages showed cardiac-restricted/cardiac-enriched expression in early stages of cardiac development for both genes. At E8.5 and E9.5, SRD5A2L2 was specifically expressed in the sinus venosus, giving rise to the inflow tract, and at lower levels in the ventricle (Fig. 5D¨C5F). SYNPO2L was detected in the heart from E8.5 and in myotomes, which give rise to skeletal muscle, at E10.5 (Fig. 5G¨C5I). In addition, whole-mount ISH on intact HFH (Fig. 5J) or hearts cut in half (Fig. 5K¨C5M) showed that the unknown genes THC2339346, SRD5A2L2, and SYNPO2L were all expressed as the positive control phospholamban in the fetal heart. As in mouse embryos, SRD5A2L2 was detected predominantly in the inflow tract of HFHs.0 C, z, d' W, u8 s+ h. J2 [6 T
7 y* _: `% j" r3 O! e- r5 Q
Figure 4. Cardiomyocyte-enriched genes at 12 days of differentiating hESCs. (A): Cluster of 15 unique genes (cluster 4a) that are upregulated at 12 days as well as in HFH of which three genes of unknown function. *, since multiple probes are present on the microarray chip, some genes are represented multiple times in one cluster. (B): Genes of cluster 4a in a time-course graph, including HFH for comparisons with 12 days with a reference line at y = 0 where the expression ratio is equal to 1. (C): Reverse transcription-polymerase chain reaction confirmation of a subset of genes (from clusters 4a and 4b) of known and unknown function that are expressed at 12 days and HFH but not or at low levels in UH. ARP was used as an internal control. (D): Quantitative real-time polymerase chain reaction of cardiac transcription factors MEF2C, TBX2, and TBX5 and cardiac marker ACTN2. (E): Schematic representation of hESCs differentiating toward the cardiac lineage compared with embryonic development. UH resemble the epiblast stage during embryonic development. Shortly after the start of hESC-END-2 coculture, (3-day) gene expression profiles indicate the onset of gastrulation with a preference for endoderm and mesoderm formation, followed by the probable formation of cardiac progenitor cells (CPCs) between 3 and 6 days, and optimal differentiation to CMs at 12 days. At the bottom, the contributions of genes from the different clusters during the different phases of differentiation are indicated. The shape of the triangles represents the relative expression of the genes within that cluster during differentiation, with two examples of genes in each triangle. The first appearance of beating areas (average starting point) is indicated by a red circle on the time line. Abbreviations: ARP, acidic ribosomal protein; CL, cluster; CM, cardiomyocyte; CPC, cardiac progenitor cell; d, day(s); hESCs, human embryonic stem cells; HFH, human fetal heart; RT, reverse transcription; UH, undifferentiated human embryonic stem cells.2 E" Q: Q2 k# v3 X6 x8 ^" b

Figure 5. In vitro and in vivo expression analysis of novel cardiac specific genes by whole-mount in situ hybridization. (A): Expression of PLN in an excised beating area at 12 days. The inset shows a nonbeating area at 12 days that shows no expression of PLN. (B): Expression of SRD5A2L2 in a 12-day beating area, whereas the nonbeating area (inset) shows no expression. (C): Expression of SYNPO2L in an excised beating area at 12 days. Nonbeating (inset) area showed no expression. (D): Expression of mouse SRD5A2L2 at day 8.5 postfertilization is restricted to the heart, with strongest expression in the ift. Dorsal to the right, anterior to the top. (E): Transverse section of an E9.5 mouse embryo with strong expression of SRD5A2L2 in both horns of the sinus venosus, lower expression in other parts of the heart. Dorsal to the top. (F): Magnification of right horn of sinus venosus showing strong expression of SRD5A2L2 at E9.5. Dorsal to the top. (G): Expression of mouse SYNPO2L at E9.5 is restricted to the heart and somites. Dorsal to the right, anterior to the top. (H): Transverse section of E8.5 mouse embryo showing SYNPO2L expression in the heart tube and developing somites. (I): Transverse section of an E10.5 mouse embryo showing SYNPO2L expression in the myotome. (J): SYNPO2L expression in human fetal heart at 11 weeks of gestation. Expression of SYNPO2L in both ventricles and stronger in both atria but absent in aorta and pulmonary vein. (K): Expression of SRD5A2L2 predominantly in the vena cava and pulmonary vein of a 14-week human fetal heart. (L): Ubiquitous expression of THC2339346 in a 17-week human fetal heart. (M): Phospholamban expression in a 20-week human fetal heart. Scale bars = 500 µm (A¨CD, G), 300 µm (E, H, I), 150 µm (F), and 2 mm (J¨CM). Abbreviations: ao, aorta; avc, atrioventricular canal; E, embryonic day; ht, heart tube; ift, inflow tract; la, left atrium; ls, left horn sinus venosus; lv, left ventricle; lvw, left ventricular wall; m, myotome; nf, neural fold; nt, neural tube; oft, outflow tract; PLN, phosholamban; pv, pulmonary vein; ra, right atrium; rs, right horn sinus venosus; rv, right ventricle; rvw, right ventricular wall; s, somite; sp, septum; SRD5A2L2, steroid reductase 5-2-like 2; st, septum transversum; svc, superior vena cava; SYNPO2L, synaptopodin 2-like; W, weeks.

DISCUSSION/ z  \- V) n2 ~

In the present study, we described the temporal gene expression of hESCs during differentiation toward the cardiac lineage. Besides information on genes that may play roles in self-maintenance or self-renewal of hESCs, genes important during the different stages of cardiac development were also identified. It is expected that genes with important roles in self-renewal of ESC would be rapidly downregulated upon differentiation. For this reason, several studies have compared the transcriptome of UH with differentiated hESCs , microarray analysis was performed on 2-, 10-, 20-, and 30-day differentiated EBs from hESCs. However, 10-day EBs had essentially unchanged levels of OCT4 expression, suggesting that differentiation by formation of EBs was relatively inefficient and slow. We demonstrated here that temporal transcriptional profiling of hESCs can be a rich source of information on the molecular mechanisms underlying human development, provided differentiation is directed and efficient, as in this study for the cardiac lineage. By cluster analysis of microarray data, we identified genes that were rapidly downregulated upon differentiation, which included OCT4 and NANOG, as well as multiple unknown genes, and confirmed their expression independently in UH colonies. Whether these genes play roles in the self-renewal of stem cells or are simply stem cell markers has to be determined in functional studies., r& s0 W) V7 S' G, v6 z# C8 a

Previously, we have shown that hESCs in coculture with END-2 cells differentiated predominantly to endoderm and mesoderm derivatives . Here, this was confirmed by microarray analysis with transient expression of the early mesodermal marker brachyury T and the endodermal markers FOX2A and SOX17, peaking after 3 days of differentiation (3¨C4 days before the appearance of the first beating areas). Other mesodermal markers or transcription factors, such as SNAI1, MSX1 MESP1, and FOXC1, were transiently upregulated, with peaks at 3 days. The detailed temporal gene profiling of differentiating hESCs presented here may lead to improvements in stepwise directed differentiation toward mesoderm and endoderm-derived cell-types or tissues, if specific signaling pathways are further activated or enhanced at the times indicated.: H+ R$ P4 v+ Z: a0 u+ O7 b' M# w

Many studies have demonstrated the importance of endoderm for the induction of cardiac differentiation . The presence of endoderm as well as mesoderm markers in differentiating hESCs suggests that differentiation to cardiomyocytes in hESC-END-2 cultures may be either direct, by the END-2 cells, and/or indirect, by the formation of endoderm derived from hESCs, inducing adjacent mesodermal cells to become cardiomyocytes.

Proteins from the WNT, TGFß, BMP, and FGF signaling pathways have been implicated in both self-renewal of stem cells as well as cardiac differentiation. Basic FGF, important for self-renewal , who demonstrated that the BMP inhibitor Noggin only induced cardiomyocyte differentiation of mESCs when applied in a very narrow time window during ESC culture.

Microarray studies have been performed on mouse P19 embryonal carcinoma cells differentiating to cardiomyocytes and during early mouse cardiogenesis, using an Nkx2.5 enhancer to select cardiac progenitor cells . The human heart, however, is not accessible at the first crucial stages of development, 4 weeks of gestation being the earliest we have obtained when the heart already has four defined chambers. Here, we have described for the first time the transcriptome of hESCs differentiated toward the cardiac lineage. Using HFH as a reference source, we were able to identify genes that were both upregulated during differentiation of hESC-CMs and in HFH. In addition to the presence of known cardiac sarcomeric genes, several unknown genes and genes not previously associated with cardiomyocytes were identified (clusters 4a and 4b). Two of these unknown genes, SRD5A2L2 and SYNPO2L, displayed cardiac expression during mouse embryonic development and in HFHs. In addition, THC2339346, THC1564329, and THC1452070, for which no mouse ortholog could be found, were expressed in HFH. These findings indicate that in vitro differentiation of hESCs as used in this study is a very powerful model for human cardiac development.
4 i2 b$ a$ d% D
Using hESCs as a model for human cardiac development will likely result in the discovery of genes not previously associated with cardiomyocyte differentiation or function and may lead to a better understanding of abnormal cardiac development. Furthermore, improved cardiomyocyte differentiation will lead to a higher production and purer populations of cardiomyocytes, facilitating entry into regenerative medicine for heart patients." `; l& s. }% `, e7 x) W; ?4 Q" k
2 O8 X- g$ p2 _/ L: k7 X* _) Y
DISCLOSURES2 K. Y/ q6 ^+ N
7 r  ^0 _6 w% N
J.K. has performed contract work for ES Cell International within the last 2 years.
, J1 ]0 A$ y+ w
ACKNOWLEDGMENTS4 l8 \6 f  g' P; k4 c
4 C2 e! {5 B! H1 p+ h; p
We thank Linda W. van Laake and Rutger J. Hassink for the isolation of the human fetal hearts and Jeroen Korving for histology. We thank ServiceXS (Leiden, The Netherlands) for help and use of Rosetta Luminator software for data analysis. Furthermore, we thank Hans C. Clevers and Chris N. Denning for critical reading of the manuscript. This study was supported by ES Cell International (J.K.) and the Interuniversity Cardiology Institute of the Netherlands (D.W.-v.O., R.P., and C.M.).
      【参考文献】
7 }+ a6 g. X; f( H8 ^1 D6 S! |* I6 Z
4 G$ ?( L% Q% c6 S/ z+ M
Thomson JA, Itskovitz-Eldor J, Shapiro SS et al. Embryonic stem cell lines derived from human blastocysts. Science 1998;282:1145¨C1147.% ]2 a8 X0 l1 O' D& e, _8 v
# p1 `3 B- ]. H
Reubinoff BE, Pera MF, Fong CY et al. Embryonic stem cell lines from human blastocysts: Somatic differentiation in vitro. Nat Biotechnol 2000;18:399¨C404.; @9 t3 X8 g# h& _, a& B6 x
, }3 ^" H# F$ O7 A+ G8 M& c
Cowan CA, Klimanskaya I, McMahon J et al. Derivation of embryonic stem-cell lines from human blastocysts. N Engl J Med 2004;350:1353¨C1356.) U7 ]0 R: R. U

Passier R, Mummery C. Origin and use of embryonic and adult stem cells in differentiation and tissue repair. Cardiovasc Res 2003;58:324¨C335.

Ramalho-Santos M, Yoon S, Matsuzaki Y et al. "Stemness": Transcriptional profiling of embryonic and adult stem cells. Science 2002;298:597¨C600.
$ J, S0 A6 C) K2 L; O; p1 N- x( q1 A
Ivanova NB, Dimos JT, Schaniel C et al. A stem cell molecular signature. Science 2002;298:601¨C604.$ Q/ b' ^5 B+ {, y- m
/ x& Z6 L3 |! p' X
Ginis I, Luo Y, Miura T et al. Differences between human and mouse embryonic stem cells. Dev Biol 2004;269:360¨C380.+ J0 w, y$ J7 w* U5 A- `3 u
. X( n4 f! G; H# g1 o. L. @+ X
Wei CL, Miura T, Robson P et al. Transcriptome profiling of human and murine ESCs identifies divergent paths required to maintain the stem cell state. STEM CELLS 2005;23:166¨C185.
& l. }& j) z1 X4 E8 P: Z4 Z: Q& g
Sperger JM, Chen X, Draper JS et al. Gene expression patterns in human embryonic stem cells and human pluripotent germ cell tumors. Proc Natl Acad Sci U S A 2003;100:13350¨C13355.1 F' F3 z  ^) z: L. F1 O* c/ X( g
3 [$ f' z: X, ]4 {# F/ o+ V
Carpenter MK, Rosler E, Rao MS. Characterization and differentiation of human embryonic stem cells. Cloning STEM CELLS 2003;5:79¨C88.
  y( E( L1 ]+ P* v7 c1 d
Zeng X, Miura T, Luo Y et al. Properties of pluripotent human embryonic stem cells BG01 and BG02. STEM CELLS 2004;22:292¨C312.' `/ {5 [0 g! a, [/ t4 m3 C) G
0 C. m; _7 u2 l" Y+ o( }3 E
Abeyta MJ, Clark AT, Rodriguez RT et al. Unique gene expression signatures of independently-derived human embryonic stem cell lines. Hum Mol Genet 2004;13:601¨C608.

Bhattacharya B, Miura T, Brandenberger R et al. Gene expression in human embryonic stem cell lines: Unique molecular signature. Blood 2004;103:2956¨C2964.

Sato N, Sanjuan IM, Heke M et al. Molecular signature of human embryonic stem cells and its comparison with the mouse. Dev Biol 2003;260:404¨C413.: j$ I: R  n, d* R0 R- P  D

Brandenberger R, Wei H, Zhang S et al. Transcriptome characterization elucidates signaling networks that control human ES cell growth and differentiation. Nat Biotechnol 2004;22:707¨C716.
2 \, P/ I" l- x6 l; }
Calhoun JD, Rao RR, Warrenfeltz S et al. Transcriptional profiling of initial differentiation events in human embryonic stem cells. Biochem Biophys Res Commun 2004;323:453¨C464.4 H3 s6 o/ T5 h4 K

Miura T, Luo Y, Khrebtukova I et al. Monitoring early differentiation events in human embryonic stem cells by massively parallel signature sequencing and expressed sequence tag scan. Stem Cells Dev 2004;13:694¨C715.& b1 t: }0 p5 L7 T: v

Dvash T, Mayshar Y, Darr H et al. Temporal gene expression during differentiation of human embryonic stem cells and embryoid bodies. Hum Reprod 2004;19:2875¨C2883.* w8 B, l2 Q4 h& o

Andrews PW, Benvenisty N, McKay R et al. The International Stem Cell Initiative: Toward benchmarks for human embryonic stem cell research. Nat Biotechnol 2005;23:795¨C797.

Mummery C, Ward D, van den Brink CE et al. Cardiomyocyte differentiation of mouse and human embryonic stem cells. J Anat 2002;200:233¨C242.% G' K+ }$ F& Q. I. _$ W
  g; Z6 F% O" B- ^% _
Mummery C, Ward-van Oostwaard D, Doevendans P et al. Differentiation of human embryonic stem cells to cardiomyocytes: Role of coculture with visceral endoderm-like cells. Circulation 2003;107:2733¨C2740.* G; t9 c! j- G" G$ {: u8 Y! m
8 C0 Y4 f0 q( F( t2 l
Xu C, Police S, Rao N et al. Characterization and enrichment of cardiomyocytes derived from human embryonic stem cells. Circ Res 2002;91:501¨C508.7 _' ]8 h6 e' ?# |
5 v9 j, X: K! L- F4 l2 q
He JQ, Ma Y, Lee Y et al. Human embryonic stem cells develop into multiple types of cardiac myocytes: Action potential characterization. Circ Res 2003;93:32¨C39.+ B9 \3 v& t$ z; U( Y! Z. W* m  @
4 x5 Y' L+ S$ b
Passier R, Oostwaard DW, Snapper J et al. Increased cardiomyocyte differentiation from human embryonic stem cells in serum-free cultures. STEM CELLS 2005;23:772¨C780.; Q* [) S% Q' r' V0 F$ C% x
2 T6 _; H3 ]& k  l9 g
Nieto MA, Bennett MF, Sargent MG et al. Cloning and developmental expression of Sna, a murine homologue of the Drosophila snail gene. Development 1992;116:227¨C237.) _6 r: }' u2 P8 K

Bettenhausen B, Hrabe de Angelis M, Simon D et al. Transient and restricted expression during mouse embryogenesis of Dll1, a murine gene closely related to Drosophila Delta. Development 1995;121:2407¨C2418.
. v7 J! u* x0 Z. \4 z0 l
Jen WC, Wettstein D, Turner D et al. The Notch ligand, X-Delta-2, mediates segmentation of the paraxial mesoderm in Xenopus embryos. Development 1997;124:1169¨C1178.$ `: S3 |6 A/ b% [5 }

Marvin MJ, Di Rocco G, Gardiner A et al. Inhibition of Wnt activity induces heart formation from posterior mesoderm. Genes Dev 2001;15:316¨C327.  h% {+ Y2 K4 ]. `) ]7 {9 v
7 |# h8 U3 V0 T% d% Z
Nakamura T, Sano M, Songyang Z et al. A Wnt- and beta-catenin-dependent pathway for mammalian cardiac myogenesis. Proc Natl Acad Sci U S A 2003;100:5834¨C5839.
5 X$ R6 k- g* @
Schneider VA, Mercola M. Wnt antagonism initiates cardiogenesis in Xenopus laevis. Genes Dev 2001;15:304¨C315.4 V/ z1 ?- u3 ~0 K* T: `& R
2 x) Z- O/ U4 P2 F' D
Laugwitz KL, Moretti A, Lam J et al. Postnatal isl1  cardioblasts enter fully differentiated cardiomyocyte lineages. Nature 2005;433:647¨C653.' y& b' R" i( t. N1 `
; g% o5 n  @0 l% t2 Y4 |# z7 `2 O
Pashmforoush M, Lu JT, Chen H et al. Nkx2¨C5 pathways and congenital heart disease; loss of ventricular myocyte lineage specification leads to progressive cardiomyopathy and complete heart block. Cell 2004;117:373¨C386.  e0 k6 w5 ~8 E. k- {2 V) Y6 A9 j
& }; k+ e4 g) g# ~8 k7 B# L# T6 ]
Schultheiss TM, Burch JB, Lassar AB. A role for bone morphogenetic proteins in the induction of cardiac myogenesis. Genes Dev 1997;11:451¨C462.9 U! z# j4 [( y, |. F: x6 H

Somi S, Buffing AA, Moorman AF et al. Dynamic patterns of expression of BMP isoforms 2, 4, 5, 6, and 7 during chicken heart development. Anat Rec A Discov Mol Cell Evol Biol 2004;279:636¨C651.
5 k$ J/ b" \4 N. t* `2 ^7 W
Lough J, Sugi Y. Endoderm and heart development. Dev Dyn 2000;217:327¨C342.5 V2 {- U( A  i. E' [7 J: z# e
  E8 V7 P" U* i, x$ P
Xu RH, Peck RM, Li DS et al. Basic FGF and suppression of BMP signaling sustain undifferentiated proliferation of human ES cells. Nat Methods 2005;2:185¨C190.9 h  c( g+ ~% T3 f! `5 H) t

Sato N, Meijer L, Skaltsounis L et al. Maintenance of pluripotency in human and mouse embryonic stem cells through activation of Wnt signaling by a pharmacological GSK-3-specific inhibitor. Nat Med 2004;10:55¨C63.2 d5 Z2 x3 F3 ]
; f1 V- s# ^. d1 u- f) g
Yuasa S, Itabashi Y, Koshimizu U et al. Transient inhibition of BMP signaling by Noggin induces cardiomyocyte differentiation of mouse embryonic stem cells. Nat Biotechnol 2005;23:607¨C611.1 e* X% D4 ]) m* |* S) |

Masino AM, Gallardo TD, Wilcox CA et al. Transcriptional regulation of cardiac progenitor cell populations. Circ Res 2004;95:389¨C397.