Skip to main content

Maternal diabetes alters transcriptional programs in the developing embryo

Abstract

Background

Maternal diabetes is a well-known risk factor for birth defects, such as heart defects and neural tube defects. The causative molecular mechanisms in the developing embryo are currently unknown, and the pathogenesis of developmental abnormalities during diabetic pregnancy is not well understood. We hypothesized that the developmental defects are due to alterations in critical developmental pathways, possibly as a result of altered gene expression. We here report results from gene expression profiling of exposed embryos from a mouse diabetes model.

Results

In comparison to normal embryos at mid-gestation, we find significantly altered gene expression levels in diabetes-exposed embryos. Independent validation of altered expression was obtained by quantitative Real Time Polymerase Chain Reaction. Sequence motifs in the promoters of diabetes-affected genes suggest potential binding of transcription factors that are involved in responses to oxidative stress and/or to hypoxia, two conditions known to be associated with diabetic pregnancies. Functional annotation shows that a sixth of the de-regulated genes have known developmental phenotypes in mouse mutants. Over 30% of the genes we have identified encode transcription factors and chromatin modifying proteins or components of signaling pathways that impinge on transcription.

Conclusion

Exposure to maternal diabetes during pregnancy alters transcriptional profiles in the developing embryo. The enrichment, within the set of de-regulated genes, of those encoding transcriptional regulatory molecules provides support for the hypothesis that maternal diabetes affects specific developmental programs.

Background

Maternal diabetes disturbs embryonic development and can cause diabetic embryopathy, with cardiovascular malformations, neural tube defects and caudal dysgenesis as the most characteristic congenital malformations [1, 2]. Diabetes-induced dysmorphologies have been ascribed to increased apoptosis [3, 4], perturbation of prostaglandin synthesis and metabolism [57], deficiency in membrane lipids [811], and generally altered metabolism in the embryo [12]. Several studies have associated oxidative stress with the maternal diabetic condition, and the administration of anti-oxidants reduced the incidence of developmental defects in experimental models of intrauterine exposure to diabetes [5, 1319]. However, it is unclear how systemic metabolic disease results in particular developmental defects that are restricted to specific tissues in diabetic embryopathy [20].

Growing evidence suggests that maternal diabetes alters expression of developmental genes in the embryo, resulting in abnormal morphogenesis. Decreased expression of Pax3, a gene involved in neural tube defects [21, 22], was found in embryos from diabetic mouse dams at gestation day 8.5, with neural tube defects evident by day 10.5 [4]. Pax3 deregulation, presumably through oxidative stress [23], is also associated with heart defects that involve neural crest cell derivatives [24]. We recently showed that Wnt signaling is affected in diabetes-exposed mouse embryos [25]. These findings support the notion that diabetic pregnancy leads to altered expression of molecules that play key roles in patterning and development of embryonic tissues, implicating altered transcriptional regulation as a potential pathogenic mechanism in diabetic embryopathy. In order to identify genes and pathways affected by maternal diabetes, we performed gene expression profiling of diabetes-exposed mouse embryos using oligonucleotide microarrays.

Results and discussion

Animal model of diabetic embryopathy

Mouse embryos were isolated from diabetic or control dams at embryonic day 10.5 (E10.5) because at this stage neural tube defects are easily detectable. The frequency of NTDs in diabetes-exposed embryos was approximately 17% (16/96 diabetes-exposed embryos) compared to 0% in normal pregnancies (0/220). Except where noted, no malformed embryos were used for gene expression studies. We found no significant differences in litter size or number of resorbed embryos between diabetic (n = 11) and non-diabetic pregnancies (n = 10; P = 0.07). We also did not detect any developmental delay in apparently unaffected diabetes-exposed embryos; their morphological appearance, i.e. features of brain development, limb development, and somite numbers were commensurate with developmental age (Kruger et al., manuscript in preparation). All experiments used the FVB inbred strain. We here report the results from two independent expression profiling experiments (Figure 1), using individual embryo samples in Experiment I and and a pooling strategy in Experiment II.

Figure 1
figure 1

Experimental approaches to determine gene expression profiles of normal and diabetes-exposed embryos. Panels A and B depict the two independent microarray experiments. Panels C and D depict expression profiles where each colored vertical line represents the expression signal for one gene and row represent individual embryos (Exp. I) or samples (Exp. II). Red represents increased expression, blue reflects decreased expression, and intermediate colors represent minor changes (the color range was chosen along an arbitrary scale). Using a hierarchical clustering algorithm (with euclidean distance metric and centroid linkage rule, as implemented in GeneSpring), these graphic representations shows that expression profiles for embryos exposed to maternal diabetes differ considerably from control unexposed embryo profiles.

Gene expression profiling in embryos exposed to diabetes

In Experiment I, we surveyed the expression profiles of 2 control and 5 diabetes-exposed embryos. In an initial comparison of expression levels between control and diabetes-exposed embryos, 302 probe sets passed the "fold-change>2" criterion, and their expression profiles were visualized using hierarchical clustering (Figure 1, Panel C). Control embryos exhibited profiles similar to each other, and differences to the expression profiles of individual diabetes-exposed embryos are visually obvious. These results support the hypothesis that maternal diabetes affects gene expression in the exposed embryos.

This survey covered about 14000 genes, and after application of the analysis criteria (see Methods for details), we identified 126 genes (~1% of the total) with expression levels that were changed in diabetes-exposed embryos by more than 2-fold relative to controls (Table 1; for full gene names and accession numbers, see Additional File 1). An additional 378 genes displayed expression differences between 1.5- and 2-fold (data not shown). The majority (83%) of the 126 genes we identified were expressed at lower levels in diabetes-exposed compared to control embryos, and this was reflected in the larger dataset as well (72.8% of the additional 378 genes were decreased in expression). This decreased gene expression was not the result of developmental delay, since the morphology of diabetes-exposed embryos was stage-appropriate.

Table 1 Genes affected by maternal diabetes classified by cellular function. see Additional file 1

Validation of microarray by quantitative real-time PCR

