PITX2 and non-canonical Wnt pathway interaction in metastatic prostate cancer

The non-canonical Wnt pathway, a regulator of cellular motility and morphology, is increasingly implicated in cancer metastasis. In a quantitative PCR array analysis of 84 Wnt pathway associated genes, both non-canonical and canonical pathways were activated in primary and metastatic tumors relative to normal prostate. Expression of the Wnt target gene PITX2 in a prostate cancer (PCa) bone metastasis was strikingly elevated over normal prostate (over 2,000-fold) and primary prostate cancer (over 200-fold). The elevation of PITX2 protein was also evident on tissue microarrays, with strong PITX2 immunostaining in PCa skeletal and, to a lesser degree, soft tissue metastases. PITX2 is associated with cell migration during normal tissue morphogenesis. In our studies, overexpression of individual PITX2A/B/C isoforms stimulated PC-3 PCa cell motility, with the PITX2A isoform imparting a specific motility advantage in the presence of non-canonical Wnt5a stimulation. Furthermore, PITX2 specific shRNA inhibited PC-3 cell migration toward bone cell derived chemoattractant. These experimental results support a pivotal role of PITX2A and non-canonical Wnt signaling in enhancement of PCa cell motility, suggest PITX2 involvement in homing of PCa to the skeleton, and are consistent with a role for PITX2 in PCa metastasis to soft and bone tissues. Our findings, which significantly expand previous evidence that PITX2 is associated with risk of PCa biochemical recurrence, indicate that variation in PITX2 expression accompanies and may promote prostate tumor progression and metastasis.

Selective regulation of p38β protein and signaling by integrin-linked kinase mediates bladder cancer cell migration

Integrin-linked kinase (ILK) and p38MAPK are protein kinases that transduce extracellular signals regulating cell migration and actin cytoskeletal organization. ILK-dependent regulation of p38MAPK is critical for mammalian kidney development and in smooth muscle cell migration, however, specific p38 isoforms has not been previously examined in ILK-regulated responses. Signaling by ILK and p38MAPK is often dysregulated in bladder cancer, and here we report a strong positive correlation between protein levels of ILK and p38β, which is the predominant isoform found in bladder cancer cells, as well as in patient-matched normal bladder and tumor samples. Knockdown by RNA interference of either p38β or ILK disrupts serum-induced, Rac1-dependent migration and actin cytoskeletal organization in bladder cancer cells. Surprisingly, ILK knockdown causes the selective reduction in p38β cellular protein level, without inhibiting p38β messenger RNA (mRNA) expression. The loss of p38β protein in ILK-depleted cells is partially rescued by the 26S proteasomal inhibitor MG132. Using co-precipitation and bimolecular fluorescent complementation assays, we find that ILK selectively forms cytoplasmic complexes with p38β. In situ proximity ligation assays further demonstrate that serum-stimulated assembly of endogenous ILK–p38β complexes is sensitive to QLT-0267, a small molecule ILK kinase inhibitor. Finally, inhibition of ILK reduces the amplitude and period of serum-induced activation of heat shock protein 27 (Hsp27), a target of p38β implicated in actin cytoskeletal reorganization. Our work identifies Hsp27 as a novel target of ILK–p38β signaling complexes, playing a key role in bladder cancer cell migration.

Assessment of gene expression of intracellular calcium channels, pumps and exchangers with epidermal growth factor-induced epithelial-mesenchymal transition in a breast cancer cell line

Background Epithelial-mesenchymal transition (EMT) is a process implicated in cancer metastasis that involves the conversion of epithelial cells to a more mesenchymal and invasive cell phenotype. In breast cancer cells EMT is associated with altered store-operated calcium influx and changes in calcium signalling mediated by activation of cell surface purinergic receptors. In this study, we investigated whether MDA-MB-468 breast cancer cells induced to undergo EMT exhibit changes in mRNA levels of calcium channels, pumps and exchangers located on intracellular calcium storing organelles, including the Golgi, mitochondria and endoplasmic reticulum (ER). Methods Epidermal growth factor (EGF) was used to induce EMT in MDA-MB-468 breast cancer cells. Serum-deprived cells were treated with EGF (50 ng/mL) for 12 h and gene expression was assessed using quantitative RT-PCR. Results and conclusions These data reveal no significant alterations in mRNA levels of the Golgi calcium pump secretory pathway calcium ATPases (SPCA1 and SPCA2), or the mitochondrial calcium uniporter (MCU) or Na+/Ca2+ exchanger (NCLX). However, EGF-induced EMT was associated with significant alterations in mRNA levels of specific ER calcium channels and pumps, including (sarco)-endoplasmic reticulum calcium ATPases (SERCAs), and inositol 1,4,5-trisphosphate receptor (IP3R) and ryanodine receptor (RYR) calcium channel isoforms. The most prominent change in gene expression between the epithelial and mesenchymal-like states was RYR2, which was enriched 45-fold in EGF-treated MDA-MB-468 cells. These findings indicate that EGF-induced EMT in breast cancer cells may be associated with major alterations in ER calcium homeostasis.

