|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
a Max-Planck Institute for Molecular Genetics (Department of Vertebrate Genomics), Berlin, Germany;
b Reproduction and Early Development Research Group, Department of Obstetrics and Gynaecology, University of Leeds, Leeds, W. Yorkshire, United Kingdom;
c Center for Reproductive Medicine & Infertility, Weill Medical College, Cornell University, New York, New York, USA
Key Words. Preimplantation development • Blastocyst • Inner cell mass • Trophectoderm • Pluripotency • Embryonic stem cells • Trophoblastic stem cells • OCT4 • CDX2 • Differentiation • Microarrays
Correspondence: James Adjaye, Ph.D., Max-Planck Institute for Molecular Genetics (Department of Vertebrate Genomics), Ihnestrasse 73, D-14195 Berlin, Germany. Telephone: 0049-30-8413-1216; Fax: 0049-30-8413-1128; e-mail: adjaye{at}molgen.mpg.de
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Several stemness genes have been identified in murine [47] and human [812] embryonic stem cell lines under culture conditions. However, the weak overlap between these genes (~25%) may be a consequence of the heterogeneous genetic backgrounds of the various cell lines or epigenetic and chromosomal instability after differentiation in vitro [13, 14]. Hence, genes for stemness may be better identified at source within the pluripotent ICM of the mammalian blastocyst.
Although specific markers for the TE and ICM, such as CDX2 and OCT4 [1517], have been identified, these are limited and underexplored. In this study, we provide the first global analysis of differential gene expression in intact human blastocysts and immunosurgically isolated ICM and TE cells for identifying novel marker genes specific to these lineages. We have uncovered gene expression patterns that suggest that the delineation of the ICM and TE from the morula may be controlled by signaling pathways, which are also crucial for determining cell fate later in embryonic development, in the ontogenesis of some cancers, and for advances in technology and the maintenance of human ES and TS cells in vitro.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Preimplantation Embryos and Generation of T7-Linked Double-Stranded cDNA and Amplified RNA
Preimplantation embryos surplus to requirement for in vitro fertilization (IVF) treatment were donated by patients for research attending the Jones Institute at Eastern Virginia Medical School (EVMS) (Norfolk, VA) and Leeds General Infirmary (LGI) (Leeds, U.K.). The IVF-derived blastocysts were similar to those used for deriving ES cells, i.e., fully expanded with a well-defined ICM and TE, lacking signs of degeneration, and at stage 3 according to the grading system described [19].
All samples were obtained after informed consent under protocols approved by the Institutional Review Board at EVMS and Research Ethics Committee of the LGI and licensed by the Human Fertilisation and Embryology Authority (U.K.). Permission was also granted by the ethics commission of the Free University Berlin, allowing the use of the generated RNA samples in Germany. Samples were washed in phosphate-buffered saline, lysed, and stored at 80°C in 50 µl Dynal lysis buffer (Dynabeads, Dynal Biotech, Bromborough, Wirral, UK, http://www.dynalbiotech.com) supplemented with 10% RNAlater (Ambion, Austin, TX, http://www.ambion.com).
cDNA samples were derived from duplicate intact blastocysts. The TE and the ICM of the human blastocyst were generated as follows: 10 pronucleate embryos were cryopreserved and, after thawing, cultured for up to 6 days in sequential medium (Irvine Scientific, Santa Anna, CA, http://www.irvinesci.com). At the blastocyst stage, the embryos underwent immunosurgery to separate the TE and ICM [20]. Unhatched blastocysts were exposed briefly to acidified Tyrodes medium (Irvine Scientific) for removing the zona pellucida. They were exposed to a rabbit antiserum raised against the BeWo trophoblastic cell line (Atlantic Antibodies, Windham, ME) and subsequently incubated with guinea pig complement (GIBCO, Grand Island, NY, http://www.invitrogen.com). Degenerated TE cells were separated from the ICM by pipetting through a finely drawn Pasteur pipette. The isolated ICM cells were washed three times in fresh medium and checked microscopically for the absence of adhering TE cells before being added to a tube containing 50 µl of RNAlater. The lysed TE cells, together with the microdroplet of medium, were also collected, washed, and transferred to separate tubes for analysis. Samples were stored for 1 month at 80°C before shipping to the U.K.
cDNA Amplification and In Vitro Transcription
T7 promoterlinked double-stranded cDNA samples derived from ICM, TE, and intact blastocysts were generated according to published protocols [21]. Briefly, mRNA was extracted from thawed lysed cells using Oligo-dT magnetic beads (Dynabeads). cDNA was generated using T7 promoterlinked oligo-dT primers for the reverse transcription (RT) step, and whole-transcriptome amplification was executed using a modified SMART amplification protocol (BD Biosciences, San Jose, CA, http://www.bdbiosciences.com). Amplified RNA was generated using the MegaScript T7 High Yield Transcription kit (Ambion). RNA purity, integrity, and concentrations were evaluated on the Agilent 2100 Bioanalyzer.
Direct Labeling of RNA and Hybridizations
Minimum Information About Microarray Experiments (MIAME) guidelines were adhered to in our experimental design [22].
Four independent labeling (including dye swaps, i.e., Cy3 and Cy5) reactions per an RNA sample (ICM, TE, and intact blastocyst) were carried out using 3 µg aRNA per reaction. Direct incorporation of Cy3 and Cy5 during RT was carried out in a 20-µl reaction volume using 1 µg of a random hexamer primer. Purification of labeled cDNAs, hybridizations, slide washes, and capturing of fluorescence images are described in detail in the online supplemental methods.
Data Analysis
Data normalization was carried out as described previously [23]. Gene expression was judged according to the comparison with a negative control sample (full details in the online supplemental methods).
For each cDNA, we performed statistical tests based on the replicate signals in experiments with ICM and TE samples. This was done for the two blastocysts separately. Three standard tests were used in parallelStudents t-test, the Welch test, and Wilcoxon rank sum test [24]. To evaluate differential expression of the genes, p values of Wilcoxon rank sum test were preferred as a reference because this test does not depend for its validity on a specific distribution (e.g., Gaussian). A recursive function was implemented for calculating the exact p values.
Pathway Analysis
Array data were used to test whether specific pathways showed differential expression. Pathways were taken from the KEGG database. Consider for each pathway i the set of related genes (xil, yil),...,(xini, yini). Here, xij, yij denotes the expression level of the j-th gene in TE and ICM, respectively. Wilcoxon matched-pairs signed rank test was used to calculate a Z-score for the differences dij = xij yij for each pathway i. These differences were ranked, and the ranks of differences with negative signs,Rneg, and those with positive signs,Rpos, were summed. The test statistic is the smaller of the two numbers, R = min{Rpos, Rneg}. If the pathway is not affected by the treatments, Rneg and Rpos will be fairly equal, but if there is a trend of underexpression or overexpression, the test statistic will be small. The Z-score is defined as
![]() |
where E is the expectation and Var is the variance of R. These were calculated as
![]() |
respectively [25].
Real-Time Reverse TranscriptionPolymerase Chain Reaction Analysis
This was carried out for a set of genes on the ABI PRISM 7900HT Sequence Detection System (Applied BioSystems, Foster City, CA, http://www.appliedbiosystems.com). Full details of the experimental set up and list of primers, annealing temperatures, and genes under investigation are shown in supplemental online Table 10 and in the online supplemental methods.
Functional Annotation of Expressed Genes
We used the gene ontology (GO) terminology related to Ensembl version 3.0 to characterize the following subsets of our expression data: All, all genes on the array; TE, all genes highly expressed in the TE (>0.9); ICM, all genes highly expressed in the ICM (>0.9),and BL, all genes highly expressed in the blastocyst (>0.9).
The GO vocabulary was taken from the Gene Ontology Web site in May 2004 (http://www.geneontology.org). This was imported into the sqlite database (http://www.sqlite.org). Data analysis was carried out using the R-Statistics software (http://www.r-project.org). We have developed a web application on our Goblet webserver (http://goblet.molgen.mpg.de/blastocyst) that allows users to view genes based on their expression pattern or on their GO annotation.
| RESULTS |
|---|
|
|
|---|
|
Expression Profiling Distinguishes ICM from TE
ICM and TE cells were isolated from two blastocysts from different individuals. Four independent hybridization experiments were performed for each biological replicate, including Cy-dye swaps. Additionally, RNA from two blastocysts was pooled to generate a reference sample. An overview of comparative gene expression in ICM and TE cells is shown in Figure 2A
and reveals as expected a high overall correlation between the data (0.90), which helps to validate them. To judge whether a given gene is expressed in these cell types, we compared its signal against a negative control sample and computed a numerical value to judge gene expression (BG-tag; see Materials and Methods). This number reflects the proportion of background noise in relation to the actual signal [23]. Typically, a BG-tag of 0.9 indicates a detectable signal for the probe (Fig. 2B
). Using this criterion, we found that 7,481 (48%) genes represented on the chip (probes) were detected in one of the three cell types (ICM, TE, or blastocyst). As demonstrated in Figure 2B
, most of these genes are either specific to the intact blastocyst (2,880) or common to all three cell types (2,031). The number of genes that are expressed exclusively in the immunosurgically isolated ICM or TE is rather low (292 and 345, respectively).
|
These newly identified genes may complement existing markers for ICM and TE, such as OCT4 and CDX2 (Figs. 1B
, 3A, 3B
) and represent diverse biological functions. For example, the TE marker SFN (stratifin) is an epithelial cell antigen, which is exclusively expressed in keratinocytes. Its role in cell proliferation and apoptosis suggests that the protein could be relevant to the regulation of growth and differentiation of multiple cell types through the protein kinase C signaling pathway [35]. Among the putative ICM markers, HMGB1 and GLTSCR2 have been identified as potential stemness genes based on expression studies in the human ES lines HSF-1, HSF-6, and H9 [10]. HMGB1 is a member of the high-mobility group of transcription factorencoding proteins that act primarily as architectural facilitators in the assembly of nucleoprotein complexes, as in the initiation of transcription of target genes. Murine Hmgb1 has been shown to be a coactivator of Oct4 [36]. GLTSCR2 is a gene of unknown function residing in the glioma tumor suppressor region of chromosome 19q. Of the proposed list of stemness genes [4, 5, 12], only ITGA6 (integrin alpha 6 chain) appeared in all the stem cell lines analyzed, but by comparing this list with our ICM markers, we have also identified ITGA5 as another candidate for stemness. These findings highlight the potential importance of integrins in cell adhesion as well as in cell surfacemediated signaling for establishing and maintaining pluripotency.
|
The full sets of data for global gene expression in the duplicate ICMs, TEs, and pooled blastocysts are presented (supplemental online Table 1). In general, the magnitude of expression recorded after RT-PCR was greater than from microarrays, an observation that is consistent with previous findings [4, 26].
Functional Annotation of Expressed Genes
The GO vocabulary provides a unified terminology for the description of genes and their products [38]. It is divided into three main categories: Molecular Function, Biological Process, and Cellular Component. A comparative analysis of these three categories at a global level within the ICM, TE, and intact blastocyst did not reveal a bias toward a particular category in these cell types (data not shown). In contrast, by repeating this analysis on Molecular Function using the 78 marker genes, 51 of which have GO annotations, we observed a slight bias toward specific molecular functions within ICM and TE cells (Fig. 4
). For example, the ribosomal proteins RPL14, RPL7A, RPL19, and RPL32, under structural molecular activity (GO:0005198), are all expressed in the ICM (Fig. 4A
). A more detailed description of these markers with respect to their chromosomal localization and ontology is presented (supplemental online Table 2). For a global overview, we combined the expression data from the ICM, TE, and intact blastocysts to create a database for searching for expression levels and related GOs (http://goblet.molgen.mpg.de/blastocyst).
|
Delineation of signaling pathways will be fundamental for understanding the mechanisms regulating pluripotency and self-renewal in cultured ES cell lines. We searched the ICM and TE data for components of these pathways by assigning p values using Wilcoxon matched-pair signed rank test, as described in Materials and Methods. This strategy is distinct from the commonly used strategy for identifying differentially expressed genes using repeated measurements with a two-sample location test, such as Students t-test or Wilcoxon rank sum test, because we directly involve groups of genes associated with particular pathways instead of conventional gene-wise analysis. The data indicate the involvement of WNT, mitogen-activated protein kinase (MAPK), transforming growth factor-ß (TGF-ß)/bone morphogenic protein (BMP), NOTCH, integrin-mediated cell adhesion, apoptosis-signaling pathways, and metabolic processes such as glycolysis, sterol biosynthesis, androgen, and estrogen metabolism. The full list of signaling and metabolic pathways identified by these analyses is in Table 1
and supplemental online Table 3. Pathway annotations were adopted from the KEGG (Kyoto Encyclopedia of Genes and Genomes; http://www.genome.jp/kegg) database.
|
Apoptosis in the Mammalian Blastocyst
Regulation of cell population size and lineage determination is mediated by cell cycle control, differentiation, and programmed cell death or apoptosis. The latter is characterized by chromatin condensation, nuclear membrane blebbing, and fragmentation in the cytoplasm and nucleus [40].
Apoptosis is evident at the blastocyst stage, if not earlier. It is mainly restricted to the ICM to regulate the size of the cell mass and perhaps to eliminate cells retaining the potential to form TE ectopically [41]. A list of expressed genes involved in the apoptosis signaling pathway is provided in supplemental online Table 4.
NOTCH Signaling in the Blastocyst
Cellcell signaling mediated by the Notch receptor determines cell fate and regulates pattern formation in many phyla. The Notch signaling pathway is operational in both the ICM and TE and is known to function in the mouse blastocyst [27]. We observed expression of several ligands and receptors in this pathway (supplemental online Table 5), consistent with a mouse model in which Wnt and Notch act sequentially to set up the initial asymmetry in the zygote and thus influencing polarity of the developing embryo [27].
TGF-ß Signaling
The TGF-ß family consists of multifunctional growth and differentiation factors regulating many cellular processes through complex signal-transduction pathways. The family members include TGF-ß isoforms, activins, and BMPs. Expression of the signaling type I and type II receptors for TGF-ß in mouse and human fertilized oocytes and blastocysts suggested a role for TGF-ß in early preimplantation development, potentially in the outgrowth of parietal endoderm [42]. Differential expression of TGF-ß isoforms, activins, BMPs, and MADHs/SMADs was also evident. In particular, BMP4, previously shown to induce the differentiation to trophoblast when overexpressed in human ES cells [43], is 2.28-fold enriched in the TE. Other components of the TGF-ß signaling cascade are shown (supplemental online Table 6).
Integrin and Cadherin-Mediated Cell Adhesion
The ICM and TE originate from the division of polar blastomeres when their cleavage furrows parallel their apical surfaces. These blastomeres polarize in response to asymmetric cell-cell contact. Pathway analysis identified signaling pathways related to integrin-mediated cell adhesion. In addition, several Na+/K+-ATPases (e.g., ATP1B3; Fig. 3B
) were overexpressed in TE, reflecting their presumptive roles in driving fluid transport across this epithelium. Other cell adhesionrelated genes were also detected, as expected where intercellular junctions are important for controlling blastocyst permeability [44]. However, there was overexpression in the TE of a subset of these genes, including Desmocollin 2 (DSC2 x1.55), Protocadherins (PCDH7 x1.67, PCDH11 x1.57, PCDHB7 x1.62), E-cadherins (CDH19 x1.9, CDH24 x1.54, CDH22 x1.82), tight junction proteins (TJP1 x1.4, TJP2 x1.8), Claudins (CLDN2 x1.4, CLDN16 x1.79, CLDN10 x2.25), and seven-pass transmembrane receptor of the cadherin superfamily (CELSR2 x1.46). For the tight junction constituents OCLN (occludins), JAM-2, (Junction adhesion molecule 2), and CGN (Cingulin), a lack of overexpression in the TE may be due to the fact that translational rather than transcriptional control is operative due to cell contact symmetry [44]. A comprehensive listing of these genes and their expression ratios is given in supplemental online Tables 1 and 7.
WNT Signaling in the Blastocyst
The WNT gene family consists of numerous conserved glycoproteins that regulate pattern formation during embryogenesis in a wide variety of tissues, including the nervous system. It has recently been shown that activation of the canonical WNT/Wnt pathway is sufficient to maintain self-renewal of both human and mouse ES cells [45] and also that this pathway is operative during human and mouse preimplantation development [21, 27]. We detected differential expression of transcripts encoding WNT ligands (WNTs), receptors of the Frizzled gene family (FZD), Frizzled-related protein family (SFRP), and intracellular signal transducers and modifiers (DVL1, AXIN). The genes encoding Casein kinase 1 alpha (CSNK1A), disheveled activator of morphogenesis 1 (DAAM1), which are agonists of the WNT pathway, are both overexpressed in the TE (Supplementary Table 8). These agonists were upregulated in differentiated ES cells [39]. In addition, we found that GSK-3B (glycogen synthase 3 kinase) expression is downregulated in the ICM, thus corroborating the reported inactivation of GSK-3B, which leads to the activation of the WNT pathway in maintaining the undifferentiated state of ES cells [45].
MAPK Signaling Pathway
The MAPK pathways regulate cell growth, differentiation, proliferation, and death. All members of the MAPK pathway have been shown to be expressed during mouse preimplantation development [46]. Differential expression of the various genes involved is shown in supplemental online Table 9.
Epigenetic Regulation of Lineage-Specific Gene Expression
De novo methylation of DNA by cytosine-5-methyltransferases is a well-characterized mechanism of epigenetic transcriptional control, and it has been shown that this mode of transcriptional control may contribute to the differentiation of the ICM and TE at the blastocyst stage [47]. Dnmt3b protein is specifically localized in the ICM of mouse blastocysts [48]. In addition, expression of DNMT3B and DNMT3A is enriched in undifferentiated human embryonic stem cells [810, 49] as well as in the ICM cell lineage (Fig. 1C
). This expression pattern suggests an important role in ICM-specific methylation in the blastocyst [47]. In contrast, transcripts of DNMT3L were expressed in both ICM and TE (Fig. 1C
). We also detected differential expression of several methyltransferases (supplemental online Table 1).
Other epigenetic regulators of X-inactivation, imprinting, maintenance of pluripotency, and the establishment of the TE lineages, including EZH2 (enhancer of zeste homologue 2), EED (embryonic ectoderm development), and CTCF (CCCTC-binding zinc finger protein), are expressed at high levels in the blastocyst and all ICM and TE samples [5052].
Imprinting
Several imprinted genes (H19, GRB10, SNURF, MEST, NAP1, UBE3A, DLX5, MAGEL2, OSBPL5/OBPH1, and ATP10A) were expressed at medium to high levels (BG-tag 70% to >90%) in the ICM, TE, and blastocyst. Our strict criteria for determining statistically significant differential expression between ICM and TE based on the microarray data may occasionally obscure more subtle differential expression patterns that are revealed when assayed by alternative methods. Real-time PCR identified a clear TE-biased (30-fold higher) expression for the imprinted gene IPL (imprinted in placenta and liver) (Fig. 3B
). Significantly, IPL/PHLDA2 is a marker of human cytotrophoblast and in the mouse Ipl restricts placental growth [53]. Thus, such imprinted genes that can act as regulators of nutrient supply at the feto-maternal interface may also influence growth and development of the early embryo.
| DISCUSSION |
|---|
|
|
|---|
The molecular profiles were consistent with expectations based on experimental mammalian embryology and were predictive of cellular physiology. For example, transcripts representing integrin- and cadherin-mediated cell adhesion, MAPK, and other gene products involved in the cell cycle and in apoptosis were identified. These components were anticipated in view of the rapid growth and dynamic morphogenesis of embryos at this transitional stage. Likewise, molecules representing pathways for glucose and sterol biosynthesis were also unsurprising. Absentees were presumably the less-abundant transcripts from these pathways.
The most interesting and potentially significant findings were among the representatives of key signaling pathways. Members of the WNT, TGF-ß/BMP, NOTCH, and phosphatidylinositol 3-kinase pathways were found. These genes are likely to be finely regulated to determine cell fate, adhesion, and migration because when expression goes awry, stem cells may undergo malignant transformation [57]. The expression of OCT4 and NANOG is already established as a marker of pluripotency, but there were other genes that were expressed in a similar temporal and spatial pattern. It will be important to investigate whether these genes are candidate markers of the ICM and their role in maintaining an undifferentiated state.
Perspectives
These data provide a reference for both normal development of the blastocyst in vivo and in vitro and for derivatives of the ICM and TE. Many human blastocysts generated in vitro fail to implant after transfer to a receptive uterus, and not all ICMs give rise to competent ES cell and trophoblastic stem cell lineages in culture. Markers of pluripotent stem cells are valuable tools for predicting cellular phenotype and are needed for advancing ES cell technology. Identification of other molecules that have roles in intracellular signaling is critical for the control of cellular differentiation and fate. This knowledge is critical for the key goal of directing the differentiation of stem cells or maintaining them as proliferating populations of pluripotent precursors, and preferably without feeder cell layers [45]. Some genes that are regulated epigenetically evidently have an unstable imprinting status, and research to safeguard the health of children born after assisted reproductive technologies must carefully assess the impact of culture conditions on the expression of these genes in preimplantation embryos [58]. Thus, there are many reasons why the molecular profiles of the ICM and TE cells in the human blastocyst are valuable, and, in the long term, they carry implications for fertility, reproductive health, and the prospects of stem cellbased technologies for treating human diseases.
| ACKNOWLEDGMENTS |
|---|
|
|
|---|
DISCLOSURES
The authors indicate no potential conflicts of interest.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
G. Cauffman, M. De Rycke, K. Sermon, I. Liebaers, and H. Van de Velde Markers that define stemness in ESC are unable to identify the totipotent cells in human preimplantation embryos Hum. Reprod., January 1, 2009; 24(1): 63 - 70. [Abstract] [Full Text] [PDF] |
||||
![]() |
I.-M. Shih and K.-T. Kuo Power of the Eternal Youth: Nanog Expression in the Gestational Choriocarcinoma Am. J. Pathol., October 1, 2008; 173(4): 911 - 914. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. M. Jones, D. S. Cram, B. Song, G. Kokkali, K. Pantos, and A. O. Trounson Novel strategy with potential to identify developmentally competent IVF blastocysts Hum. Reprod., August 1, 2008; 23(8): 1748 - 1759. [Abstract] [Full Text] [PDF] |
||||
![]() |
S J Kimber, S F Sneddon, D J Bloor, A M El-Bareg, J A Hawkhead, A D Metcalfe, F D Houghton, H J Leese, A Rutherford, B A Lieberman, et al. Expression of genes involved in early cell fate decisions in human embryos and their regulation by growth factors Reproduction, May 1, 2008; 135(5): 635 - 647. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Adjaye, R. Herwig, T. C. Brink, D. Herrmann, B. Greber, S. Sudheer, D. Groth, J. W. Carnwath, H. Lehrach, and H. Niemann Conserved molecular portraits of bovine and human blastocysts as a consequence of the transition from maternal to embryonic control of gene expression Physiol Genomics, October 19, 2007; 31(2): 315 - 327. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Sudheer and J. Adjaye Functional genomics of human pre-implantation development Brief Funct Genomic Proteomic, July 31, 2007; (2007) elm012v1. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Babaie, R. Herwig, B. Greber, T. C. Brink, W. Wruck, D. Groth, H. Lehrach, T. Burdon, and J. Adjaye Analysis of Oct4-Dependent Transcriptional Networks Regulating Self-Renewal and Pluripotency in Human Embryonic Stem Cells Stem Cells, February 1, 2007; 25(2): 500 - 510. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Cauffman, I. Liebaers, A. Van Steirteghem, and H. Van de Velde POU5F1 Isoforms Show Different Expression Patterns in Human Embryonic Stem Cells and Preimplantation Embryos Stem Cells, December 1, 2006; 24(12): 2685 - 2691. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Tabibzadeh and A. Hemmati-Brivanlou Lefty at the Crossroads of "Stemness" and Differentiative Events Stem Cells, September 1, 2006; 24(9): 1998 - 2006. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. K Riley and K. H Moley Glucose utilization and the PI3-K pathway: mechanisms for cell survival in preimplantation embryos. Reproduction, May 1, 2006; 131(5): 823 - 835. [Abstract] [Full Text] [PDF] |
||||
![]() |
|