Changes in gene expression detected by microarray were validated by quantitative real time PCR (Q-RT-PCR) for selected genes with potential relevance to diabetes or embryonic development (Table 2). Embryo samples were from different pregnancies than those employed for the microarray studies, and only embryos were used that appeared morphologically normal. Of all genes assayed, 16 exhibited differential expression in the Q-RT-PCR assay, confirming the microarray results in independent embryo samples. Three genes exhibited no differences, and seven genes were differentially expressed in both assays (P < 0.05); however, the change occurred in opposite directions. The discrepancies were traced back to (i) annotation problems: Hmga1, Lin28, Phip, (ii) different length of 3'UTR sequences where location of the microarray probe would not query all transcripts arising from the respective genes: Sema3a, Rod1, Slc2a1, or (iii) alternative splicing: Ogt, Tfrc [26], Creb1 [27]. The Q-RT-PCR assays were designed to specifically amplify a region of the transcript different from that covered by the microarray probe in order to obtain an independent measurement. For the majority of genes we tested, the independent assay confirms the initial finding that expression levels are altered in embryos exposed to maternal diabetes.

Table 2 Validation of microarray results by quantitative RT-PCR.

We cannot formally exclude altered mRNA stability as a factor causing the observed changes in mRNA levels, but it would be difficult to explain how the stability of relatively few transcripts could be altered in a selective fashion. Rather, we find that many transcription factor genes are down-regulated in their expression in diabetes-exposed embryos, and this trend is also reflected in the group of genes with 1.5- to 2-fold differences in expression; most likely therefore, the observed lower levels of gene expression are due to diminished or deregulated transcription.

Molecular classification of genes altered in response to maternal diabetes

The biological roles of many products encoded by the 126 diabetes-affected genes are known (Table 1 and Figure 2, Panel A). Most intriguingly, the largest functional category is comprised by transcription factors (15/126) and DNA-binding molecules known to affect transcriptional regulation (8/126), with both categories together comprising 18% of the identified genes. Thus, relative to the 7% of genes in the mouse genome that encode transcription factors [28], we found transcription regulatory genes highly overrepresented among our diabetes-affected genes. This is also reflected in enrichment of this category in DAVID annotation http://david.abcc.ncifcrf.gov/. A combined 22% of the genes encode cell surface receptors (10%) and signal transduction molecules (12%) that ultimately converge on transcription. The transcription factor Pax3, which was previously identified as affected by maternal diabetes [4] was not represented on the microarray, but was changed as expected when assayed by Q-RT-PCR [25]. These results indicate that transcription factors and signaling molecules are prominent targets for perturbation by maternal diabetes, and that altered transcriptional regulation plays a major role in the response of embryos to intrauterine exposure to diabetic conditions.

Figure 2
figure 2

Expression of Hif1α during mouse embryogenesis. Quantitative RT-PCR for Hif1α at various stages of development normalized to expression levels of Pole4. n = number of individual embryos tested.

Transcription regulation

If maternal diabetes deregulates cohorts of genes in the developing embryo through shared pathways of transcriptional regulation, one would expect (i) occurrence of common transcription factor binding sites (TFBS) in regulatory regions associated with multiple genes, and (ii) overrepresentation of such TFBS relative to other genes in the genome. We therefore analyzed the promoter sequences (5kb upstream of the transcription start) of our diabetes-affected genes for the presence of TFBS. As expected for genes with mostly broad expression patterns [29], we found a diverse set of conserved motifs in the upstream regions of these genes (see Additional file 2). Intriguingly, there was prominent over-representation of binding sites for the transcription factors FOXO1 and FOXO4 (-Log(p) = 12.416 and -Log(p) = 10.544, respectively), which are known to be involved in the response to oxidative stress, and for HIF1 (-Log(p) = 10.219), which is involved in the response to hypoxia. Out of 109 genes for which results were returned, FOXO1 and FOXO4 sites were enriched in 68% of the genes. NRF2 motifs were found in 76 promoters, further supporting the notion that oxidative stress may be involved in the response of diabetes-exposed embryos. HIF1 motifs were enriched in promoter regions of 22 genes (20.2%). These TFBS occur in combinations with sites for other transcription factors with a known role in responses to hypoxia, such as ATF4, E2F1 and E2F4, EGR1, ETS1, IRF1, NfkappaB, SOX9, SP3, and XBP1 (see Additional file 2 for references). Given the proposed role of oxidative stress and hypoxia in the pathogenesis of diabetic embryopathy [5, 6, 14, 17, 18, 30], it is striking that 97% of the genes affected by maternal diabetes carry in their upstream regions potential binding sites for transcription factors that are involved in responses to oxidative stress and hypoxia.

Both conditions have been reported to be associated with diabetic pregnancy [14, 31], and would predict activation of hypoxia-regulated pathways in the embryonic response to diabetic conditions. Paradoxically, the expression of Hif1α, a key regulator of embryonic responses to hypoxia [32], was reduced in our diabetes-exposed embryos at E10.5. In a post-hoc analysis of our microarray data with specific focus on HIF1 target genes [33], 22 HIF1 targets showed altered expression (fold-change >1.5) and passed one of the t-tests; 9 genes passed both statistical filters (details, see Additional file 3). Twenty HIF1-regulated genes exhibited increased expression in diabetes-exposed embryos, possibly reflecting an embryonic response to increased hypoxia. Further support for this idea comes from Hif1α message levels that are increased in diabetes-exposed embryos at E8.5 and E9.5 (Figure 3), and this increased Hif1α expression could be responsible for the increased expression of HIF target genes at E10.5. Together, these results implicate oxidative stress and hypoxia pathways in deregulated gene expression in diabetes-exposed embryos and identify the molecular targets of these pathways.

Figure 3
figure 3

Classification of diabetes-affected genes by molecular function and function in vivo. Panel A depicts the representation of molecular classes of encoded products for the diabetes-affected genes identified in Experiment I. Panel B depicts the representation of molecular classes for the diabetes-affected genes identified in Experiment II. Genes encoding products with unknown function were omitted. Panel C depicts in vivo phenotypes based upon MGI annotation for genes identified in Experiment I; Panel D depicts in vivo phenotypes for genes in Experiment II. Only knockout phenotypes (null and conditional) were included.

