Functional interactions between the LRP6 WNT co ... - Oxford Academic

3 downloads 3739 Views 389KB Size Report
Sep 15, 2010 - For Permissions, please email: [email protected]. Human Molecular ..... blast cultures. Interestingly, the level of FA ... a new application of Bayesian hierarchical selection and inter- action testing that ...
Human Molecular Genetics, 2010, Vol. 19, No. 23 doi:10.1093/hmg/ddq384 Advance Access published on September 15, 2010

4560–4572

Functional interactions between the LRP6 WNT co-receptor and folate supplementation Jason D. Gray 1,2,{, Ghunwa Nakouzi 3,{, Bozena Slowinska-Castaldo 1, Jean-Eudes Dazard 4,5, J. Sunil Rao 4,{, Joseph H. Nadeau 3,5,} and M. Elizabeth Ross 1,2,∗ 1

Laboratory of Neurogenetics and Development, 2Graduate Program in Neuroscience, Weill Medical College of Cornell University, New York, NY 10021, USA, 3Department of Genetics, 4Department of Epidemiology and Biostatistics, and 5Center for Proteomics and Bioinformatics, Case Western Reserve University, Cleveland, OH 44106, USA Received July 26, 2010; Revised and Accepted September 3, 2010

Crooked tail (Cd) mice bear a gain-of-function mutation in Lrp6, a co-receptor for canonical WNT signaling, and are a model of neural tube defects (NTDs), preventable with dietary folic acid (FA) supplementation. Whether the FA response reflects a direct influence of FA on LRP6 function was tested with prenatal supplementation in LRP6-deficient embryos. The enriched FA (10 ppm) diet reduced the occurrence of birth defects among all litters compared with the control (2 ppm FA) diet, but did so by increasing early lethality of Lrp62/2 embryos while actually increasing NTDs among nulls alive at embryonic days 10 – 13 (E10 – 13). Proliferation in cranial neural folds was reduced in homozygous Lrp62/2 mutants versus wild-type embryos at E10, and FA supplementation increased proliferation in wild-type but not mutant neuroepithelia. Canonical WNT activity was reduced in LRP6-deficient midbrain –hindbrain at E9.5, demonstrated in vivo by a TCF/LEFreporter transgene. FA levels in media modulated the canonical WNT response in NIH3T3 cells, suggesting that although FA was required for optimal WNT signaling, even modest FA elevations attenuated LRP5/6dependent canonical WNT responses. Gene expression analysis in embryos and adults showed striking interactions between targeted Lrp6 deficiency and FA supplementation, especially for mitochondrial function, folate and methionine metabolism, WNT signaling and cytoskeletal regulation that together implicate relevant signaling and metabolic pathways supporting cell proliferation, morphology and differentiation. We propose that FA supplementation rescues Lrp6Cd/Cd fetuses by normalizing hyperactive WNT activity, whereas in LRP6-deficient embryos, added FA further attenuates reduced WNT activity, thereby compromising development.

INTRODUCTION Of the more than 190 genetic mouse models of neural tube defects (NTDs) (1 –3), remarkably few have been tested for responsiveness to prenatal folic acid (FA) supplementation and fewer still show a beneficial response (1 – 3). Mice with the naturally occurring Crooked tail (Cd) mutation are an example of a folate-responsive model of rostral NTDs (exencephaly) (4). Cd results from a single-nucleotide

substitution in the low-density lipoprotein receptor-related protein 6 (Lrp6) gene (5). Lrp6 (arrow in Drosophila) encodes a single-pass transmembrane receptor that, like its paralog Lrp5, is a co-receptor with Frizzled in the canonical wingless (WNT) signaling pathway to activate b-catenin-dependent TCF/LEF-regulated transcription (6). The Lrp6Cd mutation replaces a highly conserved amino acid in the extracellular domain of LRP6 and interferes with the ability of dickkopf-1 (Dkk1) to inhibit canonical WNT



To whom correspondence should be addressed at: Weill Cornell Medical College, 1300 York Avenue, PO Box 239, New York, NY 10065, USA. Tel: +1 2127465550; Fax: +1 2127468226; Email: [email protected] These authors contributed equally to the work. ‡ Present address: Division of Biostatistics, Department of Epidemiology and Public Health, University of Miami Miller School of Medicine, Miami, FL 33136, USA. } Present address: Institute for Systems Biology, Seattle, WA 98103, USA. †

# The Author 2010. Published by Oxford University Press. All rights reserved. For Permissions, please email: [email protected]

Human Molecular Genetics, 2010, Vol. 19, No. 23

4561

signaling, resulting in sustained elevation of cytosolic b-catenin levels in the presence of WNT. Loss of Lrp6 function produces several birth defects including caudal axis truncation and limb deformities as well as exencephaly and spina bifida (7). Moreover, a hypomorphic allele of Lrp6, caused by the ringelschwanz point mutation, is associated with spina bifida in otherwise viable homozygous pups (8). Thus, both gain-of-function (GOF) and loss-of-function (LOF) Lrp6 mutations can result in NTDs. Prenatal dietary FA supplementation reduces the incidence of NTDs in Lrp6Cd/Cd mice (4). Vitamin treatments result in a true rescue of Cd mice from NTDs as there was a normal Mendelian distribution of genotypes at all supplementation levels, and phenotypes in homozygous animals were shifted to increased embryonic viability and decreased NTD occurrence with FA supplementation. Moreover, analyses of gene expression profiles and biochemical markers of FA/HCY metabolism indicate a defect in intracellular FA utilization associated with homozygosity for the Lrp6Cd mutation (9). However, it is not yet known whether this effect is specific to the Lrp6Cd allele, the Cd mouse strain background, or whether FA supplementation influences the action of LRP6 itself. In this study, the response of LRP6-deficient mice to FA supplementation was examined in fetal outcomes, in biochemical and transcriptional activity in vitro and by gene expression in vivo, in order to characterize the impact of FA on NTD incidence, LRP6 function and WNT canonical signaling.

RESULTS Folate supplementation and birth defect occurrence in Lrp6 null mice

Figure 1. Incidence of NTD and morphogenic defects in pregnancies of FA-supplemented Lrp6+/2 mice. Observed implantations were scored as resorptions, embryo dead at harvest (early lethality) or live embryo with open cranial folds (exencephalic), open caudal neuropore (spina bifida), axis truncation or limb defects and graphed as a percentage of total observed implantations. (A) Lrp62/2 embryos at E13.5 that were maintained on the 2 ppm diet demonstrate each of the phenotypes scored. The embryo at left shows tail truncation and limb defects (arrowheads) and spina bifida (arrow), while the one at right exhibits limb, trunk and exencephaly (arrow) defects. (B) The incidence of each phenotype is shown as a percentage of all conceptions, with individual embryos often displaying more than one defect. The incidence of neural tube, limb and axis defects was lower among embryos maintained on the high FA (10 ppm) diet. (C) Frequency of genotypes recovered from dams on a high or low folate diet. Litters from

Two hundred and eighty-seven embryos collected from dams that were fed diets containing 2 ppm FA (130 embryos, 29 resorptions) or 10 ppm FA (157 embryos, 29 resorptions) were examined and genotyped. Most deformities were observed in Lrp62/2 embryos with only two instances of spina bifida occurring in heterozygotes, one from each diet group. No instances of exencephaly, limb defects or axis truncation were found in any wild-type or heterozygous embryos. All Lrp6 homozygous null embryos exhibited at least two defects and several displayed three or more defects. For example, two Lrp62/2 embryos from the 2 ppm FA diet group at E12.5 –13 displayed axis truncation and limb defects, in addition to coloboma of the eye and spina bifida in one embryo and exencephaly in another (Fig. 1A). When expressed as a percent of conceptions, fewer cases of exencephaly, spina bifida, limb malformations and axis truncations were encountered in the litters maintained on the 10 ppm compared with the 2 ppm FA diet (Fig. 1B). Taken

either diet showed a significant deviation from the Mendelian ratio of 25:50:25 expected from Lrp6 +/2 × Lrp6 +/2 matings. Parallel diet studies from Lrp6+/+ matings showed a resorption rate of 10% among conceptions on the 10 ppm FA diet (Table 3). After adjusting the number of Lrp62/2 in each group by adding resorptions above the 10% baseline (green hatched bars), the 2 ppm diet approximated the predicted Mendelian distribution, but nulls were still significantly under represented in the 10 ppm FA group. ∗ P , 0.05; dP , 0.001.

4562

Human Molecular Genetics, 2010, Vol. 19, No. 23

Table 1. Genotype distribution of embryos from Lrp6+/2 × Lrp6+/2 intercrosses Dietary FA (ppm) 2 10

No. of observed +/+ +/2

D Obs/expected Lrp62/2 (%)

P-valuea

No. of litters

Average litter size

Range

2/2

37 55

20 16

246 271

0.04 0.0001

20 22

7.94 8.45

4– 12 2– 12

73 86

a

Chi-square analysis confirms the deviation from the null hypothesis of Mendelian allele frequency for heterozygous parents 1:2:1, setting the wild-type number as 1.

as a group, these deformities were found with significantly lower frequency in the higher FA diet cohort than in the 2 ppm FA-fed controls (P , 0.05). However, this was not a true rescue effect, because the number of Lrp62/2 genotypes recovered in the 2 ppm FA and especially in the 10 ppm diet cohort deviated significantly from the Mendelian expectations for a Lrp6+/2 × Lrp6+/2 intercross (Fig. 1C). Even when all resorptions above a baseline 10% (see below) are assumed to represent nulls, the recovery of Lrp62/2 embryos was significantly lower than expected in the 10 ppm FA diet group (Fig. 1C), suggesting that FA supplementation increased the severity of embryopathy occurring before implantation in Lrp62/2 mice. Tables 1 –3 further highlight these FA effects, with resorptions set aside as a separate group, as DNA could not be systematically recovered from necrotic tissue. In contrast to the expected equal numbers of homozygous wild-types and mutants among all progeny, mice displaying early embryonic or later lethality that genotyped as Lrp62/2 accounted for 46% fewer homozygous nulls than expected (P ¼ 0.04) in the 2 ppm FA diet group and more severe deviation from expected in the 10 ppm FA cohort, with 71% fewer nulls observed (P ¼ 0.0001) (Table 1). Interestingly, the deficit of recovered null embryos did not affect overall litter size, as there were no significant differences between diet groups in the average number or in the range of embryos per litter. Nevertheless, a significant shift toward early lethality was found in Lrp62/2 embryos in the FA-supplemented diet relative to control (Table 2). While on the 2 ppm FA diet, nearly all (90%) observed embryos were viable at E13, only slightly over a third (38%) of the Lrp6 null embryos on the 10 ppm FA diet were viable to that gestational age (P ¼ 0.003). Despite the lower than expected number of nulls in the 10 ppm diet group, the number of resorptions was equivalent between the 2 ppm diet (18% of implantations) and the 10 ppm diet (16% of implantations) (Table 3, P ¼ 5.1). The comparable resorption rates of wild-type embryos from Lrp6+/+ × Lrp6+/+ crosses on the 2 versus 10 ppm FA diets ruled out the possibility that the 10 ppm FA diet reduced fetal viability on the wild-type background (Table 3). The 10 ppm group from Lrp6+/2 × Lrp6+/2 parents displayed a higher rate of resorption than was observed for Lrp6+/+ × Lrp6+/+ intercrosses (9.3% at 10 ppm for Lrp6+/+ versus 15.6% for Lrp6+/2 intercrosses), suggesting that resorptions above 10% on the high folate diet were genotype-dependent. Not all resorptions could be genotyped, but those from which DNA could be recovered genotyped as nulls. These results suggest that some embryonic loss in the 10 ppm FA group occurred prior to implantation or too early in embryogenesis to be detected at E13.

Table 2. Early lethality of Lrp62/2 embryos Dietary FA (ppm)

No. of observeda Viable Non-viable embryos embryos

Viable among observed Lrp62/2 embryos (%)

2 10

18 6

90 38

2 10

P-valueb

0.003

Lrp62/2 embryos at E13 that were alive versus dead with development halted at least 24 to 48 h earlier. b Fisher’s exact test confirms the association between FA concentration and early lethality of Lrp62/2 embryos. a

Together, these data indicate that FA supplementation reduced the frequency of birth defects as a group (exencephaly and spina bifida, limb malformation and axis truncation) observed at E13 and associated with LRP6 deficiency when compared with all embryos maintained on 2 ppm FA. However, increasing dietary FA reduced these rates by shifting outcomes to earlier embryonic lethality. FA supplementation increases occurrence of NTDs in Lrp6 null mutants The initial cohort of embryos screened at E13 showed a trend suggesting a higher proportion of viable nulls display NTD when receiving the 10 ppm FA diet (6/6 nulls on 10 ppm FA had NTDs compared with 12/18 nulls on the 2 ppm FA diet). However, the shift to early embryonic lethality of Lrp62/2 embryos on the supplemented FA diet reduced the number of viable nulls in that group, impeding evaluation of statistical significance. To further explore the effect of FA on neurulation in LRP6-deficient mice, additional embryos on either diet were scored specifically for NTDs at a slightly earlier stage (E10). In these litters, all wild-type siblings displayed closed neural tubes and mutants that were scored were viable and stage appropriate. Among viable Lrp62/2 embryos from dams fed with the control diet, 68% displayed NTDs while 100% of viable null embryos on the 10 ppm FA diet had NTDs (P ¼ 0.03) (Table 4). Thus, increasing prenatal FA exacerbates embryopathy and increases NTD occurrence in LRP6-deficient mice. Effect of FA on cell proliferation and apoptosis in the LRP6-deficient neural tube A hypothesized influence of FA supplementation on NTD mechanisms is that added FA enhances cell division in the neural folds and perhaps underlying mesenchyme to promote neural tube closure (10,11). This could potentially correlate

Human Molecular Genetics, 2010, Vol. 19, No. 23

4563

Table 3. Resorption frequencies by FA diet in heterozygous (Lrp6+/2 × Lrp6+/2 ) versus wild-type (Lrp6+/+ × Lrp6+/+ ) intercrosses Dietary FA (ppm)

2 10

Lrp6+/2 × Lrp6+/2 Litters Implantations 20 22

159 186

Resorption rate (%) 18 16

P-valuea

Lrp6+/+ × Lrp6+/+ Litters Implantations

NS

9 8

71 71

Resorption rate (%) 18 10

P-valuea

NS

a

Chi-square shows no association between FA concentration and resorptions within the intercross; NS, not significant.

Table 4. Effect of maternal FA supplementation on incidence of NTDs in Lrp62/2 embryos Dietary FA (ppm)

No. of observeda No NTDs NTDs

% of Lrp62/2 embryos with NTD (%)

2 10

7 0

68 100

15 15

P-valueb

0.03

Lrp62/2 embryos at E13 or E10 that were alive and stage-appropriate. Fisher’s exact test confirms the association between FA diet and NTD in Lrp62/2 embryos.

a

b

with altered WNT signaling, because this developmentally crucial pathway is known to promote patterning and cell proliferation (12,13). We therefore examined cell proliferation in the cranial neural tube of Lrp6 wild-type, heterozygous and homozygous null embryos by quantifying the number of antiphosphohistone 3 (PH3) labeled cells in M-phase per unit length along the cranial neural folds of E10 mice (Fig. 2). In Lrp6+/+ embryos, FA supplementation significantly increased cell proliferation. Interestingly, this FA effect was abolished in heterozygous and homozygous mutant embryos, in which proliferation on either the 2 or 10 ppm FA diet was unchanged or reduced relative to wild-type embryos on the 2 ppm control diet, indicating an interaction between gene dosage and FA supplementation. Changes in WNT activity could also have an impact on cell survival (14 – 16). We therefore examined TUNEL labeling of apoptotic cells in E10 siblings (Fig. 3). There was no difference in the distribution or amounts of programmed cell death among genotype or diet conditions. While these assays sample a rapidly dynamic process at this stage making cell counts inherently inaccurate (and so counts were not performed), nevertheless it can be said that there were no qualitative genotype/diet-dependent differences in apoptosis as evaluated by TUNEL. Obviously, embryo survival overall is affected since earlier death ensues. However, increased apoptosis does not appear to be a central reason that neurulation fails in Lrp62/2 animals that are viable past E10. Canonical WNT signaling in the presence of defined FA levels We next explored whether FA supplementation modulates canonical responses to recombinant WNT3a stimulation of NIH3T3 cells in vitro. The accumulation of cytosolic b-catenin was examined in cells treated with FA

concentrations ranging from 0 to 100 mg FA/ml DMEM (Fig. 4A and B). Inclusion of FA in the defined medium did not alter the basal levels of cytosolic b-catenin. Significantly, FA depletion suppressed responses to WNT stimulation, with optimal response observed at 4 mg/ml FA. In contrast, 10 mg/ ml FA or higher concentrations attenuated WNT responses. To distinguish the cytosolic and nuclear consequences of stabilized b-catenin, WNT-stimulated transcriptional activity was examined using NIH3T3 cells transfected with the TopFlash reporter plasmid (Fig. 4C). A similar sensitivity of WNT activity to FA levels was found. These data indicate that FA levels can modulate cellular responses to WNT stimulation by influencing cytosolic accumulation of b-catenin and TCF/ LEF transcriptional activity. The importance of Lrp6 gene dosage to canonical WNT activity in the rostral neural tube was confirmed using the TCF/LEF-LacZ reporter mouse on the Lrp62/2 background, which showed a substantial decrease in the in situ b-galactosidase reaction product in the midbrain – hindbrain of LRP6-deficient embryos (Fig. 4D). Thus, continued Lrp5 gene expression does not fully compensate for LRP6 loss in the region of cranial neural tube closure point 2, which is a common location of neurulation failure in these mutants. Gene and diet effects on mRNA profiles The observation of FA influence on WNT-mediated transcription led us to examine which transcriptional pathways are associated with these FA effects on neurualtion in LRP6-deficient embryos. In particular, cranial tissue from wild-type and Lrp6+/2 E9.5 female embryos was examined to gain insight into the effects of maternal FA supplementation on molecular mechanisms leading to NTDs at a key developmental stage. In parallel, liver samples from adult female mice were used for several reasons. First, liver is the tissue where FA/homocysteine metabolism has been most comprehensively studied (17 – 19). Second, gene expression in adult heterozygous females is relevant to maternal – fetal interactions, because certain maternal metabolic features are established risk factors in the pathogenesis of NTDs (20– 22). Third, gene expression patterns identified in adult tissues from Lrp6+/2 heterozygotes may be relevant for screening assays designed to evaluate NTD risk and response to FA supplementation. Finally, Lrp6+/2 heterozygotes rather than Lrp62/2 homozygotes were used as the test group to focus on subtle primary molecular features that lead to pathological changes and to avoid secondary effects of severely compromised development in homozygous mutants (23).

4564

Human Molecular Genetics, 2010, Vol. 19, No. 23

Figure 3. FA diet supplementation does not alter the pattern of apoptosis in the E10 rostral neural tube of Lrp6-deficient embryos. Panels compare the dorsal midbrain of wild-type (WT), heterozygous (+/2) and homozygous (2/2) E10 embryos from each diet group. TUNEL labeling (arrows) in all embryos was largely confined to the dorsal tips of the neural folds where apoptosis normally occurs to facilitate intercalation of cells and promote midline fusion. Tissue was lightly counterstained with hematoxylin.

Figure 2. FA diet supplementation fails to promote proliferation in the neural tube of E10 Lrp6-deficient embryos. Panels compare the dorsal midbrain of wild-type (WT), heterozygous (+/2) and homozygous (2/2) E10 embryos from each diet group. (A) Cells in M-phase, labeled with anti-PH3 (green) and counterstained with DAPI (blue). (B) Stereological assessments of PH3labeled/total cells are compared by genotype. The labeling index of mitotic cells in Lrp6+/2 embryos is the same as wild-type but is significantly reduced in the Lrp62/2 embryos compared among siblings maintained on 2 ppm FA diet. On the 10 ppm FA diet, proliferation in the cranial neural folds was nearly doubled in WT, but was unchanged in heterozygous Lrp6+/2 embryos and was further suppressed in the homozygous nulls. n ¼ 16 WT, 21 +/2 and 14 2/2 embryos; ∗ P , 0.05, two-tailed t-test, with respect to WT on 2 ppm FA.

Of the nearly 40 000 sequences covering the mouse transcriptome that were interrogated, we focused on the curated consensus list of 1650 genes for embryos and 960 genes for liver that showed evidence for statistical significance in the genotype – diet interaction test (co-variance) and that appeared in at least one ingenuity pathway analysis (IPA) metabolic or signaling pathway. For metabolic pathways, significant overrepresentation was found for 6 of the 35 pathways in embryonic cranial tissue (Supplementary Material, Table S1) and for 36 of 67 pathways in adult liver (Supplementary Material, Table S2), whereas for signaling pathways, significant over-representation was found for 67 of 160 pathways in embryonic cranial tissues (Supplementary Material, Table S3) and for 23 of 139 pathways in adult liver (Supplementary Material, Table S4). Then we examined the ‘intersection’ set, where IPA was conducted for differentially expressed genes that were found in both embryonic cranial tissue and adult liver. Four of the 28 total metabolic pathways showed statistical significance in this intersection set (Supplementary Material, Table S5), whereas 28 of the 129 signaling pathways showed statistical significance in the same intersection set (Supplementary Material, Table S6). The ways in which selected genes co-varied in wild-type and heterozygous Lrp6+/2 mice on the 2 and 10 ppm levels of dietary FA are illustrated in Figure 5 (also see Supplementary Material, Figs S1 and S2 for graphics of all interaction effects, and Supplementary Material, Tables S7 and S8 for gene lists and gene names). We identified the three highest-ranking metabolic and signaling pathways and highlighted several other significant pathways based on prior knowledge. Several observations are relevant. At least five groups of functionally related pathways were found (Table 5). These include (with IPA pathway designations provided in parentheses): methionine and folate metabolism (purine metabolism, methionine metabolism, one carbon pool by folate, polyamine regulation), WNT and ubiquitination (protein ubiquitination, WNT/b-catenin signaling), mitochondrial dysfunction (oxidative phosphorylation, ubiquinone biosynthesis, mitochondrial dysfunction), translation control (regulation of eIF2 and p70S6K signaling, eIF2 signaling) and cytoskeletal regulation (breast cancer regulation, CCR3 signaling in eosinophils and actin cytoskeleton signaling). Collectively, these pathways affect aspects

Human Molecular Genetics, 2010, Vol. 19, No. 23

4565

DISCUSSION

Figure 4. Impact of FA on the canonical WNT signaling pathway in vitro and demonstration of reduced WNT activity in situ of E9.5 Lrp62/2 embryos. (A) Western blot of cytosolic b-catenin levels in NIH3T3 cells that were cultured for 3 days in folate-supplemented defined media prior to 2 h stimulation with recombinant WNT3a. (B) Quantification of cytosolic b-catenin in lysates from cells, using Erk2 as protein loading controls. Normalized optical density values are expressed as a percentage of the b-catenin level in unstimulated cells after 3 days in 0 mg/ml FA. (C) WNT-regulated transcriptional activity measured by TopFlash reporter assay. NIH3T3 cells were cultured for 3 days in folate-supplemented defined media prior to transfection with the TOP and FOP reporter plasmids. Normalized luminescence values are expressed as fold increase above unstimulated (Basal) cells maintained in 0 FA media over the course of the experiment. n ¼ 6 separate experiments, each run in duplicate, two-tailed t-test in both B and C. (D) Reduced canonical WNT activity in cranial folds of TCF-LacZ reporter mice on a Lrp6 2/2 background. E9.5 embryos from double heterozygous Lrp6+/2 ::TCF-LacZ+ mating pairs were fixed and processed for b-galactosidase reactions developed for 30 minutes. Lrp6+/+ ::TCF-LacZ embryos (left panel) showed high levels of activity in the midline fusion of the cranial folds (arrows) and mesencephalon. Lrp62/2 ::TCF-LacZ embryos (middle panel) showed almost no WNT activity in the midline cranial folds and much reduced activity in the mesencephalon, with residual activity attributable to remaining Lrp5 function. Lrp62/2 embryos (right panel), in which a LacZ sequence interrupts the Lrp6 locus, demonstrated no visible b-galactosidase activity at 30 minutes.

of cell proliferation, cell morphology and oxidative stress, and implicate specific molecular and biochemical functions that corroborate the in vivo and in vitro results for WNT signaling and cell proliferation. We also note that inositol metabolism was significantly altered in adult liver; inositol supplementation has been shown to reduce NTD risk in curlytail mutant mice (24).

We show that the occurrence of birth defects associated with a LOF mutation in Lrp6 is reduced among all conceptions by moderate supplementation of FA in the prenatal diet. However, in contrast to our studies in the Cd mouse line that bears a GOF mutation, the same level of FA supplementation that rescued Lrp6Cd/Cd embryos did not rescue LRP6-deficient embryos. Instead, FA supplementation of Lrp62/2 embryos shifted an already severe phenotype to earlier embryonic lethality, so that LRP6-deficient embryos responded, but with an adverse reaction, to elevated FA levels. The fact that FA levels influence WNT signaling was evident in the effect of FA supplementation on gene expression profiles in Lrp6+/2 embryos and adults. In fact, of the hundreds of pathways and networks that are traceable with IPA algorithms, it is highly significant that cell cycle regulatory and WNT pathways displayed a strong gene – diet interaction. The relevance of WNT signaling was validated using two in vitro assays of the WNT canonical pathway. In these assays, FA deficiency blunted pathway responses to WNT signaling, reflecting the importance of FA metabolism for early development (3,25). However, these assays also indicate that FA supplementation, even to a modestly elevated level, attenuates WNT signaling. These data suggest that added dietary FA likely rescued the Lrp6Cd defect by countering this net hyperactive Cd allele, whereas FA supplementation further suppressed WNT signaling to exacerbate the LOF in the Lrp6 knockout line. This indicates a direct relationship between FA supplementation and the LRP6/WNT signaling pathway. Such a direct relationship would be consistent with our previous studies indicating that the Lrp6Cd allele also impacts intracellular FA utilization (9). The mitotic index of PH3 labeling in neuroepithelial cells around the time of cranial neural tube closure (E9.5 –10) showed that raising dietary levels of FA promotes proliferation in wild-type mice. However, Lrp6 LOF mutant embryos displayed an impaired proliferative response to FA supplementation. While these results could explain the failure of added FA to rescue NTDs in Lrp62/2 embryos, they could not explain the increased loss of heterozygous mice, since there was no statistical difference in mitoses of heterozygotes compared with wild-type on the control diet or between diets in the heterozygous cranial neural folds. Moreover, Lrp6 null embryos of equivalent stage did not display a statistically significant difference in mitoses on the two diets, undermining the possibility that a proliferation defect was the sole reason for earlier embryonic loss in Lrp6 nulls receiving the 10 ppm FA diet. Similarly, TUNEL labeling ruled out an increase in apoptosis upon FA supplementation to account for the observed increase in NTD in nulls and fetal loss of heterozygous LOF mutants. Involvement of LRP6 in canonical WNT signaling suggested that FA effects on neurulation could entail transcriptional influences on embryogenesis. Indeed canonical WNT signaling in the cranial neural folds and midbrain of E9.5 embryos was diminished as assessed by the comparison of in situ b-galactosidase activity in the Lrp6+/+ ::TCF-LacZ versus Lrp62/2 ::TCF-LacZ double mutants, when sibling embryos were processed in parallel. Therefore, the loss of

4566

Human Molecular Genetics, 2010, Vol. 19, No. 23

Figure 5. Examples of gene–diet interactions. Genetic factors are shown for wild-type (WT) and mutant heterozygotes (KO), and dietary factors are shown for diets that were 2 (dashed lines) and 10 ppm (solid lines) in folate. Normalized means are indicated as arbitrary units distributed evenly on the Y-axis. Examples are shown for oxidative phosphorylation in embryos (A and B) and adults (C and D), purine, methionine and inositol metabolism in embryos (E and F) and adults (G and H), cytoskeletal regulation in embryos (which also co-varies in adults) (I and J) and WNT/b-catenin signaling in embryos (K and L). Abbreviations: ATP5D (ATP-synthase, H+ transporting, mitochondrial F1 complex, delta subunit); NDUFS8 [NADH dehydrogenase (ubiquinone) Fe-S protein 8, 23 kDa (NADH-coenzyme Q reductase)]; ATP5A1 (ATP-synthase, H+ transporting, mitochondrial F1 complex, alpha subunit 1); BACE1 (beta-site APP-cleaving enzyme); KIF1B (kinesin family member 1B); DNMT3B [DNA (cytosine-5)-methyltransferase 3 beta]; COX6C (cytochrome oxidase subunit VIc); APEX1 [APEX nuclease 1 (multifunctional DNA repair enzyme)]; PPP1CA (protein phosphatase 1 catalytic subunit alpha); MYL6 (myosin light chain 6); WNT8A (Wnt ligand 8A); FZD3 (Frizzled 3). (The complete list of examples of gene– diet interactions is provided in Supplementary Material, Figures S1 and S2.)

LRP6 in null embryo cranial folds was not fully compensated by continued LRP5 expression. In vitro assessment of WNT canonical pathway stimulation in the presence of varied FA levels revealed a bimodal relationship between WNT signaling and FA concentration, in which FA excess was as deleterious to WNT signaling as FA deficiency. Particularly striking was that relatively modest FA increases of 2.5-fold or more above the optimal level can attenuate WNT signaling in vitro. There was no effect of FA alone on the basal levels of cytosolic b-catenin or TopFlash transcriptional reporter activity in NIH3T3 fibroblast cultures. Interestingly, the level of FA that provided optimal WNT activation in vitro, 4 mg/ml, is the same concentration used in the complete Dulbecco’s modified Eagle’s medium that was previously optimized for growth and maintenance of cell lines (26). FA concentrations of 10 mg/ml or higher blunted the cellular response to recombinant WNT3a in both b-catenin stabilization and TopFlash reporter assays. Therefore, FA supplementation has an impact on transcriptional activation by the canonical WNT signaling pathway. We examined gene – diet interactions because of our interest in the impact of dietary folate supplementation on NTDs and expression profiles in wild-type and Lrp6 heterozygous mutants in two key tissues at different life stages. Our expression studies were focused on E9.5 crania because

NTDs arise in embryonic cranial folds where key aspects of folate metabolism in mitochondria take place (27) and on adult liver where most of the methylation reactions and 50% of all methionine metabolism occurs (28). In addition, liver function in both embryos and adults depends on WNT signaling (29,30), and methionine metabolism in liver can be adversely affected in mice that have mutations in WNT signaling pathway genes (9,23). The statistical analysis used here is a new application of Bayesian hierarchical selection and interaction testing that provides a robust and highly reproducible result. Our interrogation of gene – diet interactions revealed several metabolic and signaling pathways that are significantly affected by Lrp6 mutation and FA supplementation in cranial development during embryogenesis, in liver function in adults or in both. Understanding the gene- and diet-dependent interplay among these physiological pathways may lead to general strategies for preventing NTDs, regardless of their responsiveness to folate. We propose that mutations in WNT signaling genes such as Lrp6 not only adversely affect methionine metabolism, but also secondarily affect methylation reactions, DNA synthesis (purine and pyrimidine metabolism), translation control and, through polyamine metabolism, the balance of cell proliferation, mitochondrial functions and perhaps metabolism of alternative methyl donors (Table 5). This proposal nicely integrates our molecular and pathway studies with known functions of these pathways.

Human Molecular Genetics, 2010, Vol. 19, No. 23

4567

Table 5. Pathways most significantly associated with diet–genotype interaction Class and pathway

P-valuea

Genesb

Embryo: metabolism 1. Oxidative phosphorylation

4.7 × 10210

NDUFA4, ATP6V0B, COX6A1, ATP5D, ATP6V1C1, NDUFB5, COX8A, COX5B, UQCRB, NDUFB8, NDUFA1, NDUFB9, NDUFAB1, ATP5H , ATP5J2, UQCRFS1, ATP6V1G1, ATP5F1, ATP5I, ATP5G3, COX4I1, ATP5J, NDUFC1, COX7A2, COX17, COX6B1, ATP5A1, COX6C, NDUFV3, NDUFA13, PPA1, ATP5C1, NDUFV2, ATP5B, NDUFS8, UQCRC2, COX5A, NDUFA3, UQCRC1, NDUFA7 ATP5D, KIF1B, NME2, POLRMT, ATP11B, APRT, IMPDH2, POLR3A, DDX3Y, ATP5H, PSMC2, ATP5F1, ATP5G3, ATP5J, NME4, RP2, POLR1B, ATP5A1, RRM2, SMARCA5, PSMD6, RALBP1, NT5C2, ATP5C1, INO80, ATP5B, PDE6D, NME1, ATP6V0B, DDX39, DHX16, VCP, ATP5J2, ATP6V1G1, ACTC1, DHX15, ATP5I, PKM2, ABCD3, PRIM2, POLG, ADCY6, PAICS, DHX8, CASK, HSPD1, ATIC, BAT1, GMPS, PDE5A ALDH2, MAOB, MIF, DHCR24, PRDX1, ALDH1A3, ELOVL2, GOT2, PRDX6, ELOVL6, PRDX2 NDUFA4, NDUFC1, NDUFV3, NDUFB5, NDUFB8, NDUFA1, NDUFA13, NDUFB9, UFSP2, NDUFV2, NDUFS8, NDUFAB1, NDUFA3, NDUFA7 DNMT3B, SRM, AMD1, MAT2A, MTAP, AHCY

2. Purine metabolism

6.4 × 1025

3. Phenylalanine metabolism 4. Ubiquinone biosynthesis

2.6 × 1024 7.7 × 1023

8. Methionine metabolism

0.06

Embryo: signaling 1. Mitochondrial dysfunction

3.6 × 1029

2. Protein ubiquitination

1.4 × 1028

3. Regulation of eIF4 and p70S6K signaling

4.2 × 1028

4. eIF2 signaling

5.6 × 1028

37. Wnt/b-catenin signaling

8.0 × 1023

64. Polyamine regulation in colon cancer

4.9 × 1022

Adult liver: metabolism 1. Oxidative phosphorylation

7.6 × 10212

2. Fatty acid metabolism

3.7 × 1027

3. Tryptophan metabolism

1.5 × 1026

7. Ubiquinone biosynthesis

1.5 × 1024

12. Inositol metabolism

2.2 × 1023

36. One carbon pool by folate

3.8 × 1022

Adult liver: signaling 1. Mitochondrial dysfunction

5.3 × 10210

2. LPS/IL-1 mediated inhibition of RXR function

2.0 × 1029

3. Acute phase response signaling

1.4 × 1028

4. Polyamine regulation

1.4 × 1022

NDUFA4, COX6A1, NDUFB5, COX8A, COX5B, NDUFB8, UQCRB, NDUFB9, MAOB, CASP9, NDUFAB1, UQCRFS1, HTRA2, GPX4, COX4I1, ATP5J, COX17, COX6B1, COX7A2, COX6C, ATP5A1, NDUFV3, MAPK12, GPX7, NDUFA13, GSR, ATP5C1, NDUFS8, NDUFV2, ATP5B, UQCRC2, COX5A, NDUFA3, UQCRC1, NDUFA7 USP24, UBE2H, USP45, PSMA3, PSMA7, UBE2L3, UBE2N, FBXW7, HSPA5, PSMB6, USP48, TCEB2, PSMD14, PSMA2, PSMC2, PSMB4, USP28, UBB, PSMA6, USP9X, PSMD6, UBE2S, NEDD4, UBE2D1, FZR1, HSPA8, PSMD11, CBL, PSMB2, ZBTB12, PSMA5, PSMA4, ANAPC5, UBE2E2, HSP90AA1, USP46, PSMD4, UBC, USP25, UBE2E1, BIRC2, UBE2C, UBE2I EIF2S3, EIF2S2, PPP2CA, SOS2, EIF1, RPS6, PAIP2, EIF2S1, EIF4G1, EIF4E, EIF4EBP1, EIF3B, EIF4G2, EIF1AX, PPM1L, PABPC1, RPS6KB1, NRAS, EIF3F, GRB2, RAC1, EIF3E, MAPK12, PPP2R1A, EIF3C, EIF4A3, EIF4A1, EIF3I, PPP2R5E EIF2S3, EIF2S2, PPP1CC, SOS2, EIF1, PPP1CB, RPS6, EIF2S1, EIF4G1, EIF4E, EIF3B, EIF4G2, EIF1AX, EIF5, PPP1CA, PABPC1, NRAS, EIF3F, GRB2, RAC1, EIF3E, EIF3C, EIF4A3, EIF4A1, EIF3I PPP2CA, CSNK1A1, TLE1, WNT8B, CCND1, CSNK1E, SOX2, PPM1L, SOX18, SOX4, UBB, GJA1, CREBBP, ACVR1, RAC1, SOX11, APC, CSNK2A2, CDH2, PPP2R1A, FZD3, SFRP1, PPP2R5E, UBC, WNT5A SAT1, OAZ1, ODC1, APC

SDHB, ATP6V1D, ATP6V0C, ATP5D, COX6A1, NDUFB5, COX8A, COX5B, UQCRB, NDUFA1, NDUFA2, ATP5G1, NDUFC2, NDUFAB1, ATP5H (includes EG:10476), ATP5I, NDUFV1, COX7A2, COX6C, ATP5A1, ATP5C1, UCRC, ATP5B, NDUFV2, NDUFS8, COX5A, NDUFA12, SDHD, NDUFA3, UQCRC1, NDUFA7, COX15 CYP2D26, ECHS1, CYP2B9, CYP2C70, CYP2C44, ACAD9, GCDH, CYP2D13, ACAA1B, HADHB, ALDH2, SLC27A5, CYP2B13, CYP1A2, CYP3A4, ALDH1A7, ACSL5, ACAT1, CYP2C19, CYP2B6, SDS, HADH, ALDH7A1 CYP2D26, ECHS1, CYP2B9, CYP2C70, CYP2C44, METTL7B, NEDD4, GCDH, CYP2D13, TMLHE, HADHB, ALDH2, CYP2B13, MAOB, CYP1A2, CYP3A4, ALDH1A7, CAT, ACAT1, CYP2C19, CYP2B6, SDS, HADH, ALDH7A1 NDUFV1, NDUFB5, METTL7B, NDUFA1, NDUFA2, NDUFV2, NDUFS8, NDUFC2, NDUFAB1, NDUFA12, NDUFA3, ALDH6A1, NDUFA7 HSD17B13, COX6C, TPI1, ALDOB, LSR, NDUFA12, CYP2B6, PTGR2, APEX1, ALDH6A1, IDH3B, NDUFA7, FADS1 SARDH, ALDH1L1, SHMT1, SHMT2 SDHB, COX6A1, COX5B, NDUFB5, COX8A, PSEN2, UQCRB, NDUFA2, MAOB, NDUFAB1, GPX4, NDUFV1, COX7A2, ATP5A1, COX6C, BACE1, ATP5C1, ATP5B, NDUFV2, NDUFS8, CAT, COX5A, SDHD, NDUFA12, NDUFA3, UQCRC1, NDUFA7, COX15 APOE, SLC10A1, LIPC, GSTM5, ALDH1L1, HMGCS2, FMO5, SLC27A5, MAOB, ACSL5, GSTM4, ALAS1, XPO1, FMO4, HMGCS1, SLCO1B3, ALDH6A1, ALDH7A1, MGST1, FMO2, ABCB11, ALDH8A1, GSTO1, PAPSS2, APOC1, ABCC2, CYP3A4, SULT1A1, CAT, NR5A2, CYP2B6, PLTP, CYP2C19, NDST1, SULT1B1 TCF4, ITIH3, FTL, HAMP, FN1, APOA1, APOA2, SOS2, AMBP, SOCS2, PDPK1, KRAS, CP, NR3C1, RBP1, HRG, ITIH2, SERPINA1, FGB, AGT, TTR, C1S, AHSG, IL6R, SERPINF1, C5, SERPINF2, KLKB1, HP, RRAS2, TF, RBP4 TCF4, OAZ1, KRAS, APC

Embryo –adult intersection: metabolism (four significant pathways total) 1. Oxidative phosphorylation 2.4 × 10210 COX6A1, ATP5D, NDUFB5, COX8A, COX5B, NDUFA1, ATP5C1, ATP5B, NDUFV2, NDUFAB1, COX5A, NDUFA3, UQCRC1, ATP5I Continued

4568

Human Molecular Genetics, 2010, Vol. 19, No. 23

Table 5. Continued Class and pathway 2. Ubiquinone biosynthesis 3. Glycolysis and gluconeogenesis 8. Purine metabolism

P-valuea

Genesb

8.9 × 1024 1.9 × 1022 2.5 × 1022

NDUFV2, NDUFAB1, NDUFB5, NDUFA3, NDUFA1 PGK1, AKR1A1, GPI, GAPDH ATP5C1, ATP5B, ATP5D, NME2, PSMC2, ATP11B, ATP5I, NT5C2

Embryo– adult intersection: signaling (28 significant pathways total) 1. Mitochondrial dysfunction 9.3 × 10210 COX6A1, NDUFB5, COX8A, COX5B, ATP5C1, MAOB, ATP5B, NDUFV2, NDUFAB1, COX5A, NDUFA3, GPX4, UQCRC1 ROCK1, GNB1, GNAS, SOS2, GNB2L1, TUBA4A, PPP1CA, PPP1R14B 2. Breast cancer regulation by 1.4 × 1023 Stathmin1 ROCK1, GNB1, GNAS, CFL1, GNB2L1, PPP1CA 3. CCR3 signaling in eosinophils 1.8 × 1023 ROCK1, CFL1, MYL6, ACTB, SOS2, PPP1CA, APC, TMSB4X 5. Actin cytoskeleton signaling 2.9 × 1023 OAZ1, APC 13. Polyamine regulation in colon 1.4 × 1022 cancer Numbers refer to the IPA ranking for the appropriate tissue (cranial tissue, embryo versus liver, adult) and pathway (metabolic versus signaling) (see Supplementary Materials for the complete listing of pathways and genes). Significance of association between the genes in the data set and the canonical pathway. b Genes within the interaction data set that are significantly associated with the canonical pathways. a

Pathway analysis of gene–diet interactions suggests a central role for methionine and polyamine metabolism in these developmental and physiological processes. In particular, Sadenosylmethionine (AdoMet, also known as SAM) is used for methylating a wide variety of molecules such as DNA, RNA, lipids and proteins, for polyamine synthesis, a source of AdoMet radicals and a component of riboswitches that control transcription and translation (31). Availability of AdoMet for these reactions depends heavily on dietary folate. Polyamines (putrescine, cadaverine, spermidine and spermine), which are DNA-binding organic molecules with at least two amino groups, are ubiquitous in rapidly dividing cells (32). AdoMet decarboxylase converts AdoMet to dcAdoMet, which then is used to make spermine and spermidine. Although cell growth is retarded or blocked when polyamine pools are depleted, cells rigorously maintain polyamine levels at the expense of methionine metabolism, which in turn leads to imbalanced AdoMet and nucleotide pools, thereby adversely affecting methylation as well as purine and pyrimidine metabolism (33). The interdependency between methionine metabolism and polyamine biosynthesis exacerbates the functional consequences of folate deficiency, with ramifications for DNA synthesis, cell proliferation and mitochondrial functions. Degradation of polyamines links methionine metabolism with mitochondrial function and alternative methyl donors. Degradation of spermine and spermidine uses acetyl-CoA and, under the action of succinate dehydrogenase (SDH1), leads to production of succinate, which is essential for complex II functions in mitochondria. Thus, demands of polyamine metabolism impact methionine metabolism in two ways, first by preferentially using AdoMet for polyamine synthesis rather than for methylation reactions, and second by using acetyl-CoA for polyamine degradation rather than for choline metabolism, which could affect availability of betaine, an alternative methyl donor. Together, these results and observations suggest that disruption of methionine metabolism in WNT mutant mice leads to diverse metabolic consequences that are remediated or exacerbated, depending on genetic conditions, with supplementary folate. These studies provide a novel perspective to the gene– environment interactions that contribute to NTD formation.

With respect to developmental signaling pathways, not only FA deficit but also elevated FA levels can impair WNT pathway function. Depending on individual genotype, FA supplementation can be deleterious to neurulation and embryogenesis. In the Lrp6 LOF mutants, the effects of supplementation on proliferation did not explain the negative impact of added FA on the heterozygous and homozygous mutants. Similarly, we did not detect an enhancement of apoptotic cell death in mutant embryos compared with wild-type. Thus, the impact of the detected pathway changes related to genotype and diet likely rests not in cell proliferation or cell death but in cell function and possibly cell morphology (oxidative stress and cytoskeletal regulation). Finally, these data in Lrp6 mutant mice also demonstrate that the mechanism of birth defect prevention by FA supplementation (reduction of NTD occurrence by amelioration versus increased embryonic loss) can vary not only with the genes associated with increased risk, but also with the type of mutation within those genes (i.e. LOF versus GOF). Clearly, a detailed understanding of genetic risk and interactions between folate metabolism and other developmentally important pathways will be required to optimize birth defect prevention strategies for individual families.

MATERIALS AND METHODS Animals All procedures involving animals were carried out in accordance with the National Institutes of Health Guide for the Care and Use of Laboratory Animals and were approved by the Institutional Animal Care and Use Committee at Weill Medical College of Cornell University. Gene-trap mice in which the Lrp6 locus was inactivated (7) have been backcrossed more than 12 generations to the C3H/HeJ background. Mice were housed in climate-controlled Thoren units with a 12 h light – dark cycle. NTD phenotyping Mating pairs of Lrp6 +/2 mice were maintained on a defined diet containing 10 or 2 ppm FA (Research Diets Inc., New

Human Molecular Genetics, 2010, Vol. 19, No. 23

Brunswick, NJ, USA) for two generations prior to tissue or embryo collection. Embryos from timed pregnant females were harvested at or near E13 and fixed in 4% paraformaldehyde to preserve embryos for phenotyping. Implantations were scored by visual inspection as undergoing resorption, early lethality (dead at harvest with development halted at least 24– 48 h earlier) or as live embryos that displayed developmental staging comparable with wild-type littermates and signs of active cardiovascular circulation at the time of collection. Embryos were further scored as having one or more defects including spina bifida (open caudal neural tube), exencephaly (open cranial folds), caudal axis truncation and/or limb deformities. Genotyping of embryos and postnatal pups was performed by PCR as previously described using tissue from the embryonic yolk sac or tail snips, respectively (7). Statistical comparisons using chi-square analysis and Fisher’s exact or Student’s t-tests (two-tailed) were performed with Microsoft Excel and GraphPad QuickCalcs (http://grap hpad.com/quickcalcs/chisquared1.cfm; http://graphpad.com/ quickcalcs/contingency1.cfm, accessed April 2009). The threshold for significance was P , 0.05. Immunohistochemistry and TUNEL labeling Heterozygous pairs were placed in a mating cage in the evening and were separated the following morning, which was designated embryonic day 0.5 (E0.5). Embryos from Lrp6 +/2 × Lrp6 +/2 intercrosses within each diet group were harvested on E10 and yolk sacs were collected for genotyping. Embryos were fixed in 4% paraformaldehyde overnight at 48C, and then transferred to 0.1 M PBS prior to paraffin processing (Tissue Tek 2000, Miles Laboratories). Embedded tissues were sectioned coronally at 6 mm, mounted on adhesive-coated slides (Fisher Scientific), deparaffinized and rehydrated. Antigen retrieval was performed with Reveal (BioCare Medical) and quenched in 3% hydrogen peroxide. Tissues were blocked in SNIPER (BioCare Medical) for 30 min at room temperature prior to incubation in anti-PH3 antibody (Upstate Biotechnology, 16-189, 1:1000) overnight at 48C, to specifically label dividing cells in M-phase. Immunolabeling was visualized using a secondary antibody conjugated to a fluorophore (Alexa Fluor 488, Molecular Probes, A-11070, 1:500). Slides were cover slipped with DAPI mounting media (Vectastatin) and photographed on a Nikon Optiphot-2 compound microscope fitted with a Spot Insight CCD camera (Diagnostic Instruments). For TUNEL labeling, rehydrated sections were permeabilized with 20 mg/ml proteinase K for 15 min at room temperature and quenched with 3% hydrogen peroxide. The ApoTag Plus Apoptosis Detection kit (Millipore, S7101) was used per the manufacturer’s instructions to label dying cells using the immunoperoxidase reaction. Slides were briefly (two dips) counterstained with hemotoxylin before they were dehydrated and cover slipped for imaging. b-catenin stabilization assay The levels of cytosolic b-catenin in WNT-stimulated cultured cells were measured as previously detailed (5,34), with the following modifications. Briefly, FA powder (Sigma) was

4569

dissolved in DMEM without FA (Specialty Media, Phillipsburg, NJ, USA) and culture media were prepared with 10% calf serum + DMEM with specific FA concentrations (0, 1, 4, 10, 50 or 100 mg FA/ml DMEM). Two hours after NIH3T3 cells were seeded in a six-well plate, growth media containing 10% calf serum+DMEM were replaced with culture media of specific FA concentrations. Cells were cultured for 72 h in folate-specific media prior to assays, which were run as previously described (34). Cultures were incubated in vehicle with or without 40 ng recombinant WNT3a (R&D Systems, 1324-WN) for 2 h. Following collection of the cytosolic fraction, cell lysates were normalized for protein concentration using the BCA Assay (Pierce) and western transferred to nitrocellulose membranes. Blots were blocked at room temperature with 1% non-fat dairy milk for 1 h before incubation at 48C overnight in primary antibody, either anti-b-catenin (Santa Cruz, sc-7963, 1:1000) or anti-Erk2 (Santa Cruz, sc-154, 1:5000) diluted in block solution. Blots were washed 2 × 15 min in TBS w/0.05% Tween-20 before incubation with horseradish peroxidaseconjugated secondary antibodies for 1 h at room temperature. Chemiluminescence reagents (Pierce, pico-34078 or femto-34096) were use for immunolabeled band detection. X-ray film exposures were scanned with an optical densitometer and analyzed using Quantity One software (Bio-Rad). Band densities of b-catenin were normalized for protein loading to the amount of Erk2 present in each lane. Transfection and TopFlash reporter assay A TCF/LEF in vitro reporter assay was used to examine the effect of FA on WNT-induced transcriptional activity. NIH3T3 cells (100 000 cells per well of a 24-well plate) were plated in DMEM containing 0, 1, 4, 10, 50 or 100 mg/ml FA. After 3 days, cells were transfected with 80 ng of TopFlash or FopFlash reporter plasmids using Lipofectamine according to the manufacturer’s protocol (BioRad, #170-3350). Cells were allowed to recover 36 h in the FA-defined medium after which, either 0 or 200 ng of recombinant WNT3a (R&D Systems, 1324-WN) was added per well. Cultures were incubated another 8 h prior to collection and lysis of cells for luciferase reporter assay using a Microlumat Plus luminometer with automated sample injection (EG&G Berthold, LB96V). WNT activity (TopFlash) was normalized to a WNT-insensitive luciferase control (FopFlash) plasmid. In situ assessment of canonical WNT activity A TCF/LEF-LacZ reporter mouse line (35) was crossed with Lrp6 mutants and the progeny of Lpr6+/2 ::TCF-LacZ doubleheterozygous parents were evaluated for canonical WNT activity in the cranial folds at E9.5. Embryos were fixed with 4% paraformaldehyde for 30 min at 48C, then washed three times for 20 min in 0.1 M PBS at 48C before incubation in X-gal wash buffer (0.02% Igepal, 0.01% Na-deoxycholate, 2 mM MgCl2, 100 mM Na – phosphate) for 5 min at room temperature, transfered to X-gal staining solution [1 mg/ml X-Gal, 5 mM K3Fe(CN)6, 5 mM K4Fe(CN)6 in X-Gal wash buffer] and incubated at 378C. After 30 min, embryos were rinsed in 0.1 M PBS, 3 × 5 min then 2 × 1 h and post-fixed with 4%

4570

Human Molecular Genetics, 2010, Vol. 19, No. 23

paraformaldehyde overnight at 48C. Images were collected with Leica Applications Suite software running a Leica DFC310 FX camera fitted to a Leica M165 FC stereomicroscope. Tissue collection and microarray For expression studies in adult tissues, approximately 0.5 cm3 of liver was excised from each of the experimental female mice and placed in RNA-Later (Ambion) solution at 48C overnight to prevent RNA degradation. Embryonic tissues for comparative expression studies were collected from pregnant dams at E9.5, so that cranial tissues from the otic vesicle rostrally were placed in RNA-later solution and frozen at 2808C. Remaining tissues were genotyped and RNA was extracted from female Lrp6+/+ and Lrp6+/2 littermates. The RNeasy Mini Kit (Qiagen) was used to isolate total RNA from liver and embryo tissues. This method included a DNAse treatment to digest genomic DNA followed by column purification. The RNA samples were used to make cRNA probes that were hybridized to the Affymetrix Mouse 430 2.0 expression array, which includes over 39 000 transcripts. Procedures were carried out according to Affymetrix protocols for single-round amplifications. The design consisted of two groups; the first included four of each of the Lrp6 wild-type and heterozygous animals on control 2 ppm FA, while the other group was fed a 10 ppm FA diet. This design was treated as a two-way factorial arrangement. In this design, main effects as well as the interaction effect could be evaluated. However, the interaction effect was of primary interest here as this evaluates the differential impact of FA supplementation in Lrp6 wild-type and heterozygous animals. Background correction, normalization and summarization of probe-level data The ‘robust multichip analysis’ method was used for background adjustment of perfect match (PM) probe intensities (36). With this method, model-based background-corrected PM signals were computed assuming a normal–exponential convolution model for the distribution of probe intensities. Probe intensities were corrected for global standardization, variance-stabilization and normality. The sample intensities were all brought on the same scale through standardization (scaling and centering) of each feature. To ensure that the assumptions of normality and homoscedasticity across the samples between experimental groups were met in our statistical model (e.g. when detecting differential expression), the data were transformed using the generalized log (also known as arcsinh(x)) (37,38). After transformation, the data were returned on an approximate base-2 log scale, where variances were approximately independent of the mean intensities. We applied a probe-specific probe set summarization algorithm to convert the previous normalized probe level signals into expression values. We used the ‘PM only’ probe-specific method, along with the ‘Medianpolish’ summarization method (36). Statistics on genotype – diet interaction effects Bayesian hierarchical model selection. The interaction effect of FA supplementation with Lrp6 variants was evaluated

with a specialized Bayesian hierarchical model known as ‘spike and slab’ regression, which was introduced by Lempers (39) and Mitchell and Beauchamp (40) as a Bayesian approach to model selection in linear regression models. The expression ‘spike and slab’ refers to the prior for the regression coefficients used in their Bayesian hierarchy which assumed the regression coefficients were mutually independent with a two-point mixture distribution composed of a uniform flat distribution (the slab) and a degenerate distribution at zero (the spike). George and McCulloch (41) used a scale (variance) mixture of two normal distributions, which was highly advantageous and led to an efficient Gibbs sampling algorithm that heavily popularized the spike and slab approach. These were extended to the class of rescaled spike and slab models by Ishwaran and Rao (42). Rescaling the response by the square root of the sample size divided by the root mean square error led to a non-vanishing penalization effect, which, when used in tandem with a continuous bimodal prior, ensures many useful model selection properties for the posterior mean of the regression coefficients (42,43). Recently, Ishwaran and Rao (43) looked in depth at the geometry of generalized ridge regression (GRR), a method introduced by Hoerl and Kennard (44) as a way to overcome correlated and ill-conditioned regression settings. This analysis showed that GRR possesses unique advantages in high dimensional correlated settings, i.e. the P bigger than n problem, but more interestingly, that weighted GRR (WGRR) regression could potentially be even more effective. Noting that the posterior mean from a rescaled spike and slab model is a type of WGRR estimator, they showed that these estimators, when coupled with dimension reduction, could be highly effective for prediction in high dimensional large P, small n problems. Actual variable selection is achieved with a three-step process that includes: (i) pre-filtering of variables to a manageable dimension size (because P ? n), followed by (ii) spike and slab regression estimation, which is in turn followed by (iii) a non-negative least squares step. The spike and slab estimates have a selective shrinkage property so that only variables whose effects are truly null are estimated in a way that they are shrunk towards zero, i.e. the null, and away from true signal, which is not shrunken. Here true signal is estimated with a specialized (non-shrunk) generalized ridge estimator known as the minimum least squares estimator (43). The application of non-negative least squares shrinks spike and slab-regression GRR estimates to those that correspond to noise down to zero (43). This step simply allows one to make a formal decision about which GRR effects are truly noise versus signal, without the need to set an a priori significance level. The selective shrinkage property of spike and slab regression GRR estimates allows this last step to be particularly effective. As shown by Ishwaran and Rao (43), this threestep process also has optimal variable selection properties including variable selection consistency, i.e. with probability going to 1, only true signal is identified—i.e. no false positives or negatives. Although this is an asymptotic result, it also translates into good finite sample properties. Interaction effect. An interaction effect was deemed important not only if it survived the above thresholding procedure, but

Human Molecular Genetics, 2010, Vol. 19, No. 23

also if at least one of its corresponding main effects were also found to be significant. Analysis of genotype – diet interactions identified genes that were differentially expressed between the wild-type and heterozygous states of Lrp6 and that differed between the 2 and the 10 ppm FA diets. All calculations were performed using a new R library called ‘SpikeSlab’ available from the R CRAN repository (www.cran.r-project. com). Building consensus lists of genes. In our Bayesian hierarchical model selection approach, our interest is in generating a posterior distribution of the regression coefficients, from which we focus on the posterior mean, in order to determine which genes enter the model. Because the joint posterior distributions that are generated do not have closed form solutions, the coefficient estimation process requires approximating (sampling from) these joint posterior distributions by means of the Gibbs sampler and Markov chain Monte Carlo (MCMC) algorithm. Due to the stochastic nature of the algorithm used, we obtain slightly different lists of genes for each run of the MCMC algorithm. One way around this is to build a consensus list of all the genes that appear at a certain frequency over several outputted models. In this way, we generated a consensus list of 2195 genes for the embryo data in which all were found to co-vary 100% reproducibly (20 times out of 20 MCMC runs). Similarly, we built a consensus list of 1112 genes for the liver. Curation. All interaction plots and functional pathway analyses were done with default IPA filters and after curation of the consensus lists. For the curation KIAA, RIKEN, UNKNOWN, PSEUDOGENES, PREDICTED and duplicated probes from the same gene with opposite coefficient signs were removed. Otherwise the median coefficient was taken for consistent duplicated probes. Final curated consensus lists consisted of 1650 genes for the embryo and 960 for the liver. Pathway and network analysis. IPA software (Ingenuity Systemsw, www.ingenuity.com) was used to identify biological pathways and networks defined by the interaction-significant genes. Affymetrix identifiers for these genes were uploaded into the software. Each gene identifier was mapped to its corresponding gene object in the IPA knowledge base. Differentially expressed genes were then mapped to the global molecular network that is developed from information in the Ingenuity Knowledge Base. Many of the genes that mapped to the global molecular network were associated with a canonical pathway in the Ingenuity Knowledge Base and were thus eligible for canonical pathway analysis. Significance of the association between the data set and the canonical pathway was measured in two ways: (i) a ratio of the number of differentially expressed genes from the data set that map to the pathway divided by the total number of all the genes that exist in the canonical pathway, and (ii) Fischer’s exact test was used to calculate a P-value for the probability that the association between the genes in the data set and the canonical pathway is explained by chance alone. Pathways passing the P-value threshold of 0.05 were considered significantly associated with our data

4571

set. IPA also associated many differentially expressed genes with networks in the Ingenuity Knowledge Base. Networks for these genes were then algorithmically generated based on their connectivity. A score that was the negative log of the P-value of the right-tailed Fisher’s exact test was assigned for each network. This score takes into account the number of eligible genes in our data set and the size of the network to calculate the fit between each network and the genes in the data set.

SUPPLEMENTARY MATERIAL Supplementary Material is available at HMG online. Conflict of Interest statement. None declared.

FUNDING This work was supported by the NIH (NS058979), The WCMC-TMHRI initiative to M.E.R., NIH (NRSA NS059562) to J.D.G. and NIH (HD037485) to J.H.N.

REFERENCES 1. Harris, M.J. and Juriloff, D.M. (2007) Mouse mutants with neural tube closure defects and their role in understanding human neural tube defects. Birth Defects Res. A Clin. Mol. Teratol., 79, 187– 210. 2. Harris, M.J. and Juriloff, D.M. (2010) An update to the list of mouse mutants with neural tube closure defects and advances toward a complete genetic perspective of neural tube closure. Birth Defects Res. A Clin. Mol. Teratol., 88, 653– 669. 3. Gray, J.D. and Ross, M.E. (2009) Mechanistic insights into folate supplementation from Crooked tail and other NTD-prone mutant mice. Birth Defects Res. A Clin. Mol. Teratol., 85, 314– 321. 4. Carter, M., Ulrich, S., Oofuji, Y., Williams, D.A. and Ross, M.E. (1999) Crooked tail (Cd) models human folate-responsive neural tube defects. Hum. Mol. Genet., 8, 2199–2204. 5. Carter, M., Chen, X., Slowinska, B., Minnerath, S., Glickstein, S., Shi, L., Campagne, F., Weinstein, H. and Ross, M.E. (2005) Crooked tail (Cd) model of human folate-responsive neural tube defects is mutated in Wnt coreceptor lipoprotein receptor-related protein 6. Proc. Natl Acad. Sci. USA, 102, 12843– 12848. 6. Schweizer, L. and Varmus, H. (2003) Wnt/Wingless signaling through beta-catenin requires the function of both LRP/Arrow and frizzled classes of receptors. BMC Cell Biol., 4, 4. 7. Pinson, K.I., Brennan, J., Monkley, S., Avery, B.J. and Skarnes, W.C. (2000) An LDL-receptor-related protein mediates Wnt signalling in mice. Nature, 407, 535– 538. 8. Kokubu, C., Heinzmann, U., Kokubu, T., Sakai, N., Kubota, T., Kawai, M., Wahl, M.B., Galceran, J., Grosschedl, R., Ozono, K. et al. (2004) Skeletal defects in ringelschwanz mutant mice reveal that Lrp6 is required for proper somitogenesis and osteogenesis. Development, 131, 5469– 5480. 9. Ernest, S., Carter, M., Shao, H., Hosack, A., Lerner, N., Colmenares, C., Rosenblatt, D.S., Pao, Y.H., Ross, M.E. and Nadeau, J.H. (2006) Parallel changes in metabolite and expression profiles in crooked-tail mutant and folate-reduced wild-type mice. Hum. Mol. Genet., 15, 3387–3393. 10. Copp, A.J., Greene, N.D. and Murdoch, J.N. (2003) The genetic basis of mammalian neurulation. Nat. Rev. Genet., 4, 784 –793. 11. Beaudin, A.E. and Stover, P.J. (2007) Folate-mediated one-carbon metabolism and neural tube defects: balancing genome synthesis and gene expression. Birth Defects Res. C Embryo Today, 81, 183– 203. 12. Logan, C.Y. and Nusse, R. (2004) The Wnt signaling pathway in development and disease. Annu. Rev. Cell Dev. Biol., 20, 781– 810. 13. Yang, Y. (2009) Growth and patterning in the limb: signaling gradients make the decision. Sci. Signal, 2, pe3.

4572

Human Molecular Genetics, 2010, Vol. 19, No. 23

14. Zerlin, M., Julius, M.A. and Kitajewski, J. (2008) Wnt/Frizzled signaling in angiogenesis. Angiogenesis, 11, 63–69. 15. Kubota, T., Michigami, T. and Ozono, K. (2009) Wnt signaling in bone metabolism. J. Bone Miner. Metab., 27, 265– 271. 16. Hackam, A.S. (2005) The Wnt signaling pathway in retinal degenerations. IUBMB Life, 57, 381– 388. 17. Finkelstein, J.E., Hauser, E.R., Leonard, C.O. and Brusilow, S.W. (1990) Late-onset ornithine transcarbamylase deficiency in male patients. J. Pediatr., 117, 897– 902. 18. Nijhout, H.F., Reed, M.C., Lam, S.L., Shane, B., Gregory, J.F. 3rd and Ulrich, C.M. (2006) In silico experimentation with a model of hepatic mitochondrial folate metabolism. Theor. Biol. Med. Model., 3, 40. 19. MacFarlane, A.J., Perry, C.A., Girnary, H.H., Gao, D., Allen, R.H., Stabler, S.P., Shane, B. and Stover, P.J. (2009) Mthfd1 is an essential gene in mice and alters biomarkers of impaired one-carbon metabolism. J. Biol. Chem., 284, 1533– 1539. 20. Detrait, E.R., George, T.M., Etchevers, H.C., Gilbert, J.R., Vekemans, M. and Speer, M.C. (2005) Human neural tube defects: developmental biology, epidemiology, and genetics. Neurotoxicol. Teratol., 27, 515– 524. 21. Hague, W.M. (2003) Homocysteine and pregnancy. Best Pract. Res. Clin. Obstet. Gynaecol., 17, 459–469. 22. van Guldener, C. and Stehouwer, C.D. (2001) Homocysteine-lowering treatment: an overview. Expert Opin. Pharmacother., 2, 1449– 1460. 23. Ernest, S., Christensen, B., Gilfix, B.M., Mamer, O.A., Hosack, A., Rodier, M., Colmenares, C., McGrath, J., Bale, A., Balling, R. et al. (2002) Genetic and molecular control of folate-homocysteine metabolism in mutant mice. Mamm. Genome, 13, 259– 267. 24. Greene, N.D. and Copp, A.J. (1997) Inositol prevents folate-resistant neural tube defects in the mouse. Nat. Med., 3, 60– 66. 25. Ross, M.E. (2010) Gene-environment interactions: folate metabolism and the embryonic nervous system. Wiley Interdiscip. Rev. Syst. Biol. Med., in press. 26. Dulbecco, R. and Freeman, G. (1959) Plaque production by the polyoma virus. Virology, 8, 396– 397. 27. Pike, S.T., Rajendra, R., Artzt, K. and Appling, D.R. (2010) Mitochondrial C1-tetrahydrofolate synthase (MTHFD1L) supports the flow of mitochondrial one-carbon units into the methyl cycle in embryos. J. Biol. Chem., 285, 4612–4620. 28. Mato, J.M., Corrales, F.J., Lu, S.C. and Avila, M.A. (2002) S-Adenosylmethionine: a control switch that regulates liver function. FASEB J., 16, 15– 26.

29. Thompson, M.D. and Monga, S.P. (2007) WNT/beta-catenin signaling in liver health and disease. Hepatology, 45, 1298–1305. 30. Nejak-Bowen, K. and Monga, S.P. (2008) Wnt/beta-catenin signaling in hepatic organogenesis. Organogenesis, 4, 92– 99. 31. Loenen, W.A. (2006) S-adenosylmethionine: jack of all trades and master of everything? Biochem. Soc. Trans., 34, 330–333. 32. Wallace, H.M., Fraser, A.V. and Hughes, A. (2003) A perspective of polyamine metabolism. Biochem. J., 376, 1 –14. 33. Bistulfi, G., Diegelman, P., Foster, B.A., Kramer, D.L., Porter, C.W. and Smiraglia, D.J. (2009) Polyamine biosynthesis impacts cellular folate requirements necessary to maintain S-adenosylmethionine and nucleotide pools. FASEB J., 23, 2888– 2897. 34. Giarre, M., Semenov, M.V. and Brown, A.M. (1998) Wnt signaling stabilizes the dual-function protein beta-catenin in diverse cell types. Ann. N. Y. Acad. Sci., 857, 43– 55. 35. Mohamed, O.A., Clarke, H.J. and Dufort, D. (2004) Beta-catenin signaling marks the prospective site of primitive streak formation in the mouse embryo. Dev. Dyn., 231, 416–424. 36. Irizarry, R.A., Bolstad, B.M., Collin, F., Cope, L.M., Hobbs, B. and Speed, T.P. (2003) Summaries of Affymetrix GeneChip probe level data. Nucleic Acids Res., 31, e15. 37. Huber, W., von Heydebreck, A., Sultmann, H., Poustka, A. and Vingron, M. (2002) Variance stabilization applied to microarray data calibration and to the quantification of differential expression. Bioinformatics, 18(Suppl. 1), S96–S104. 38. Durbin, B.P., Hardin, J.S., Hawkins, D.M. and Rocke, D.M. (2002) A variance-stabilizing transformation for gene-expression microarray data. Bioinformatics, 18(Suppl. 1), S105– S110. 39. Lempers, F.B. (1971) Posterior Probabilities of Alternative Linear Models. Rotterdam University Press, Rotterdam. 40. Mitchell, T.J. and Beauchamp, J.J. (1988) Bayesian variable selection in linear regression. J. Am. Stat. Assoc., 83, 1023–1036. 41. George, E.I. and McCullouch, R.E. (1993) Variable selection via Gibbs sampling. J. Am. Stat. Assoc., 88, 881 –889. 42. Ishwaran, H. and Rao, J.S. (2005) Spike and slab gene selection for multigroup microarray data. J. Am. Stat. Assoc., 100, 764 –780. 43. Ishwaran, H. and Rao, J.S. (under review) Generalized ridge regression: geometry and computational solutions when p is larger than n. Roy. Stat. Soc. Ser. B Theor. Method. 44. Hoerl, A.E. and Kennard, R.W. (1970) Ridge regression: Biased estimation for non-orthogonal problems. Technometrics, 12, 55–67.