Alteration of gene expression by alcohol exposure at early neurulation
© Zhou et al; licensee BioMed Central Ltd. 2011
Received: 6 September 2010
Accepted: 21 February 2011
Published: 21 February 2011
Skip to main content
© Zhou et al; licensee BioMed Central Ltd. 2011
Received: 6 September 2010
Accepted: 21 February 2011
Published: 21 February 2011
We have previously demonstrated that alcohol exposure at early neurulation induces growth retardation, neural tube abnormalities, and alteration of DNA methylation. To explore the global gene expression changes which may underline these developmental defects, microarray analyses were performed in a whole embryo mouse culture model that allows control over alcohol and embryonic variables.
Alcohol caused teratogenesis in brain, heart, forelimb, and optic vesicle; a subset of the embryos also showed cranial neural tube defects. In microarray analysis (accession number GSM9545), adopting hypothesis-driven Gene Set Enrichment Analysis (GSEA) informatics and intersection analysis of two independent experiments, we found that there was a collective reduction in expression of neural specification genes (neurogenin, Sox5, Bhlhe22), neural growth factor genes [Igf1, Efemp1, Klf10 (Tieg), and Edil3], and alteration of genes involved in cell growth, apoptosis, histone variants, eye and heart development. There was also a reduction of retinol binding protein 1 (Rbp1), and de novo expression of aldehyde dehydrogenase 1B1 (Aldh1B1). Remarkably, four key hematopoiesis genes (glycophorin A, adducin 2, beta-2 microglobulin, and ceruloplasmin) were absent after alcohol treatment, and histone variant genes were reduced. The down-regulation of the neurospecification and the neurotrophic genes were further confirmed by quantitative RT-PCR. Furthermore, the gene expression profile demonstrated distinct subgroups which corresponded with two distinct alcohol-related neural tube phenotypes: an open (ALC-NTO) and a closed neural tube (ALC-NTC). Further, the epidermal growth factor signaling pathway and histone variants were specifically altered in ALC-NTO, and a greater number of neurotrophic/growth factor genes were down-regulated in the ALC-NTO than in the ALC-NTC embryos.
This study revealed a set of genes vulnerable to alcohol exposure and genes that were associated with neural tube defects during early neurulation.
Children born to women who drink heavily during pregnancy are at risk for various developmental disorders, collectively called Fetal Alcohol Spectrum Disorder (FASD). Fetal Alcohol Syndrome (FAS) is a severe form of FASD in which the affected child is diagnosed with growth retardation, abnormal central nervous system development (typically including microencephaly), and a characteristic pattern of abnormal facial features [1–4]; organ dysmorphology, particularly of the eye and heart, may be evident in FAS cases as well [5, 6]. Disruption of complex molecular cascades that regulate embryonic morphogenesis likely are responsible for the teratogenic effects of alcohol. Potential mechanisms include metabolic stress, reduced signaling by transcription factors, retinoic acid or growth factors, disrupted cell-cell interactions, impaired cell proliferation, and apoptosis [7–16]. Several of these mechanisms may have direct roles in causing the cell death and growth retardation in multiple systems, including brain and head (for review see ).
Expression of a number of genes during development was reported to be affected by alcohol in different experimental paradigms, including homeobox genes such as Msx2  and sonic hedgehog [19, 20], neurotrophic molecules (e.g. ADNP gene ), fetal liver kinase 1 (Flk1) ), retinol-related genes (e.g. Crabp1 and Fabp4; ), nucleotide excision repair gene, (Ercc6l) , stress-related genes (e.g. heat shock protein 47 ), and differentiation and apoptosis genes such as Timp4, Bmp15, Rnf25, Akt1, Tulp4, Dexras1 . These altered genes suggest potential mechanisms for the abnormal development in FASD. However, the wide-ranging developmental abnormalities in FASD are likely a consequence of the interaction of multiple genes. Examination of global gene expression provides a holistic view of genes that potentially interact and collaboratively contribute to the abnormal development. Alcohol exposure induced changes in a group of cellular adhesion genes (e.g. L1cam and integrin) in neuroblastoma cells . A brief ethanol exposure (3 h) at gestation day 8 (E8) in mouse embryos altered expression of genes of metabolic, cell programming and cytoskeletal signaling pathways . An earlier alcohol exposure at E6-E8 also altered a set of genes related to PLUNC, neurofilament, and pale ear .
In animal models of prenatal alcohol exposure, sources of variability include the pattern, concentration, amount, and developmental stage of alcohol exposure, maternal stress, embryonic growth and maturation of embryos between litters and even within a given litter and within inbred strains of mice . Control of all these variables in rapidly developing embryos is virtually unattainable in vivo. To limit these variables, a whole embryonic culture [30, 31] was adopted, including stage alignment based on somite number, in which the pattern, amount and concentration of alcohol and embryonic staging were controlled. Inbred C57BL/6 mice, with known susceptibility to ethanol teratogenesis [32, 33], were used for this study.
Differences in the dose and timing of alcohol exposure are known contributors to variation in the phenotypic spectrum in FASD. Understanding the pattern of gene alterations that co-vary with different outcomes produced by different alcohol doses or developmental timing of exposure would provide valuable insights into mechanisms underlying this phenotypic variability. As development is highly dynamic throughout gestation, we asked how alcohol exposure might affect genome-wide gene expression at the critical stage of neurulation (E8-10), when the nervous system (and other major organs) are actively forming in mouse. We have shown that at this key stage, neural tube formation was highly sensitive to the alcohol insult . DNA methylation was altered, with the degree of change commensurate with severity of neural tube defect . In the current study, in an initial experiment, cluster analysis indicated distinct differences in gene expression not only between control- and alcohol-treated embryos, but also between two phenotypic subsets of alcohol-treated embryos discernable at the end of alcohol treatment, one group which had a closed neural tube (ALC-NTC) and the other group with an open neural tube (ALC-NTO). A second study with a larger set of arrays was then performed in which alcohol-treated embryos of both neural tube phenotypes were specifically compared. We report here the correlation of alcohol-induced embryonic growth retardation and neural tube abnormalities with changes in expression in networks of genes known to regulate embryonic growth, organ development, and neural specification processes.
Embryonic dysmorphology after alcohol exposure, scored according to Maele-Fabry et al,1992.
3 ± 0
2.80 ± 0.08
2.86 ± 0.10
2.70 ± 0.14
2.77 ± 0.09
2.15 ± 0.21
2.17 ± 0.26
2.11 ± 0.39
4.76 ± 0.10
3.81 ± 0.27*
4.57 ± 0.14
2.29 ± 0.29** ^^
4.52 ± 0.11
3.71 ± 0.27
4.50 ± 0.14
2.14 ± 0.14** ^^
4.71 ± 0.10
3.86 ± 0.24*
4.50 ± 0.14
2.57 ± 0.30** ^^
Caudal Neural Tube
4.76 ± 0.12
4.11 ± 0.19*
4.09 ± 0.26*
4.14 ± 0.26*
4.80 ± 0.09
4.33 ± 0.19
4.59 ± 0.19
3.81 ± 0.36*
4.80 ± 0.10
4.10 ± 0.16**
4.15 ± 0.19*
4.00 ± 0.31*
2.01 ± 0.06
1.51 ± 0.13**
1.48 ± 0.18*
1.57 ± 0.20
0.53 ± 0.10
0.20 ± 0.08*
0.21 ± 0.09
0.19 ± 0.14
2.08 ± 0.11
1.99 ± 0.09
2.12 ± 0.08
1.71 ± 0.18
2.41 ± 0.14
2.06 ± 0.16
2.21 ± 0.18
1.76 ± 0.30
0.47 ± 0.08
0.26 ± 0.08
0.29 ± 0.11
0.20 ± 0.13
3.59 ± 0.14
2.87 ± 0.14**
3.02 ± 0.17*
2.57 ± 0.20**
3.95 ± 0.12
3.68 ± 0.11
3.88 ± 0.10
3.29 ± 0.18* ^
4.81 ± 0.09
4.38 ± 0.16
4.50 ± 0.17
4.14 ± 0.34
53.97 ± 0.66
45.83 ± 1.54**
49.14 ± 1.54**
39.23 ± 1.72** ^^
Among 127 samples of alcohol-treated embryos, 34 (27%) had various degrees of incomplete neural tube closing (Figure 1); this compares to 3 (2%) out of the 139 controls. These openings in the neural tube mostly occurred in the head fold, although delayed or incomplete neural tube closure in midbrain and hindbrain was also seen. The abnormalities and developmental delays are clearly more severe in ALC-NTO than in ALC-NTC subgroups, particularly in development of the neural axis including hindbrain, midbrain, forebrain, otic vesicle.
Genes with changed expression†
Exp. 1 Fold Change
Exp. 2 Fold Change
expressed sequence AI415282
ATPase, H+ transporting, lysosomal protein 2
CDNA sequence BC008163
calcium/calmodulin-dependent serine protein kinase
CDC-like kinase 1
CDC like kinase 4
CREBBP/EP300 inhibitory protein 1
Cysteine rich protein 61
Dachshund 2 (Drosophila)
early B-cell factor 1
early B-cell factor 2
early B-cell factor 3
EGF-like repeats and discoidin I-like domains 3
EGF-containing fibulin-like extracellular matrix protein 1
forkhead box D1
histone 1, H3a
Histone 1, H4i
histone 3, H2a
insulin-like growth factor 1
Lectin, galactose binding, soluble 1
Melanoma antigen, family H, 1
myc target 1
N-ethylmaleimide sensitive fusion protein beta
N-myc downstream regulated gene 1
proviral integration site 1
Protein phosphatase 1, regulatory subunit 14A
pentaxin related gene
RAB11a, member RAS oncogene family
retinol binding protein 1, cellular
ribosomal protein L13a
ribosomal protein L17
Synapse associated protein 1
Tissue inhibitor of metalloproteinase 3
ubiquitin-conjugating enzyme E2B, RAD6 homology
Vascular cell adhesion molecule 1
RIKEN cDNA 1110008H02 gene
RIKEN cDNA 2010011I20 gene
RIKEN cDNA 2310034L04 gene
RIKEN cDNA 5033414D02 gene
RIKEN cDNA 5230400G24 gene
RIKEN cDNA 5730420B22 gene
RIKEN cDNA A630082K20 gene
acyl-CoA synthetase long-chain family member 6
ATPase, Na+/K+ transporting, alpha 1 polypeptide
expressed sequence AW547365
Expressed sequence C78212
Amino acid metabolism
carbamoyl-phosphate synthetase 2, aspartate transcarbamylase, and dihydroorotase
carnitine deficiency-associated gene expressed in ventricle 3
RIKEN cDNA E130306I01 gene
exosome component 2
high mobility group AT-hook 2
Hydroxysteroid 11-beta dehydrogenase 2
Insulin degrading enzyme
interferon alpha responsive gene
Nuclear Protein Transport
integrin alpha 6
Kruppel-like factor 16
NADH dehydrogenase (ubiquinone) Fe-S protein 1
PHD finger protein 13
Proteasome (prosome, macropain) 26 S subunit, non-ATPase, 3
pentatricopeptide repeat domain 1
ras homolog gene family, member U
RNA polymerase 1-4
serum amyloid A 2
solute carrier family 27 (fatty acid transporter), member 4
transformation related protein 53 binding protein 1
ubiquitin-conjugating enzyme E2, J1
RIKEN cDNA 0610007A15 gene
Adult male corpora quadrigemina cDNA, RIKEN full-length enriched library, clone:B230210C03 product:u
RIKEN cDNA 1300001I01 gene
RIKEN cDNA 1700017B05 gene
RIKEN cDNA 1700054N08 gene
RIKEN cDNA 4632417K18 gene
RIKEN cDNA 4930485D02 gene
RIKEN cDNA 5930416I19 gene
Genes in Experiment 2 that are turned on or off by alcohol treatment.
Change in Alcohol-treated
Adducin 2 (beta)
Complement component factor i
F-box only protein 2
GTP cyclohydrolase 1
Growth factor independent 1B
Natriuretic peptide precursor type B
Phosphatidylinositol membrane-associated 1
Protective protein for beta-Galactosidase
Tachykinin receptor 2
Acyl-Coenzyme A binding domain containing 5
F-box only protein 2
FERM domain containing 3
Programmed cell death 4
Aldehyde dehydrogenase 1 family, member B1
Four GSEA analyses were conducted within each experiment: control versus all alcohol-treated (ALC), control versus ALC-NTC, control versus ALC-NTO, and ALC-NTC versus ALC-NTO. As 415 GO gene sets and 191 stem cell related gene set were pre-selected, there were totally 4 × (415+191) = 2424 GSEA tests. We found 15 gene sets that were significant at 5% and shared the same enrichment direction in both experiments. By chance, one would expect only 2424 × (0.05 × 0.05 × 0.5) = 3; therefore, the FDR is 3/15 = 20%. The significant gene sets common to the two experiments are outlined below.
GSEA for Early Developmental Biology GO sets.
Gene Set Description
Control vs ALC-NTO/ALC-NTC (see legends)
(Ctgf, Igfbp2, Emp1, Osm, Cyr61, Gap43, Crim1, Tgfb3, Igfbp7, Nov, Emp3), Gpc3, Csf1, Socs2, Bmp6, Bmp4, Inhbb, Lepre1, Wrn, Wig1, Cish
Regulation of cell growth
Regulation of growth
( Ctgf, Igfbp2, Osm, Cyr61, Gap43, Crim1, Igfbp7, Nov), Gpc3, Csf1, Socs2 ,
Insulin-like growth factor binding
Ctgf, Igfbp2, Cyr61, Crim1, Igfbp7, Nov
Ctgf, Anxa2, Cyr61, Thbs1, Vegfa, Tie1, Elk3, Flt1, Crhr2, Vegfc, Kdr, Bmp4, Adra2b, Tnfrsf12a
Mab21l1, Neurod1, Neurod4, Ntrk2, Fkbp8, Bmpr1b, Crb1, Stat3, Tspan5, Pax6, Bmp4, Map3k1
Hist3h2b, a; Hist1h3f; Hist1h1c; Hist1h2b, c; Hist1h3a; H1f0; Smarca2; Nap1l3
ALC-NTO vs ALC-NTC
Growth, Growth retardation
Epidermal growth factor receptor (EGFR) signaling pathway
Pde6g, Egfr, Hbegf
(Enriched in ALC-NTC)
The growth-related genes represented the largest group of affected genes. There were 5 GO sets of growth-associated genes (Table 4). Many of these genes, identified by GSEA in both experiments, were also identified in Experiment 2 at the single gene level; e.g. the Growth gene set (GO:0040007): Ctgf, Igfbp2, Emp1, Osm, Cyr61, Gap43, Crim1, Tgfb3, Nov, Socs2, and Wrn were significantly reduced in Experiment 2, and Igfbp7, Emp3, Bmp4, Bmp6, Inhbb, Wig1, and Cish were reduced but did not reach the criteria for significance. The additional growth genes in Epidermal growth factor receptor (EGFR) signaling pathway GO group appear to be reduced to a greater extent in ALC-NTO than in ALC-NTC (Table 4).
GSEA for Stem Cell Related Gene Sets
Control vs ALC
Other ECM Molecules~^
Ctgf, Thbs2, Tgfbi, Ecm1
ECM Protease Inhibitors~,#
TGF-β Activin-responsiv~ e
Junb, Fos, Tgfbi, Pdgfb, Tgfb1i1, Igf1
Other Regulators of Cell Differentiation (Neural Specification) #
Elavl3, Neurod1, Neurod4, Nhlh1, Neurog1, Nhlh2, Neurog3, Spock2, Neurog2
Other Related growth Factor^
Other Related growth Factor#
Ctgf Hgf Igf1
RT-PCR confirmation of differences in gene expression: Neural specification genes from Experiment 1.
basic helix-loop-helix domain, class B5
SRY-box containing gene 5
myosin light chain
RT-PCR confirmation of differences in gene expression: Growth/neurotrophic factor genes from Experiment 2.
The abnormal embryonic development resulting from the alcohol treatment at this specific stage of development (Figure 1; Table 1) was consistent with our previous report  and those of others [41, 42]. Two different facets of abnormal development could be identified: growth delay and frank teratogenesis. Delays in growth were also evident by the significant reductions in the total RNA per embryo and in the delayed morphological staging (Table 1). The affected structures were derived from each of the three germ layers, i.e., neural tube and brain vesicles (ectoderm), somites and cardiovascular system (mesoderm/endoderm), and involved a wide range of tissues and organs (e.g., heart, head, limbs). Alterations in all of these have been observed in FAS cases. The teratogenic consequences were evident as dysmorphology of various organs (central nervous system, eye, and heart) that involved pathogenic effects beyond just the observed delay of the normal course of development. Examples include enlarged heart primordium and abnormally enlarged ventricular chambers, detached pericardial sac, small forebrain, flat telencephalic vesicle, failure in neural tube closure, and small and irregularly shaped eyes.
We observed in Experiment 1 that gene expression profiles from alcohol treatment of embryos in this controlled culture system yielded two distinguishable patterns; comparison to the morphological data revealed that these were correlated with two different phenotypes: open (ALC-NTO) and closed neural tubes (ALC-NTC). The phenotypes and correlated gene expression differences were reproduced in Experiment 2. The embryos with open neural tubes (ALC-NTO) had more severe delays in brain and otic development than those with closed neural tubes (ALC-NTC) (Table 1). These different phenotypes are consistent with our previous in vivo observation in a liquid diet model of prenatal alcohol exposure in C57BL/6 mice, which resulted in partial penetration of incomplete neural tube closure (as late as embryonic day 15) and a cascade of deficits in midline structural development . Finding this difference in development in experimentally controlled culture conditions indicates either a stochastic event or that an extremely sensitive gene-environment interaction is involved, e.g. different outcomes based on small differences in developmental stage at the time of exposure or small differences in tissue concentrations of alcohol across embryos. We have recently found greater DNA hypermethylation in ALC-NTO than in ALC-NTC embryos, particularly in genes on chromosomes 7, 10, and X. Remarkably, there was a >10 fold increase in the number of hypermethlyated genes on chromosomes 10 and X in ALC-NTO than ALC-NTC .
Both the ALC-NTC and the ALC-NTO embryos demonstrated lower expression of genes in sets related to cell growth, growth factors, heart (angiogenesis), and eye (in NTC vs. Control) (Table 4 ; Table 7). The ALC-NTC and ALC-NTO embryos also differed in other sets of functionally related genes. The histone gene set was selectively reduced in ALC-NTO compared to controls. The epidermal growth factor signaling pathway genes were lower in ALC-NTO than ALC-NTC (Table 4). At the single gene analysis level, Experiment 2 showed a greater number of neurotrophic/growth factor genes were down-regulated in ALC-NTO than in ALC-NTC groups, particularly in the TGFβ, NTF3, S100, and EGF families. These differences in gene expression between the ALC-NTO and ALC-NTC embryos appear to be correlated with the more severe teratogenic trajectory of the ALC-NTO group, but causal relationships have yet to be established.
The neural tube abnormality may either be a delay in neural tube closure or a neural tube defect. In either case, a delay in closing of the neural tube is associated with deficits in midline brain development due to disruption of the timing of critical events of early brain development. At more mature stages, such midline deficits include craniofacial abnormalities, corpus callosum, olfactory bulb, cerebellum, and raphe neuron formation [43–50].
Green and colleagues  reported that a 3 to 4 h binge-like alcohol exposure, with blood alcohol concentration 300 to 400 mg/dL at E8, produced a major abnormality in craniofacial and eye development in C57BL/6 mice at E15 or E17 (effects in the C57BL/6J substrain were greater than in the C57BL/6N substrain). Alterations of gene expression were reported to occur within hours of alcohol exposure at E8; these genes included metabolic and cellular gene, down-regulated ribosome and proteasome pathways; upregulated glycolysis and pentose phosphate, tight junction, and Wnt signaling pathways, as well as other cellular profile genes. In another study, a comparable high dose of alcohol exposure at an earlier stage, E6-E8, produced growth retardation, abnormal tail torsion, open neural tube, reduction of somite number, and other malformations . The altered gene expression at E10 included cytoskeletal (Neurofilament), signal transduction (Zinc finger protein, MAP kinase related, Transcription factor Nf2l2), and metabolic genes (lactate dehydrogenase, Aldolase 1). In the current study, a similar dose of alcohol exposure at the stage of neurulation (E8-10) produced a major neural and cardiovascular retardation and other organ system abnormalities. The trends of gene expression are consistent with the observed developmental delay and growth retardation in FASD. Among the genes with reduced expression in the alcohol-treated embryos were those involved in growth retardation, neural development, heart and hematopoiesis, and epigenetics. Among the identified functionally related gene sets, the most notable effect was the down regulation of growth-related genes, which represented the largest group of affected genes (Table 4). These genes provide plausible candidates for mechanistic links to the observed embryonic growth retardation.
Expression of neural specification genes (Table 5 and 7) and neurotrophic/growth factor genes (Table 4 and 7) was also reduced by the ethanol exposure. These participate in neuronal specification, neural stem cell differentiation, and neural fate determination [51–55]. Suppression of these genes predicts a downstream reduction in the early formation of neural cells. Null neurog 1 (Ngn1) or neurog 2 (Ngn2) leads to sensory abnormality [56, 57]). These differential expression of neuronal specification/patterning genes together with neurotrophic genes supports the dysmorphism and developmental delay of neural tube and fore-to mid-brain formation. The Igf1 and EGF genes were also identified by a microarray study with 3 h alcohol treatment  indicating they are altered early after ethanol exposure. The down-regulation of these neural specification and neural trophic/growth factor genes may play a major role in the neurodevelopmental deficit observed in the current study and featured in FASD.
Although heterogeneity of tissue arising from use of whole embryos might have masked some changes in specific tissues, two functional gene sets, optic vesicle and the heart (Table 4), were identified and specifically linked to our observed developmental delay and abnormalities. Also, the collective down-regulation of key hematopoiesis genes that were either absent (Table 3) or reduced (Table 2) is consistent with the reduced blood circulation observed in the embryos.
Many histone genes related to epigenetic regulation of transcription were affected by ethanol (Table 4). The reduction of many histone variants would alter chromatin organization, affecting transcription at a global level [58, 59]; this may be an important effect of the alcohol that leads to the reduction of total RNA and induced growth retardation. Modification of epigenetic processes is a potential mechanism by which alcohol may alter gene expression during development, and may be an important candidate mechanism for the pathophysiology of fetal alcohol syndrome.
Other genes that were present in the control group but absent in the alcohol-treated group (Table 3) likely reflect a delay in onset or a strong inhibition of normal expression at this stage of development. Among them, four hematopoiesis genes [glycophorin A (Gypa), adducin 2 (Add2), beta-2 microglobulin (B2m), and ceruloplasmin (Cp)] associated with blood cell formation were absent in the alcohol-treated groups; these genes are key components in the pathway of white and red blood cell formation [36, 38, 60–62]. The absence of these genes is in agreement with the low circulating blood cells seen in alcohol treated embryos (Figure 2). The expression of aldehyde dehydrogenase 1B1 (Aldh1b1) was induced in both of our experiments by alcohol treatment during this period of early neurulation (Table 2 last row). Because Aldh1b1 encodes an efficient enzyme for breakdown of acetaldehyde formed during metabolism of ethanol, this up-regulation is likely a detoxification response to the high level of ethanol in the environment. However, the metabolism of other substrates of this enzyme (e.g., retinoic acid, corticosteroids, biogenic amines, neurotransmitters, and lipids) that are required for normal development may be adversely affected by this increase in Aldh1b1 expression [63, 64].
In summary, alcohol exposure during the period of early neurulation at ~E8-E10, is predominantly inhibitory to gene expression, particularly the neural developmental genes. We found major reductions in gene sets involved in neurospecification, neural growth factors, cell growth and hematopoiesis. These effects on gene expression parallel the growth delay and developmental abnormalities including brain, neural tube, eye, heart, blood cells, and embryonic vascularization which are major targets in FASD. Our study, in conjunction with others that use different developmental periods of alcohol exposure, provides an important portfolio of alcohol-induced changes in gene expression associated with altered development. Together, these gene profiles should contribute to the generation of testable new hypotheses concerning the mechanistic path from gene expression changes to embryonic structural deficits, and for causal mechanisms of alcohol-induced teratogenesis (e.g., brain growth retardation, neural tube midline deficit, craniofacial dysmorphology) in fetal alcohol spectrum disorder. Two such hypotheses emerge from the current study. The first is that alcohol causes a delay in development of the nervous system by inhibiting specific sets of genes involved in neural development (Ngn, Nhlh, Sox, Igf, Ntf, and Egf). The second is that neural tube defects are mediated by the inhibition of genes in the epidermal growth factor signaling pathway and genes encoding histone variants.
All experimental procedures were approved by the Institutional Animal Care and Use Committee of the Indiana University School of Medicine (Indianapolis, IN) and are in accordance with the guidelines of the Institutional Animal Care and Use Committee of the National Institute on Drug Abuse, National Institutes of Health, and the Guide for the Care and Use of Laboratory Animals . Two-month-old C57BL/6 mice (~20 g) were purchased from Harlan, Inc. (Indianapolis, IN). Upon arrival, breeder mice were individually housed and acclimated for at least one week before mating began. The mice were maintained on a reverse 12 h light-dark cycle (lights on: 19:00 - 07:00) and provided with laboratory chow and water ad libitum. Two females were placed with one male for two hours between 08:00 and 10:00. When a vaginal plug was detected after the mating period, it was designated as embryonic day 0 (E0). On E8.25 at 15:00, dams were sacrificed using CO2 gas. The embryos were treated at this stage, which is the beginning of neurulation. The window of 46 hrs treatment covered the stages of the formation of the major organs, neural specification and patterning. These stages are known to be vulnerable to alcohol .
The technique for whole embryo culture was based on the methods described by New . The gravid uterus was removed and placed in sterile PBS (0.1 M phosphate buffer containing saline) at 37°C. The embryo in the visceral yolk sac along with a small piece of the ectoplacental cone (hereafter called embryo, unless otherwise stated) was carefully removed from the deciduas tissues and the Reichert's membrane in PBS containing 4% fetal bovine serum (Sigma, St Louise. MO). After removal, three embryos bearing 3-5 somites (E8.25) were incubated in a culture bottle in 20 mL of medium which consisted of 70% immediately centrifuged heat-inactivated rat serum (Harlan Sprague-Dawley, Inc, Indianapolis, IN) and 30% phosphate buffered saline (137 mM NaCl, 2.7 mM KCl, 0.5 mM MgCl2, 8 mM Na2HPO4, 1.47 mM KH2PO4, 0.9 mM CaCl2, 5.6 mM glucose, 0.33 mM sodium pyruvate, pH7.4) supplemented with 20 units/ml penicillin and 20 units/ml streptomycin (Sigma, St. Louis, MO), and gassed with 5% O2, 5% CO2, and 90% N2 in a rotating culture system (B.T.C. Precision Incubator Unit, B.T.C. engineering, Cambridge, England, 36 rpm) for 2 h. After 2 h, treatment was initiated by transferring embryos into the same medium with or without 88 mM ethanol in isotonic buffer. The bottles were gassed for an additional 20 h with 5% O2, 5% CO2, and 90% N2, and then between 22 h and 46 h with 20% O2, 5% CO2, and 75% N2. The culture medium in alcohol and control cultures was replaced with fresh medium (with or without ethanol, respectively) 22 h after the start of the treatment. In this culture system, it was previously determined that the media alcohol concentration declined from 88 mM to 44 mM over the course of the experiment. Alcohol concentrations in this range (44-88 mM) have been commonly used in whole embryo cultures to generate FAS-related structural malformations [41, 42, 67] in multiple strains of mice , and are comparable to blood alcohol concentrations produced by in vivo doses of acute ethanol injections that produce teratogenic effects in mice during this embryonic period . This level, though high, is within the range attained by human alcoholics [69, 70].
All cultures were terminated 46 hrs from the beginning of treatment. The concentration of ethanol in the medium was assayed at three time points on each day (0 [initial], 12, and 22 hours on the first day; at 0 [after media change], 12, and 24 hours on the second day) in a separate group of embryos not used for the analyses, to avoid the potential confounding effects of drawing samples from the cultures. Media samples from alcohol- or vehicle-treated cultures were assayed in duplicate for alcohol concentrations using an Analox alcohol analyzer (Analox Instruments USA, Lunenburg, MA).
At the end of culture, viability was confirmed by observing the blood circulation of the yolk sac and the beating heart. Cultured embryos were quickly immersed in 0.7 ml TRIzol (Invitrogen, Carlsbad, CA) and homogenized for extracting total RNA for the RT-PCR and microarray processes (see microarray section, below), or fixed in 4% paraformaldehyde in PBS for the evaluation of the developmental status.
Whole embryos were used because the dysmorphology is observed throughout tissue derived from the three germ layers and in various developing organs (e.g., head fold, caudal neural tube, heart, lung bud, somites, and limbs). Also, dissection of the millimeter size embryos would unavoidably introduce another source of variability: whole embryos yield sufficient total RNA for single embryo analysis, whereas dissected tissues yield too little RNA and would require pooling or amplification for microarray analysis. Although this limits the resolution of genes contributing to different topographic changes, we thought that obtaining a complete gene expression profile in parallel with this widespread alcohol-induced teratogenesis in the embryo would be informative.
The analysis of embryo dysmorphology was performed as described by van Maele-Fabry et al.  and in our previous report . The morphological features of the developing embryo, including the allantois, flexion, heart, caudal neural tube, hind-brain, midbrain, forebrain, otic system, optic system, branchial bars, maxillary process, mandibular process, forelimb, hindlimb, and somites, were examined and scored for any malformations using the ordinal scales of our previous report . Scores for each of the above features were typically not normally distributed, so they were analyzed statistically by the non-parametric Mann-Whitney U test. The number of somites was normally distributed, so those data were analyzed by Student's t-test, using StatView software (SAS Institute, Inc. Cary, NC).
Two microarray experiments were performed. In Experiment 1, total RNA was isolated from individual whole embryos (4 vehicle control, 4 alcohol treated). Each embryo was immediately immersed in 700 ml TRIzol (Invitrogen) and homogenized using a Polytron. Extraction followed the TRIzol protocol. Ethanol precipitated RNA was resuspended in DEPC water. RNA was cleaned up using RNeasy mini-kit (Qiagen, Valencia, CA) The quality of RNA was assessed by the Agilent Bioanalyzer (Agilent Technologies, Waldbronn, Germany)and by spectrophotometry from 220 nm to 350 nm; concentration was determined from A260. Typical total RNA yields were 5-10 μg/embryo. Microarray analysis was performed at the Center for Medical Genomics at the Indiana University School of Medicine. Labeling and hybridization to Affymetrix Mouse Genome 430A GeneChips® (Affymetrix, Santa Clara, CA) were carried out following the manufacturer's suggested procedure. Fragmented biotinylated RNA from each embryo was separately hybridized to its own GeneChip for 17 hours at 42°C. The microarray analysis revealed striking differences among the 4 alcohol treated samples, which segregated as two separate pairs rather than one set of four; subsequently, it was noted that one pair of embryos had an open neural tube (ALC-NTO) and the other pair had the neural tube closed (ALC-NTC). All 4 control embryos had closed neural tubes.
Experiment 2 was designed to follow-up these initial results and provide an independent test of the gene expression correlations with the two neural tube phenotypes. Total RNA was isolated from individual embryos (4 vehicle control, 7 alcohol treated: 4 ALC-NTO, 3 ALC-NTC). RNA extraction and microarray analysis was as described above, except that Affymetrix Mouse Genome 430 2.0 GeneChips® (Affymetrix, Santa Clara, CA) were used.
The Mouse Genome 430A chip contains over 22,600 probe sets representing transcripts and variants from over 14,000 well-characterized mouse genes. The newer Mouse Genome 430 2.0 Array contains all of the probe sets present on the earlier 430A chip plus additional probe sets for a total of approximately 45,000 probe sets that analyze the expression of over 39,000 transcripts and variants from over 34,000 well characterized mouse genes. The differences in feature size and probe set content make direct comparisons inappropriate, due to scanning and scaling issues, but because the probe sets on the 430A are present on the 430 2.0 array, those can be compared at the level of gene lists.
The data from independent arrays (each with RNA from a single embryo) for each of the treatments were extracted using the Affymetrix Microarray Suite 5.0 (MAS5) algorithm. Data for both experiments have been deposited in GEO/NCBI and have been assigned series accession number GSE9545 and sample numbers GSM241642 through GSM241660.
To minimize false positive results, only genes detected ("present" by the MAS5 algorithm) on at least half of all individual arrays in at least one experimental condition were retained for further analysis. This avoids data that primarily represent "noise" [72, 73].
To detect differentially expressed genes, control samples were compared to ALC-NTC samples, or ALC-NTO samples, or their combination, using a Welch's t-test on the log-transformed signals. To see genes that were similarly affected in both experiments, we intersected the gene lists. To avoid missing genes that met a stringent significance threshold in one experiment but were just beyond that threshold in the second, we chose p ≤ 0.05 as the threshold for each experiment. Given that the two experiments were independent, the probability that a gene overlaps by chance and differs in the same up/down direction in both experiments is (0.05)*(0.05)/2 = 0.00125. False discovery rate (FDR) was calculated based on the number of genes expected to be significant and in the same direction in both experiments under the null hypothesis/the number of such genes actually found.
Hierarchical clustering with average linkage function was used to construct a dendrogram based upon all genes that were present on at least half of the arrays in an experimental group.
Gene Set Enrichment Analysis (GSEA) [74, 75] was carried out to identify groups of related genes that were differentially expressed. GSEA analyses were conducted for 4 different comparisons: control vs. ALC, control vs. ALC/NTC, control vs. ALC/NTO, and ALC/NTC vs. ALC/NTO. The top ranked genes in a significant gene set, in the region up to the maximum score, were considered significant. To reduce multiple testing issues, the GSEA in this study was conducted using two gene set databases designed to test the hypotheses that groups of genes related to Early Development or Stem Cells were differentially affected by alcohol.
(b) Stem Cell Related Gene Sets: 191 GO categories related to stem cells, neurogenesis, osteogenesis, extracellular matrix, developmental signal transduction pathway, cell cycle, growth factor, TGFβ/BMP signaling, Wnt signaling, and notch signaling were developed by Superarray Bioscience http://www.superarray.com. The gene set information is listed in Additional file 3 (shown with consent of Superarray Bioscience, Frederick, MD).
Primers for qRT-PCR
After Experiment 2, we decided to test the three groups (control, ALC/NTO, ALC/NTC) as pools, and chose growth/neurotrophic genes. A separate experiment was carried out with embryonic treatments identical to those used in Experiment 1. Whole embryos were homogenized in TRIzol (Invitrogen) using a Mini-Bead-Beater-8 (Bipspec products, INC, Bartlesville, OK), and total RNA isolation was as described above. Two different pools were created for each condition: Control1 (n = 12), ALC/NTC1 (n = 16), ALC/NTO1 (n = 5), Control2 (n = 5), ALC/NTC2 (n = 9), ALC/NTO2 (n = 6). The relative quantification of expression of each RNA pool was performed using the ABI Prism 7700 Sequence Detection System and calculated using the standard curve method (Applied Biosystems, User Bulletin #2; http:////www.appliedbiosystems.com). In each experiment, a relative expression level was determined for the two pools from each group in triplicate; 3-4 repeat experiments were performed, resulting in 18-24 values from each group. The treatment groups were compared with one way ANOVA followed by Student's t test.
This study was supported by NIAAA P50 AA07611 to the Alcohol Research Center at Indiana University School of Medicine (Center PI, David W. Crabb; FAS project, Feng C. Zhou), and in part by AA016698 to FCZ. Microarray studies were carried out in the Center for Medical Genomics at Indiana University School of Medicine, which is partially supported by the Indiana Genomic Initiative at Indiana University (INGEN®); INGEN® is supported in part by the Lilly Endowment, Inc. We thank Ms. Li-jun Ni for her technical assistance in embryonic culture, and Ronald Jerome and Chunxiao Zhu for technical assistance in preparation of microarray.
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.