Functional roles of genes deregulated by maternal diabetes

Our identification of genes whose expression is affected by exposure to maternal diabetes suggests that those genes could be involved in the developmental defects in diabetic pregnancies. This notion is supported by qualitative expression information (MGI) for 69 of these genes, with 62 detected in the embryonic CNS, and 35 in the embryonic cardiovascular system. With both CNS and heart frequently being affected in diabetic embryopathy, genes with abnormal expression in these tissues might contribute to pathogenesis of birth defects.

Functional data also support this hypothesis: most remarkably, about one fourth of the diabetes-affected genes are known genes for which a functional role in embryonic development was established experimentally: in knockout mutant mice, 35 of the genes we identified have been shown, by genetics, to be required for mouse embryonic development (Table 3 and Figure 3, Panel C; references are given in Supplementary Table 4). This implies that these genes, under conditions of maternal diabetes but in the absence of genetic alterations, are subject to gene-environment interactions and respond to the intra-uterine environment of a diabetic pregnancy. Deficiencies in 15 of our genes (Agtr2, Cxadr, Dysf, Hif1a, Il6st, Itgav, Pdgfra, Pxn, Sema3a, Tgfbr1, Vcl, Adam10, Epha3, Efnb2, and Sfrs2) have been shown to cause cardiovascular defects or disease in mouse models. In humans, the risk for congenital heart disease is 2.8 times higher in infants born to diabetic mothers compared to the offspring of non-diabetic mothers [2] with higher odds ratios for specific malformations [34, 35]. In light of this, it is possible that down-regulation by maternal diabetes of one or more of the genes we discovered could contribute to abnormalities of the heart [36].

Table 3 In vivo function of genes affected in diabetes exposed embryos.
Table 4 Wnt-pathway genes affected by maternal diabetes.

Central nervous system malformations occur in about 5% of children born to diabetic mothers [35], which represents an up to 15-fold higher risk of over pregnancies unaffected by diabetes. Intriguingly, we found seven genes affected by maternal diabetes that previously have been associated with neural tube defects (NTD): Hif1a, Pdgfra, Twsg1, Adam10, Tgfbr1, Tfrc, and Vcl (for references, see Additional file 4). Thus, de-regulated expression of these genes in diabetes-exposed embryos might predispose embryos to neural tube defects. We also analyzed our dataset from Experiment I for differences between the 5 embryos that were exposed to maternal diabetes, of which two exhibited defective closure of the neural tube. In this comparison, we identified only two genes (Etnk1, Gmfb) that passed the criteria filter of >1.5-fold change and both statistical significance tests. Both genes were detected in the initial diabetes-exposed versus normal comparison; we did not discover any genes that were uniquely altered in NTD embryos, implying that differences between NTD-affected and -unaffected individuals with regard to gene expression are mostly quantitative. Our results are consistent with the hypothesis that diabetes of the mother alters expression of specific known heart-defect and neural-tube-defect genes, and that these genes may be responsible for the birth defects in diabetic pregnancies.

Confirmation of major findings by a separate profiling experiment

Our intial survey employed embryo samples that came from two pregnant females: one STZ-treated diabetic and a control untreated dam. This presents the theoretical possibility that any differences between progeny of the two dams reflect differences between pregnancies in addition to diabetic state. Also, we used individual embryo samples, and this approach is likely to incur substantial variability in the data and thus understimation of molecular changes. To address both concerns, we conducted Experiment II, in which equal amounts of RNA was pooled from three embryos of same gestational age into one sample, with each embryo derived from a different dam; for the diabetic as well as the control condition, we prepared three such pools, respectively (Figure 1, Panel B). Taking advantage of technical advances, these samples were processed and hybridized to the Affymetrix Mouse 430 2.0 chip, which surveys 39000 transcripts. Using the same analysis criteria as before, we identified 2231 transcripts of which 276 (12.37%) showed increased levels of expression in the diabetes-exposed samples, and 1955 (87.63%) exhibited decreased expression compared to the controls (Figure 1, Panel D). Thus, we confirm the general trend in the results from Experiment I. Of the differentially expressed transcripts, 179 lacked identifying features, such as a name, RefSeq or ENSEMBLE IDs, or Unigene Mm. cluster number; they also lacked any annotation information. This left us with 2052 annotatable genes. Classification by molecular functions revealed a distribution of molecular properties (Figure 3, Panel B) highly similar to that in Experiment I (Figure 3, Panel A). Again, genes encoding transcription factor and DNA-binding regulatory molecules were significantly enriched, accounting for 16.3% of the deregulated genes; strong enrichment was confirmed by DAVID annotation. Annotation for function in vivo identified 1836 gene entries in MGI; for 1095 of those, phenotype information was not available. However, 747 genes were associated with documented phenotypes in mouse mutants, of which 388 are developmental phenotypes by virtue of embryonic, neonatal, or perinatal death of homozygous mutant offspring. Again, the distribution of particular phenotypes in Experiment II (Figure 3, Panel D) was very similar to that of Experiment I (Figure 3, Panel C). Metabolic abnormalities were reported for mutants of 46 genes, and evidence for abnormal growth (pre- and post-natal) was obtained for 279 genes. Most notably, 161 genes are known to be associated with heart defects when mutated, and 112 genes are known to play causal roles in neural tube defects. This is only a fraction (35%) of the more than 300 NTD genes contained in the MGI database (as of October 1, 2008). Similarly, from the published collection of 170 mouse mutants with neural tube defects [37] for which the underlying molecular defect is known, 55 genes (32% of 170) were identified in Experiment II. Taken together, these results indicate that maternal diabetes affects specific pathogenic pathways leading to NTDs. Except for two genes, all NTD genes exhibited decreased expression on the arrays. In summary, the main findings of the initial microarray experiment were confirmed.

