A molecular scheme for improved characterization of human embryonic stem cell lines
BMC Biology volume 4, Article number: 28 (2006)
Human embryonic stem cells (hESC) offer a renewable source of a wide range of cell types for use in research and cell-based therapies to treat disease. Inspection of protein markers provides important information about the current state of the cells and data for subsequent manipulations. However, hESC must be routinely analyzed at the genomic level to guard against deleterious changes during extensive propagation, expansion, and manipulation in vitro.
We found that short tandem repeat (STR) analysis, human leukocyte antigen (HLA) typing, single nucleotide polymorphism (SNP) genomic analysis, mitochondrial DNA sequencing, and gene expression analysis by microarray can be used to fully describe any hESC culture in terms of its identity, stability, and undifferentiated state.
Here we describe, using molecular biology alone, a comprehensive characterization of 17 different hESC lines. The use of amplified nucleic acids means that for the first time full characterization of hESC lines can be performed with little time investment and a minimum of material. The information thus gained will facilitate comparison of lines and replication of results between laboratories.
Human embryonic stem cells (hESC) are a potentially limitless, albeit controversial, source of therapeutic cells for numerous diseases and injuries. It is likely that different hESC lines are best suited to different uses, but at present, it is rare for any laboratory to work with more than a few lines. One reason is the expense of these lines; another is that in the USA hESC are governed by a dual-track policy. Cell lines derived prior to 9 August 2001 (currently about 20 available lines) can be examined using federal funds, and currently much of the available information on hESC biology has been generated using these funds and cell lines . However, cell lines derived after this date are far more numerous, but while it is legal to work on these lines using non-federal funds, information on their properties remains sparse. Government-funded researchers are reluctant to use these lines given the difficulties in accounting for federal and non-federal funds.
The lack of comparative analysis of hESC lines matters, because the properties and behavior of each line are uniquely shaped by their histories. It has become clear that different derivations produce hESC lines that are similar overall, but with inherent differences in gene expression, methylation status, X-chromosome inactivation, rate of self-renewal, and ability to differentiate [2–4]. More importantly, the behavior of cells and their phenotypic state changes as culture conditions and the stress to which they are subjected is altered, and permanent genomic changes frequently occur as passage numbers increase [5–7]. This has led to great difficulty in comparing results from one laboratory with another and even comparing results with different passages of the same cell line.
Therefore, thorough and routine characterization of hESC lines is essential to avoid compromising the validity of results. The most common characterization method for hESC is immunocytochemical analysis of a handful of markers, including SSEA-3, SSEA-4, TRA-1-60, TRA-1-80, and OCT-3/4 . The next most frequent is reverse transcription PCR, which is used for those group of genes whose expression is involved in maintenance of the undifferentiated state [9, 10]. While these assays certainly give indications of the undifferentiated state of the cells, they do not address other issues such as pluripotentiality or the degree of culture adaptation and genomic instability.
To facilitate comparisons among lines, the hESC research community has begun to develop a number of tools. Work is proceeding toward conditions that support the propagation of all lines , sets of markers that truly define the undifferentiated and unadapted state of the cells [7, 12–14], and markers predictive of the differentiation capacity of the cells . The work presented here is part of efforts to create a database of the properties of each line and to identify a reference standard for comparisons between laboratories. To this end, we have assembled a set of molecular tests for hESC lines that assess identity, stability of the nuclear and mitochondrial genomes, histocompatibility profile, and the undifferentiated state of the cells. Some of these assays have been previously performed on individual lines, but to our knowledge, no single group has used all of these tests on any single line, and few comparisons between lines have been made publicly available.
In this paper, we describe the analysis of multiple lines and show that this entire set of tests can be performed with a minimal sample size and over a short time period, and that these tests allow comparison of datasets across cell lines (a critical requirement to permit rapid progress in the field). We suggest that an internet database of hESC characterization data and standard reference materials will permit the research community to readily compare cell lines without the necessity of growing all lines in their own laboratory or tracking pre- and post-9 August cell lines.
The 17 cell lines analyzed are listed in Table 1. The lines were grown on feeders, or feeder-free on fibronectin, using bovine, human, or artificial serum according to the providers' protocols (see Methods).
Measures of identity
The results of short tandem repeat (STR) analysis of these lines are shown in Figure 1. The eight STR markers plus amelogenin used here are capable of identifying individual human genomes with a discrimination of 1 in hundreds of millions . Each of the early passage lines studied here was shown to be genetically distinct from the others. Each late passage line or karyotypically abnormal subline (BG01V, SA02) is identical to the early passage line, showing no evidence of overgrowth by contaminating cells. Only two instances of possibly heterogeneous cultures are observed. Cultures of the I6 line were grown in two separate laboratories, and one of these cultures shows a reproducible abnormality in the chromosome 13 microsatellite marker D13S317. One culture at passage 44 shows a third allele (with 12 repeats) for this marker; the low intensity of the signal suggests either an abnormal subpopulation or lessened PCR efficiency. The second culture, grown independently to passage 47, shows two normal alleles at all loci. The second abnormal observation was also found at the D13S317 locus in the line SA02. This line shows three distinct alleles at passage 29, consistent with trisomy 13 as shown in the karyotype . In contrast, the karyotypically normal subline SA02.5, which was subcloned from SA02, carries only two alleles on chromosome 13. The lines H7 and H9 show no changes visible by STR between lower and higher passage numbers.
The histocompatibility profile of each line was determined by sequence-specific oligonucleotide probes, and is summarized in Table 2. Each line has a unique histocompatibility profile, while the two sublines SA02.5 and BG01V are identical to their parental lines. Only the NTERA-2 embryonal carcinoma line shows an atypical histocompatibility profile, being homozygous at all tested loci.
Measures of stability
Finer resolution analysis of genomic stability is afforded by single nucleotide polymorphism (SNP) assays. In this analysis, over 100,000 SNPs were assayed by microarray. The average SNP spacing in the genome is 26 kb.
Figure 2 shows examples of hESC lines compared at low and high passage (H7 and H9), and comparison of karyotypically normal lines (BG01 and SA02.5) with karyotypically abnormal derivative lines (BG01V and SA02). SNP analysis can detect loss of heterozygosity (LOH) at any locus where two alleles were originally present, and is frequently used in analysis of tumors for discovery of potential oncogenes and anti-oncogenes . In addition, comparison of intensity signals can be used to calculate copy number at all loci [19, 20], and this has been recently used to detect genomic amplifications and deletions in hESC .
Comparison of the signal intensity and separation angle of all SNPs shows distinct discontinuities over certain chromosomes (Figure 2). This is the same approach used with Affymetrix SNP arrays to determine copy number ; however, reference BeadArray data from known diploid loci is required to translate intensities into copy number. Nevertheless, local changes in the allele intensity ratio (theta) are seen in chromosomes 12, 17, and a portion of chromosome X of BG01V, chromosome 13 of SA02, and within the q arm of chromosome 1 in H9. All of these areas have previously been shown to be duplicated in these cell lines . No other significant discontinuities were found for any other chromosome of any line.
Sequencing of the entire mitochondrial genome (16,544bp) was performed using an oligonucleotide sequencing microarray (Affymetrix MitoChip v2.0), as previously described . Briefly, 'early' and 'late' passages of nine paired hESC lines (BG01, BG02, BG03, HES2, HES3, SA01, SA02, H7, and H9) were sequenced for alterations in mtDNA in the course of in vitro passage. The results of this study have been published elsewhere , and six heteroplasmic sequence alterations occurring in two of nine (22%) later passage hESC lines were identified and confirmed by conventional dideoxy sequencing. The raw MitoChip data on the nine hESC lines is available on request (from AM).
Measures of pluripotency
Gene-expression analysis was performed by BeadArray (Illumina) to assay the transcriptome of each line. RNA expression profiles of several lines were directly compared with RNA from a pool of three undifferentiated lines previously used to compare expression by EST enumeration and MPSS [15, 22] (Figure 3). The normal lines tested (H9, I6, BG01, BG02, and BG03) show an average correlation coefficient (r 2) of 0.93 for genes detected with high confidence (>0.99). Karyotypically abnormal ES and EC lines (BG01V and NTERA-2) produce r 2 values of 0.92 and 0.88, respectively. This suggests that overall gene expression is more sensitive to the difference between normal and culture-adapted lines than is the use of immunocytochemical markers.
Not all genes contribute equally to the pluripotency of hESC, of course. Recent work comparing the transcriptomes of multiple hESC lines has identified core sets of 105 ES-specific genes  and 194 genes specifically upregulated in differentiating embryoid bodies , which have been used to create focused arrays for hESC testing. Of the ES-specific genes, 96 could be identified and detected at high confidence in our pooled hESC lines (Table 4). Nearly all of these were detected at similar levels in all of the individual hESC lines tested (H9, BG01, BG01V, BG02, BG03 and I6) and in the EC line NTERA-2. Every line but BG01 under- or over-expressed by 2.5-fold at least one ES marker, especially I6, which underexpressed 11 genes including TDGF-1, REX-1, and UTF-1. Few ES-specific genes were overexpressed, except for the variant line BG01V, which overexpressed four genes, including ABCG2, and NTERA-2, which overexpressed two genes, including the germ cell tumor marker GDF3.
Of the 192 differentiation-specific genes, 168 could be detected at a significant level in the pooled hESC lines. Approximately 89–96% of these were also detected in any individual line, except for BG01V, in which 134 genes were significantly expressed. This agrees with previous findings of low-level expression of diverse differentiation markers in undifferentiated hESC . NTERA-2 and I6 underexpressed the greatest number of these markers, while H9 and NTERA-2 also overexpressed six EB-specific genes each. In H9, three of these were collagen genes, suggesting the presence of human fibroblasts possibly differentiated from the hESC grown under feeder-free conditions.
The technical demands of culturing human embryonic stem cells are already well known in the field. It is now becoming clearer that frequent characterization of the cells in culture must be added to the burden. The present data confirm the general finding that hESC lines in culture accumulate genomic, epigenetic, and mitochondrial changes within as few as 22 passages . As a result, laboratories attempting to expand hESC lines must subject them to frequent characterization tests sensitive to aberrations below the resolution of G-banded karyotypes. The number and frequency of required tests also increases the value of rapid and convenient test methods, such as the many array-based molecular tests currently on the market or in development. This manuscript presents a molecular characterization scheme for hESC lines that helps to define the state of the cells. The tests include measures of identity, stability, and self-renewal, and the undifferentiated state of the cells. In addition, we have provided data on potential reference standards for hESC research. We show by examples that each of these tests offers important insights into the state of the hESC line and that the combination of these tests provides a unique characterization profile of a given cell line that is sensitive to the deleterious changes that can occur in culture.
This initiative, led by the American Type Culture Collection (ATCC) and supported by most of the providers of hESC lines, represents a complementary effort that uses cell lines approved by the National Institutes of Health to develop a basis for comparison across different laboratories. Equally important, by providing data in the same format on a new post-9 August line, it facilitates the direct comparison of lines eligible and ineligible for federal funding to determine which cell line is best suited to the individual researcher's needs. Such a comparison has not to our knowledge been previously performed, as it is largely impossible for most researchers in the USA.
Our estimate of testing costs for this characterization scheme is that these can be completed for a single line for approximately $2500. The identity, stability, and undifferentiated state of lines can be assessed with less than 5 μg of genomic DNA and 100 ng of total RNA. This can be obtained from 2–3 million cells, comparable with the amount required to run a single karyotype analysis using Giemsa-banded metaphase spreads. We further note that it will not be necessary for every researcher to repeatedly perform every one of these tests. We have identified commercial service providers both within and outside the USA, which can perform these analyses, and ATCC and other banks have offered to similarly characterize the lines they maintain.
The characterization scheme is composed of tests chosen on the basis of comprehensiveness, wide availability, and cost efficiency. For example, STR analysis as a measure of identity is a simple PCR-based analysis that is available in a kit from multiple providers, and can be readily run in any laboratory. One advantage of this assay is that it also readily distinguishes male and female samples, and in theory can be used to examine relatedness between samples and the larger database sets that have been generated for other purposes. ATCC is developing a comprehensive database of DNA profiles of hESC and other human cell lines based on STR loci . This test allows for discrimination of at least 1 in 100 million individuals .
HLA typing of hESC lines is essential, both for future use in transplantation and for the immediate goal of estimating how many hESC lines are required to meet the needs of all potential patients . HLA typing at the allelic level included HLA-A, HLA-B, HLA-C, HLA-DRB, and HLA-DQB loci.
SNP arrays are commonly used for mapping new markers of genetic disease and for detecting LOH in cancer. Recently, algorithms have been developed that allow oligonucleotide SNP arrays to measure chromosomal copy number at high resolution [19, 20]. This expands the utility of SNPs to detect non-reciprocal translocations, aneuploidy, and partial amplifications or deletions of chromosomes. Maitra et al have used this method to show that hESC in culture can acquire amplifications or deletions of small chromosomal regions . The SNP array method is of great value in characterizing hESC lines, as it improves the resolution for detecting genetic alterations. Conventional Giemsa-stained karyotypes or spectral karyotypes have a resolution of about 1–5 Mb. Comparative genomic hybridization (CGH) of metaphase spreads has a resolution of approximately 20 Mb, and array-based CGH with genomic or cDNA clones has improved this to <1 Mb. Bacterial artificial chromosome clone arrays have now improved the resolution to <100 kb. In contrast, currently available arrays of 100,000 SNPs have a potential resolution of about 25 kb. With over 2 million SNPs already mapped in the human genome, increased density arrays could drop this even further.
Like SNP genotyping, sequencing of the mitochondrial DNA offers the dual benefits of measuring identity (across 16,544 bases comprising a large number of polymorphic haplotypes) as well as stability over passages. Mitochondrial genes accumulate mutations in the somatic cells of aging animals. Mice in which mtDNA damage is accelerated show numerous symptoms of premature aging . The connections between tissue stem cells and aging are still controversial, but it has been demonstrated that changes in mitochondrial properties correlate with reduced competence of adult stem cells . Thus, although mtDNA mutations are quite common in somatic cells, the functional consequences may be greater in stem cells. Array-based sequencing of mtDNA is cost effective, and the reduced number of reactions needed in long-range PCR reactions and the automated analysis of genotype data using the RA Tools software (see Methods) significantly decreases analysis time.
Measures of identity and stability are generally of less concern to researchers growing hESCs than are the potency and undifferentiated state of the cells. These qualities are generally assessed by antibody staining for protein biomarkers. However, with the ability to rapidly assess expression of all human genes, it is possible to measure the expression of undifferentiated markers and any conceivable differentiation marker in a single step. The gene-expression profile of pluripotent hESCs has not been completely defined, but much work has been carried out toward defining common sets in undifferentiated and differentiating hESCs [12–14]. The use of this microarray method to distinguish undifferentiated and differentiated hESC cultures has been demonstrated . This gives the possibility of assigning a quality score to any hESC culture, which will be of value in improving culture methods.
The main drawback of these molecular methods is that the cells are assayed as an amalgam, so the ability to detect differentiated cells or aberrant subpopulations must be demonstrated. For example, we have estimated that at least 20% of the population must carry a genetic abnormality before it can be detected by a SNP array (data not shown). Currently, this method is often used to screen for LOH in cancer. Normal cells within the tumor severely degrade the sensitivity of LOH detection when as few as 10% are present, but copy-number detection is less sensitive to mixed populations . Similarly, the presence of low levels of differentiated cells is likely to escape notice when surveying global gene expression. The data presented here shows the overexpression of collagen genes in a sample of hESC cultured in feeder-free conditions, which suggests that the microarray method can detect differentiating subpopulations, but there is no way to retrospectively determine the level of differentiation.
Owing to the demands of culturing hESCs, it is rare that any one laboratory will bring together numerous lines to compare them under identical conditions. Thus, comparisons must be facilitated between laboratories by the promulgation of standard reference materials. The present work has thus compared multiple hESC lines to several potential reference materials, each having advantages in different areas of hESC research.
The greatest need for standards is in the area of gene-expression studies, as whole genome arrays are still expensive enough to discourage broad comparisons in any one laboratory, and the results are not easily compared between laboratories. For gene expression studies, we have created a pool of RNA from three independent hESC lines. This pool shows great concordance of expression with individual normal hESC lines, as the correlation coefficients (r 2) with normal lines average 0.87, or 0.93 for genes expressed at high confidence levels. Even if the line H9, one of the lines in the pooled samples, is excluded, the r 2 of highly confident expression for normal ES lines averages 0.93. Thus, the pooled sample shows high colinearity with hESC lines of diverse origin. As the RNA pool can be captured in the form of a cDNA library relatively inexpensively, it can easily be widely disseminated.
Alternatively, we have also provided data on embryonal carcinoma and karyotypically abnormal hESC lines, which are potential reference materials. As these lines are freely available, relatively inexpensive, and not subject to licensing and patent permissions, it is feasible to run these lines as internal/normalization controls. Internal reference standards also allow normalization of data to posted datasets available worldwide. NTERA-2 is an embryonal carcinoma line that expresses the same early embryo antigens now used to identify undifferentiated hESC. In global gene-expression studies, such as the data presented here, NTERA-2 closely mimics hESC lines, showing an average r 2 with the hESC pool of 0.84 for all genes and 0.88 for genes at high confidence. Reference material in the form of RNA or cDNA could easily be produced in large lots from this cell line, for distribution around the world. The NTERA-2 line also has advantages for antibody-based comparisons, as the line is well adapted to culture and grows quickly without feeder layers, attachment matrices, or the addition of any specialized growth factors other than serum. If the phenotype of NTERA-2 were compared between laboratories, this last would probably be the major source of variation. Thus, an effort could be made to introduce standardized serum-free conditions for the growth of NTERA-2 and validate the reproducibility of its gene expression across laboratories in order to create an unlimited source of reference material for hESC.
While the undifferentiated phenotype of NTERA-2 is quite similar to that of hESCs, its behavior in differentiation paradigms is quite different. It shows little spontaneous differentiation in culture without induction by retinoic acid or other chemicals [29, 30]. Therefore, it may be more advantageous to have a reference material that is grown in conditions identical to normal hESC and shows a similar capacity for differentiation in vitro  and in vivo . The aneuploid line BG01V has these characteristics. Compared with NTERA-2, the overall gene-expression profile of BG01V is more similar to undifferentiated hESC lines (r 2 = 0.87 for all genes, 0.92 for genes at high confidence). The gene expression of differentiating cultures is being examined, but is much more likely to mimic differentiating hESCs than NTERA-2 would. For whole-cell studies of differentiation, BG01V is available at low cost. The growth conditions for BG01V, as for hESC in general, have not yet been standardized to ensure reproducible results across laboratories.
In summary, we have found that the molecular characterization of hESC can be performed with minimal material, time, and expense. The series of tests presented here provides a method for frequent monitoring of many of the properties thought to be involved in stem-cell pluripotency. We believe that such characterization provides a useful set of results that will aid the cooperative international stem-cell effort, and will allow progress to be more rapid than in the past. We also hope that as cooperative standards become more widely accepted, commercial organizations will see fit to offer these as tools or kits to make the work of individual researchers easier.
Establishment and culture of human embryonic stem cell lines
All hESC lines were grown under similar conditions, with variations dictated by the providers' protocols. The standard culture medium was Dulbecco minimal essential medium (DMEM)/F12-Glutamax, 20% KnockOut Serum Replacement (KSR), 2 mM nonessential amino acids, 100 μM beta-mercaptoethanol (all from Invitrogen, Carlsbad, CA, USA), 50 μg/ml penicillin/streptomycin, and 4 ng/ml human recombinant basic fibroblast growth factor (FGF2; PeproTech Inc., Rocky Hill, NJ, USA) Feeder-free cultures were supplemented with FGF2 at 20 ng/ml.
TE06 (I6): Cells were cultured for 30 passages on mouse embryonic fibroblast (MEF) layers in KnockOut DMEM supplemented as above, except for the use of 5% fetal calf serum and 15% KSR (Invitrogen).
BG01, BG02, BG03 and BG01V: Cells were maintained under feeder-free condition on fibronectin-coated plates in medium that had been conditioned by mouse-embryo fibroblasts for 24 hours.
WA09 (H9): Human ES cell line H9 (WiCell, Madison, WI, USA) were cultured on feeder layers derived from mitotically inactivated WA09-derived mesenchyme or under feeder-free conditions on Matrigel- or human laminin-coated plates for at least 10 passages.
HUES-7: cells were maintained on a feeder layer of mitomycin C-treated CF-1 mouse embryonic fibroblasts (MEF) (ATCC, SCRC-1040.2). Growth medium was ES-DMEM (ATCC) supplemented with plasmanate (10%) (Bayer), KSR (10%; Invitrogen), L-alanyl-L-glutamine (2.0 mM; ATCC), non-essential amino acids (1X; ATCC), β-mercaptoethanol (0.055 mM; Invitrogen), penicillin (100 IU/mL) and streptomycin (100 μg/mL) (both ATCC), and FGF2 (12 ng/mL; R&D Systems, Inc.).
Cells were passaged every 4–5 days using collagenase IV (200 U/mL; Invitrogen) except HUES-7, which were passaged using 0.05% trypsin/EDTA (Invitrogen).
Frozen cell pellets were resuspended in 1× phosphate-buffered saline (PBS). A 20-mL aliquot was spotted onto a labeled FTA® card (Whatman) and allowed to dry. The FTA card lyses the cells on contact and binds the DNA to the paper surface. Prior to PCR, a portion of the dried spot was removed with a Harris punch, washed three times with purification reagent (Whatman), washed once with Tris-EDTA (TE) buffer (pH 8.0), and allowed to dry. STR analysis was conducted using the multiplex-PCR-based PowerPlex 1.2 kit (Promega Corporation). Loci analyzed include D5S818, D13S317, D7S820, D16S539, vWA, TH01, amelogenin, TP0X and CSF1P0. Electropherogram data were collected on an ABI 310 genetic analyzer (Applied Biosystems). Data was analyzed using Genescan 3.1 and Genotyper 2.0 software (Applied Biosystems). The resulting profile was imported into an in-house database and screened against all other baseline profiles of all samples tested by the ATCC. STR analysis of the NTERA-2 cell line was performed similarly using isolated genomic DNAs.
Genomic DNA was isolated from cells using the GenElute™ mammalian genomic DNA miniprep kit (Sigma-Aldrich, Inc.). HLA DNA typing was performed by utilizing hybridization of PCR-amplified DNA with sequence-specific oligonucleotide probes (SSOP) (Tepnel Lifecodes Corporation). The target DNA is amplified by PCR and then allowed to denature and rehybridize to complementary DNA probes conjugated to fluorescently coded microspheres. A flow analyzer identifies the fluorescent intensity on each microsphere, and the determined HLA type is based on the reaction pattern compared with patterns associated with public HLA gene sequences. Assays were performed to determine the HLA-A, HLA-B, HLA-C, HLA-DRB, and HLA-DQB loci.
Genomic DNA isolated as above was amplified following the Infinium™ whole genome genotyping assay  and hybridized to a prototype Sentrix BeadArray (Illumina Inc., San Diego, CA, USA). Allele calls were made by GenCall software (Illumina Inc.).
Gene expression analysis
Feeder-free cultures were prepared for gene expression analysis by manually harvesting individual colonies with uniform typical undifferentiated ES cell morphology. Feeder-containing cultures were harvested by scraping both ES and feeder cells. WA01 (H1), WA07 (H2), and WA09 (H9) cells were cultured under feeder-free conditions, and pooled prior to RNA isolation.
RNA was isolated from cultured cells using the Qiagen RNEasy kit (Qiagen, Inc, Valencia, CA, USA). Sample amplification was performed using 100 ng of total RNA as input material by the method of Van Gelder et al  using the Illumina® RNA amplification kit (Ambion Inc., Austin, TX, USA) following the manufacturer's instructions. Labeling was achieved by incorporating biotin-16-UTP (Perkin-Elmer Life and Analytical Sciences, Boston, MA, USA) present at a ratio of 1:1 with unlabeled UTP. Labeled, amplified material (700 ng per array) was hybridized to a pilot version of the Illumina Ref-8 BeadChip according to the manufacturer's instructions (Illumina Inc.). Arrays were washed then stained with Amersham FluoroLink streptavidin-Cy3 (GE Healthcare Bio-Sciences, Little Chalfont, Buckinghamshire, UK) following the BeadChip manual. Arrays were scanned with an Illumina BeadArray Reader confocal scanner according to the manufacturer's instructions. Array data processing and analysis was performed using Illumina BeadStudio software.
Mitochondrial DNA sequencing
Genomic DNA was extracted from paired hESC lines, and PCR amplified in three long-range PCR reactions as previously described . Amplified DNA was pooled, fragmented. and labeled as described in the Affymetrix CustomSeq™ resequencing protocol. MitoChips were hybridized overnight, washed on the Affymetrix fluidics station using the pre-programmed CustomSeq™ resequencing wash protocols, and scanned. Data analysis was performed using RA Tools, a modified version of the previously described adaptive background genotype-calling scheme (ABACUS) ; the open-source software is available from the Drosophila Population Genomics Project . Briefly, RA Tools uses an objective statistical framework to assign each genotype call a 'quality score', which is the difference between the log10 likelihood of the best-fitting and the second-best-fitting statistical model for assigning a genotype at any position on the sequencing array. The total quality-score threshold (totThresh) is the quality score that a given base has to exceed in order to be called. Increasing this value requires increased support for base calls, and consequently, fewer bases are called. Bases that fail to reach this threshold are called 'N.' The optimum total threshold quality score is determined empirically to be 12, which yields the highest base-call rate with the lowest discrepancy between genotypes for replicate samples. Confirmation of array-based sequencing was performed with a dye-terminator platform using the ABI Big Dye cycle sequencing kit (Applied Biosystems).
National Institutes of Health: NOT-OD-02-005. Notice of criteria for federal funding of research on existing human embryonic stem cells and establishment of NIH human embryonic stem cell registry. Released 7 November. 2001, [http://grants.nih.gov/grants/guide/notice-files/NOT-OD-02-005.html]
Abeyta MJ, Clark AT, Rodriguez RT, Bodnar MS, Pera RA, Firpo MT: Unique gene expression signatures of independently-derived human embryonic stem cell lines. Hum Mol Genet. 2004, 13: 601-608. 10.1093/hmg/ddh068.
Hoffman LM, Hall L, Batten JL, Young H, Pardasani D, Baetge EE, Lawrence J, Carpenter MK: X-inactivation status varies in human embryonic stem cell lines. Stem Cells. 2005, 23: 1468-1478. 10.1634/stemcells.2004-0371.
Skottman H, Mikkola M, Lundin K, Olsson C, Stromberg AM, Tuuri T, Otonkoski T, Hovatta O, Lahesmaa R: Gene expression signatures of seven individual human embryonic stem cell lines. Stem Cells. 2005, 23: 1343-1356. 10.1634/stemcells.2004-0341.
Draper JS, Smith K, Gokhale P, Moore HD, Maltby E, Johnson J, Meisner L, Zwaka TP, Thomson JA, Andrews PW: Recurrent gain of chromosomes 17q and 12 in cultured human embryonic stem cells. Nat Biotechnol. 2004, 22: 53-54. 10.1038/nbt922.
Maitra A, Arking DE, Shivapurkar N, Ikeda M, Stastny V, Kassauei K, Sui G, Cutler DJ, Liu Y, Brimble SN, Noaksson K, Hyllner J, Schulz TC, Zeng X, Freed WJ, Crook J, Abraham S, Colman A, Sartipy P, Matsui SI, Carpenter M, Gazdar AF, Rao M, Chakravarti A: Genomic alterations in cultured human embryonic stem cells. Nat Genet. 2005, 37: 1099-1103. 10.1038/ng1631.
Enver T, Soneji S, Joshi C, Brown J, Iborra F, Orntoft T, Thykjaer T, Maltby E, Smith K, Dawud RA, Jones M, Matin M, Gokhale P, Draper J, Andrews PW: Cellular differentiation hierarchies in normal and culture-adapted human embryonic stem cells. Hum Mol Genet. 2005, 14: 3129-3140. 10.1093/hmg/ddi345.
Draper JS, Pigott C, Thomson JA, Andrews PW: Surface antigens of human embryonic stem cells: changes upon differentiation in culture. J Anat. 2002, 200: 249-258. 10.1046/j.1469-7580.2002.00030.x.
Niwa H, Miyazaki J, Smith AG: Quantitative expression of Oct-3/4 defines differentiation, dedifferentiation or self-renewal of ES cells. Nat Genet. 2000, 24: 372-376. 10.1038/74199.
Noaksson K, Zoric N, Zeng X, Rao MS, Hyllner J, Semb H, Kubista M, Sartipy P: Monitoring differentiation of human embryonic stem cells using real-time PCR. Stem Cells. 2005, 23: 1460-1467. 10.1634/stemcells.2005-0093.
Rao BM, Zandstra PW: Culture development for human embryonic stem cell propagation: molecular aspects and challenges. Curr Opin Biotechnol. 2005, 16: 568-576. 10.1016/j.copbio.2005.08.001.
Yang AX, Mejido J, Luo Y, Zeng X, Schwartz C, Wu T, Thies RS, Bhattacharya B, Han J, Freed B, Rao M, Puri RK: Development of a focused microarray to assess human embryonic stem cell differentiation. Stem Cells Dev. 2005, 14: 270-284. 10.1089/scd.2005.14.270.
Bhattacharya B, Miura T, Brandenberger R, Mejido J, Luo Y, Yang AX, Joshi BH, Ginis I, Thies RS, Amit M, Lyons I, Condie BG, Itskovitz-Eldor J, Rao MS, Puri RK: Gene expression in human embryonic stem cell lines: unique molecular signature. Blood. 2004, 103: 2956-2964. 10.1182/blood-2003-09-3314.
Bhattacharya B, Cai J, Luo Y, Miura T, Mejido J, Brimble SN, Zeng X, Schulz TC, Rao MS, Puri RK: Comparison of the gene expression profile of undifferentiated human embryonic stem cell lines and differentiating embryoid bodies. BMC Dev Biol. 2005, 5: 22-10.1186/1471-213X-5-22.
Miura T, Luo Y, Khrebtukova I, Brandenberger R, Zhou D, Thies RS, Vasicek T, Young H, Lebkowski J, Carpenter MK, Rao MS: 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-715. 10.1089/scd.2004.13.694.
Lins AM, Micka KA, Sprecher CJ, Taylor JA, Bacher JW, Rabbach DR, Bever RA, Creacy SD, Schumm JW: Development and population study of an eight-locus short tandem repeat (STR) multiplex system. J Forensic Sci. 1998, 43: 1168-1180.
NIH Human Embryonic Stem Cell Registry. [http://stemcells.nih.gov/research/registry/]
Zhou X, Rao NP, Cole SW, Mok SC, Chen Z, Wong DT: Progress in concurrent analysis of loss of heterozygosity and comparative genomic hybridization utilizing high density single nucleotide polymorphism arrays. Cancer Genet Cytogenet. 2005, 159: 53-57. 10.1016/j.cancergencyto.2004.09.014.
Zhao X, Li C, Paez JG, Chin K, Janne PA, Chen TH, Girard L, Minna J, Christiani D, Leo C, Gray JW, Sellers WR, Meyerson M: An integrated view of copy number and allelic alterations in the cancer genome using single nucleotide polymorphism arrays. Cancer Res. 2004, 64: 3060-3071. 10.1158/0008-5472.CAN-03-3308.
Nannya Y, Sanada M, Nakazaki K, Hosoya N, Wang L, Hangaishi A, Kurokawa M, Chiba S, Bailey DK, Kennedy GC, Ogawa S: A robust algorithm for copy number detection using high-density oligonucleotide single nucleotide polymorphism genotyping arrays. Cancer Res. 2005, 65: 6071-6079. 10.1158/0008-5472.CAN-05-0465.
Huang J, Wei W, Zhang J, Liu G, Bignell GR, Stratton MR, Futreal PA, Wooster R, Jones KW, Shapero MH: Whole genome DNA copy number changes identified by high density oligonucleotide arrays. Hum Genomics. 2004, 1: 287-299.
Brandenberger R, Khrebtukova I, Thies RS, Miura T, Jingli C, Puri R, Vasicek T, Lebkowski J, Rao M: MPSS profiling of human embryonic stem cells. BMC Dev Biol. 2004, 4: 10-10.1186/1471-213X-4-10.
Zipori D: The nature of stem cells: state rather than entity. Nat Rev Genet. 2004, 5: 873-878. 10.1038/nrg1475.
ATCC: STR Database for Human Cell Lines. [http://www.atcc.org/common/cultures/str.cfm]
Taylor CJ, Bolton EM, Pocock S, Sharples LD, Pedersen RA, Bradley JA: Banking on human embryonic stem cells: estimating the number of donor cell lines needed for HLA matching. Lancet. 2005, 366: 2019-2025. 10.1016/S0140-6736(05)67813-0.
Trifunovic A, Wredenberg A, Falkenberg M, Spelbrink JN, Rovio AT, Bruder CE, Bohlooly YM, Gidlof S, Oldfors A, Wibom R, Tornell J, Jacobs HT, Larsson NG: Premature ageing in mice expressing defective mitochondrial DNA polymerase. Nature. 2004, 429: 417-423. 10.1038/nature02517.
Lonergan T, Brenner C, Bavister B: Differentiation-related changes in mitochondrial properties as indicators of stem cell competence. J Cell Physiol. 2006
LaFramboise T, Weir BA, Zhao X, Beroukhim R, Li C, Harrington D, Sellers WR, Meyerson M: Allele-specific amplification in cancer revealed by SNP array analysis. PLoS Comput Biol. 2005, 1: e65-10.1371/journal.pcbi.0010065.
Andrews PW: Retinoic acid induces neuronal differentiation of a cloned human embryonal carcinoma cell line in vitro. Dev Biol. 1984, 103: 285-293. 10.1016/0012-1606(84)90316-6.
Andrews PW, Gonczol E, Plotkin SA, Dignazio M, Oosterhuis JW: Differentiation of TERA-2 human embryonal carcinoma cells into neurons and HCMV permissive cells. Induction by agents other than retinoic acid. Differentiation. 1986, 31: 119-126. 10.1111/j.1432-0436.1986.tb00392.x.
Zeng X, Chen J, Liu Y, Luo Y, Schulz TC, Robins AJ, Rao MS, Freed WJ: BG01V: a variant human embryonic stem cell line which exhibits rapid growth after passaging and reliable dopaminergic differentiation. Restor Neurol Neurosci. 2004, 22: 421-428.
Gunderson KL, Steemers FJ, Lee G, Mendoza LG, Chee MS: A genome-wide scalable SNP genotyping assay using microarray technology. Nat Genet. 2005, 37: 549-554. 10.1038/ng1547.
Van Gelder RN, von Zastrow ME, Yool A, Dement WC, Barchas JD, Eberwine JH: Amplified RNA synthesized from limited quantities of heterogeneous cDNA. Proc Natl Acad Sci USA. 1990, 87: 1663-1667. 10.1073/pnas.87.5.1663.
Maitra A, Cohen Y, Gillespie SE, Mambo E, Fukushima N, Hoque MO, Shah N, Goggins M, Califano J, Sidransky D, Chakravarti A: The Human MitoChip: a high-throughput sequencing microarray for mitochondrial mutation detection. Genome Res. 2004, 14: 812-819. 10.1101/gr.2228504.
Cutler DJ, Zwick ME, Carrasquillo MM, Yohn CT, Tobin KP, Kashuk C, Mathews DJ, Shah NA, Eichler EE, Warrington JA, Chakravarti A: High-throughput variation detection and genotyping using microarrays. Genome Res. 2001, 11: 1913-1925.
Drosophila Population Genomics Project. [http://www.dpgp.org/]
Brimble SN, Zeng X, Weiler DA, Luo Y, Liu Y, Lyons IG, Freed WJ, Robins AJ, Rao MS, Schulz TC: Karyotypic stability, genotyping, differentiation, feeder-free maintenance, and gene expression sampling in three human embryonic stem cell lines derived prior to August 9, 2001. Stem Cells Dev. 2004, 13: 585-597. 10.1089/scd.2004.13.585.
Amit M, Itskovitz-Eldor J: Derivation and spontaneous differentiation of human embryonic stem cells. J Anat. 2002, 200: 225-232. 10.1046/j.1469-7580.2002.00032.x.
Thomson JA, Itskovitz-Eldor J, Shapiro SS, Waknitz MA, Swiergiel JJ, Marshall VS, Jones JM: Embryonic stem cell lines derived from human blastocysts. Science. 1998, 282: 1145-1147. 10.1126/science.282.5391.1145.
Reubinoff BE, Pera MF, Fong C, Trounson A, Bongso A: Embryonic stem cell lines from human blastocysts: somatic differentiation in vitro. Nature Biotechnology. 2000, 18: 399-404. 10.1038/74447.
Heins N, Englund MC, Sjoblom C, Dahl U, Tonning A, Bergh C, Lindahl A, Hanson C, Semb H: Derivation, characterization, and differentiation of human embryonic stem cells. Stem Cells. 2004, 22: 367-376. 10.1634/stemcells.22-3-367.
Mandal A, Tipnis S, Pal R, Ravindran G, Bose B, Patki A, Rao MS, Khanna A: Characterization and In Vitro Differentiation Potential of a New Human Embryonic Stem Cell Line, Relicell®hES1. Differentiation. 2006, 74: 81-90. 10.1111/j.1432-0436.2006.00051.x.
Cowan CA, Klimanskaya I, McMahon J, Atienza J, Witmyer J, Zucker JP, Wang S, Morton CC, McMahon AP, Powers D, Melton DA: Derivation of embryonic stem-cell lines from human blastocysts. N Engl J Med. 2004, 350: 1353-1356. 10.1056/NEJMsr040330.
This study was supported by a National Institute on Aging contract (N01AG40002; Characterization of Human Embryonic Stem Cell Lines). We thank Dr Doug Melton for providing the HUES-7 line and Dr David Barker and Illumina Inc. for providing whole-genome genotyping data. AM received support from the Sol Goldman Pancreatic Cancer Research Center and the Maryland Cigarette Restitution Fund.
Support for the analysis for non-federally approved lines was obtained by ATCC Manassas, VA, USA. NIH employees listed as contributors were responsible solely for the supply of federally approved lines, and their contribution was limited to the analysis of lines derived pre- 9 August 2001.
GS performed the STR analysis of all lines. YL, SS, and JL performed gene expression analysis. CO, WX, and XZ grew hESC lines for analysis and provided materials to outside testers. RJ analyzed the genotyping data. AM developed and performed the mitochondrial-DNA sequencing. MSR and JMA developed the characterization scheme and coordinated all materials and tests. RJ and MSR wrote the manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Josephson, R., Sykes, G., Liu, Y. et al. A molecular scheme for improved characterization of human embryonic stem cell lines . BMC Biol 4, 28 (2006). https://doi.org/10.1186/1741-7007-4-28
- Human Leukocyte Antigen
- Short Tandem Repeat
- Undifferentiated State
- Single Nucleotide Polymorphism Array
- hESC Line