Characterisation of a glutathione reductase gene and its genetic locus from pea (Pisum sativum L.)

A cDNA encoding the chloroplast/mitochondrial form of glutathione reductase (GR:EC 1,6,4,2) from pea (Pisum sativum L.) was used to map a single GR locus, named GORI. In two domesticated genotypes of pea (cv, Birte and JI 399) it is likely that the GORI locus contains a single gene. However, in a semi-domesticated land race of pea sequences were detected but closely related sets of GR gene sequences were in JI 281 represent either a second intact gene or a partial or pseudogene copy. A GR gene was cloned from ev. Birte, sequenced and its structure analysed. No features of the transcription or structure of the gene suggested a mechanism for generating any more than one form of . From these data plus previously published biochemical evidence was suggested a second, distinct gene encoding for the cytosolic form of GR should be present in peas. The GORI-encoded GR mRNA can be detected in all main organs of the plant and no alternative spliced species was present which could perhaps account for the generation of multiple isoforms of GR. The mismatch between the number of charge-separable isoforms in pea and the proposed number suggests that different GR isoforms arise by some form of post-transnational modification.

Modulation of breast cancer progression and differentiation by the gp30/neregulin

In the last decade we have come to understand that the growth of cancer cells in general and of breast cancer in particular depends, in many cases, upon growth factors that will bind to and activate their receptors. One of these growth factor receptors is the erbB-2 protein which plays an important role in the prognosis of breast cancer and is overexpressed in nearly 30% of human breast cancer patients. While evidence accumulates to support the relationship between erbB-2 overexpression and poor overall survival in breast cancer, understanding of the biological consequence(s) of erbB-2 overexpression remains elusive. Our recent discovery of the gp30 has allowed us to identify a number of related but distinct biological endpoints which appear responsive to signal transduction through the erbB-2 receptor. These endpoints of growth, invasiveness, and differentiation have clear implications for the emergence, maintenance and/or control of malignancy, and represent established endpoints in the assessment of malignant progression in breast cancer. We have shown that gp30 induces a biphasic growth effect on cells with erbB-2 over-expression. We have recently determined the protein sequence of gp30 and cloned its full length cDNA sequence. We have also cloned two additional forms to the ligand, that are believed to be different isoforms. We are currently expressing the different forms in order to determine their biological effects. To elucidate the cellular mechanisms underlying cell growth inhibition by gp30, we tested the effect of this ligand on cell growth and differentiation of the human breast cancer cells which overexpress erbB-2 and cells which express low levels of this protooncogene. High concentrations of ligand induced differentiation of cells overexpressing erbB-2, as measured by inhibition of cell growth, and increased synthesis of milk components, and modulation of E-cadherin and up- regulation of c-jun and c-fos. These findings indicate that ligand-induced growth inhibition in cells overexpressing erbB-2 is associated with an apparent induction of differentiation. The availability of gp30 derived synthetic peptides and its full cDNAs provides tools necessary to acquire a better understanding of the mechanism of action of the this ligands and the erbB-2 receptor in breast cancer.

Antisense-mediated suppression of hyaluronan synthase 2 inhibits the tumorigenesis and progression of breast cancer