Indeed, of the 126 genes whose expression was altered by more than 2-fold in Experiment I, 67 were also recovered above the 2-fold change cut-off in Experiment II. Of the 378 probe sets with expression level altered between 1.5-fold and 2-fold, 187 were shared. Thus, of the 504 probe sets with altered expression in Experiment I, more than half (254) were confirmed in Experiment II, providing independent validation for the major results of the first experiment. This 50% confirmation rate for independent microarray experiments in the same biological paradigm agrees well with similar findings for independent yeast microarray results [38]. Thus, employing individual embryo samples as well as a pooling strategy, we have identified molecular targets in the embryo that respond to maternal diabetes. Also noteworthy is that Experiment II confirmed our earlier candidate gene studies that showed components of the Wnt pathway altered in diabetes-exposed embryos [25]. In fact, 43 genes with roles in Wnt signaling are affected by maternal diabetes (Table 4); with exception of Cank2d, Rbx1 and Tcf7l2 (which are upregulated), the expression levels of all of these genes are decreased in diabetes-exposed embryos compared to controls. This finding provides further support for our hypothesis that maternal diabetes affects specific developmental programs.

Using mRNA from whole individual embryos allowed us to survey a broad range of embryonic tissues that are potentially affected by maternal diabetes. This approach might have "missed" effects on genes that are expressed only in small cell populations of the embryo. However, in Experiment II, we identify 19 of the 47 published genes that were found altered more than 1.5-fold in microarray analysis of cranial neural tube tissue from diabetes-exposed embryos with neural tube defects at E 11.5 [39], and four of those genes are shared with Experiment I. Concordance was found for increased expression of Bnip3, and for decreased levels of En2, Hes6, Ina, Map3k7, Med1, Msx1, Mtap1B, Ngn2, Notch1, TgfβII, Doublecortin, Protocadherin18, Tgfβρεχεπτορ1, TopoIIβ, with the latter four genes confirmed also in Experiment I. Notch3, Nr2f2, Shh, and Tial1, were increased in dissected neural tube [39] but decreased in whole embryos, indicating that they may be deregulated in multiple tissues. Nonetheless, the overlap between results from different laboratories, despite differences in experimental design, provides additional validation to our findings.

We cannot currently distinguish which of the changes in gene expression are in direct response to the diabetic milieu, and which are indirect changes downstream of altered transcription factor expression, potentially increased hypoxia [14] or alterations in yolk sac [40] or placenta [Salbaum, Kruger, Pavlinkova, Zhang, and Kappen, manuscript submitted]. In this regard, it is interesting to note that we find no congruence to genes reported as affected by maternal diabetes in yolk sac of E12 rat embryos [41]. This indicates not only that both yolk sac and embryo gene expression are affected by maternal diabetes, but that extra-embryonic tissues respond differently than the embryo proper. It is noteworthy that among the deregulated genes with known phenotypes in mouse mutants, over 100 have been reported to be associated with placental alterations. Even though we have only surveyed the embryo proper, this is suggestive evidence that placental gene expression may also be altered in diabetic pregnancies. Our findings are consistent with the idea that altered gene expression in the embryo, as de-regulated by maternal diabetes, plays an important role in the pathogenesis of diabetes-induced birth defects [2, 42].

Implications for prevention of adverse outcomes from diabetic pregnancies

High glucose levels during critical periods of morphogenesis appear to be the major teratogen in diabetic pregnancy. In experimental animals, excess glucose is sufficient to cause dysmorphogenesis of embryos in glucose-injected dams or in whole embryo culture [11, 4345]. The precise mechanism(s) by which hyperglycemia induces diabetic embryopathy is(are) not clear, although involvement of the Glut2 (Slc2a2) transporter has been demonstrated [45]. Several studies report increased oxidative stress in embryos in a diabetic environment, and the administration of antioxidants, such as vitamins C or E, can reduce the occurrence of developmental defects [13, 17, 46, 47]. Which genes are functionally involved in these responses in diabetes-exposed embryos, and which mechanisms provide for the protective effect of anti-oxidant treatment in diabetic embryopathy remains to be investigated, but it is likely that one or more of the genes we have identified constitute targets in the antioxidant response. Similarly, folate supplementation has been shown to be protective against NTDs in diabetic pregnancies [46, 48]. Interestingly, the gene encoding platelet derived growth factor receptor α (Pdgfrα), mutants of which exhibit neural tube defects [49], is folate-responsive in mice [50]. Genes whose expression is altered in diabetes-exposed embryos thus represent excellent candidates for folate-responsive genes, and may mediate the beneficial effect of folate in the prevention of neural tube and other developmental defects.

Methods

Animals

Diabetes was induced in 7–9 week old female FVB mice by two intraperitoneal injections of 100 mg/kg body weight Streptozotocin in 50 mM sodium citrate buffer at pH4.5 (STZ; Sigma, St. Louis, MO) within a one-week interval. The dams were set up for mating no earlier than 7 days after the last injection, and the day of detection of a vaginal plug was counted as day 0.5 of gestation. We used embryos only from dams (n = 11) whose blood glucose levels exceeded 250 mg/dl; average glucose levels were 148 mg/dl(± 18) before STZ treatment, 337 mg/dl(± 79) on the day of mating, and 528 mg/dl(± 70) on the day of embryo harvest.

Microarray Analysis

Total RNA was isolated from embryos at embryonic day 10.5 (E10.5) using Trizol® (Invitrogen, Carlsbad, CA). We processed 2 controls and five diabetes-exposed embryos; the latter group included two specimen with neural tube defects (NTD) so as to capture the full phenotype spectrum of diabetes-exposure in pregnancy. Individual RNA samples (5 μg) from whole embryos were reverse transcribed (Invitrogen) and labeled (Affymetrix, Santa Clara, CA). In Experiment I, samples were individually hybridized to 7 Affymetrix430A2.0 chips, which were scanned using a GeneChip3000 scanner; Affymetrix GCOS imaging software was used for quality control. In Experiment II, equal amounts of RNA prepared from 3 individual embryos were pooled into one sample; each embryo was from a different pregnancy and three such pools were constructed for a total of 9 control embryos, and independently, for 9 diabetes-exposed embryos; all embryos were morphologically normal. Expression levels and "Present", "Marginal", "Absent" flags were determined with default parameters through comparison of matched and mismatched oligonucleotides for the respective gene sequence.

