Population analysis of vitamin D receptor polymorphisms and the role of genetic ancestry in an admixed population.

The vitamin D receptor (VDR) is an essential protein related to bone metabolism. Some VDR alleles are differentially distributed among ethnic populations and display variable patterns of linkage disequilibrium (LD). In this study, 200 unrelated Brazilians were genotyped using 21 VDR single nucleotide polymorphisms (SNPs) and 28 ancestry informative markers. The patterns of LD and haplotype distribution were compared among Brazilian and the HapMap populations of African (YRI), European (CEU) and Asian (JPT+CHB) origins. Conditional regression and haplotype-specific analysis were performed using estimates of individual genetic ancestry in Brazilians as a quantitative trait. Similar patterns of LD were observed in the 5' and 3' gene regions. However, the frequency distribution of haplotype blocks varied among populations. Conditional regression analysis identified haplotypes associated with European and Amerindian ancestry, but not with the proportion of African ancestry. Individual ancestry estimates were associated with VDR haplotypes. These findings reinforce the need to correct for population stratification when performing genetic association studies in admixed populations.

Introduction

The vitamin D receptor (VDR) is a member of the superfamily of nuclear receptors for steroid hormones that functions as a ligand-activated transcription factor (Dusso ). The VDR associated with the secosteroid hormone 1,25-dihydroxyvitamin D3 (1,25(OH)2 Vitamin D3) and heterodimerized with the retinoid X receptor (RXR) binds to vitamin D3 response elements in the promoter region of responsive genes (Dusso ). The genes that are up- or down-regulated by the complex of vitamin D3, VDR, RXR and other recruited proteins are associated with calcium homeostasis, bone metabolism, cell cycle, immunomodulation and other hormonal systems (Dusso ; Lips, 2007). The broad range of vitamin D functions has focused attention on the VDR gene as an important candidate gene that could explain variations in specific phenotypes possibly connected with vitamin D metabolism (Valdivielso and Fernandez, 2006).

Much information has been generated since the first description of VDR polymorphism (Morrison ) and has led to intense investigation of the allelic variation in the VDR gene in different ethnic populations (Nejentsev ; Thakkinstian ; Fang ). The conflicting data on the association of VDR polymorphisms with specific phenotypes is sometimes confusing. The reasons commonly given to explain the difficulty in reproducing many results include uncontrolled environmental factors, population stratification, locus heterogeneity and different linkage disequilibrium (LD) patterns (Nejentsev ; Thakkinstian ; Fang ). The VDR region consists essentially of three haplotype blocks located in the intergenic region of the VDR and COL1A1 gene, the 5′ promoter region and the 3′ region encompassing the untranslated region, with the frequency distributions of LD and haplotypes varying among European, African and Asian populations (Nejentsev ; Fang ).

Specific variations in the allelic frequencies of VDR polymorphisms among Europeans, Africans, Amerindians and Asians could increase the risk of spurious associations in studies of recently admixed populations such as Brazilians (Rosenberg and Nordborg, 2006). The current Brazilian population is one of the most heterogeneous in the world, descending from an admixture of Europeans, Amerindians and Africans during the last five centuries. The use of ancestry informative markers (AIMs) has revealed ample genetic heterogeneity in the Brazilian population (Callegari-Jacques ; Parra ; Marrero ; Lins ) and this characteristic may be used to control population stratification in association studies (Suarez-Kurtz ).

Association studies in admixed populations relying only on self-reported ancestry or physical features to arrange volunteers in homogenous groups may produce spurious associations because of stratification generated by admixture (Cardon and Palmer, 2003; Ziv and Burchard, 2003; Suarez-Kurtz ). This is particularly important when DNA markers used to conduct association studies and the phenotype investigated display different frequency distributions among the reference groups that gave rise to the admixed population (Pritchard and Donnelly, 2001; Rosenberg and Nordborg, 2006). Bone phenotypes differ among Africans and Europeans (Gilsanz ; Jones ) and several polymorphisms, including those of the VDR gene, have been investigated as candidates to explain its quantitative variation. The purpose of this study was to perform a greater in-depth analysis of variability in the VDR gene in admixed Brazilians and correlate this variability with individual genetic ancestry estimates in order to identify possible pitfalls when performing association studies in an admixed population.