The progression of several cancers is correlated with the increased synthesis of the glycosaminoglycan, hyaluronan. Hyaluronan is synthesized at the plasma membrane by various isoforms of hyaluronan synthases (HAS). The importance of HAS2 expression in highly invasive breast cancer was characterized by the antisense inhibition of HAS2 (ASHAS2). The effect of HAS2 inhibition on cell proliferation, migration, hyaluronan metabolism, and receptor status was characterized in vitro, whereas the effect on tumorigenicity and metastasis was established in vivo. HAS2 inhibition resulted in a 24-hour lag in proliferation that was concomitant to transient arrest of 79% of the cell population in G 0-G1. Inhibition of HAS2 did not alter the expression of the other HAS isoforms, whereas hyaluronidase (HYAL2) and the hyaluronan receptor, CD44, were significantly down-regulated. ASHAS2 cells accumulated greater amounts of high molecular weight hyaluronan (>10,000 kDa) in the culture medium, whereas mock and parental cells liberated less hyaluronan of three distinct molecular weights (100, 400, and 3,000 kDa). The inhibition of HAS2 in the highly invasive MDA-MB-231 breast cancer cell line inhibited the initiation and progression of primary and secondary tumor formation following s.c. and intracardiac inoculation into nude mice, whereas controls readily established both primary and secondary tumors. The lack of primary and secondary tumor formation was manifested by increased survival times where ASHAS2 animals survived 172% longer than the control animals. Collectively, these unique results strongly implicate the central role of HAS2 in the initiation and progression of breast cancer, potentially highlighting the codependency between HAS2, CD44, and HYAL2 expression. ©2005 American Association for Cancer Research.

Interaction between Arabidopsis Brca2 and its partners Rad51, Dmc1, and Dss1

The Arabidopsis (Arabidopsis thaliana) orthologs of Brca2, a protein whose mutations are involved in breast cancer in humans, were previously shown to be essential at meiosis. In an attempt to better understand the Brca2-interacting properties, we examined four partners of the two isoforms of Brca2 identified in Arabidopsis (AtRad51, AtDmc1, and two AtDss1 isoforms). The two Brca2 and the two Dss1 isoforms are named AtBrca2(IV), AtBrca2(V), AtDss1(I), and AtDss1(V) after their chromosomal localization. We first show that both AtBrca2 proteins can interact with either AtRad51 or AtDmc1 in vitro, and that the N-terminal region of AtBrca2 is responsible for these interactions. More specifically, the BRC motifs (so called because iterated in the Brca2 protein) in Brca2 are involved in these interactions: BRC motif number 2 (BRC2) alone can interact with AtDmc1, whereas BRC motif number 4 (BRC4) recognizes AtRad51. The human Rad51 and Dmc1 proteins themselves can interact with either the complete (HsRad51) or a shorter version of AtBrca2 (HsRad51 or HsDmc1) that comprises all four BRC motifs. We also identified two Arabidopsis isoforms of Dss1, another known partner of Brca2 in other organisms. Although all four Brca2 and Dss1 proteins are much conserved, AtBrca2(IV) interacts with only one of these AtDss1 proteins, whereas AtBrca2(V) interacts with both of them. Finally, we show for the first time that an AtBrca2 protein could bind two different partners at the same time: AtRad51 and AtDss1(I), or AtDmc1 and AtDss1(I).

The impact of prolonged hyperinsulinaemia on glucose transport in equine skeletal muscle and digital lamellae

REASONS FOR PERFORMING STUDY An increased incidence of metabolic disease in horses has led to heightened recognition of the pathological consequences of insulin resistance (IR). Laminitis, failure of the weight-bearing digital lamellae, is an important consequence. Altered trafficking of specialised glucose transporters (GLUTs) responsible for glucose uptake, are central to the dysregulation of glucose metabolism and may play a role in laminitis pathophysiology. OBJECTIVES We hypothesised that prolonged hyperinsulinaemia alters the regulation of glucose transport in insulin-sensitive tissue and digital lamellae. Our objectives were to compare the relative protein expression of major GLUT isoforms in striated muscle and digital lamellae in healthy horses and during hyperinsulinaemia. STUDY DESIGN Randomised, controlled study. METHODS Prolonged hyperinsulinaemia and lamellar damage were induced by a prolonged-euglycaemic hyperinsulinaemic clamp (p-EHC) or a prolonged-glucose infusion (p-GI) and results were compared to electrolyte-treated controls. GLUT protein expression was examined with immunoblotting. RESULTS Lamellar tissue contained more GLUT1 protein than skeletal muscle (p = 0.002) and less GLUT4 than the heart (p = 0.037). During marked hyperinsulinaemia and acute laminitis (induced by the p-EHC), GLUT1 protein expression was decreased in skeletal muscle (p = 0.029) but unchanged in the lamellae, while novel GLUTs (8; 12) were increased in the lamellae (p = 0.03), but not skeletal muscle. However, moderate hyperinsulinaemia and subclinical laminitis (induced by the p-GI) did not cause differential GLUT protein expression in the lamellae vs. control horses. CONCLUSIONS The results suggest that lamellar tissue functions independently of insulin and that IR may not be an essential component of laminitis aetiology. Marked differences in GLUT expression exist between insulin-sensitive and insulin-independent tissues during metabolic dysfunction in horses. The different expression profiles of novel GLUTs during acute and subclinical laminitis may be important to disease pathophysiology and require further investigation.