Statistical analyses were performed using GeneSpring7 (Silicon Genetics, Redwood City, CA) and CyberT [51]http://cybert.microarray.ics.uci.edu/. We grouped the data for control embryos and those for diabetes-exposed embryosm respectively, and filtered results in three steps: (i) "expression", i.e. "present" or "marginal" in at least one of seven samples; (ii) "statistical significance" between control and experimental samples of P < 0.05 in both CyberT and the t-test in GeneSpring; and (iii) "fold change", i.e. difference between control and diabetes-exposed samples of beyond either two-fold or 1.5-fold. The rationale for employing complementary data analysis packages and details of data transformation have been described elsewhere [52].

Of 22690 probe sets present on the arrays http://www.affymetrix.com/products_services/arrays/specific/mouse430a_2.affx, 15364 probes exhibited signals in at least one of the 7 samples, and 302 probe sets differed by more than two fold between these samples. Of these, 180 probe sets passed the t-test in GeneSpring (P < 0.05; Welch's test assuming unequal variances; false discovery rate set at 0.05), and 174 probes yielded P-values below P < 0.05 in Cyber-T. Permutation of the order of tests (significance first, fold-change second) identified differential signals from 2262 probe sets (Cyber-T), with 575 probes displaying differences between 1.5 and 2-fold, and 174 probes with differences greater than 2-fold between controls and exposed embryos. Regardless of order of filtering criteria, identical sets of probes were recovered, thus validating the analysis process. Between Cyber-T and GeneSpring, 145 sets passed both statistical filters, and after removal of duplicates, 126 genes were found to be differentially expressed above the 2-fold cut-off criterion.

For the second experiment, Affymetrix Mouse 430 2.0 arrays were used, which contain 45101 probe sets http://www.affymetrix.com/products_services/arrays/specific/mouse430_2.affx. 29687 probe sets exhibited signals in at least one of 6 samples; differences reached statistical significance at p < 0.05 for 9835 probes in the t-test (P < 0.05; Welch's test assuming unequal variances; false discovery rate set at 0.05) implemented in GeneSpring (GX version 9). 5915 probe sets exhibited differences greater than 1.5-fold, with 2796 differentially expressed greater than 2-fold. Cyber-T identified 13770 probe sets with statistical significance, of which 3992 exceeded the 1.5-fold change level, and an additional 4601 exceeded the greater than 2-fold criterion. After removal of internal controls, 5688 probe sets passed the filtering criteria for statistical significance in both Cyber-T and GeneSpring and exhibited >1.5-fold change between experimental and control samples, of which 2634 probe sets were identified with greater than 2-fold differential expression. Reduction of duplicates for a given gene was done by judgement call factoring in signal intensity, P-value, distribution of calls ("Present" was judged as more reliable than "Marginal") and fold-change; only one entry per gene was retained for a total of 2231 transcripts with differential expression greater than 2-fold.

The primary data files are available at the NCBI Gene Expression Omnibus repository (Accession number GSE9675).

Quantitative Real-time PCR

Quantitative Real-Time PCR (Q-RT-PCR) using an ABI Prism7000 instrument was performed as described [53] on cDNA samples from individual diabetes-exposed embryos (5 litters) and controls (4 litters), or pools of 4–5 control embryos from the same litter (4 litters) isolated at E10.5 (for details, see legend to Table 2). At E9.5, 6 control and 9 diabetes-exposed embryos were selected from 3 litters each, respectively, and E8.5 embryos were from 4 litters (10 controls) and 5 litters (9 diabetes-exposed embryos). All embryos used for Q-RT-PCR were morphologically normal. Normalization was done to Polymerase epsilon 4 (Pole4) cDNA in the same sample; Pole4 levels were unaffected by maternal diabetes on Experiment I and Experiment II arrays. Differences between samples (n = individual embryos except where noted otherwise) were evaluated for statistical significance using an unpaired two-tailed t-test. Primers (Additional file 5) were positioned across exon-exon junctions to exclude amplification of potentially contaminating DNA. The amplification products were designed to originate from a different region of the mRNA than that detected by probes on the microarray, in order to provide independent confirmation of expression measurements.

Annotation for tissue expression, function and mutant mouse phenotypes

Information on gene expression in embryos, where available, was collected from MGI http://www.informatics.jax.org. Molecular function attributes were based on GO-annotation (NetAffx™ https://www.affymetrix.com/analysis/netaffx/index.affx, updated as of July 21, 2008), supplemented with information from ENSEMBL and UCSC genome browsers and PubMed. Information on mutant phenotypes was obtained from MGI (as of October 21, 2008) for null and conditional alleles.

Transcription factor binding site prediction

Whole Genome rVISTA http://genome.lbl.gov/vista/index.shtml was used to identify transcription factor binding sites that are conserved between mouse and human and are over-represented in the 5 Kb upstream regions of our maternal diabetes affected genes relative to all 5 Kb upstream regions in the human genome (P-value < 0.006).