Material and Methods

Population sample

The Brazilian population sample (BRZ) consisted of 200 unrelated healthy subjects randomly chosen from individuals involved in no-cost paternity investigations from 2003 to 2005. All subjects signed an informed consent form that allowed the use of their DNA samples for paternity testing and further population genetics research. To avoid bias during analysis no attempt was made to classify the subjects according to morphological or social traits. The subjects were allocated to one of five groups (n = 40 each) based on their birthplace in one of the five geopolitical regions (Midwest, Northeast, North, Southeast and South) of Brazil. The research protocol was approved by the university Ethics Committee.

HapMap data and genotyping

The VDR genotypes of the HapMap population samples were retrieved from an online database (Data Rel 21a/phaseII Jan07, on NCBI B35 assembly, dbSNP b125). The total sample consisted of 89 unrelated East-Asian individuals (ASN) comprising 45 Han Chinese from Beijing (CHB) and 44 Japanese from Tokyo (JPT), 60 unrelated individuals from northern and western European origin (CEU) and 60 unrelated Yoruba individuals (YRI) from Ibadan, Nigeria.

The choice of VDR SNPs was based on markers of HapMap phase I and phase II data that were polymorphic in at least one population and dispersed with average intervening distances of 5 kb; Haploview software (Barrett ) was used to establish the LD patterns. Subsequently, a minimum set of SNPs representing the original LD blocks was selected with a 90% prediction coverage (Lins ). Sixteen SNPs were selected in addition to the five most studied VDR polymorphisms (rs11568820-Cdx2, rs10735810-FokI, rs1544410-BsmI, rs7975232-ApaI and rs731236-TaqI) used in our previous studies (Gentil , 2009; Lins , 2009; Moreno Lima ).

Estimates of genetic admixture in the Brazilian samples were calculated using a set of 28 autosomal ancestry informative markers (AIMs) selected from previous studies that reported large differences in allele frequency among European, African and Native American populations. Detailed procedures for calculating the ancestry estimates are described by Lins who used the same population as the present study.

PCR primers and single base extension primers were designed using Primer3 based on recommendations of the SNaPshot Multiplex Kit protocol (Applied Biosystems). The SNPs were assembled into three multiplex panels and then genotyped by a modified single base extension methodology described elsewhere (Lins ).

Statistical analysis

Estimates of allele frequency, deviations from Hardy-Weinberg equilibrium and pairwise genetic distance estimates based on Wright’s Fst statistics were calculated using Arlequin v. 3.01 (Excoffier ).

The linkage disequilibrium analyses were done by estimating the parameters D’ and r. The structure of haplotype blocks in each population was defined by solid spine of LD algorithm in Haploview version 3.32 (Barrett ). This criterion defines a block when the first and last markers are in strong LD with all intermediate markers, thereby providing more robust block boundaries. After block definition, the haplotypes of the HapMap population samples were estimated for the blocks established in BRZ in order to compare the haplotype distributions. Haplotype estimates were then calculated using Whap software (Purcell ) which is based on the expectation-maximization algorithm and uses the estimates of posterior probabilities to account for the ambiguity of haplotype phase in subsequent association tests. This package was developed to handle quantitative traits and covariates for regression analysis. In this case, the individual ancestry estimated for each ethnic group was set as a quantitative trait in BRZ. Initially, a likelihood ratio test (LRT) on the omnibus conditional regression test indicated whether there was a significant influence of haplotypes on the trait. Then, a haplotype-specific regression-based test comparing the haplotype effect against all other haplotypes indicated whether the effect observed on a specific haplotype was significant (Purcell ). For quantitative traits, the conditional analysis approach is more robust and was chosen to model genotype conditional on trait, instead of trait conditional on genotype (the usual approach in such analyses) (Purcell ). In addition, a standard approach run in Phase version 2.0.2 software (Stephens and Donnelly, 2003) was used to estimate recombination hotspots by comparing the median of recombination parameters among SNP pairs with the background rate assumed for the general human population (Crawford ).