Cloning and tissue distribution of novel splice variants of the ovine ghrelin gene

Background The ghrelin axis is involved in the regulation of metabolism, energy balance, and the immune, cardiovascular and reproductive systems. The manipulation of this axis has potential for improving economically valuable traits in production animals, and polymorphisms in the ghrelin (GHRL) and ghrelin receptor (GHSR) genes have been associated with growth and carcass traits. Here we investigate the structure and expression of the ghrelin gene (GHRL) in sheep, Ovis aries. Results We identify two ghrelin mRNA isoforms, which we have designated Δex2 preproghrelin and Δex2,3 preproghrelin. Expression of Δex2,3 preproghrelin is likely to be restricted to ruminants, and would encode truncated ghrelin and a novel C-terminal peptide. Both Δex2 preproghrelin and canonical preproghrelin mRNA isoforms were expressed in a range of tissues. Expression of the Δex2,3 preproghrelin isoform, however, was restricted to white blood cells (WBC; where the wild-type preproghrelin isoform is not co-expressed), and gastrointestinal tissues. Expression of Δex2 preproghrelin and Δex2,3 preproghrelin mRNA was elevated in white blood cells in response to parasitic worm (helminth) infection in genetically susceptible sheep, but not in resistant sheep. Conclusions The restricted expression of the novel preproghrelin variants and their distinct WBC expression pattern during parasite infection may indicate a novel link between the ghrelin axis and metabolic and immune function in ruminants.

Analysis, characterisation and expression of gill-expressed carbonic anhydrase genes in the freshwater crayfish Cherax quadricarinatus

Changes in water quality parameters such as pH and salinity can have a significant effect on productivity of aquaculture species. Similarly, relative osmotic pressure influences various physiological processes and regulates expression of a number of osmoregulatory genes. Among those, carbonic anhydrase (CA) plays a key role in systemic acid–base balance and ion regulation. Redclaw crayfish (Cherax quadricarinatus) are unique in their ability to thrive in environments with naturally varied pH levels, suggesting unique adaptation to pH stress. To date, however, no studies have focused on identification and characterisation of CA or other osmoregulatory genes in C. quadricarinatus. Here, we analysed the redclaw gill transcriptome and characterized CA genes along with a number of other key osmoregulatory genes that were identified in the transcriptome. We also examined patterns of gene expression of these CA genes when exposed to three pH treatments. In total, 72,382,710 paired end Illumina reads were assembled into 36,128 contigs with an average length of 800 bp. Approximately 37% of contigs received significant BLAST hits and 22% were assigned gene ontology terms. Three full length CA isoforms; cytoplasmic CA (ChqCAc), glycosyl-phosphatidylinositol-linked CA (ChqCAg), and β-CA (ChqCA-beta) as well as two partial CA gene sequences were identified. Both partial CA genes showed high similarity to ChqCAg and appeared to be duplicated from the ChqCAg. Full length coding sequences of Na+/K+-ATPase, V-type H+-ATPase, sarcoplasmic Ca+-ATPase, arginine kinase, calreticulin and Cl− channel protein 2 were also identified. Only the ChqCAc gene showed significant differences in expression across the three pH treatments. These data provide valuable information on the gill expressed CA genes and their expression patterns in freshwater crayfish. Overall our data suggest an important role for the ChqCAc gene in response to changes in pH and in systemic acid–base balance in freshwater crayfish.

Total transcriptome, proteome, and allergome of Johnson grass pollen, which is important for allergic rhinitis in subtropical regions

