|Home | About | Journals | Submit | Contact Us | Français|
Low dietary calcium (Ca) intake during growth limits peak bone mass but physiological adaptation can prevent this adverse effect. To assess the genetic control on the physiologic response to dietary Ca restriction (RCR) we conducted a study in 51 BXD lines fed either 0.5% (basal) or 0.25% (low) Ca diets from 4–12 wks of age (n=8/line/diet). Ca absorption (CaAbs), femur bone mineral density (BMD), and bone mineral content (BMC) were examined. ANCOVA with body size as covariate was used to detect significant line and diet main effects, and line-by-diet interactions. Body size-corrected residuals were used for linkage mapping and to estimate heritability (h2). Loci controlling the phenotypes were identified using composite interval mapping on each diet and for the RCR. h2 of basal phenotypes (0.37– 0.43) and their RCR (0.32–0.38) was moderate. For each phenotype we identified multiple QTL on each diet and for the RCR. Several loci affected multiple traits: Chr 1 (88.3–90.6 cM, CaAbs, BMC), Chr 4 (45.8–49.2 cM, CaAbs, BMD, BMC), Chr 8 (28.6–31.6 cM, CaAbs, BMD RCR), and Chr 15 (13.6–24 cM, BMD, BMC), and (32.3–36 cM, CaAbs RCR, BMD). This suggests that gene clusters may regulate interdependent bone-related phenotypes. Using in silico expression QTL (eQTL) mapping and bioinformatic tools we identified novel candidates for the regulation of bone under Ca stress (Ext1, Deptor), and for the first time, we report genes modulating Ca absorption (Inadl, Sc4mol, Sh3rf1 and Dennd3), and both Ca and bone metabolism (Tceanc2, Tll1 and Aadat). Our data reveal gene-by-diet interactions and the existence of novel relationships between bone and Ca metabolism during growth. This article is protected by copyright. All rights reserved
Attaining peak bone mass (PBM) during growth is a critical determinant of osteoporotic fracture in humans.(1) While 60% or more of the total variance in PBM is influenced by genetics,(2) there is growing evidence that environmental factors interact with genetics to define the extent of bone accrual during growth.(3,4) Dietary calcium (Ca) is an essential environmental factor required for skeletal mineralization. Inadequate Ca intake during growth compromises attainment of PBM and increases the risk of fracture later in life.(1) Fortunately, physiological mechanisms are in place to increase intestinal Ca absorption efficiency and renal Ca reabsorption to protect bone in response to dietary Ca stress.(5)
Several investigators have reported that the impact of dietary Ca intake on bone and Ca homeostasis is different across racial groups.(6–8) Similarly, our group has reported that bone parameters, Ca absorption efficiency, and their response to dietary Ca stress were significantly affected by genetic background in a diverse population of 11 inbred mouse lines.(9) Overall, these studies demonstrate that gene-by-environment (GxE) interactions can affect bone and Ca metabolism. Although the impact of genetic variants on bone phenotypes has been examined in many genome-wide association (GWA) and quantitative trait loci (QTL) mapping studies, very few studies have focused on dissecting the GxE interactions influencing bone and none of them have studied Ca absorption. Moreover, despite the importance of dietary Ca intake for attaining and maintaining bone mass, no studies have been conducted to identify the genetic loci regulating bone in different dietary Ca environments. Hence, our understanding of the gene-by diet (GxD) interactions affecting the attainment of PBM is limited.
The purpose of this study was to identify the genetic loci that regulate Ca and bone metabolism on a controlled, optimal Ca diet and to determine how genetics influences the response of mice to dietary Ca stress. For this study we used 51 BXD recombinant inbred (RI) lines(10) and well-defined diets. This is the largest genetic mapping study for bone phenotypes in one RI panel, it is the first mapping study to evaluate intestinal Ca absorption efficiency in any mouse model, and it is the first to evaluate GxD interactions on bone and Ca metabolism.
BXD recombinant inbred (RI) mouse lines are each defined by a fixed recombination pattern of alleles from the C57BL/6J (B6) and DBA/2J (DBA) inbred mouse lines.(11) Because each BXD line is inbred, individuals within a line are genetically identical and the panel can be used for experiments that include biological replicates and environmental interventions that test for the existence of gene-by-environment interactions.(11)
Male mice from 51 BXD RI lines were obtained at 4 wks of age from The Jackson Labs (Bar Harbor, ME) (the lines used are listed in Suppl. Table 1). Four to eight mice per line were obtained at a time and at arrival, an equal number of mice from each line were randomly assigned to either a 0.5% Ca (adequate) or 0.25% Ca (low) diet (AIN93G base with 200 IU vitamin D3/kg diet, Research Diets, New Brunswick, NJ) (n=8/diet/line). Dietary Ca levels were chosen to meet the rodent dietary Ca requirement (0.5% Ca) or to elicit a response to Ca restriction in target organs.(9) Food and water were provided ad libitum. Mice were group-housed (2–4 mice/cage) on TEK-fresh bedding (Harlan Laboratories) in conventional shoebox cages in the Purdue University animal facilities and maintained in an UVB light-free environment on a 12 h-light/dark cycle under standard conditions of temperature and humidity. At 12 wks of age mice were deprived of food 12 h prior to sacrifice. The day of the study, mice were anesthetized with an intraperitoneal injection of ketamine (22 mg/mL) and xylazine (33 mg/mL), and Ca absorption was measured by 45Ca appearance in the serum 10 min after an oral gavage test.(9) Individuals conducting absorption tests and sample collection were blinded to genotype and dietary treatment. Femora were placed in 10% neutral buffered formalin for 48 h and stored in 70% ethanol afterwards. Following removal of muscle from the fixed femora, femur length was recorded using a digital caliper (Mitutoyo America Corporation, Aurora, IL) and BMC (g) and BMD (g/cm2) were determined using a PIXImus II densitometer (Lunar; GE-Healthcare, Madison, WI) at the Indiana University School of Medicine. Scans were conducted in air on 10 bones at a time placed on the company-supplied plastic imaging plate (0.25 × 0.25 mm focal spot size; energy=80 kV, current=400 μA, fixed threshold of 1320). Samples were randomized and positioned uniformly along a grid to minimize potential scan position effects. Regions of interest were manually set for each bone and analysis was done by one investigator (RMR). The coefficient of variation for measurements of femora is 3–5%.(12) All of the experiments were approved by the Purdue University Animal Care and Use Committee.
Statistical analyses were conducted using SAS Enterprise Guide 4.2 (SAS Institute Inc., Cary, NC). Phenotype values with a z-score in the extreme 2.5% of either end of a line/diet group distribution were removed as outliers. Adherence to a normal distribution was assessed by Anderson-Darling tests. For non-normally distributed data, the Box-Cox transformation was performed. Ca absorption values required a natural log transformation to normalize the data. A new phenotype, the response to dietary Ca restriction (RCR), was calculated for each endpoint as we have previously reported.(9) ANCOVA was used to test for significant main effects of genetic background (i.e. line), diet, and line-by-diet interactions; femur length (FL) and body weight (BW) were used as covariates in the model.(13) ANCOVA was also used to determine the effect of genetic background on the RCR parameter. Significant covariate effects were observed on both diets and the RCR for BMD and BMC (BW, FL) and on both diets for Ca absorption (FL). Data were expressed as ANCOVA adjusted least square means (LSmean) ± SEM. Differences were considered significant when p<0.05. The Tukey-Kramer post-hoc test was used to determine differences between LSmeans.
Body size corrected residuals were used for linkage mapping and to calculate heritability estimates. For each phenotype the covariate effect of BW and FL was determined by Pearson’s correlation and removed by linear regression.(13) BMD and BMC were affected by both BW and FL whereas Ca absorption was only affected by FL. Narrow-sense heritability (h2) of bone and Ca absorption phenotypes was calculated for each diet as well as for the RCR parameter using the r2 of a one-way ANOVA (main effect=line). Line means (n=51) for each phenotype on the 0.5 % and 0.25% Ca diets and for their RCR were used for genetic mapping. Using 51 BXD RI strains gave us the power to reliably detect quantitative trait loci (QTL) accounting for 15% of the variance observed in the population with a power of ~ 0.80 and p<0.05.(14)
Marker information and BXD genotypes were downloaded from The GeneNetwork (www.genenetwork.org/genotypes/BXD.geno) and the sex-averaged genetic location of each marker was updated using the New Standard Genetic Map for the Laboratory Mouse(15) (http://cgd.jax.org/mousemapconverter). Markers with duplicated genetic locations or perfectly correlated genotypes in our panel of BXD lines were removed. The final genetic map contained 1558 markers. Composite interval mapping (CIM) was conducted on line means of raw (including BW and FL) and co-variate corrected phenotype data using Windows QTL Cartographer v2.5_011 (http://statgen.ncsu.edu/qtlcart/WQTLCart.htm). Forward selection identified 5 significant background markers. CIM was carried out using a Haldane map function, 2 cM walking speed, and a 10 cM window. Each diet (0.5% or 0.25% Ca) group and the RCR were mapped separately. Permutations (n=500) were used to determine significance for each analysis. A peak with a logarithm (base 10) of odds (LOD) ≥ 2 was considered putative. QTL plots for raw data are provided in Supplemental Fig. 1.
The QTL candidate region was defined using 1-LOD drop support intervals which approximate 95% confidence intervals.(16) Genetic locations of peaks (cM) were converted to base pairs positions (GRCm38) using the Mouse Map Converter tool. If a QTL was found within ≤ 5cM distance of another, the two loci were reported together and the location expressed as a range. QTL candidate regions were populated with genome features including protein-coding genes, non-coding RNA genes, gene fragments, and unclassified genes from the Mouse Genome Informatics (MGI) database (informatics.jax.org/).(17) Genes were identified as the region from the first to last exon plus 5000 bp upstream of the gene coordinates for exon 1. Regions that were identical by descent (IBD) between B6 and DBA were identified using the Mouse Phylogeny Viewer (http://msub.csbio.unc.edu/)(18) and 100% IBD regions were removed. Genes that were less than 100% IBD and the remaining non-IBD regions were used to query for polymorphisms between B6 and DBA using the SNP/variation query tool at the Mouse Phenome Database (MPD).(19) The Sanger1 and Sanger2 mouse data sets were merged to retrieve variation between the B6 and DBA lines. MPD annotations (dbSNP 138) were used to categorize polymorphisms by gene attribute: intronic, mRNA un-translated region (5′ and 3′ UTR), promoter region (5000 bp upstream from Exon 1), and exon-associated (i.e. synonymous and non-synonymous codons, stop and gain codons, splice sites, or frameshift mutations). Effects of non-synonymous amino acid changes and insertion/deletions were examined for potential functional effects using PROVEAN v1.1 (cutoff = −2.5).(20) Genome-wide DNase1 hypersensitive site (DNase1 HSS) data from the mouse ENCODE project was obtained from the UCSC genome browser and used to identify potential regulatory regions within the QTL regions; DNase1 HSS peak information from 10 mouse adult tissues was merged (fat pad, genital fat pad, heart, kidney, large intestine, liver, lung, skeletal muscle, spleen, and brain) and then overlaid with the polymorphisms in the QTL region.
Local (cis) and distal (trans) expression QTL (eQTL) analysis was conducted in silico using microarray databases from the BXD panel at the GeneNetwork (www.genenetwork.org/webqtl/main.py). The GeneNetwork does not contain small intestine microarray data so only bone (GN accession: GN414) and kidney mRNA (GN Accession: GN240) data were used for cis-eQTL detection and only the bone database was used for trans-eQTL. Significance was defined as a Likelihood Ratio Statistic (LRS) score ≥ 15 for cis-eQTL and LRS ≥ 20 for trans-eQTL. The NCBI37/mm9 location of the QTL (Mb) plus an inclusion buffer of 10 Mb was used for the search. Significant associations between transcript abundance of eQTL and our phenotype(s) of interest were determined by Pearson’s Correlation coefficient.
For candidate genes within QTL regions, we searched for functional annotation using gene ontology (GO) terms at MGI. Information regarding functional association of the genes with skeletal and/or Ca metabolism, and their role in biological pathways was obtained from the Coremine Medical Explorer (www.coremine.com) and the Ingenuity Target Explorer (https://targetexplorer.ingenuity.com). Finally, information of tissue gene expression patterns was obtained using BioGPS (http://biogps.org/).
All mice were healthy throughout the study. Body weight (BW) was significantly different among the lines (p<0.0001) but it was not affected by dietary Ca intake. BW variation was observed in the BXD population (mean 25.6 ± 0.7 g, range: 19.7 to 31.5 g) so phenotype data was adjusted for the effect of body size (see methods). Variability in Ca absorption, BMD and BMC was observed among the 51 RI lines for each diet group and for their response to dietary Ca restriction (RCR) (Z scores shown in Suppl. Fig. 2 and Suppl. Fig. 3; LSMeans in Suppl. Table 2). Ca absorption, and femoral BMD and BMC were each significantly affected by line (p<0.0001) and diet (p<0.0001) but only Ca absorption was affected by a significant line-by-diet interaction (p=0.0006). The average effect of dietary Ca restriction in the BXD panel was a 25.5% increase in Ca absorption efficiency (p<0.0001), a 3.5% decline in BMD (p<0.0001), and a 3.9% reduction in BMC (p<0.0001).
The narrow sense heritability (h2) for BMD and BMC was 43% and 40% in mice on the basal diet and 51% and 48% in mice on the low Ca diet groups, respectively. Heritability of Ca absorption efficiency was 37% in the basal Ca diet group and 40% in the low Ca diet group. The responses to dietary Ca restriction (RCR) for BMD, BMC and Ca absorption were all significantly affected by genetic background (p<0.0001) and had a heritability of 32%, 38% and 36%, respectively. The RCR for each phenotype ranged from −37.3% to 108.9% for Ca absorption, −15 to 11.4% for BMD, and −27.3 to 23.6% for BMC across the BXD panel.
We evaluated associations between bone parameters and Ca absorption efficiency (Suppl. Table 3). BMD was highly correlated to BMC on either diet, and their responses to Ca restriction were also positively and significantly correlated with each other (r > 0.9, p<0.0001). The RCR for Ca absorption was positively correlated to Ca absorption on the low Ca diet (r= 0.49, p<0.001) but negatively correlated to Ca absorption on the basal Ca diet (r = −0.38, p=0.005). This indicates that lines with high basal Ca absorption were less likely to increase Ca absorption efficiency when fed a low Ca diet. The RCR for BMD or BMC were not correlated to their phenotypes on the basal or on the low Ca diets (r = 0.22, p=0.1 and r = 0.32, p = 0.22, respectively). However, the responses of BMD and BMC to dietary Ca restriction were significantly correlated to Ca absorption efficiency on both the low (r=0.52, p < 0.0001 and r=0.36, p=0.010) and basal Ca diets (r=0.40, p = 0.003 and r=0.37, p=0.007) (Fig. 1).
We identified multiple significant QTL for each phenotype in each diet group and for the RCR (co-variate adjusted, Fig. 2; unadjusted raw, Supplmental Fig. 1). When QTL mapping results across analyses were considered together (3 phenotypes, 2 diets, and the RCR), 15 significant QTLs were identified (Table 1, Suppl. Fig. 4). We identified significant QTL controlling the BMD RCR on chromosome (Chr) 8 (QTL8b, 26.9–32.7 cM); the BMC RCR in Chr 1 (QTL1a 66.8–69.03 cM), Chr 4 (QTL4, 44.8–50.7 cM) and Chr 9 (QTL9b, 57–59.74 cM); and for the Ca absorption RCR in Chr 3 (QTL3, 27.9–31.7 cM), Chr 8 (QTL8a, 22.7–26.9 cM), and Chr 15 (QTL15b, 30.87–37.6 cM). Ten of the 15 QTL we identified were significant for more than one phenotype (QTL1b, 3, 4, 7, 8b, 9b, 10, 12, 15a, 15b). Within each phenotype, none of the QTL for the RCR overlapped with those identified at baseline. Loci that modulated traits in both diet groups were observed for BMC on Chr1 and Chr15 and for Ca absorption on Chr12 (Table 1). A list of genome features after IBD filtering as well as the number of polymorphisms associated with non-IBD genes for each these loci are given in Suppl. Table 4.
7 QTL regions were selected for characterization based on: the amount of variation explained by a single QTL (r2), mapping of multiple phenotypes to a loci, and identification of a loci controlling a phenotype in two or more conditions (Table 1). The detailed classification of functional polymorphisms and insertion/deletions (Indels) for each of these loci is in Suppl. Table 5. Genes with stop codon, frameshift, splice site mutations, deleterious non-synonymous polymorphisms, or indels were automatically considered as candidate genes. A frameshift mutation was detected for QTL4 (Chr 4, rs263625036, 3110021N24Rik) and 2 splice site donor (SSD) polymorphisms were identified for QTL8b (Chr8, rs33461690, AW046200, G/C and rs31742488, Tll1, G/A). Potentially deleterious non-synonymous polymorphisms were identified in QTL1b (Chr1, rs32873773, Lyplal1, I58T); QTL4 (Chr 4, rs31915781, Gm12689, S50F, and rs31774248, Gm12689, E 77A); QTL8b (Chr 8, rs33014006, Adam29, S163I, and rs238747469, Gm4975, R18L); and QTL15a (Chr15, rs32502839, Pkhd1l1, N3877T, and rs47743747, Abra, L374V).
Cis-eQTL overlapping with our QTL regions whose mRNA levels were significantly correlated with the phenotype and whose gene was not identical by descent between B6 and DBA were considered candidate genes.(21) Cis-eQTL meeting our criterion were detected in: QTL1b (Illumina probe: scl000933.1_43 lacking annotation), QTL4 (Tceanc2, Acot11 and Inadl), QTL8b (Apela, Aadat, Sh3rf1 and Sc4mol), QTL15a (Ext1, Ennp2, Deptor), and QTL15b (Dennd3) (Table 2).
The candidate gene list for each QTL included polymorphisms affecting protein coding and cis-eQTL genes. Bioinformatic analysis identified genes whose polymorphisms might affect our phenotypes based on their reported function, their tissue gene expression pattern, and their reported interaction with other genes or proteins. No deleterious mutations or significant cis-eQTL were detected for QTL2 or QTL3. Based on our analysis we identified candidate genes for loci regulating basal Ca absorption (Inadl, Sc4mol, Sh3rf1), the Ca absorption RCR (Dennd3); bone mass on one of the two diets (Acot11, Ext1, Enpp2 and Deptor); and both Ca absorption and basal bone mass (Lyplal1) or bone RCR (Tceanc2, Aadat and Tll1) (Suppl. Table 6).
Previous studies have shown that the change in Ca absorption efficiency in response to changes in dietary Ca intake has a strong genetic component.(6,8,9) Consistent with these findings, we reported that basal Ca absorption efficiency and its response to dietary Ca restriction (RCR) are heterogeneous in the BXD RI panel and that they have heritability estimates (37 % and 36% respectively) that are similar to those reported for other Ca and bone metabolism phenotypes in mice and humans.(22,23) We also found that the genetic regulation of the adaptive response (RCR) is independent from the regulation of Ca absorption efficiency or bone mass at baseline indicating that the loci for the RCR contain genes sensitive to the dietary Ca environment.
Several genes have been proposed to contribute to basal and vitamin D-regulated intestinal Ca absorption(24) including candidates from a recent genome-wide study that identified forty-five new 1,25(OH)2D-target genes in the intestine.(25) However, these genes were either located in regions that are identical by descent between B6 and DBA, located outside the QTL we identified, or contained polymorphisms whose variants have no consequence on protein function or mRNA level. Thus, our QTL mapping study has revealed novel variation controlling intestinal Ca absorption and it has the potential to reveal new genes that control Ca absorption efficiency and its response to dietary Ca changes.
Within our list of prioritized loci, several influenced either basal Ca absorption efficiency (QTL4, 8b, 1b) or its RCR (QTL3,15b). Our bioinformatics analysis identified several candidate genes within each loci (Suppl. Table 6). However, the functional link between our candidate genes and Ca metabolism is not straight forward because there is limited information on the role of the genes (Apela) or because their function has not been linked to Ca homeostasis: Lyplal1 (protein depalmitoylation), Sc4mol (lipid/cholesterol metabolism), Sh3rf1 (Rac/JNK signaling), Dennd3 (regulation of GTPase activity). Nevertheless, further studies on these genes could reveal novel roles for them in the regulation of Ca metabolism.
One candidate gene of particular interest is the Inadl gene in QTL4. Inadl helps organize multimeric protein complexes in the plasma membrane(26) and contributes to the formation of tight junctions (TJ).(27) Inadl could influence Ca absorption by affecting the formation of Ca transporter complexes including either TRPV6 (at the apical membrane) or PMCA1b (at the basolateral membrane). It could also affect the TJ proteins that mediate the passive diffusion of small ions across the epithelia (e.g. claudins 2 and 12(28)). Since Inadl expression level is positively correlated to Ca absorption efficiency, we hypothesize that elevated levels of Inadl could enhance the proper formation of transmembrane Ca transport complexes and ultimately increase intestinal Ca absorption.
Many studies have examined the independent effects of genetics and diet on bone but very few of these studies are focused on GXD interactions. Our work is the first carefully controlled study to identify loci underlying the bone responses to dietary Ca stress. We found that genetic factors contribute up to 43% of the variation observed in BMD and BMC at baseline and up to 38% of the variation in RCR. Several of our bone loci were located within the 95% confidence intervals of loci from other studies: i.e. QTL1b (87.6 cM, whole body aBMD), QTL10 (62.5 cM, femur aBMD), and QTL15b (33.8 cM, femur aBMD; and 38.7 cM, whole-body aBMD).(29) Within our loci, we found a number of genes known to regulate bone metabolism in humans: osteonectin (Spock 3/SPARC)(30) within QTL8b, osteoprotegerin (TNFRSF11B)(31) in QTL15a, and the leptin receptor (Lepr) gene (32) in QTL4. However, the Lepr gene is 100% IBD between B6 and DBA mice and neither of the other two genes contain deleterious polymorphisms affecting protein function or gene expression. Hence, it is unlikely that these three genes are contributors to the phenotype variation observed in BMD and BMC in the BXD population.
Wang et al.(33) recently used a panel of male mice from 46 BXD RI lines and identified a significant loci controlling femur and tibia BMD in the same genomic region where our QTL15a is located. They identified a non-synonymous coding polymorphism (V813L) in the Trichorhino-phalangeal Syndrome Type I Protein (TRPS1) gene (rs32398060, B6=C, DBA=G) as the variant controlling the phenotype. However, we found that the polymorphism was classified as neutral using the PROVEAN software (score = −0.01), suggesting that the sequence variation does not affect protein function. In addition, while TRPS1 mRNA is a cis-eQTL at the QTL15a locus, TRPS1 mRNA levels were not correlated with either BMD or BMC in our BXD population. Therefore, we do not consider TRPS1 to be a strong candidate gene for QTL15a.
The most promising candidate gene in QTL15a is exostoses multiple 1 (Ext1). Ext1 has glycosyltransferase activity required for heparan sulfate (HS) biosynthesis,(34) a component of the extracellular matrix (ECM) critical for endochondral development.(35) Ext1 expression is very high in osteoblasts and osteoclasts, and mice with a hypomorphic mutation in Ext1 have impaired chondrocyte and bone differentiation, and defects in HS production.(35) In humans, a splice variant of Ext1 has been reported as a causative locus for hereditary multiple exostosis (HME), a disease characterized by a low PBM and osteoporosis,(36) and a polymorphism near the Ext1 gene was recently identified as responsible for height variation in a population of Korean adults.(37) In our study, bone Ext1 mRNA level was positively correlated with BMD and BMC in mice on the low Ca diet. We hypothesize that reduced Ext1 transcript levels will reduce production of HS, reduce signals for endochondral ossification, and lower BMD and BMC, especially under dietary Ca stress.
In addition to the loci regulating basal BMD and BMC, we identified loci for their RCR in QTL8b, QTL1a, QTL4, and QTL9b. These loci did not co-localize with loci controlling basal bone mass. This indicates that there are independent genetic controls on basal bone mass and its RCR. As a result, future bone genetics studies in humans should control for habitual dietary Ca intake. QTL4 and QTL8b also controlled basal Ca absorption efficiency. Others have shown that higher Ca absorption efficiency during growth is associated with higher bone deposition in children(38) and mice(9) while reduced Ca absorption efficiency in adulthood is associated with increased hip fracture risk.(39,40) Consistent with a functional relationship between bone and Ca absorption efficiency, we found that the RCR for BMC and BMD were significantly, positively correlated to intestinal Ca absorption efficiency (Fig. 1, Suppl. Table 3). These findings suggest that genetic factors leading to high Ca absorption efficiency protect bone under dietary Ca stress. In addition, the co-localization of these complementary phenotypes within a locus suggests the existence of novel relationships between bone and Ca metabolism during growth.
QTL4 is a particularly interesting locus controlling BMC and BMD on the basal Ca diet, the BMC RCR, and basal Ca absorption efficiency (Table 1). Others have mapped a number of other bone-related traits to this locus including femur mechanical properties,(41) spine and femur BMD, and spine and femur μCT phenotypes.(42,43) Several potential candidate genes are in this region but there is not an obvious link between their biological roles and Ca/bone metabolism. However, one gene, Tceanc2, encodes a protein whose GO annotation identifies it as involved in DNA binding, regulation of mRNA processing, and transcription elongation by RNA polymerase II. Tceanc2 mRNA expression was positively correlated with basal BMD and negatively correlated with basal Ca absorption and the BMC RCR. Although, very little is known regarding the role of this gene in mice or humans, its fundamental biological function in regulatory events suggests it could impact many physiologic processes, including bone and Ca metabolism.
QTL8b controls both basal Ca absorption efficiency and the BMD RCR. Two interesting candidate genes, Tll1 and Aadat, are within this locus. Tll1 is highly expressed in osteoblasts and is a member of the mammalian bone morphogenetic protein-1 (BMP1)-like proteinases that regulate the formation, organization, and mineralization of extracellular matrix (ECM).(44) In BXD mice, a G/A polymorphism affects a splice donor site in Tll1. Since the polymorphism is associated with a positive impact on the BMD RCR and basal Ca absorption, we hypothesize that alternative splicing results in a more active isoform of the Tll1 protein that increases the formation of ECM matrix, increases bone mass, and protects bone from dietary Ca restriction. However, while this hypothesis explains the role of Tll1 in the regulation of bone mass, the role of this gene in the modulation of Ca absorption is less clear and may be indirect.
Aadat mRNA levels were negatively associated with both basal Ca absorption and the RCR for BMD and this suggests Aadat is a promising co-regulator of Ca and bone metabolism. Aadat catalyzes the production of kynurenic acid (KYNA), a non-competitive glutamate receptor antagonist.(45) Glutamate-mediated excitation of neurotransmitters can stimulate thyroid function.(46) In rats, glutamate receptor agonists can increase TSH and thyroid hormone concentrations whereas glutamate receptor antagonists decrease them.(47) A recent human GWAS identified Aadat as a novel locus associated with the variation in free T4 levels.(48) We hypothesize that decreased Aadat levels increase thyroid hormone levels which others have shown can enhance intestinal Ca absorption,(49) chondrocytes maturation, osteoblast differentiation, and osteoclasts activation.(50)
Our study has several strengths and weaknesses that should be noted. One strength of our approach was that by using the BXD RILs, with numerous recombination events during inbreeding, we had higher resolution mapping than is possible with F2 populations. In addition, since the RI have genotypes that are fixed as homozygotes they can be used for the study of GXE interactions. As a result, our study is unique in that it has determined the genetic architecture of the response of bone and Ca metabolism to the stress environment of low dietary Ca intake. Some limitations in our study also exist. First, the BXD panel was built with two founders so captures just a portion of the genetic diversity available in the mouse genome. Second, although we used more lines than any other bone RI mapping study, the relatively small size of our population (51 lines) limits our power to detect small effect size QTLs. Third, we studied only male mice, and so we are likely missing important loci that are influenced by the female sex hormones. Finally, the assignment of arbitrary gene boundaries for our QTL characterization may have missed other plausible candidate genes.
In summary, our study is novel in that it is the first to examine either the genetic architecture of intestinal Ca absorption efficiency or the response of bone mass and Ca absorption to dietary Ca restriction. Our data revealed strong genetic effects on basal Ca and bone metabolism and their responses to Ca stress. Importantly, the genetic variation controlling the response to low Ca intake is independent of the genetic variation controlling baseline levels of these phenotypes. This demonstrates that there are independent genetic effects that should be respected in the design of future bone genetic studies. Moreover, our study confirms the existence of GxD interactions influencing bone accrual during growth and shows that accrual of bone during growth is dependent on high Ca absorption efficiency. A critical finding from our study was the evidence for novel regulation of Ca and bone homeostasis under conditions of Ca stress, and the possibility of common loci controlling Ca and bone phenotypes during growth.
This work was supported by NIH grant ES019103 to JCF. CONACyT, Mexico provided a partial graduate scholarship to PRF.
†This article has been accepted for publication and undergone full peer review but has not been through the copyediting, typesetting, pagination and proofreading process, which may lead to differences between this version and the Version of Record. Please cite this article as doi: [10.1002/jbmr.2760]
JCF is on the Scientific Advisory Board for Innophos, Inc. All other authors have no conflicts of interest to report.
Authors’ roles: Study design: JCF. Study execution: RAR, LW and PRF. Data collection: RAR. Data analysis: PRF, LW, MZ. Data interpretation: JCF, PRF, RAR. Drafting manuscript: PRF and JCF. Approval of final manuscript: JCF, PRF, RAR, LW and MZ. JCF takes responsibility for the integrity of the data analysis.
Additional Supporting Information may be found in the online version of this article.