Results

The allele frequency distribution for the VDR gene was similar in the five geopolitical regions (Table 1) and a pairwise Fst test identified no significant difference among them (all p-values were > 0.050). However, a significant difference was found when the Brazilian subgroups were combined to form one group (BRZ) and compared with the HapMap populations in a pairwise Fst test. In this case, the Brazilian population was genetically more distant from the HapMap African derived population (BRZ-YRI Fst = 0.154; p < 0.001) than from the HapMap population with European background (BRZ-CEU Fst = 0.012; p = 0.009). The analysis of individual loci showed that only four out of 19 loci were significantly different (p < 0.05) in the pair BRZ-CEU (Table 2). In contrast, in the other population pairs, only a few loci did not differ significantly in their allele frequencies (Table 2).

Table 1

SNPs genotyped in the VDR gene and their chromosomal position (chromosome 12 genomic contig, NT_029419.11) and allelic frequencies in the samples from the HapMap populations, the total Brazilian population and the five geopolitical regional samples.

SNP (locus)	Position	Allele	CEU	YRI	ASN	BRZ	N	NE	MW	SE	S
rs4077869	46591911	A	N.A.	0.400	0.978	0.667	0.538	0.650	0.650	0.731	0.769
rs11568820 (Cdx2)	46588812	C	0.792	0.017	0.506	0.684	0.663	0.705	0.637	0.688	0.731
rs4516035	46586093	A	0.575	0.992	0.989	0.788	0.778	0.786	0.789	0.684	0.958
rs10783219	46581755	A	0.331	0.000	0.416	0.306	0.244	0.325	0.333	0.282	0.346
rs7302235	46579105	A	0.741	0.425	N.A.	0.703	0.637	0.776	0.724	0.643	0.731
rs3890734	46575622	C	0.675	0.877	0.994	0.684	0.679	0.688	0.700	0.675	0.676
rs2853559	46569072	A	0.424	0.169	0.421	0.449	0.414	0.544	0.545	0.395	0.361
rs2853564	46564754	A	0.583	0.908	0.579	0.372	0.262	0.338	0.410	0.488	0.363
rs2254210	46559981	C	0.633	0.658	0.624	0.688	0.700	0.712	0.675	0.613	0.738
rs10735810 (FokI)	46559162	C	0.525	0.833	0.646	0.674	0.675	0.744	0.667	0.628	0.658
rs886441	46549231	A	0.808	0.583	0.978	0.767	0.800	0.718	0.868	0.705	0.744
rs2239179	46544033	C	0.417	0.292	0.247	0.443	0.462	0.408	0.500	0.463	0.371
rs2248098	46539623	A	0.425	0.381	0.663	0.465	0.387	0.551	0.425	0.450	0.512
rs2239185	46530826	A	N.A.	0.542	0.347	0.568	0.552	0.574	0.600	0.615	0.500
rs1544410 (BsmI)	46526102	C	0.525	0.712	0.921	0.605	0.637	0.590	0.618	0.526	0.654
rs7975232 (ApaI)	46525104	G	0.424	0.375	0.645	0.460	0.475	0.500	0.425	0.385	0.512
rs731236 (TaqI)	46525024	A	0.526	0.750	0.933	0.627	0.625	0.603	0.700	0.551	0.654
rs9729	46522890	G	0.414	0.337	0.657	0.390	0.375	0.395	0.412	0.295	0.474
rs7968585	46518360	A	0.592	0.633	0.326	0.525	0.500	0.483	0.433	0.645	0.528
rs11608702	46515035	A	0.358	0.183	0.618	0.347	0.275	0.350	0.395	0.295	0.423
rs2544040	46509213	C	1.000	0.885	1.000	0.985	1.000	1.000	1.000	1.000	0.925