Background Genomic data are lacking for many allergen sources. To circumvent this limitation, we implemented a strategy to reveal the repertoire of pollen allergens of a grass with clinical importance in subtropical regions, where an increasing proportion of the world's population resides. Objective We sought to identify and immunologically characterize the allergenic components of the Panicoideae Johnson grass pollen (JGP; Sorghum halepense). Methods The total pollen transcriptome, proteome, and allergome of JGP were documented. Serum IgE reactivities with pollen and purified allergens were assessed in 64 patients with grass pollen allergy from a subtropical region. Results Purified Sor h 1 and Sor h 13 were identified as clinically important allergen components of JGP with serum IgE reactivity in 49 (76%) and 28 (43.8%), respectively, of patients with grass pollen allergy. Within whole JGP, multiple cDNA transcripts and peptide spectra belonging to grass pollen allergen families 1, 2, 4, 7, 11, 12, 13, and 25 were identified. Pollen allergens restricted to subtropical grasses (groups 22-24) were also present within the JGP transcriptome and proteome. Mass spectrometry confirmed the IgE-reactive components of JGP included isoforms of Sor h 1, Sor h 2, Sor h 13, and Sor h 23. Conclusion Our integrated molecular approach revealed qualitative differences between the allergenic components of JGP and temperate grass pollens. Knowledge of these newly identified allergens has the potential to improve specific diagnosis and allergen immunotherapy treatment for patients with grass pollen allergy in subtropical regions and reduce the burden of allergic respiratory disease globally.

Purification of the major group 1 allergen from bahia grass pollen, Pas n 1

Background Group 1 grass pollen allergens are glycoproteins of the β-expansin family. They are a predominant component of pollen and are potent allergens with a high frequency of serum IgE reactivity in grass pollen-allergic patients. Bahia grass is distinct from temperate grasses and has a prolonged pollination period and wide distribution in warmer climates. Here we describe the purification of the group 1 pollen allergen, Pas n 1, from Bahia grass (Paspalum notatum), an important subtropical aeroallergen source. Methods Pas n 1 was purified from an aqueous Bahia grass pollen extract by ammonium sulphate precipitation, hydrophobic interaction and size exclusion chromatography, and assessed by one- and two-dimensional gel electrophoresis, immunoblotting and ELISA. Results Pas n 1 was purified to a single 29-kDa protein band containing two dominant isoforms detected by an allergen-specific monoclonal antibody and serum IgE of a Bahia grass pollen-allergic donor. The frequency of serum IgE reactivity with purified Pas n 1 in 51 Bahia grass pollen-allergic patients was 90.6%. Serum IgE reactivity with purified Pas n 1 was highly correlated with serum IgE reactivity with Bahia grass pollen extract and recombinant Pas n 1 (r = 0.821 and 0.913, respectively). Conclusions Pas n 1 is a major allergen reactive at high frequency with serum IgE of Bahia grass pollen-allergic patients. Purified natural Pas n 1 has utility for improved specific diagnosis and immunotherapy for Bahia grass pollen allergy.

The dominant 55kDa allergen of the subtropical Bahia grass (Paspalum notatum) pollen is a group 13 pollen allergen, Pas n 13

Bahia grass, Paspalum notatum, is an important pollen allergen source with a long season of pollination and wide distribution in subtropical and temperate regions. We aimed to characterize the 55. kDa allergen of Bahia grass pollen (BaGP) and ascertain its clinical importance. BaGP extract was separated by 2D-PAGE and immunoblotted with serum IgE of a grass pollen-allergic patient. The amino-terminal protein sequence of the predominant allergen isoform at 55. kDa had similarity with the group 13 allergens of Timothy grass and maize pollen, Phl p 13 and Zea m 13. Four sequences obtained by rapid amplification of the allergen cDNA ends represented multiple isoforms of Pas n 13. The predicted full length cDNA for Pas n 13 encoded a 423 amino acid glycoprotein including a signal peptide of 28 residues and with a predicted pI of 7.0. Tandem mass spectrometry of tryptic peptides of 2D gel spots identified peptides specific to the deduced amino acid sequence for each of the four Pas n 13 cDNA, representing 47% of the predicted mature protein sequence of Pas n 13. There was 80.6% and 72.6% amino acid identity with Zea m 13 and Phl p 13, respectively. Reactivity with a Phl p 13-specific monoclonal antibody AF6 supported designation of this allergen as Pas n 13. The allergen was purified from BaGP extract by ammonium sulphate precipitation, hydrophobic interaction and size exclusion chromatography. Purified Pas n 13 reacted with serum IgE of 34 of 71 (48%) grass pollen-allergic patients and specifically inhibited IgE reactivity with the 55. kDa band of BaGP for two grass pollen-allergic donors. Four isoforms of Pas n 13 from pI 6.3-7.8 had IgE-reactivity with grass pollen allergic sera. The allergenic activity of purified Pas n 13 was demonstrated by activation of basophils from whole blood of three grass pollen-allergic donors tested but not control donors. Pas n 13 is thus a clinically relevant pollen allergen of the subtropical Bahia grass likely to be important in eliciting seasonal allergic rhinitis and asthma in grass pollen-allergic patients.