References

  1. Kucera J: Rate and type of congenital anomalies among offspring of diabetic women. J Reprod Med. 1971, 7 (2): 73-82.

    CAS  PubMed  Google Scholar 

  2. Martinez-Frias ML: Epidemiological analysis of outcomes of pregnancy in diabetic mothers: identification of the most characteristic and most frequent congenital anomalies. Am J Med Genet. 1994, 51: 108-113. 10.1002/ajmg.1320510206.

    Article  CAS  PubMed  Google Scholar 

  3. Reece EA, Ma XD, Zhao Z, Wu YK, Dhanasekaran D: Aberrant patterns of cellular communication in diabetes-induced embryopathy in rats: II, apoptotic pathways. Am J Obstet Gynecol. 2005, 192 (3): 967-972. 10.1016/j.ajog.2004.10.592.

    Article  PubMed  Google Scholar 

  4. Phelan SA, Ito M, Loeken MR: Neural tube defects in embryos of diabetic mice: role of the Pax-3 gene and apoptosis. Diabetes. 1997, 46 (7): 1189-1197. 10.2337/diabetes.46.7.1189.

    Article  CAS  PubMed  Google Scholar 

  5. Wentzel P, Eriksson UJ: A diabetes-like environment increases malformation rate and diminishes prostaglandin E(2) in rat embryos: reversal by administration of vitamin E and folic acid. Birth Defects Res A Clin Mol Teratol. 2005, 73 (7): 506-511. 10.1002/bdra.20145.

    Article  CAS  PubMed  Google Scholar 

  6. Wentzel P, Welsh N, Eriksson UJ: Developmental damage, increased lipid peroxidation, diminished cyclooxygenase-2 gene expression, and lowered prostaglandin E2 levels in rat embryos exposed to a diabetic environment. Diabetes. 1999, 48 (4): 813-820. 10.2337/diabetes.48.4.813.

    Article  CAS  PubMed  Google Scholar 

  7. Piddington R, Joyce J, Dhanasekaran P, Baker L: Diabetes mellitus affects prostaglandin E2 levels in mouse embryos during neurulation. Diabetologia. 1996, 39 (8): 915-920. 10.1007/BF00403910.

    Article  CAS  PubMed  Google Scholar 

  8. Goldman AS, Baker L, Piddington R, Marx B, Herold R, Egler J: Hyperglycemia-induced teratogenesis is mediated by a functional deficiency of arachidonic acid. Proc Natl Acad Sci USA. 1985, 82 (23): 8227-8231. 10.1073/pnas.82.23.8227.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Sussman I, Matschinsky FM: Diabetes affects sorbitol and myo-inositol levels of neuroectodermal tissue during embryogenesis in rat. Diabetes. 1988, 37 (7): 974-981. 10.2337/diabetes.37.7.974.

    Article  CAS  PubMed  Google Scholar 

  10. Khandelwal M, Reece EA, Wu YK, Borenstein M: Dietary myo-inositol therapy in hyperglycemia-induced embryopathy. Teratology. 1998, 57 (2): 79-84. 10.1002/(SICI)1096-9926(199802)57:2<79::AID-TERA6>3.0.CO;2-1.

    Article  CAS  PubMed  Google Scholar 

  11. Wentzel P, Wentzel CR, Gareskog MB, Eriksson UJ: Induction of embryonic dysmorphogenesis by high glucose concentration, disturbed inositol metabolism, and inhibited protein kinase C activity. Teratology. 2001, 63 (5): 193-201. 10.1002/tera.1034.

    Article  CAS  PubMed  Google Scholar 

  12. Yang X, Borg LA, Eriksson UJ: Altered metabolism and superoxide generation in neural tissue of rat embryos exposed to high glucose. Am J Physiol. 1997, 272 (1 Pt 1): E173-180.

    CAS  PubMed  Google Scholar 

  13. Reece EA, Wu YK, Zhao Z, Dhanasekaran D: Dietary vitamin and lipid therapy rescues aberrant signaling and apoptosis and prevents hyperglycemia-induced diabetic embryopathy in rats. Am J Obstet Gynecol. 2006, 194 (2): 580-585. 10.1016/j.ajog.2005.08.052.

    Article  CAS  PubMed  Google Scholar 

  14. Li R, Chase M, Jung SK, Smith PJ, Loeken MR: Hypoxic stress in diabetic pregnancy contributes to impaired embryo gene expression and defective development by inducing oxidative stress. Am J Physiol Endocrinol Metab. 2005, 289 (4): E591-599. 10.1152/ajpendo.00441.2004.

    Article  CAS  PubMed  Google Scholar 

  15. Sakamaki H, Akazawa S, Ishibashi M, Izumino K, Takino H, Yamasaki H, Yamaguchi Y, Goto S, Urata Y, Kondo T, et al: Significance of glutathione-dependent antioxidant system in diabetes-induced embryonic malformations. Diabetes. 1999, 48 (5): 1138-1144. 10.2337/diabetes.48.5.1138.

    Article  CAS  PubMed  Google Scholar 

  16. Sivan E, Lee YC, Wu YK, Reece EA: Free radical scavenging enzymes in fetal dysmorphogenesis among offspring of diabetic rats. Teratology. 1997, 56 (6): 343-349. 10.1002/(SICI)1096-9926(199712)56:6<343::AID-TERA1>3.0.CO;2-X.

    Article  CAS  PubMed  Google Scholar 

  17. Cederberg J, Siman CM, Eriksson UJ: Combined treatment with vitamin E and vitamin C decreases oxidative stress and improves fetal outcome in experimental diabetic pregnancy. Pediatr Res. 2001, 49 (6): 755-762. 10.1203/00006450-200106000-00007.

    Article  CAS  PubMed  Google Scholar 

  18. Chang TI, Horal M, Jain SK, Wang F, Patel R, Loeken MR: Oxidant regulation of gene expression and neural tube development: Insights gained from diabetic pregnancy on molecular causes of neural tube defects. Diabetologia. 2003, 46 (4): 538-545.

    CAS  PubMed  Google Scholar 

  19. Zangen SW, Ryu S, Ornoy A: Alterations in the expression of antioxidant genes and the levels of transcription factor NF-Kappa B in relation to diabetic embryopathy in the Cohen Diabetic rat model. Birth Defects Res A Clin Mol Teratol. 2006, 76 (2): 107-114. 10.1002/bdra.20227.

    Article  CAS  PubMed  Google Scholar 

  20. Greene MF: Diabetic embryopathy 2001: moving beyond the "diabetic milieu". Teratology. 2001, 63: 116-118. 10.1002/tera.1021.

    Article  CAS  PubMed  Google Scholar 

  21. Epstein DJ, Vekemans M, Gros P: Splotch (Sp2H), a mutation affecting development of the mouse neural tube, shows a deletion within the paired homeodomain of Pax-3. Cell. 1991, 67 (4): 767-774. 10.1016/0092-8674(91)90071-6.

    Article  CAS  PubMed  Google Scholar 

  22. Epstein DJ, Vogan KJ, Trasler DG, Gros P: A mutation within intron 3 of the Pax-3 gene produces aberrantly spliced mRNA transcripts in the splotch (Sp) mouse mutant. Proc Natl Acad Sci USA. 1993, 90 (2): 532-536. 10.1073/pnas.90.2.532.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  23. Morgan SC, Relaix F, Sandell LL, Loeken MR: Oxidative stress during diabetic pregnancy disrupts cardiac neural crest migration and causes outflow tract defects. Birth Defects Res A Clin Mol Teratol. 2008, 82 (6): 453-463. 10.1002/bdra.20457.

    Article  CAS  PubMed  Google Scholar 

  24. Morgan SC, Lee HY, Relaix F, Sandell LL, Levorse JM, Loeken MR: Cardiac outflow tract septation failure in Pax3-deficient embryos is due to p53-dependent regulation of migrating cardiac neural crest. Mech Dev. 2008, 125 (9–10): 757-767. 10.1016/j.mod.2008.07.003.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Pavlinkova G, Salbaum JM, Kappen C: Wnt signaling in caudal dysgenesis and diabetic embryopathy. Birth Defects Res A Clin Mol Teratol. 2008, 82: 710-719. 10.1002/bdra.20495.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. Stearne PA, Pietersz GA, Goding JW: cDNA cloning of the murine transferrin receptor: sequence of trans-membrane and adjacent regions. J Immunol. 1985, 134 (5): 3474-3479.

    CAS  PubMed  Google Scholar 

  27. Yang L, Lanier ER, Kraig E: Identification of a novel, spliced variant of CREB that is preferentially expressed in the thymus. J Immunol. 1997, 158 (6): 2522-2525.

    CAS  PubMed  Google Scholar 

  28. Gray PA, Fu H, Luo P, Zhao Q, Yu J, Ferrari A, Tenzen T, Yuk DI, Tsung EF, Cai Z, et al: Mouse brain organization revealed through direct genome-scale TF expression analysis. Science. 2004, 306 (5705): 2255-2257. 10.1126/science.1104935.

    Article  CAS  PubMed  Google Scholar 

  29. Schug J, Schuller W-P, Kappen C, Salbaum JM, Bucan M, Stoeckert CJ: Promoter Features Related to Tissue Specificity as Measured by Shannon Entropy. Genome Biology. 2005, 6 (4): R33-10.1186/gb-2005-6-4-r33.

    Article  PubMed Central  PubMed  Google Scholar 

  30. Reece EA, Homko CJ, Wu YK, Wiznitzer A: The role of free radicals and membrane lipids in diabetes-induced congenital malformations. J Soc Gynecol Investig. 1998, 5 (4): 178-187. 10.1016/S1071-5576(98)00008-2.

    Article  CAS  PubMed  Google Scholar 

  31. Sakamaki H, Akazawa S, Ishibashi M, Izumino K, Takino H, Yamasaki H, Yamaguchi Y, Goto S, Urata Y, Kondo T, et al: Significance of glutathione-dependent antioxidant system in diabetes-induced embryonic malformations. Diabetes. 1999, 48 (5): 1138-1144. 10.2337/diabetes.48.5.1138.

    Article  CAS  PubMed  Google Scholar 

  32. Iyer NV, Kotch LE, Agani F, Leung SW, Laughner E, Wenger RH, Gassmann M, Gearhart JD, Lawler AM, Yu AY, et al: Cellular and developmental control of O2 homeostasis by hypoxia-inducible factor 1 alpha. Genes Dev. 1998, 12 (2): 149-162. 10.1101/gad.12.2.149.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. Semenza GL: Targeting HIF-1 for cancer therapy. Nat Rev Cancer. 2003, 3 (10): 721-732. 10.1038/nrc1187.

    Article  CAS  PubMed  Google Scholar 

  34. Loffredo CA, Wilson PD, Ferencz C: Maternal diabetes: an independent risk factor for major cardiovascular malformations with increased mortality of affected infants. Teratology. 2001, 64 (2): 98-106. 10.1002/tera.1051.

    Article  CAS  PubMed  Google Scholar 

  35. Becerra JE, Khoury MJ, Cordero JF, Erickson JD: Diabetes mellitus during pregnancy and the risks for specific birth defects: a population-based case-control study. Pediatrics. 1990, 85 (1): 1-9.

    CAS  PubMed  Google Scholar 

  36. Jovanovic L, Knopp RH, Kim H, Cefalu WT, Zhu XD, Lee YJ, Simpson JL, Mills JL: Elevated pregnancy losses at high and low extremes of maternal glucose in early normal and diabetic pregnancy: evidence for a protective adaptation in diabetes. Diabetes Care. 2005, 28 (5): 1113-1117. 10.2337/diacare.28.5.1113.

    Article  PubMed  Google Scholar 

  37. Harris MJ, Juriloff DM: Mouse mutants with neural tube closure defects and their role in understanding human neural tube defects. Birth Defects Res A Clin Mol Teratol. 2007, 79 (3): 187-210. 10.1002/bdra.20333.

    Article  CAS  PubMed  Google Scholar 

  38. Marguerat S, Jensen TS, de Lichtenberg U, Wilhelm BT, Jensen LJ, Bahler J: The more the merrier: comparative analysis of microarray studies on cell cycle-regulated genes in fission yeast. Yeast. 2006, 23 (4): 261-277. 10.1002/yea.1351.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  39. Jiang B, Kumar SD, Loh WT, Manikandan J, Ling EA, Tay SS, Dheen ST: Global gene expression analysis of cranial neural tubes in embryos of diabetic mice. J Neurosci Res. 2008, 86 (16): 3481-3493. 10.1002/jnr.21800.

    Article  CAS  PubMed  Google Scholar 

  40. Reece EA, Pinter E, Homko C, Wu Y-K, Naftolin F: The Yolk Sac Theory: Closing the Circle on Why Diabetes-Associated Malformations Occur. J Soc Gynecol Investig. 1994, 1 (1): 3-13.

    CAS  PubMed  Google Scholar 

  41. Reece EA, Ji I, Wu YK, Zhao Z: Characterization of differential gene expression profiles in diabetic embryopathy using DNA microarray analysis. Am J Obstet Gynecol. 2006, 195 (4): 1075-1080. 10.1016/j.ajog.2006.05.054.

    Article  CAS  PubMed  Google Scholar 

  42. Goto MP, Goldman AS: Diabetic embryopathy. Curr Opin Pediatr. 1994, 6 (4): 486-491. 10.1097/00008480-199408000-00023.

    Article  CAS  PubMed  Google Scholar 

  43. Fine EL, Horal M, Chang TI, Fortin G, Loeken MR: Evidence that elevated glucose causes altered gene expression, apoptosis, and neural tube defects in a mouse model of diabetic pregnancy. Diabetes. 1999, 48 (12): 2454-2462. 10.2337/diabetes.48.12.2454.

    Article  CAS  PubMed  Google Scholar 

  44. Kumar SD, Dheen ST, Tay SS: Maternal diabetes induces congenital heart defects in mice by altering the expression of genes involved in cardiovascular development. Cardiovasc Diabetol. 2007, 6 (1): 34-10.1186/1475-2840-6-34.

    Article  PubMed Central  PubMed  Google Scholar 

  45. Li R, Thorens B, Loeken MR: Expression of the gene encoding the high-Km glucose transporter 2 by the early postimplantation mouse embryo is essential for neural tube defects associated with diabetic embryopathy. Diabetologia. 2007, 50 (3): 682-689. 10.1007/s00125-006-0579-7.

    Article  CAS  PubMed  Google Scholar 

  46. Gareskog M, Eriksson UJ, Wentzel P: Combined supplementation of folic acid and vitamin E diminishes diabetes-induced embryotoxicity in rats. Birth Defects Res A Clin Mol Teratol. 2006, 76 (6): 483-490. 10.1002/bdra.20278.

    Article  PubMed  Google Scholar 

  47. Reece EA, Wu YK: Prevention of diabetic embryopathy in offspring of diabetic rats with use of a cocktail of deficient substrates and an antioxidant. Am J Obstet Gynecol. 1997, 176 (4): 790-797. 10.1016/S0002-9378(97)70602-1.

    Article  CAS  PubMed  Google Scholar 

  48. Wentzel P, Gareskog M, Eriksson UJ: Folic acid supplementation diminishes diabetes- and glucose-induced dysmorphogenesis in rat embryos in vivo and in vitro. Diabetes. 2005, 54 (2): 546-553. 10.2337/diabetes.54.2.546.

    Article  CAS  PubMed  Google Scholar 

  49. Pickett EA, Olsen GS, Tallquist MD: Disruption of PDGFRalpha-initiated PI3K activation and migration of somite derivatives leads to spina bifida. Development. 2008, 135 (3): 589-598. 10.1242/dev.013763.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  50. Spiegelstein O, Cabrera RM, Bozinov D, Wlodarczyk B, Finnell RH: Folate-regulated changes in gene expression in the anterior neural tube of folate binding protein-1 (Folbp1)-deficient murine embryos. Neurochem Res. 2004, 29 (6): 1105-1112. 10.1023/B:NERE.0000023597.37698.13.

    Article  CAS  PubMed  Google Scholar 

  51. Baldi P, Long AD: A Bayesian framework for the analysis of microarray expression data: regularized t-test and statistical inferences of gene changes. Bioinformatics. 2001, 17 (6): 509-519. 10.1093/bioinformatics/17.6.509.

    Article  CAS  PubMed  Google Scholar 

  52. Kappen C, Pavlinkova G, Kruger C, Salbaum JM: Analysis of altered gene expression in diabetic embryopathy. Comprehensive Toxicology. Edited by: McQueen CA. 2008, Oxford, United Kingdom: Elsevier, 2

    Google Scholar 

  53. Kruger C, Talmadge C, Kappen C: Expression of folate pathway genes in the cartilage of Hoxd4 and Hoxc8 transgenic mice. Birth Defects Res A Clin Mol Teratol. 2006, 76 (4): 216-229. 10.1002/bdra.20245.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We are grateful for technical assistance by Diane Costanzo, Dana S'aulis, and the UNMC microarray core facility, which received support from the NCRR through P20RR016469 and P20RR018788. We are also grateful to Drs. Claudia Kruger and Daniel Geschwind (UCLA) for advice on Q-RT-PCR and microarray interpretation, respectively. G.P. was funded through a supplement to RO1-HD34706 to C.K., and J.M.S. was funded through RO1-HD055528. All authors have read and approved the final version of the manuscript.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Claudia Kappen.