HapMap populations: CEU = European ancestry, YRI = African ancestry, ASN = Asian ancestry; Total Brazilian population (BRZ) and regional Brazilian samples: N = North, NE = Northeast, MW = Midwest, SE = Southeast, S = South. N.A. = Data not available.

Table 2

Pairwise population differences (Fst) for locus-by-locus analysis for the HapMap samples and total Brazilian population.

refSNP (locus)	YRI-BRZ		ASN-BRZ		CEU-BRZ		CEU-YRI		CEU-ASN		YRI-ASN
	Fst	p	Fst	p	Fst	p	Fst	p	Fst	p	Fst	p
rs4077869	0.133	< 0.001	0.230	< 0.001	N.A.	N.A.	N.A.	N.A.	N.A.	N.A.	0.601	< 0.001
rs11568820 (Cdx2)	0.566	< 0.001	0.062	< 0.001	0.023	0.022	0.766	< 0.001	0.153	< 0.001	0.430	< 0.001
rs4516035	0.176	< 0.001	0.199	< 0.001	0.092	0.001	0.403	< 0.001	0.445	< 0.001	0.007	1.000
rs10783219	0.219	< 0.001	0.022	0.017	0.004	1.000	0.327	< 0.001	0.009	0.183	0.369	< 0.001
rs7302235	0.146	< 0.001	N.A.	N.A.	0.003	0.814	0.178	< 0.001	N.A.	N.A.	N.A.	N.A.
rs3890734	0.083	< 0.001	0.237	< 0.001	0.005	1.000	0.104	< 0.001	0.347	< 0.001	0.122	< 0.001
rs2853559	0.146	< 0.001	0.002	1.000	0.004	1.000	0.136	< 0.001	0.007	1.000	0.129	< 0.001
rs2853564	0.417	< 0.001	0.079	< 0.001	0.082	0.001	0.238	< 0.001	0.007	1.000	0.227	< 0.001
rs2254210	0.003	1.000	0.005	0.204	0.001	0.431	0.007	1.000	0.008	1.000	0.004	1.000
rs10735810 (FokI)	0.054	< 0.001	0.002	1.000	0.042	0.003	0.190	< 0.001	0.023	0.062	0.077	0.001
rs886441	0.074	< 0.001	0.143	< 0.001	0.001	0.562	0.105	0.001	0.154	< 0.001	0.407	< 0.001
rs2239179	0.041	0.002	0.074	< 0.001	0.004	1.000	0.024	0.074	0.057	0.007	0.003	0.701
rs2248098	0.008	0.147	0.071	< 0.001	0.003	0.838	0.006	1.000	0.101	< 0.001	0.141	< 0.001
rs2239185	0.004	1.000	0.089	< 0.001	N.A.	N.A.	N.A.	N.A.	N.A.	N.A.	0.068	0.004
rs1544410 (BsmI)	0.018	0.067	0.208	< 0.001	0.007	0.203	0.063	0.783	0.349	< 0.001	0.142	< 0.001
rs7975232 (ApaI)	0.009	0.134	0.062	< 0.001	0.003	1.000	0.004	0.001	0.088	< 0.001	0.131	< 0.001
rs731236 (TaqI)	0.027	0.021	0.202	< 0.001	0.015	0.093	0.095	0.309	0.373	< 0.001	0.123	< 0.001
rs9729	0.001	0.533	0.128	< 0.001	0.005	1.000	0.004	1.000	0.106	< 0.001	0.181	< 0.001
rs7968585	0.017	0.064	0.072	< 0.001	0.002	0.386	0.004	0.003	0.128	< 0.001	0.169	< 0.001
rs11608702	0.055	0.003	0.134	< 0.001	0.006	1.000	0.067	< 0.001	0.120	< 0.001	0.315	< 0.001
rs2544040	0.128	< 0.001	0.007	0.170	0.005	0.349	0.116	< 0.001	N.A.	N.A.	0.144	< 0.001