Mosaic PPM1D mutations are associated with predisposition to breast and ovarian cancer

Improved sequencing technologies offer unprecedented opportunities for investigating the role of rare genetic variation in common disease. However, there are considerable challenges with respect to study design, data analysis and replication. Using pooled next-generation sequencing of 507 genes implicated in the repair of DNA in 1,150 samples, an analytical strategy focused on protein-truncating variants (PTVs) and a large-scale sequencing case-control replication experiment in 13,642 individuals, here we show that rare PTVs in the p53-inducible protein phosphatase PPM1D are associated with predisposition to breast cancer and ovarian cancer. PPM1D PTV mutations were present in 25 out of 7,781 cases versus 1 out of 5,861 controls (P = 1.12 × 10-5), including 18 mutations in 6,912 individuals with breast cancer (P = 2.42 × 10-4) and 12 mutations in 1,121 individuals with ovarian cancer (P = 3.10 × 10-9). Notably, all of the identified PPM1D PTVs were mosaic in lymphocyte DNA and clustered within a 370-base-pair region in the final exon of the gene, carboxy-terminal to the phosphatase catalytic domain. Functional studies demonstrate that the mutations result in enhanced suppression of p53 in response to ionizing radiation exposure, suggesting that the mutant alleles encode hyperactive PPM1D isoforms. Thus, although the mutations cause premature protein truncation, they do not result in the simple loss-of-function effect typically associated with this class of variant, but instead probably have a gain-of-function effect. Our results have implications for the detection and management of breast and ovarian cancer risk. More generally, these data provide new insights into the role of rare and of mosaic genetic variants in common conditions, and the use of sequencing in their identification.

Genetic control of bone density and turnover: Role of the collagen 1α1, estrogen receptor, and vitamin D receptor genes

Genetic factors are known to influence both the peak bone mass and probably the rate of change in bone density. A range of regulatory and structural genes has been proposed to be involved including collagen 1α1 (COL1A1), the estrogen receptor (ER), and the vitamin D receptor (VDR), but the actual genes involved are uncertain. We therefore studied the role of the COL1A1 and VDR loci in control of bone density by linkage in 45 dizygotic twin pairs and 29 nuclear families comprising 120 individuals. The influences on bone density of polymorphisms of COL1A1, VDR, and ER were studied by association both cross-sectionally and longitudinally in 193 elderly postmenopausal women (average age, 69 years) over a mean follow-up time of 6.3 years. Weak linkage of the COL1A1 locus with bone density was observed in both twins and families (p = 0.02 in both data sets), confirming previous observations of linkage of this locus with bone density. Association between the MscI polymorphism of COL1A1 and rate of lumbar spine bone loss was observed with significant gene-environment interaction related to dietary calcium intake (p = 0.0006). In the lowest tertile of dietary calcium intake, carriers of "s" alleles lost more bone than "SS" homozygotes (p = 0.01), whereas the opposite was observed in the highest dietary calcium intake (p = 0.003). Association also was observed between rate of bone loss at both the femoral neck and the lumbar spine and the TaqI VDR polymorphism (p = 0.03). This association was strongest in those in the lowest tertile of calcium intake, also suggesting the presence of gene-environment interaction involving dietary calcium and VDR, influencing bone turnover. No significant association was observed between the PvuII ER polymorphism alone or in combination with VDR or COL1A1 genotypes, with either bone density or its rate of change. These data support the involvement of COL1A1 in determination of bone density and the interaction of both COL1A1 and VDR with calcium intake in regulation of change of bone density over time.