Human Molecular Genetics Advance Access originally published online on October 18, 2006
Human Molecular Genetics 2006 15(23):3387-3393; doi:10.1093/hmg/ddl415
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Parallel changes in metabolite and expression profiles in crooked-tail mutant and folate-reduced wild-type mice
1 Department of Genetics, 2 Center for Computational Genomics and 3 Department of Electrical Engineering and Computer Science, Case Western Reserve University School of Medicine, 10900 Euclid Ave Cleveland, OH 44106 USA, 4 Department of Neurology and Neuroscience, Weill Medical College, Cornell University New York, NY, USA, 5 Departments of Human Genetics, Medicine, Pediatrics and Biology, McGill University, Montreal, Quebec, Canada and 6 Department of Cancer Biology, Lerner Research Institute, Cleveland Clinic Foundation, Cleveland, OH, USA
* To whom correspondence should be addressed. Tel: +1 2163680581; Fax: +1 2163683832; Email: emailjhn4{at}case.edu
Received July 24, 2006; Revised September 12, 2006; Accepted October 10, 2006
| ABSTRACT |
|---|
|
|
|---|
Anomalies in homocysteine (HCY) and folate metabolism are associated with common birth defects and adult diseases, several of which can be suppressed with dietary folate supplementation. Although supplementation reduces the occurrence and severity of neural tube defects (NTDs), many cases are resistant to these beneficial effects. The basis for variable response and biomarkers that predict responsiveness are unknown. Crooked-tail (Cd) mutant mice are an important model of folate-responsive NTDs. To identify features that are diagnostic for responsiveness versus resistance to dietary folate supplementation, we surveyed metabolite and expression levels in liver samples from folate-supplemented, folate-reduced and control diets in Cd mutant and wild-type adult females. Cd homozygotes had normal total homocysteine (tHcy) levels suggesting that folate suppresses NTDs through a mechanism that does not involve modulating serum tHcy levels. Instead, parallel changes in metabolite and expression profiles in folate-supplemented Cd/Cd homozygotes and folate-reduced +/+ and Cd/+ mice suggest that Crooked-tail homozygotes have a defect in the utilization of intracellular folate. Then, by combining these expression and metabolite profile results with published results for other models and their controls, two clusters were found, one of which included several folate-responsive NTD models and the other previously untested and presumably folate-resistant models. The predictive value of these profiles was verified by demonstrating that NTDs of Ski/ mutant mice, whose profile suggested resistance to folate supplementation, were not suppressed with dietary folate supplementation. These results raise the possibility of using metabolite and expression profiles to distinguish folate-responsive and resistance adult females who are at risk for bearing fetuses with an NTD.
| INTRODUCTION |
|---|
|
|
|---|
Birth defects such as neural tube defects (NTDs) and adult diseases such as vascular disease, end-stage renal disease, cognitive dysfunction and certain cancers are often associated with anomalies in homocysteine (HCY) and folate metabolism (1,2). Disease pathogenesis is thought to result directly from the detrimental effects of elevated total homocysteine (tHcy) levels, folate deficiency or perhaps indirectly from redox reactions that are mediated through HCY, cysteine and glutathione metabolism (3). tHcy and folate levels are usually inversely related (1) and dietary supplementation with folate in many, but not all cases, lowers tHCY levels, restores normal folate levels and reduces the risk for NTDs (4), colon cancer (5) and ischemic heart disease (6). Responsiveness to folate supplementation varies considerably for poorly understood reasons. With predicted responses based on genetic, molecular and metabolic features, dietary supplementation might be targeted to individuals who might benefit most and alternative treatments identified for folate-resistant cases.
NTD models in mice vary remarkably in their response to various dietary supplementations and nutrient treatments. Several models involving genes such as CART1 (7), CITED2 (8), Crooked-tail (9,10), FOLBP1 (11) and RFC1 (12) show reduced incidence and severity of NTDs after folate treatments, whereas other mutants such as axial-degeneration (13,14), bent-tail (15) and curly-tail (16) are resistant to the beneficial effects of folate supplementation. Responsiveness of some, but not all, of these mutants to dietary folate supplementation raises the possibility of identifying metabolite and expression biomarkers for differential response. These and other markers in functionally related metabolic and signaling pathways could then be tested for their relevance to NTDs in humans.
Crooked-tail (Cd) mutant mice are an important model of folate-responsive NTDs (9,10). Cd homozygous mutant mice show early lethality, exencephaly and small body weight in rare survivors; a crooked tail is found in Cd/+ heterozygotes (9). A diet supplemented with folic acid reduces early embryonic lethality and incidence of exencephaly in homozygous Cd/Cd mice (9). The Crooked-tail phenotype results from a mutation in the low density lipoprotein receptor-related protein 6 (Lrp6) gene, a co-receptor with Frizzled that is required for WNT signaling (10).
To identify diagnostic features that distinguish folate-responsiveness and -resistance in adult females, we undertook a systematic study of expression and metabolite levels in liver samples from Cd and wild-type female mice that were raised on a folate-supplemented, a folate-reduced diet or a control diet. We focused on liver samples from adult females partly because folateHCY metabolism in liver has been extensively characterized in this organ and because risk for NTDs has in some cases been attributed to maternal metabolic anomalies (2). We found that profiles for folate-supplemented Cd/Cd homozygotes were indistinguishable and that folate levels were normal in Cd/Cd homozygotes, arguing that these mice, are functionally folate-reduced. Combining in a joint analysis, these results with those of a previous study for other NTD models (17), we discovered two main clusters of strains and mutants, one that was composed of mutants that respond to folate-supplementation and another whose responsiveness was unknown. We predicted that the second cluster included folate-resistant mice. Tests with one of these mutants, which involved an engineered deficiency of the SKI protein, verified these profile predictions.
| RESULTS |
|---|
|
|
|---|
We surveyed serum levels of tHcy and folate as key metabolic indicators of the HCY and folate pathways in wild-type (+/+), heterozygous (Cd/+) and homozygous (Cd/Cd) females that were fed a control, folate-reduced or folate-supplemented diet. As expected (1,2), tHcy levels were elevated only in mice that were raised on the folate-reduced diet (Fig. 1A). Contrary to expectations, however, Cd mice had normal tHcy levels as well as normal levels of extracellular and intracellular folate compared with wild-types and heterozygotes, regardless of diet (Figs 1B and C). Therefore, the beneficial effects of folate treatments in Cd females were not mediated through a reduction in tHcy levels.
|
Two metabolites of HCY metabolism, total glutathione (tGsh) and total cysteine (tCys), were surveyed in mice of the various genotypes that were raised on the alternative diets. Glutathione (GSH) is a major cellular antioxidant and maternal glutathione may play an important role in protecting the embryo during early development by scavenging endogenous reactive oxygen species (18). tGsh levels were consistently elevated in Cd/Cd mice and tended to increase with decreasing amounts of dietary folate (Supplementary Material, Fig. S1), suggesting that the Cd mutation may adversely affect antioxidant functions in the liver. Cysteine plays numerous roles such as substrate for protein synthesis and probably as an extracellular antioxidant (19). However, no obvious patterns were observed for tCys levels among genotypes and diets (Supplementary Material, Fig. S2).
Hepatic expression profiles were used to identify genotypes and diets (as perturbations) that clustered together and to identify the patterns of metabolites and gene expression (as consequences of perturbations) that distinguished each cluster. Given the focus on folate and HCY metabolism, the miniarrays included most of the genes that are directly involved in folate and HCY metabolism as well as an equal number of unrelated control genes. Two main clusters of strains (A and B) and two main clusters of gene expression and metabolite profiles (1 and 2) were found (Fig. 2A). Cluster A included wild-type +/+ and Cd/+ mice on the control and folate-supplemented diets (4 and 10 mg FA/kg, respectively), whereas Cluster B included mutant Cd/Cd homozygotes on the control and folate-supplemented diets and wild-type +/+ and Cd/+ mice on the folate-reduced diet (0 mg FA/kg). Wild-type mice that were raised on the folate-reduced diet and Cd homozygous mutant mice that were fed the folate-supplemented shared several striking molecular features. For example, both showed significantly increased MTHFR and SHMT RNA levels and significantly reduced MTR RNA levels (Fig. 2A). Cluster 1 included genes in various metabolic pathways (HCY and folate metabolism, inflammatory response, lipid transport, phosphoinositide signaling and glucose transport and metabolism), whereas Cluster 2 included the remainder of the genes except those involved in phosphoinositide signaling and in glucose transport and metabolism. Jack-knifing, a resampling method that consists of iteratively deleting each unit of a data set (20), was used to test whether particular genes or pathways were critical for defining the clusters. Sequential elimination of each gene and separately each pathway did not affect the fundamental structure or composition of the various clusters, suggesting that clustering did not depend on results for particular genes or pathways (data not shown).
|
To test whether similarities among expression profiles might be diagnostic for responsiveness to folate supplementation, we compared expression profiles for the various Cd mice with those for Apob, Gli3, Pax3, Ptch and Ski mutant mice (Fig. 2B). Heterozygosity for these mutations adversely affects HCY levels and expression profiles for genes involved in HCY and folate metabolism (17). Interestingly, mice with phenotypic defects that are folate-responsive grouped together: one cluster contained Cd (9,10), Apc (21) and Pax3 (22,23) mutant mice, whereas the other cluster contained Apob, Ptch and Ski, whose responsiveness to folate is unknown. These results predict that, despite having anomalies in HCYfolate metabolism (17), Apob, Ptch and Ski mutant mice would be folate-resistant, or perhaps respond to folate treatments in a manner that is functionally distinct from the response of mutants in the first group.
To test this hypothesis, we treated Ski mutant mice with folic acid during the first half of gestation. Ski mutants were tested in both the C57BL/6J (N9N10 backcross generations) and 129S6 (N4 backcross generation) backgrounds to determine whether folic acid treatment reduced the incidence of either the clefting or exencephaly phenotypes observed in these strains (24). As predicted, folic acid had no effect on the incidence of either strain-dependent phenotype in Ski mutant homozygotes (Table 1).
|
| DISCUSSION |
|---|
|
|
|---|
Mutations in several genes in WNT signaling (Lrp6 and Apc) and the related sonic hedgehog (SHH) signal transduction pathway (Ski, Ptch1, Gli3 and Pax3) cause NTDs and other birth defects as well as adult diseases associated with anomalies in folate and HCY pathways (25). Previous work showed that metabolite and expression levels for HCY and folate metabolism are adversely affected in mice that are heterozygous for mutations in the WNT and SHH pathways (17). Moreover, the ability of dietary folate supplementation to suppress phenotypes in many of these mutant mice and similarly folate deficiency to adversely affect SHH signal transduction suggests a reciprocal interdependency between these metabolic and signaling pathways. These interdependencies link purine and pyrimidine synthesis, amino acid metabolism, DNA, protein and lipid methylation in HCY and folate metabolism with cell adhesion, cell death, growth factors, differentiation, proliferation, positional information and transcriptional regulation in WNT and SHH signaling (1,2,25).
The striking similarities among the metabolite and gene expression profiles in folate-supplemented Cd/Cd mice and in folate-reduced +/+ and Cd/+ mice suggest that Cd homozygotes have a defect in folate utilization despite high folate levels in the diet, serum and liver. Cd/Cd cells either do not have adequate access to intracellular folate or are not able to utilize the available folate effectively. Perhaps, by increasing folate levels with dietary supplementation, a threshold is reached that restores normal neural tube development. Alternatively, perhaps folate utilization differs between fetal metabolism, when correction occurs, and adults, when diagnostic measurements are usually made.
NTDs are a classic example of multifactorial traits involving interactions between both genetic and environmental factors. Unraveling these complexities in humans is challenging. Because both genetic and environmental factors are readily controlled mouse models enable discovery of genes, molecules, pathways and dietary supplements that can be prioritized for testing in humans. The present study and published work (17) illustrates these complexities. For example, the cluster to which the Gli3-Xt mutation belonged depends on the inbred background on which it occurred (Fig. 2B). Similarly, the cluster to which the C57BL/6J (B6) strain belonged depended on which if any mutant was congenic on that background, with B6, B6-Apc*min and B6-Pax3*Sp included in cluster 1, but B6-Apob and B6-Gli3*xt included in cluster 2 (Fig. 2B). Thus, both the identity of the mutated gene and the genetic background influence gene expression and metabolite profiles. Moreover, folate deficiency may be a necessary but not sufficient condition to cause NTDs because folate-reduced wild-type mice do not show NTDs. Thus, mutations in Lrp6 and other genes in SHH and WNT signaling may be sensitizing mutations, whose phenotypic manifestations depend in part on genetic background and in part on folate levels in the diet.
The metabolic and gene expression profiles presented here identified two clusters of strains and mutants that may be generally predictive of folic acid responsiveness of neurulation on NTD-prone backgrounds. Combined with previous observations (17), adverse phenotypes associated with mutant genes in WNT signaling (9,10,21) and SHH signaling (22,23) can be suppressed with folate treatment. These pathways provide a promising system for genetic, metabolite and expression screening for NTD risk assessments in humans. Precisely, how the unique features of these profiles affect WNT and sonic hedgehog related protein functions and the general applicability of these profiles to other NTD defects and adult diseases associated with anomalies in HCY and folate metabolism remain to be determined.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Crooked-tail mutants: study population, diets and metabolite measures
Female mice were used because in humans maternal HCY as well as fetal metabolism are risk factors for fetal NTDs (14). Virgin females were used to control for possible effects of pregnancy on HCY levels. An additional reason for focusing on females is that among human anencephalies and certain mouse exencephalies, such as Cd, females are more frequently affected than males (25). Mice were raised from conception on control [4 mg folic acid (FA)/kg chow], folate-supplemented (10 mg FA/kg chow) or folate-reduced (0 mg FA/kg chow) defined diets (P.J. Noyes, NH). Homozygous Cd/Cd mice on the folate-reduced diet were not analyzed because of their severely reduced viability (9).
Cd genotyping
The polymorphism detection protocol was described previously (9,10).
Blood and tissue samples
After deep anesthetization of 68-week-old virgin female mice, blood samples were obtained via cardiac puncture and were collected in non-heparinized tubes. After centrifugation, serum samples were stored at 80°C. At autopsy, liver samples were placed immediately on dry ice.
Metabolite measurements
The HPLC method of Ubbink and Vermaak (26) was used to measure total serum HCY, glutathione and cysteine. Total folate serum and hepatic levels were measured with a microbiological enzymatic assay (27).
Ski mutants: study population and folate treatment
Pregnant females were treated daily with folic acid (5 mg/kg body weight, by i.p. injection) from E0.5 (day of plug) to E9.5. Fetuses were collected at E16.5E18.5, examined for craniofacial phenotypes and genotyped by Southern blotting (28).
Expression arrays
Arrays of gene-specific PCR products were used to evaluate the impact of mutation and diet on the RNA abundance of genes involved in HCYfolate metabolism as well as a panel of genetically and biochemically independent genes as controls. The methods that were used for making the arrays, gene-specific reagents, hybridization conditions and data analysis were described previously (17). Supplementary Material, Table S1 lists the genes that were included on the array.
Data analysis
We used two methods to analyze the expression data: hierarchical clustering (29) and mean-shift algorithm clustering (29,30). For the hierarchical clustering, the Multiple Viewer Experiment (MeV) software package developed by The Institute of Genome Research was used. The hierarchical analysis was statistically supported by jack-knife resampling with 1000 iterations, with resampling conducted on both strains and genes separately.
The data were also analyzed using the mean-shift algorithm (30,31), an unsupervised and non-parametric data density estimation procedure. Clustering is a natural application for the mean-shift algorithm. Each data point was assigned to the closest mode (cluster) center, starting with an arbitrary assumed mode center and then estimating the gradient of the data density function. At every step, we discovered that the assumed mode center is not really the center of the distribution at that location and adjusted by shifting the temporary position of the mean in the direction of the gradient.
| SUPPLEMENTARY MATERIAL |
|---|
|
|
|---|
Supplementary Material is available at HMG Online.
| ACKNOWLEDGEMENTS |
|---|
This work was supported by grants from the NIH (HL58982 to J.H.N., HD30728 to C.C. and RO1NS34998 to M.E.R.), and by gift from the Charles B. Wang Foundation to Y.H.P. and J.H.N.
Conflict of Interest statement. There are no conflicts of interest to report.
| REFERENCES |
|---|
|
|
|---|
- Mudd S.H., Levy H.L., Skraus J.P. (2001) Disorders of transsulfuration. In Scriver C.R., Beaudet A.L., Sly W.S., Valle D. (Eds.). The Metabolic and Molecular Bases of Inherited Disease 8th edn. (McGraw Hill, New York) Vol. II: pp. 20072056.
- Rosenblatt D.S. and Fenton W.A. (2001) Inherited disorders of folate and cobalamin transport and metabolism. In Scriver C.R., Beaudet A.L., Sly W.S., Valle D. (Eds.). The Metabolic and Molecular Bases of Inherited Disease 8th edn. (McGraw Hill, New York) Vol. III: pp. 38973933.
- Heinecke J.W. (2001) Unique aspects of sulfur chemistry: homocysteine and lipid oxidation. In Carmel R. and Jacobsen D.W. (Eds.). Homocysteine in Health and Disease(Cambridge University Press, Cambridge) pp. 3238.
- MRC Vitamin Study Research Group. (1991) Prevention of neural tube defects: results of the Medical Research Council Vitamin Study. Lancet 338:131137.[CrossRef][ISI][Medline]
-
Giovannucci E., Stampfer M.J., Colditz G.A., Hunter D.J., Fuchs C., Rosner B.A., Speizer F.E., Willett W.C. (1998) Multivitamin use, folate, and colon cancer in women in the Nurses Health Study. Ann. Intern. Med. 129:517524.
[Abstract/Free Full Text] -
Schnyder G., Roffi M., Pin R., Flammer Y., Lange H., Eberli F.R., Meier B., Turi Z.G., Hess O.M. (2001) Decreased rate of coronary restenosis after lowering of plasma homocysteine levels. N. Engl. J. Med. 345:15931600.
[Abstract/Free Full Text] - Zhao Q., Behringer R.R., de Crombrugghe B. (1996) Prenatal folic acid treatment suppresses acrania and meroanencephaly in mice mutant for the Cart1 homeobox gene. Nat. Genet. 13:275283.[CrossRef][ISI][Medline]
-
Barbera J.P., Rodriguez T.A., Greene N.D., Weninger W.J., Simeone A., Copp A.J., Beddington R.S., Dunwoodie S. (2002) Folic acid prevents exencephaly in Cited2 deficient mice. Hum. Mol. Genet. 11:283293.
[Abstract/Free Full Text] -
Carter M., Ulrich S., Oofuji Y., Williams D.A., Ross M.E. (1999) Crooked tail (Cd) models human folate-responsive neural tube defects. Hum. Mol. Genet. 8:21992204.
[Abstract/Free Full Text] -
Carter M., Chen X., Minnerath S., Ross M.E. (2005) Crooked tail (Cd) model of human folate-responsive neural defects is mutated in Wnt coreceptor lipoprotein receptor-related protein 6. Proc. Natl Acad. Sci. USA 102:1284312848.
[Abstract/Free Full Text] - Piedrahita J.A., Oetama B., Bennett G.D., van Waes J., Kamen B.A., Richardson B.A., Lacey J., Anderson R.G., Finnell R.H. (1999) Mice lacking the folic acid-binding protein Folbp1 are defective in early embryonic development. Nat. Genet. 23:228232.[CrossRef][ISI][Medline]
-
Zhao R., Russell R.G., Wang Y., Liu L., Gao F., Kneitz B., Edelmann W., Goldman I.D. (2001) Rescue of embryonic lethality in reduced folate carrier-deficient mice by maternal folic acid supplementation reveals early neonatal failure of hematopoietic organs. J. Biol. Chem. 276:1022410228.
[Abstract/Free Full Text] - Essien F.B. (1992) Maternal methionine supplementation promotes the remediation of axial defects in Axd mouse neural tube mutants. Teratology 45:205212.[CrossRef][ISI][Medline]
-
Essien F.B. and Wannberg S.L. (1993) Methionine but not folinic acid or vitamin B-12 alters the frequency of neural tube defects in Axd mutant mice. J. Nutr. 123:2734.
[Abstract/Free Full Text] - Franke B., Klootwijk R., Lemmers B., de Kovel C.G., Steegers-Theunissen R.P., Mariman E.C. (2003) Phenotype of the neural tube defect mouse model bent tail is not sensitive to maternal folinic acid, myo-inositol, or zinc supplementation. Birth Defects Res. A. Clin. Mol. Teratol. 67:979984.[CrossRef][ISI][Medline]
- Greene N.D. and Copp A.J. (1997) Inositol prevents folate-resistant neural tube defects in the mouse. Nat. Med. 3:6066.[CrossRef][ISI][Medline]
- 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 folatehomocysteine metabolism in mutant mice. Mamm. Genome 13:259267.[CrossRef][ISI][Medline]
-
Legge M. and Sellens M.H. (1991) Free radical scavengers ameliorate the 2-cell block in mouse embryo culture. Hum. Reprod. 6:867871.
[Abstract/Free Full Text] -
Ueland P.M., Mansoor M.A., Guttormsen A.B., Muller F., Aukrust P., Refsum H., Svardal A.M. (1996) Reduced, oxidized and protein-bound forms of homocysteine and other aminothiols in plasma comprise the redox thiol statusa possible element of the extracellular antioxidant defense system. J. Nutr. 126:1281S1284S.
[Abstract/Free Full Text] - Lynch M. and Walsh B. (1998) Sib analysis. Genetics and Analysis of Quantitative Traits(Sinauer Associates, Sunderland) pp. 553580.
-
Song J., Medline A., Mason J.B., Gallinger S., Kim Y.I. (2000) Effects of dietary folate on intestinal tumorigenesis in the apcMin mouse. Cancer Res. 60:54345440.
[Abstract/Free Full Text] -
Fleming A. and Copp A.J. (1998) Embryonic folate metabolism and mouse neural tube defects. Science 280:21072109.
[Abstract/Free Full Text] - Li D., Pickell L., Liu L., Rosen R. (2006) Impact of methylenetetrahydrofolate reductase deficiency and low dietary folate on the development of neural tube defects in splotch mice. Birth Defects Res. A. Clin. Mol. Teratol. 76:5559.[CrossRef][ISI][Medline]
- Colmenares C., Heilstedt H.A., Shaffer L.G., Schwartz S., Berk M., Murray J.C., Stavnezer E. (2002) Loss of the SKI proto-oncogene in individuals affected with 1p36 deletion syndrome is predicted by strain-dependent defects in Ski/ mice. Nat. Genet. 30:106109.[CrossRef][ISI][Medline]
- Copp A.J., Greene N.D., Murdoch J.N. (2003) The genetic basis of mammalian neurulation. Nat. Rev. Genet. 4:784793.[CrossRef][ISI][Medline]
- Ubbink J.B., Hayward Vermaak W.J., Bissbort S. (1991) Rapid high-performance liquid chromatographic assay for total homocysteine levels in human serum. J. Chromatogr. 565:441446.[ISI][Medline]
- Molloy A.M. and Scott J.M. (1997) Microbiological assay for serum, plasma, and red cell folate using cryopreserved, microtiter plate method. Meth. Enzymol. 281:4353.[CrossRef][ISI][Medline]
-
Berk M., Desai S.Y., Heyman H.C., Colmenares C. (1997) Mice lacking the ski proto-oncogene have defects in neurulation, craniofacial, patterning, and skeletal muscle development. Genes Dev. 11:20292039.
[Abstract/Free Full Text] -
Eisen M.B., Spellman P.T., Brown P.O., Botstein D. (1998) Cluster analysis and display of genome-wide expression patterns. Proc. Natl Acad. Sci. USA 95:1486314868.
[Abstract/Free Full Text] - Fukunaga K. and Hostetler L.D. (1975) The estimation of the gradient of a density function with applications in pattern recognition. IEEE Trans. Inf. Theory IT-21:3240.[CrossRef]
-
Cheng Y. (1995) Mean shift, mode seeking, and clustering. IEEE Trans. Inf. Theory 17:790799.[CrossRef]
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