HapMap populations: CEU = European derived ancestry, YRI = African derived ancestry, ASN = Asian derived ancestry; BRZ = total Brazilian population. N.A. = Data not available (includes monomorphic loci and missing data). p < 0.05 indicates a significant difference.

Information on the VDR haplotype structure in the HapMap data showed that haplotype extension was greater in CEU, followed by ASN and YRI, which had more blocks of lower extension, compared to the others (Figure 1). The CEU and ASN populations had similar LD patterns, with two blocks in the 5′ region, one of them identical, and one block in the 3′ region. A difference was observed only in the length of the first 5′ and 3′ haplotype blocks.

Figure 1

Distribution of the haplotype blocks among populations (ASN – East-Asian individuals, CEU – northern and western European origin, YRI – Yoruba individuals from Ibadan, Nigeria) and estimates of the pairwise recombination rates in the Brazilian (BRZ) population. Horizontal bars represent the extent of the haplotype blocks in each population. Values on the y-axis are the median of the factor by which the recombination rate between loci (x-axis) exceeds the background recombination rate. The SNPs rs4077869, rs2239185, rs7302235 and rs4516053 are shown to facilitate location since they were excluded from population analysis.

For a comparative inter-population analysis of haplotype block diversity, four SNPs were excluded from the BRZ dataset either because they lacked genotypes in HapMap (e.g., rs4077869 and rs2239185 missing in the CEU population and rs7302235 missing in the ASN population) or deviated from Hardy-Weinberg expectations in BRZ (rs4516035 p = 0.001). Overall, two blocks were observed in the Brazilian population, one in the 5′ gene region and another at the end of the transcription region and the 3′ UTR region of the VDR gene (Figure 1). The 5′ haplotype block contained the Cdx2, rs10783219 and rs3890734 SNPs and extended 13 kb, with mean linkage disequilibrium measures of D’ = 0.924 and r2 = 0.175. The 3′ haplotype block consisted of rs2248098, BsmI, ApaI, TaqI, rs9729, rs7968585 SNPs and spanned 21 kb, with mean linkage disequilibrium measures of D’ = 0.849 and r2 = 0.475. The estimated haplotypes for each block identified five haplotypes in the 5′ region (Table 3) and 19 in the 3′ region of the gene (Table 4) when all populations were considered.

Table 3

Distribution of the 5′ haplotype frequencies among different populations and the association between Brazilian populations and European ancestry.

		Haplotype frequency				Haplotype-specific test (BRZ)
HAP ID	5′ HAP	YRI	ASN	CEU	BRZ	EUR-ancestry
						B	p
5′H01	CTT	0.017	−	0.325	0.314	0.048	0.036 *
5′H02	TTC	0.883	0.494	0.208	0.309	− 0.059	0.013 *
5′H03	CAC	−	0.410	0.342	0.292	0.020	0.395
5′H04	CTC	−	0.090	0.125	0.086	− 0.039	0.281
5′H05	TTT	0.100	−	−	−	−	−

HAP ID = Haplotype identity; 5′ HAP = Haplotypes in the 5′ gene region consisting of Cdx2, rs10783219 and rs3890734 SNPs. HapMap populations: ASN = Asian derived ancestry, CEU = European derived ancestry, YRI = African derived ancestry; BRZ = total Brazilian population. EUR = European genetic ancestry. B = raw regression coefficient. p < 0.05 indicates a significant correlation.

Table 4

Distribution of the 3′ haplotype frequencies among different populations and the association between Brazilian populations and European (EUR) and Amerindian (AMR) ancestry.

		Haplotype frequency				Haplotype-specific test (BRZ)