Additional information

Authors' contributions

GP performed microarray analyses and PCR assays, collected information for annotations and drafted a first version of the manuscript, JMS oversaw the statistical analysis by Cyber-T, designed experiment II and led the annotation effort, CK conceived of the study, performed the annotation for experiment II and wrote the manuscript.

Electronic supplementary material

12864_2009_2158_MOESM1_ESM.doc

Additional file 1: Genes with altered expression in diabetes-exposed embryos. The file contains a list of genes identified by microarray analysis. (DOC 170 KB)

12864_2009_2158_MOESM2_ESM.doc

Additional file 2: Transcription factor binding sites overrepresented in promoters of genes affected by maternal diabetes. The file contains a list of putative transcription factor binding sites and respective references. (DOC 116 KB)

12864_2009_2158_MOESM3_ESM.doc

Additional file 3: Known HIF1 target genes that exhibit altered expression in diabetes-exposed embryos. The file contains a list of HIF1 target genes and respective references. (DOC 174 KB)

12864_2009_2158_MOESM4_ESM.doc

Additional file 4: References for in vivo function of genes affected in diabetes exposed embryos. The file contains a list of references for the in vivo function of particular genes. (DOC 489 KB)

12864_2009_2158_MOESM5_ESM.doc

Additional file 5: Primer sequences for quantitative RT-PCR. The file contains a list of RefSeqIDs, position information and sequences for primers to amplify particular genes. The amplification rate for each primer pair is also listed. (DOC 60 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Pavlinkova, G., Salbaum, J.M. & Kappen, C. Maternal diabetes alters transcriptional programs in the developing embryo. BMC Genomics 10, 274 (2009). https://doi.org/10.1186/1471-2164-10-274

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-10-274

Keywords