|
|
||||||||
OPEN ACCESS ARTICLE
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
EMBRYONIC STEM CELLS |
aDivision of Developmental and Regenerative Medicine, Department of Obstetrics, Gynecology, and Reproductive Sciences, University of Pittsburgh School of Medicine, Pittsburgh Development Center, Magee-Womens Research Institute and Foundation, Pittsburgh, Pennsylvania;
bUniversity of Pittsburgh School of Medicine, Dorothy P. and Richard P. Simmons Center for Interstitial Lung Disease, Pulmonary, Allergy, and Critical Care Medicine, University of Pittsburgh Medical Center, Pittsburgh, Pennsylvania
Key Words. Embryonic stem cells • Pluripotency • Differentiation • Gene expression • Primates
Correspondence: Christopher S. Navara, Ph.D., 204 Craft Avenue, Pittsburgh, Pennsylvania 15213, USA. Telephone: (412) 641-2430; Fax: (412) 641-2410; e-mail: cnavara{at}pdc.magee.edu
Received April 19, 2007;
accepted for publication July 2, 2007.
First published online in STEM CELLS EXPRESS July 19, 2007.
| ABSTRACT |
|---|
|
|
|---|
1 year with stable diploid karyotype (except for one stable trisomy) and expresses in vitro pluripotency markers, and eight have already formed teratomas. Unlike the heterogeneous gene expression profiles found among hESCs, these nhpESCs display remarkably homogeneous profiles (>97%), with full-sibling lines nearly identical (>98.2%). Female nhpESCs express genes distinct from their brother lines; these sensitive analyses are enabled because of the very low background differences. Experimental comparisons among these primate ESCs may prove more reliable than currently available hESCs, since they are akin to inbred mouse strains in which genetic variables are also nearly eliminated. Finally, contrasting the biological similarities among these lines with the heterogeneous hESCs might suggest that additional, more uniform hESC lines are justified. Taken together, pedigreed primate ESCs display homogeneous and reliable expression profiles. These similarities to mouse ESCs suggest that heterogeneities found among hESCs likely result from their disparate origins rather than intrinsic biological limitations with primate embryonic stem cells. Disclosure of potential conflicts of interest is found at the end of this article.
| INTRODUCTION |
|---|
|
|
|---|
ESCs from mice remain unique in their ability to generate chimeric fetuses and offspring after reaggregation with diploid or tetraploid mouse embryos, and mESCs are capable of transmission through the germ line to form either fertile sperm or eggs [10, 11]. Nuclear reprogramming with either abnormal mouse eggs or in vitro has just been reported using mouse cells [12–15]. Embryos from various mammals differ markedly in their biological capacities, with domestic species succeeding in somatic cell nuclear transfer [16] before mice [17], and yet ESC lines from other species, including rats, have yet to be characterized unequivocally [18, 19]. Mice are unrivalled for their experimental utilities, and mouse embryos are also capable of arresting development in utero at the blastocyst stage, a process known as diapause [20]. These properties may also contribute to the remarkable features of mESCs to remain stable in the undifferentiated state while continuing to proliferate. This also raises questions as to whether ESC cells from other mammals, including primates, will share some or all of these properties.
Notwithstanding the invaluable research contributions of mESCs, there are still some limitations for modeling human diseases with mESCs, prompting primate ESC derivations first in monkeys in 1995 (from in vivo embryos flushed after mating) [21] and then in humans in 1998 [22]. Since the original five hESC lines published by Thomson et al. [22], dozens of new hESC lines have been established [5, 23–26] from clinically discarded embryos generated by Assisted Reproductive Technologies (ART) clinics around the world and donated anonymously by infertile patients/couples. The journal Nature reports >275 hESC lines worldwide [27], and the International Society for Stem Cell Research's registry cites well over 100 hESC lines (http://isscr.org/science/sclines.htm). Nonhuman primates afford research opportunities beyond those available from human ART specimens, including the deliberate production of embryos with defined genetics using pedigreed fertile gametes as well as chimera and germ-line transmission testing for pluripotency. The sequencing of the rhesus genome [28] now permits direct genomic comparison among monkeys, apes, and humans. American and Japanese researchers have derived several macaque nonhuman primate ESC (nhpESC) lines [29–31], including a parthenogenetic ESC line [32]. Pluripotency has been demonstrated in nhpESCs using hESC criteria, that is, ability to form tissues representative of all three germ layers in vitro and in teratomas. To address whether primate ESCs are biologically heterogeneous or whether variations among hESCs might be due to their varying origins, we established pedigreed nhpESCs using gametes from fertile rhesus and compared expression profiles and pluripotency within and among these family groups.
Here, we describe the derivation of ten new nonhuman primate embryonic stem cell lines of known pedigree and demonstrate that they are pluripotent using in vitro and teratoma assays. Derivation success was directly related to embryo quality as measured by rapid preimplantation development where rate of development correlated with likelihood of success in ESC line establishment. The gene expression patterns of these 10 lines are remarkably homogeneous but with still-detectable differences among families and between sexes. Taken together, this suggests that the variability among hESCs is likely the result of their derivation from genetically diverse embryos in disparate research facilities around the world, and that the observed heterogeneity among hESCs is not due to intrinsic biological variability, as compared with mice, but rather is due to the restricted numbers and qualities of the hESC lines widely available.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Establishment, Characterizations, and Confirmations of Pedigreed nhpESC Lines
Pluripotency Markers Detected by Immunocytochemistry.
nhpESCs were assayed for characteristic pluripotency markers Oct-4, Nanog, stage-specific embryonic antigen (SSEA)-4, TRA 1-81, and TRA 1-60 as well as the negative hESC-nhpESC marker SSEA-1 (Oct-4, TRA 1-81, and TRA 1-60 from [Santa Cruz Biotechnology Inc., Santa Cruz, CA, http://www.scbt.com]; Nanog from [R&D Systems Inc., Minneapolis, http://www.rndsystems.com]; SSEA-4 and SSEA-1 from [Developmental Studies Hybridoma Bank, Iowa City, IA, http://www.uiowa.edu/
dshbwww]). Immunocytochemistry was performed on undifferentiated colonies after fixation by 2% paraformaldehyde in phosphate-buffered saline (PBS) for 40 minutes. After washes in PBS + 1% Triton X-100 (PBS-Tx; Sigma) for 15 minutes, nonspecific binding of the primary, excluding Nanog, and secondary antibodies was blocked by a 30-minute incubation in PBS containing 5% goat serum and 0.3% bovine serum albumin. Nanog staining was carried out without the blocking step. Primary antibodies were applied for 40 minutes at 37°C in a humidified chamber, washed extensively in PBS-Tx, and then incubated in appropriate fluorescently labeled secondary antibodies (40 minutes). DNA was detected with 5 µM TOTO-3 (Molecular Probes, Eugene, OR, http://probes.invitrogen.com) added at room temperature 20 minutes prior to mounting coverslips in Vectashield antifade medium (Vector Laboratories, Burlingame, CA, http://www.vectorlabs.com) to retard photobleaching. Slides were examined using laser scanning confocal microscopy [37].
Pluripotency Markers Detected by Reverse Transcription-Polymerase Chain Reaction. Pluripotent nhpESCs were collected by scraping and then pelleted by centrifugation at 200g for 5 minutes. RNA was isolated using TRIzol (Invitrogen). Briefly, 100 µl of TRIzol was added per 100–1,000 cells and vortexed to lyse cells to homogeneity; 200 µl of chloroform (Sigma) was added followed by mixing and centrifugation for 5 minutes at 25,000g. The RNA containing supernatant was removed and the RNA pelleted using 600 µl of 100% isopropanol added to the supernatant and incubated at –20°C for at least 4 hours. The sample was centrifuged twice at 13,000g for 30 minutes at 4°C, with an ethanol wash in between, followed by air-drying. The RNA pellet was reconstituted in nuclease-free water and treated with 1 µl of DNase I for 30 minutes at 37°C. cDNA was prepared using the ImProm-II Reverse Transcription System (Progen, Heidelberg, Germany, http://www.progen.de) according to manufacturer's directions. Primers used were Oct-4, forward cgaccatctgccgctttgag and reverse ccccctgtcccccattccta, Nanog forward ctgtgatttgtgggcctgaa and reverse tgtttgcctttgggactggt, Rex1 forward gcgtacgcaaattaaagtccaga and reverse cagcatcctaaacagctcgcagaat, and Sox-2 forward cccccggcggcaatagca and reverse tcggcgccggggagatacat.
Pluripotency Demonstrated in Teratomas. Please see supplemental online Methods.
Independent External Genetic Testing for Pedigreed Confirmations. Pedigree of each cell line described was confirmed independently by DNA genotyping by Dr. Cecilia Penedo, Veterinary Genetics Laboratory, University of California, Davis (supplemental online Table 3).
Cytogenetic Analysis to Assay Normal Karyotype. Independent cytogenetic investigations were performed by the clinical cytogenetics facility of the University of Pittsburgh Cancer Institute under the supervision of Dr. Susanne Gollin and also within our lab [35].
Gene Expression Profiling
RNA Extraction. Total RNA was extracted using the TRIzol protocol (100 microliters per 104–105 cells) [38] followed by the Qiagen RNeasy Micro Kit (Hilden, Germany, http://www.qiagen.com) according to the manufacturer's recommendations to "clean up" the RNA isolated using TRIzol. RNA quality and quantity were determined using a NanoDrop spectrophotometer (NanoDrop, Wilmington, DE, http://www.nanodrop.com) and Agilent Bioanalyzer (Agilent Technologies, Palo Alto, CA, http://www.agilent.com).
Preparation of Labeled cRNA.
One microgram of total RNA was used to start the manual target preparation using the CodeLink Expression Bioarray System (Amersham Biosciences, Piscataway, NJ, http://www.amersham.com). Briefly, double-stranded cDNA synthesis was performed with a T7 oligo(dT) primer followed by purification. This cDNA was used as a template for in vitro transcription with biotin labeled nucleotides. Fifteen micrograms of the labeled cRNA were hybridized to Affymetrix rhesus macaque genome 49 format arrays (GeneChip Rhesus Macaque Genome Array, catalog number 900656; Affymetrix, Santa Clara, CA, http://www.affymetrix.com), followed by washing and staining with streptavidin phycoerythrin as recommended by the manufacturer. Arrays were scanned on an Affymetrix GeneChip 3000 Scanner. The arrays contain 52,303 probe sets that represent
30,000 human orthologs and expressed sequence tags (ESTs). Affymetrix GCOS software was used for the scanning of the probe arrays, and the probe intensity analysis and normalization were performed using RMAExpress [39].
Microarray Data Analysis. The gene expression analysis protocol has been described previously [40]. Statistical analysis was performed using the ScoreGene gene expression package (http://www.cs.huji.ac.il/labs/compbio/scoregenes), and data visualization was performed using Genomica (http://genomica.weizmann.ac.il) [41], Spotfire Decision Site 8.0 (Spotfire Inc., Göteborg, Sweden, http://spotfire.tibco.com), and TreeView (http://jtreeview.sourceforge.net). The RMA output for every gene was divided by the geometric mean of all the values for the same gene and was log base 2 transformed. In the analysis, we only included transcripts with locus link numbers. To determine the differentially expressed genes, we used t test or the nonparametric threshold number of misclassifications [42]. In the clustering, we included only genes that had a p value <.01 in both scoring methods and a ratio more than twofold in any of the pairwise comparisons. False discovery rate (FDR) analysis was carried out as described previously [43].
| RESULTS |
|---|
|
|
|---|
|
|
|
Growth characteristics of the isolated lines are variable but seemingly independent of sex, family relatedness, or passage number. Each line has a cell cycle time of approximately 20 hours, and individual cultures are passaged weekly. Cell lines isolated most recently (2406, 2706, 2906, 3006, 3106, and 3806) grow the most robustly regardless of passage number, suggesting perhaps that greater experience deriving ESC lines not only improves derivation efficiency but also quality of the newest lines.
The overall derivation success rate of these nhpESCs was 24.4% of the total number of blastocysts (Table 1) and 34.3% of those ICMs plated successfully (10 of 35). Importantly, the success rate was strongly dependent on the age of the blastocyst when immunosurgery was performed. We chose only fully expanded blastocysts for derivations. Embryos reached full expansion between days 9 and 11 of culture, counting the day of ICSI as day 1; 46.7% of the blastocysts that reached full expansion on day 9 successfully established nhpESC lines, whereas those that reached full expansion on day 10 resulted in stable nhpESC lines only 10.5% of the time, although no difference was observed in the plating efficiency between these two groups. These results were especially apparent when embryos were compared within the same clutch. Supplemental online Table 1 cites results comparing nhpESC lines derived from the same ovulation event within a single experiment where nhpESCs were established from day 9 blastocysts but in none of these cases from blastocysts at day 10. We also noted that Indian origin rhesus macaques (captive bred in the U.S. to other Indian rhesus) were more successful in nhpESC establishment compared with rhesus of Chinese origin (first generation imports from China).
|
31,000 annotations that represent all genes with a known gene ID or human orthologs and including ESTs (all raw data can be found at Geo National Center for Biotechnology Information accession number GSE7534). Each analysis was performed in triplicate except nhp206 (only two samples) and nhp2906 (single sample). To avoid the risk of introducing artifactual variability into the analysis, we chose to minimize variability among samples by preparing all samples using the same reagents at the same time. This precluded repeating nhp206 and nhp2906 in separate experiments. Since nhp2906 was represented only by a single sample, it was excluded from the average cluster analysis. However, all samples are included and compared in supplemental online Figure 1. We performed hierarchical clustering of the average gene expression using optimal leaf order analysis [44] and found that all of these nhpESC lines were extremely similar to each other in their gene expression (Fig. 3, black indicates no difference). As an indication of the consistency of this analysis, similar relationships were also observed when individual samples were clustered (supplemental online Fig. 1). To further quantitate the similarity, we calculated the correlation values among individual samples. The lowest value was 94.4% between samples 1a and 31c and was as high as 99% between many of the samples (supplemental online Table 3), indicating very little variability in the individual samples and confirming the similarity among these 10 ESC lines. Even among these highly similar lines, it was observed that full siblings had expression patterns more similar to each other than to their half siblings, indicating a familial influence on expression (Fig. 3).
This relationship can be seen more clearly in the heat map depicted in Figure 4. This figure was generated by using Genomica software and calculating the log ratio of the correlation among samples (supplemental online Table 2). As expected, samples from the same line are closely related (darker colors). Moreover, the similarity is observed not only among individual samples of the same line, but also among lines that are full siblings. Box A and Box C represent the comparisons of samples within a family; note the outlined dark area and the surrounding lighter areas indicating similarity within the family. Box B highlights the dissimilarity when comparing samples between the two families; note the light color within this box, clearly distinguishing full siblings from half siblings.
|
|
| DISCUSSION |
|---|
|
|
|---|
Considerable heterogeneity among the available hESC lines including growth rate and genetic and epigenetic stability and gene expression profiles has been reported [2]. These gene expression comparisons have primarily focused on "stemness" genes, that is, those expressed in all pluripotent or multipotent lines but not in differentiated or somatic cells [3, 49–54]. The first two studies [51, 52] identified approximately 250 putative genes involved in mESC pluripotency, and similar genes were identified in hESCs. Abeyta et al. [1] compared two female lines, HSF-6 and H9 (derived in different laboratories), with one male line, HSF-1, derived in the same laboratory as HSF-6. HSF-6 expression was more similar to HSF-1 than to the female line derived elsewhere, suggesting that the variability among lines might be due to derivation and culture differences. However, these interpretations are complicated by the unknown pedigrees even when comparing cell lines isolated within the same laboratory as well as causes of the underlying infertility of the patient/couple.
More recent comprehensive expression studies confirm variability [4, 5, 55] among hESC lines, including among lines generated within the same lab by the same personnel. Li et al. [25] generated five new lines and found that
10% of the probes were expressed in all five lines, although large variations among lines precluded direct line-by-line comparisons. Skottman et al. derived seven new hESC lines and, by comparing
10,000 genes, found
300 genes unique to each line. Interestingly, the gene expression profiles among four lines derived in Finland clustered more closely together than to their three lines derived in Sweden, raising questions regarding lab-to-lab variations as well as heredity. It will be interesting to determine whether similar relatedness in gene expression patterns is observed in related preimplantation embryos. Although not conducted in this report, methods exist to address this question in both humans [56, 57] and nonhuman primates [58]. This experiment would also be useful for determining whether the great similarity observed among these lines is related to extended in vitro culture under similar conditions. In mice, strain-specific differences in mESC gene profiles have been demonstrated [6, 7]. Taken together, these results suggest that knowledge of strain or family group origin might be essential for fully accurate interpretations. Access to infertility therapies in most countries is skewed by socioeconomic parameters, and it is likely that the donated embryos and resultant hESC lines are not representative of full demographic diversity. Consequently, even without direct knowledge of genetic origins, additional attention should be devoted to ensuring hESC research resources that are both well characterized and also inclusive of the population diversity.
Here, microarray comparisons within and between these nhpESC families show remarkably homogeneous gene expression profiles. It is important to highlight that these studies compared pluripotent ESCs with each other. Transcriptomic analysis of their progeny differentiated in vitro, in vivo in teratomas, and in utero in chimera is underway. The consistent and homogeneous gene expression of these nhpESCs (>96% identity) enables subtle analyses among ESCs that might otherwise be impossible, since the intrinsic variability described in hESCs would swamp minor differences. Recent transcriptional profiles of rhesus ESCs analyzed against human annotations reported an 85% [59] concordance. Our results achieved with rhesus ESCs are analyzed using the recently available rhesus gene annotations (http://www.affymetrix.com/products/arrays/specific/rhesus_macaque.affx). We found similar results when using the human annotations (not shown), reinforcing the reliability of our rhesus/rhesus results as well as the previous human/rhesus transcriptomics analysis [59]. Furthermore, analysis of average expression produced nearly identical results to comparisons of each of the individual triplicate samples, demonstrating the consistency of the Affymetrix rhesus chips and of the sample preparation. Additionally, the optimal leaf order analysis clustering precisely conforms with the experimentally designed parental origins and sex of each ESC line, providing another separate assurance of experimental accuracy beyond the independent genetic testing by the Veterinary Genetics Laboratory, University of California, Davis, on the genetic fidelity of each line reported here.
Sex-specific differences between the male and female cell lines are observed as a subset of genes overexpressed in males and underexpressed in females and a different subset overexpressed in female and underexpressed in male nhpESCs. Many but not all genes overexpressed in male lines are Y-chromosome linked, just as X-linked genes are overexpressed in the female lines. X inactivation is variable in pluripotent hESCs [60], perhaps explaining why not all sex linked genes are overexpressed in lines of their respective sex. Additionally, we analyzed gene expression using Affymetrix monkey arrays, but the annotations analyzed are from human homologs and there probably does not exist a 1:1 ratio between these two species. Finally, we do not expect all genes to be expressed at this stage of development, and nonexpressed genes would be excluded from the analysis. However, the expression profile of nhp906 (trisomy 15) is not strikingly different from the euploid ones, suggesting that there are stringent mechanisms to compensate for gene dosages in ESCs. With the establishment of primate ESC lines that have consistent genetic expression comes the ability to use these lines for detailed studies of infertility and development. ESC lines established using pedigreed primates could be used as bioassays to probe pluripotency and developmental differences among blastocysts of varying qualities, those generated with differing ART methods (e.g., in vitro fertilization, ICSI, or spermatid injections), or after environmental exposures during gametogenesis or preimplantation development. These abilities to extract small but biologically meaningful differences in gene expression profiles are possible because of the commonality of the background profile. Parallel, complementary studies with hESCs might possibly uncover novel pathways for the donors' underlying infertility, improved procedures for derivations, strategies to retain the fidelity of genomic imprints, and optimal growth and differentiation protocols.
Nearly identical expression profiles comparing different lines of either pluripotent mESCs or, as shown here, nhpESCs reinforce the reliability and biological utility of these research resources. Identifying the source or sources responsible for the newly recognized discordant profiling results with hESCs is important both for basic biology and also due to the clinical implications of hESC research findings. The existing federally approved hESC lines were derived from embryos discarded as suboptimal by the collaborating ART clinics and donated anonymously by the infertile patient/couple; therefore, no information is known regarding the relatedness of the existing lines. Also, the couple's infertility might have resulted from underlying genetic or epigenetic problems, further confounding the hESC line's status. Additionally, intrinsic defects might have resulted in the embryo's developmental delays, rendering them of limited clinical utility and therefore donated for fundamental studies. Finally, the lines were derived in labs distributed around the world and, thus, interlaboratory variability may play a role. hESC heterogeneity might be attributed to genetic variation due to parentage or sex, underlying causes of the couple's infertility or developmental delay of the embryo, variations in derivation and ESC establishment procedures and labs, and/or biological variations intrinsic to primate ESCs. Due to regulatory oversight policies, determination of these sources of hESC heterogeneity is not possible. In our study, embryos derived from rhesus monkeys of Indian origin more successfully produced ESC lines than did those from Chinese origin monkeys. These results suggest that hESC lines should be closely examined for ethnicity and utility parameters such as differentiation and transplantation, a set of experiments not possible with federal funding under the current guidelines.
Here, using embryos generated exclusively from gametes obtained from fertile rhesus, we find that nhpESC lines are more successfully generated from prime embryos versus subprime ones. This is somewhat surprising since the establishment of two-dimensional cell cultures would have been predicted to demand less stringency on the embryo versus the spatiotemporal challenges of implantation and pregnancy. Although neither pluripotency nor transcriptomic comparisons among the ESCs from the best embryos (lines 106, 206, 306, 2706, 2906, 3006, and 3106) and the subprime ones (lines 906, 2406, and 3806) show important differences, it is noteworthy that the one aneuploid line is from a subprime embryo (line 906). Furthermore, our most recent lines grow most consistently and reliably, suggesting that increased ESC experience improves the quality of the derived line and, by extrapolation, that newer lines are superior to older ones.
Questions remain as to exactly how similar primate ESCs are to mESCs, and both limitations in the availability of primate embryos and ESC lines as well as appropriate ethical restrictions on certain types of experiments makes it challenging to answer these questions with certainly. However, nhpESC studies continue to offer opportunities for bridging the gap between mouse and responsible human ESC research in several ways, including the following: pluripotency and germ line transmission assayed in intra- and interspecific chimeras; derivations from prime versus delayed embryos; lines generated with deliberate and defined epigenetic and genetic characteristics including mitochondrial DNA pedigrees [61, 62]; families, including inbred ones, of related ESCs generated and derived simultaneously versus sibling ESCs generated at differing times to analyze environmental influences; aging or transformation of the cells in vitro as studied by changes between batches of varying passage numbers; gametogenesis potentials and functional tests of gametes produced in vitro or after intra- and interspecific chimeras; genetic manipulation of nhp gametes or embryos such as transgenic, knockout/knockin, or small interfering RNA to produce more meaningful animal models of human disease; and transplantation investigations in which differentiated progeny from these cells are transplanted into full-sibling offspring generated by the identical rhesus monkeys or other related primates. These familial-related nhpESCs now also afford opportunities to explore changes with passaging and other environmental factors as well as sex specificity. However, with only a single decade of hESC and nhpESC research, and with understandable research constraints on some types of hESC investigations, comprehensive and rigorous comparisons between and among the primate ESC lines is only now underway. We suggest that ready access to more homogeneous and reliable primate ESC lines—human and nonhuman—will accelerate reliable findings of biomedical importance.
| CONCLUSION |
|---|
|
|
|---|
| DISCLOSURE OF POTENTIAL CONFLICTS OF INTEREST |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
|
|
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
F. Wianny, A. Bernat, C. Huissoud, G. Marcy, S. Markossian, V. Cortay, P. Giroud, V. Leviel, H. Kennedy, P. Savatier, et al. Derivation and Cloning of a Novel Rhesus Embryonic Stem Cell Line Stably Expressing Tau-Green Fluorescent Protein Stem Cells, June 1, 2008; 26(6): 1444 - 1453. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||