HAP ID	3′ HAP	YRI	ASN	CEU	BRZ	AMR-ancestry		EUR-ancestry
						B	p	B	p
3′H01	ACGAGG	0.083	0.635	0.408	0.347	0.001	0.906	− 0.013	0.546
3′H02	GTTGTA	0.142	0.056	0.475	0.321	− 0.048	0.001*	0.053	0.020*
3′H03	GCTATA	0.125	0.219	0.100	0.101	0.041	0.018*	− 0.060	0.078
3′H04	GCTATG	−	0.022	−	0.035	0.045	0.099	− 0.105	0.041*
3′H05	ACGAGA	−	−	−	0.028	− 0.054	0.518	0.113	0.241
3′H06	GTGGTA	−	−	−	0.027	0.043	0.187	− 0.040	0.553
3′H07	ACTATA	0.117	0.011	−	0.023	0.029	0.487	0.023	0.845
3′H08	GCGATA	0.058	−	−	0.022	0.057	0.093	− 0.030	0.737
3′H09	ACGATG	−	−	−	0.020	− 0.045	0.378	0.088	0.262
3′H10	GTTGTG	−	−	−	0.018	0.034	0.452	− 0.003	0.979
3′H11	ACTATG	−	−	−	0.016	0.049	0.209	− 0.119	0.156
3′H12	GTTATA	0.017	0.011	−	0.016	0.003	0.944	− 0.090	0.181
3′H13	ATTGTA	−	−	−	0.013	− 0.134	0.161	0.280	0.082
3′H14	ACTGTA	0.092	−	−	0.012	− 0.122	0.164	0.199	0.188
3′H15	GCGAGG	0.167	0.017	−	−	−	−	−	−
3′H16	GTTAGG	0.108	−	−	−	−	−	−	−
3′H17	ACGATA	0.050	−	−	−	−	−	−	−
3′H18	ATTATA	0.017	0.011	−	−	−	−	−	−
3′H19	GCTGTA	−	0.011	−	−	−	−	−	−

HAP ID = Haplotype identity; 3′ HAP = Haplotypes in the 5′ gene region consisting of rs2248098, BsmI, ApaI, TaqI, rs9729 and rs7968585 SNPs. HapMap populations: ASN = Asian derived ancestry, CEU = European derived ancestry, YRI = African derived ancestry; BRZ = total Brazilian population. AMR = Amerindian genetic ancestry; EUR = European genetic ancestry. B = raw regression coefficient.

(p < 0.05) indicates a significant correlation.

Regions with no haplotype blocks and a low LD were found between SNPs rs3890734 and rs2853559 and between SNPs rs2254210, FokI and rs886441, for which no block structure was found in any of the populations studied (Figure 1). The test to identify possible recombination hotspots showed two regions of greater intensity relative to the background recombination rate. In Brazilians, the region at FokI and rs886441 had a recombination rate 60 times higher than the background rate, whereas between rs3890734 and rs2853559 the rate was 18 times higher (Figure 1).

The test for population structure in the Brazilian population using autosomal AIMs identified a higher probability for a three-hybrid population, and assigned estimates of contributions as 0.771 ± 0.044 for European, 0.143 ± 0.019 for African and 0.085 ± 0.015 for Amerindian. The individual estimates of ancestry proportion showed that most individuals had a widely distributed three-hybrid pattern of variation with a trend towards a higher contribution by European ancestry. The individual estimates were later used as a quantitative trait in conditional regression tests.

The conditional regression omnibus test indicated a significant correlation with the 5′ haplotypes only when estimates of European ancestry were used (p = 0.027) and for 3′ haplotypes there was a significant correlation when using estimates of Amerindian and European ancestry (p = .018 and 0.041, respectively); there was no significance when using African ancestry proportions (p > 0.111). Haplotype-specific tests in the 5′ region identified two haplotypes with opposing raw regression coefficients (B) indicating that the first haplotype (5′H01) was positively related to European ancestry while the second (5′H02) had a negative correlation (Table 3). In the 3′ region, haplotypes with opposing raw regression coefficient signs (B) were significant for Amerindian ancestry (3′H02 and 3′H03) and two haplotypes (3′H02 and 3′H04) were significant for European ancestry (Table 4). When the 3′ haplotype block was reduced to only the Bsm-Apa-Taq haplotype, a similar effect of ancestry was observed, but in this case the haplotype with a positive regression in Amerindian ancestry had a negative regression in European ancestry and vice-versa (Table 5).

Table 5

Distribution of the Bsm-Apa-Taq haplotype frequencies among different populations and the association between Brazilian populations and European (EUR) and Amerindian (AMR) ancestry.

Bsm-Apa-Taq		Haplotype frequency				Haplotype-specific test (BRZ)
HAP ID	HAP	YRI	ASN	CEU	BRZ	AMR-ancestry		EUR-ancestry
						B	p	B	p
baT	CGA	0.367	0.657	0.417	0.418	0.013	0.279	− 0.006	0.766
BAt	TTG	0.158	0.062	0.475	0.332	− 0.044	0.001*	0.058	0.008*
bAT	CTA	0.258	0.258	0.108	0.175	0.037	0.005*	− 0.065	0.007*
Bat	TGG	−	−	−	0.023	− 0.090	0.189	0.071	0.422
BAT	TTA	0.125	0.017	−	0.022	− 0.008	0.858	− 0.089	0.155
BaT	TGA	−	−	−	0.017	− 0.022	0.760	− 0.022	0.800
bAt	CTG	0.083	−	−	0.013	− 0.073	0.286	0.121	0.242

HAP ID = Haplotype identity; HAP = Haplotypes in the gene region consisting of BsmI, ApaI and TaqI SNPs. HapMap populations: ASN = Asian derived ancestry, CEU = European derived ancestry, YRI = African derived ancestry; BRZ = total Brazilian population. AMR = Amerindian genetic ancestry; EUR = European genetic ancestry. B = raw regression coefficient.

(p < 0.05) indicates a significant correlation.

Discussion

In this work, a panel of allelic diversity at the VDR gene locus was generated and the effect of genetic ancestry on haplotype distribution in an admixed population was evaluated. The results described here extend the information about VDR genotypes and provide the first map of VDR haplotypes for a Brazilian population.

The measures of population admixture evaluated with AIMs revealed that the biogeographical ancestral structure of the Brazilian population was European, African and Amerindian (in this order), as previously described (Lins ). However, Fst analysis of VDR SNPs revealed a greater distinction between the Brazilian population and the HapMap population of African origin than with the HapMap population of European origin. These results indicate that, in an admixed population, recent admixture cannot always eliminate ancestral LD block structures along chromosomes (Sawyer ; Tang ). Consequently, complex levels of population admixture can create analytical risks in research involving loci associated with susceptibility to disease and in populations with different profiles, especially in those undergoing admixture. This is likely the case of the VDR gene and bone metabolism phenotypes, such as bone mineral density and osteoporosis (Thakkinstian ; Uitterlinden ).

Some studies have clearly demonstrated significant associations between European ancestry and body composition traits in admixed populations (Bonilla ; Shaffer ), while others have shown the relationship between VDR haplotypes and fracture risk in Whites or osteoporosis in European and Asian populations (Thakkinstian ; Fang ). The present work improves our understanding of the effect of genetic ancestry on the distribution of VDR gene haplotypes in admixed populations and reinforces the importance of correcting for population admixture in genetic association studies that encompass bone-related phenotypes.

The comparative analysis of VDR revealed genetic heterogeneity involving different haplotype blocks in the four populations studied. Considering the recent admixture of Latin American people, especially Brazilians, the variation in the patterns of LD seen here is not surprising in view of demographic events and genetic factors such as drift and recombination during the process of admixture (Gabriel ; Liu ; Sawyer ). Haplotype structure analysis in other admixed populations has also revealed the importance of genetic heterogeneity since linkage disequilibrium increases or breaks down differently in different populations (Moraes ; Boldt ; Lohmueller ; Nakamoto ).

Random genetic drift generates large diversity among populations with the same continental origin, e.g., African, European, Amerindian or Asian (Rosenberg ). In the present case, the HapMap populations used for comparison were not the most representative sources for the Brazilian parental population studied here and differed from those used to estimate individual autossomic ancestry (which also do not represent the Brazilian parental populations). This limitation is extremely important since the groups used here represent more general continental populations and their characteristics should therefore not be extrapolated to specific populations such as the many that constituted the admixture in current Brazilians, particularly when studying disease-related polymorphisms. As with considerations about genetic ancestry and disease-associated SNPs, the patterns of LD and identification of genetic ancestry blocks should also be addressed prior to the selection of tagSNPs for association studies in admixed populations (Lins ).

The distribution of the 5′ haplotype for the VDR gene in BRZ was similar to that of the CEU population. The regression coefficient confirmed this relationship. In this case, the haplotype that showed a positive association with European ancestry in Brazilians (5′H01) was absent in the ASN population and had a low frequency in YRI (Table 3). In contrast, the haplotype that correlated negatively with European ancestry (5′H02) had a higher frequency in YRI and ASN than in CEU. These results indicate that both European and non-European ancestries contributed to haplotype variation in our admixed population. Notably, this same phenomenon occurred in the 3′ region haplotypes, but in this case there was also a significant association with the Amerindian contribution (Table 4).

The African ancestry showed no correlation with any of the haplotypes examined. Indeed, variation in the extent and amount of LD in the VDR gene is lower in African populations than in European or East-Asians (Nejentsev ; Fang ), suggesting that extensive recombination precluded the extension of LD in African populations. In agreement with this, the 3′ haplotypes defined by the BRZ population exhibited high diversity and low frequencies in YRI (Table 4) when compared to Bsm-Apa-Taq haplotypes, which had a higher LD (Table 5). Moreover, the FokI SNP was not in linkage disequilibrium with any of the SNPs in the populations tested. This observation corroborates previous findings (Nejentsev ; Fang ) and suggests a major site of haplotype breakage and recombination in the VDR gene that is independent of ethnicity.

Previous studies of the VDR gene in Brazilians only sampled a few loci but some complex phenotypes (Lazaretti-Castro ; Hauache ; de Brito Junior ; Maistro ; Goulart ; Gentil , 2009; Moreno Lima ; Rezende ). Some studies have tried to correct the genetic heterogeneity of the Brazilian population by using self-reported ancestry or physical traits as proxy for different ethnic groups. For instance, Rezende reported no difference in the distribution of VDR haplotypes in Brazilian self-reported Blacks and Whites. This may be explained by the fact that dissociation between physical appearance and genetic ancestry in Brazilians (Parra ; Marrero ) may have created spurious similarity in genotype and haplotype frequencies among these groups because of the population substructure. As shown here, admixture in the Brazilian population provides the opportunity to segregate the contribution of individual ancestry from genetic ancestry blocks at the VDR gene locus.

In conclusion, the results of this investigation provide a large map of haplotypes for the entire VDR gene and intragenic regions in a carefully sampled Brazilian population. Comparison with the HapMap data was essential for understanding the patterns of LD and haplotype variation among populations and for elucidating the effects of admixture on this diversity. Our findings also support studies that involve tagSNP selection between different ethnic populations (Nejentsev ; Lins ). Genetic polymorphisms such as those observed here may partly explain ethnic differences in vitamin D3 status and their relationship to bone phenotype and hormonal homeostasis, particularly in elderly women, as well as the role of environmental factors such as diet, lifestyle and sun exposure (Uitterlinden ; Dusso ; Valdivielso and Fernandez, 2006).