Current limitations of SNP data from the public domain for studies of complex disorders: a test for ten candidate genes for obesity and osteoporosis.

BACKGROUND: Public SNP databases are frequently used to choose SNPs for candidate genes in the association and linkage studies of complex disorders. However, their utility for such studies of diseases with ethnic-dependent background has never been evaluated.
RESULTS: To estimate the accuracy and completeness of SNP public databases, we analyzed the allele frequencies of 41 SNPs in 10 candidate genes for obesity and/or osteoporosis in a large American-Caucasian sample (1,873 individuals from 405 nuclear families) by PCR-invader assay. We compared our results with those from the databases and other published studies. Of the 41 SNPs, 8 were monomorphic in our sample. Twelve were reported for the first time for Caucasians and the other 29 SNPs in our sample essentially confirmed the respective allele frequencies for Caucasians in the databases and previous studies. The comparison of our data with other ethnic groups showed significant differentiation between the three major world ethnic groups at some SNPs (Caucasians and Africans differed at 3 of the 18 shared SNPs, and Caucasians and Asians differed at 13 of the 22 shared SNPs). This genetic differentiation may have an important implication for studying the well-known ethnic differences in the prevalence of obesity and osteoporosis, and complex disorders in general.
CONCLUSION: A comparative analysis of the SNP data of the candidate genes obtained in the present study, as well as those retrieved from the public domain, suggests that the databases may currently have serious limitations for studying complex disorders with an ethnic-dependent background due to the incomplete and uneven representation of the candidate SNPs in the databases for the major ethnic groups. This conclusion attests to the imperative necessity of large-scale and accurate characterization of these SNPs in different ethnic groups.

Background

A single nucleotide polymorphism (SNP) is generally defined as a stable substitution of a single base with a frequency of more than 0.01 in at least one population [1]. In human genetic studies, SNPs are simply referred to as bi-allelic markers since the other types (tri-allelic and tetra-allelic SNPs) are very rare in the human genome [2]. SNPs have been recognized as an important tool in human genetics and medicine [3,4]. Because SNPs are abundant and scattered throughout the whole human genome, they have been widely used in the genetic association studies of various complex diseases such as obesity, osteoporosis, asthma, hypertension (see, e.g., [5-7]). Various analyses of SNPs across the human genome have been also conducted to determine haplotype patterns in human populations [8-12]. These data are very useful to study the genetic basis of common complex diseases [13].

At present there are several SNP public databases. The largest are dbSNP and HGVbase containing together several million SNPs [14]. There also exist other relatively small or specific SNP databases, e.g., TSC [15], JSNP [16], HOWDY [17], GeneSNPs [18]. While the databases have been continuously expanded, the quality and completeness of the deposited SNP data remain to be of particular importance and to be assessed. Recent evaluation of quality and comprehensiveness of about four million candidate SNPs from some public and Celera databases showed that about 6–12% SNPs could not be validated [19]. They represent rare variants, population-specific SNPs and sequencing errors. The other studies reported that only about 50% SNPs are common in any given population. Therefore, measurement of allele frequency and linkage disequilibrium for SNPs in databases should be done to efficiently select a minimal subset of SNPs for population association studies of complex diseases [20,21].

Obesity and osteoporosis are common complex disorders with continuously growing burden and cost for their prevention and treatment. For example, the most recent data for the United States, derived from the third National Health and Nutrition Examination Survey (1988–94), showed ~20% of US men and ~25% of US women are obese and these proportions are increasing [22]. The same tendency has been observed in other human populations [23]. Osteoporosis is another major health problem, particularly in the elderly. More than 40% of postmenopausal women, on average, will suffer at least one osteoporotic fracture [24]. This disorder incurred an estimated direct cost of ~14 billion dollars in the USA alone in 1995 [25].

SNPs of candidate genes for obesity and osteoporosis may be of particular importance in genetic studies of these disorders. The allele frequencies are important in the selection of SNPs for studying complex diseases [26]. For example, association studies of osteoporosis usually generate inconsistent results in different ethnic groups [27]. One important reason is that polymorphisms associated significantly with certain osteoporotic phenotypes in a given ethnic group may be absent or rare in another ethnic group [28-31]. Moreover, racial differences in the prevalence of certain allele could account for certain proportion of disease trait variation between different ethnicities [32]. Finally, a comparison of SNP allele frequencies among different ethnicities can provide valuable information for mapping by admixture linkage disequilibrium (MALD) [33].

Despite the importance of the above data, no attempts have been made to check the quality and completeness of SNPs of candidate genes for obesity and osteoporosis in the public databases. In fact, most allele frequencies in the databases were obtained from studies of relative small samples (dozens or so), which could yield large sampling errors. In addition to the databases, there are abundant SNP data from many association/linkage studies of these two common complex diseases. These individually published data have not yet been assessed in reference to each other, or to those in the databases. In most association studies, the regular sample sizes used to detect allele frequencies have only been of the order of a hundred or so. However, to obtain reliable estimates of allele frequency distributions, relatively larger sample sizes are needed. For example, for an SNP, to let allele frequency estimate's |error| = 0.05, the required sample sizes for frequency distribution (0.5, 0.5), (0.7, 0.3), and (0.9, 0.1) should be at least 1150, 1000, and 400, respectively [34]. Otherwise, the estimates of allele frequencies may be biased to a large extent.

In the present study, we determined allele frequencies of 41 SNPs in 10 candidate genes for obesity and/or osteoporosis in a large American-Caucasian sample (1,873 subjects from 405 nuclear families). We then compared them with the corresponding data for the major ethnic groups obtained from the public databases and various individual studies to check their accuracy and completeness, and to compare the SNP allele frequencies in different ethnic groups.

Results

SNP polymorphism in the studied Caucasian sample

Of the 41 candidate SNPs for obesity and/or osteoporosis identified in the studied Caucasian sample (Table 1), 29 (70.7%) had minor allele frequency ≥0.1. Eight SNPs (19.5%) were monomorphic. Allele frequencies of twelve SNPs in Caucasians were reported for the first time (Table 2). Five of these SNPs were monomorphic and the other seven were polymorphic with minor allele frequencies ranged 0.019–0.398.

Table 1

Summary information about the studied SNPs

SNPa	Geneb	Namec	dbSNP Accession	Mutationd	Domain	Minor allele frequency	Forward primer	Reverse primer
1	APOE		ss12568587	G-C	5' UTR	0.357	TCCCCAGGAGCCGGTGA	CCCCAAGCCCGACCCC
2	APOE		ss12568609	G-A	Intron 2	0.399	CCTCAGGTGATCTGCCCGTTTC	ACTCCTGGGCTCAAGTGATCCTC
3	APOE		ss12568607	T-C	Exon 4	0.149	CGGGCACGGCTGTCCAA	CGAGCATGGCCTGCACCTC
4	APOE		ss12568612	C-T	Exon 4	0.087	GCTGCGTAAGCGGCTCC	GCGGCCCTGTTCCACC
5	COL1α1		ss12568606	G-T	5' UTR	0.154	GCACCCTGCCCTAGACCAC	CCTAGTGCCAGCGACTGCA
6	COL1α1	Sp1	ss12568597	G-T	Intron 1	0.188	CCAATCAGCCGCTCCCATTC	CATCGGGAGGGCAGGCTC
7	COL1α1		ss12568598	G-T	Exon 8	0.000	GGAAGACTGGGATGAGGGCA	GGCTCGCCAGGCTCACC
8	COL1α1		ss12568584	G-A	Exon 45	0.019	CTCAGCCTTCCCTGGCCAA	AGGCGGAAGTTCCATTGGCATC
9	ER-α		ss12568579	A-G	Exon 1	0.484	TTGAGCTGCGGACGGTTCA	CGCCGGTTTCTGAGCCTTC
10	ER-α	PvuII	ss12568596	T-C	Intron 1	0.449	TGGGATTCCAGGCATGAACCAC	TGGCGTCGATTATCTGAATTTGGCC
11	ER-α		ss12568619	G-A	Intron 3	0.257	CCCAGAAACAAGTCATCTGCTATTGACA	TGTAACAAAAGGTTAACAATGGTTAGCCC
12	ER-α		ss12568618	G-C	Exon 4	0.218	ACAGCCTGGCCTTGTCCC	CAGGTTGGTCAGTAAGCCCATCA
13	ER-α		ss12568585	G-A	Intron 4	0.098	GATCAATGAAGTGGGTCTTGAAAAACCAA	GGTGACAAGCTGGAAATCTAAGCTTCA
14	ER-α		ss12568605	G-A	Intron 6	0.119	GGAACGGCCCTTGGAAATTGTAAA	CTGCCTACAGAATACAGTCAGCCA
15	ER-α		ss12568617	G-A	Exon 8	0.203	TCGCATTCCTTGCAAAAGTATTACATCAC	CAAGCAAATGAATGGCCACTCATCTAGAAA
16	IL-6		ss12568616	C-G	5' UTR	0.434	GGGCAGAATGAGCCTCAGACATC	GACATGCCAAAGTGCTGAGTCACTAATA
17	IL-6		ss12568586	A-T	Exon 4	0.005	CCTCCACTGCAAAGGATTTATTCAACA	CATGTCCTGCAGCCACTGGTTC
18	LEPR		ss12568615	A-G	Exon 4	0.290	AGATTTAAGTTGTCTTGCATGCCACC	TTAAGCCCAGCATCCATTAGCTATTCTTTC
19	LEPR		ss12568604	G-C	Exon 14	0.182	GAGTAATTGGAGCAATCCAGCCTACA	GCTTCAGCCACTGTACATCTTAGCTC
20	LEPR		ss12568614	G-A	Exon 20	0.353	GCCACGCTGATCAGCAACTC	CCCTTGACTTGTCAGTCAAAAGCAC
21	PTHR1		ss12568589	G-A	Intron 1	0.001	GACTTACATTAGGATTCAAGGTTACTGCCA	GGGACGCAAGCCTGAGTCC
22	PTHR1		ss12568592	A-G	Intron 2	0.398	GCAGAACCCTAAGGGCTTGTCA	GGCGGGACCCAGGATACA
23	PTHR1		ss12568591	C-T	Intron 8	0.374	CGAGCCTCAATTCAGGTGAATCTAACC	CCCGCCCCAAGTGGAACA
24	PTHR1		ss12568588	G-A	Intron 10	0.397	CCTTGAGCCCTTGGTTTTCCTTTC	GCTCCGGGAACAAAAAGTGGATCA
25	PTHR1		ss12568590	G-A	Exon 13	0.380	CTACAAGGCTCAAATTGCCCCAAA	TTGGCGTCCACTACATTGTCTTCA
26	TGF-β1		ss12568613	C-T	5' UTR	0.311	GGGCCCAGTTTCCCTATCTGTAAA	CTGGGCCACCGTCCTCATC
27	TGF-β1		ss12568603	+C/-C	Intron 4	0.021	CCACGCCCCACTTATCTATCCC	GGAAAGGCCGGTTCATGCCA
28	TGF-β1		ss12568593	C-T	Exon 5	0.008	CAGGCTACAAGGCTCACCTGAA	GGTTCACTACCGGCCGC
29	TGF-β1		ss12568602	C-T	Intron 5	0.274	GGCTTGTCTTAAGCATTGCGTGAAATTAA	GTACAGCTGCCGCACGC
30	TNFR2		ss12568594	T-C	5' UTR	0.000	TGCACTCGGCCTGTTTAGACTC	CTGTTCATCTGCCCCCTGCC
31	TNFR2		ss12568611	T-G	Exon 6	0.198	AGCCACCCCAGCCACTC	GCTTGGAGCAGTGCTGGGTTC
32	UCP3		ss12568601	C-T	5' UTR	0.265	CACTGCCCTCACCAGCCA	GTGAGTCCTGCCACGGCA
33	UCP3		ss12568600	G-A	Exon 2	0.000	GCCCTAAAGGGACTGGGCA	GAAAGGTAACGAGGTCAGCAAAACA
34	UCP3		ss12568595	T-C	Exon 3	0.258	TGATTCCCGTAACATCTGGACTTTCATC	CTGCCTAAATCCCCTTAGCAGAAAAAAAAA
35	UCP3		ss12568580	T-C	Exon 5	0.447	CCTAACAGGAACTTTGCCCAACATCA	TCCACGGAGTTCTGGGTTCC
36	UCP3		ss12568599	C-T	Exon 7	0.000	CCTAACAGGAACTTTGCCCAACATCA	TCCACGGAGTTCTGGGTTCC
37	VDR		ss12568583	G-A	5' UTR	0.281	CAGCATGCCTGTCCTCAGC	CCAGTACTGCCAGCTCCCA
38	VDR	FokI	ss12568581	C-T	Exon 2	0.373	TGGCCCTGGCACTGACTC	GGCACGTTCCGGTCAAAGTC
39	VDR		ss12568582	C-T	Exon 4	0.000	GGACAGTCTGCGGCCCA	CCCTACTCCCTGGGCCC
40	VDR	BsmI	ss12568610	G-A	Intron 8	0.419	GTGCCCCTCACTGCCCTTA	CCTCAAATAACAGGAATGTTGAGCCCA
41	VDR	TaqI	ss12568608	T-C	Exon 9	0.408	GGGCCAGGCAGTGGTATCAC	AGGTCGGCTAGCTTCTGGATCA

aDesignation in the present study. bAbbreviations: APOE, apolipoprotein E; COL1α1, collagen type I α1; ER-α, estrogen receptor-α; IL-6, interleukin-6; LEPR, leptin receptor; PTHR1, parathyroid hormone (PTH)/PTH-related peptide receptor type 1; TGF-α1, transforming growth factor-ß1; TNFR2, tumor necrosis factor receptor 2; UCP3, uncoupling protein 3; VDR, vitamin D (1,25-dihydroxyvitamin D3) receptor. cCommon name. dMinor alleles are given in bold.

Table 2

Minor allele frequencies of SNPs firstly reported for Caucasians by this study

Gene	SNP	Database or referencea	Sample size	Minor allele frequency
COL1α1	7	[52]	292	0.000
COL1α1	8	dbSNP	1,861	0.019
ER-α	14	JSNP	1,861	0.119
IL-6	17	dbSNP	190	0.005
PTHR1	21	JSNP	380	0.001
PTHR1	22	JSNP	1,861	0.398
PTHR1	23	JSNP	1,863	0.374
PTHR1	24	JSNP	1,855	0.397
PTHR1	25	JSNP	1,851	0.380
TNFR2	30	JSNP	281	0.000
VDR	37	JSNP	1,861	0.281
VDR	39	dbSNP	190	0.000

For gene abbreviations and SNP designations see Table 1. aThe source of information about the given SNP.

Data on the candidate SNPs for obesity and/or osteoporosis from major ethnic groups: an analysis of literature and databases

Caucasians

In available literature and public SNP databases, we found data about allele frequencies of 29 SNPs (representing 9 genes) in Caucasians, which were shared with the present study. The corresponding data are given in Table 3. The allele frequencies of only 14 SNPs had been previously reported for Caucasians to the public databases (dbSNP and/or HGVbase). No significant differences were found between the allele frequencies in our study and in the public databases, except for SNP31 in Australians (P = 0.05, Table 3).

Table 3

Comparison of SNP allele frequencies of the studied Caucasian sample with the data on Caucasians from databases and literature

Gene	SNP	Reference population	Database or reference	Data from a database or reference		Caucasiansa
				Frequencyb	n	Frequency	n
APOE	1	American Caucasians	[53]	0.33	220	0.357	1,862	0.48
	2	American Caucasians	dbSNP	0.33	24	0.399	1,839	0.64
		Europeans from North Karelia, Finland	dbSNP	0.4	24	0.399	1,839	0.84
	3	European-Americans	dbSNP	0.12	24	0.149	1,858	0.91
		Caucasiansc	HGVbase	0.16	152	0.149	1,858	0.81
		Europeans from North Karelia, Finland	dbSNP	0.23	24	0.149	1,858	0.37
	4	European-Americans	dbSNP	0.19	24	0.087	1,783	0.16
		Caucasiansc	HGVbase	0.08	152	0.087	1,783	0.87
		Europeans from North Karelia, Finland	dbSNP	0.04	24	0.087	1,783	0.25
COL1 α1	5	White women with Spanish ancestors	[54]	0.13	256	0.154	1,856	0.36
	6	American Caucasians	[31]	0.187	637	0.188	1,855	0.98
ER-α	9	Australian females	[55]	0.44	125	0.484	1,861	0.39
		Greek males	[56]	0.49	50	0.484	1,861	0.95
	10	American females	[57]	0.45	253	0.449	1,861	0.96
		Finnish females	[58]	0.42	322	0.449	1,861	0.38
	11	SANGER 12 DNAs of Caucasian origin	dbSNP	0.08	12	0.257	1,857	0.32
	12	Australian females	[55]	0.23	120	0.218	1,859	0.84
		Slovenian females	[59]	0.3	85	0.218	1,859	0.10
		Caucasiansc	dbSNP	0.27	31	0.218	1,859	0.59
	13	Caucasiansc	HGVbase	0.08	42	0.098	1,863	1.00
	15	Caucasiansc	dbSNP	0.18	31	0.203	1,863	0.90
IL-6	16	Caucasians	dbSNP	0.5	31	0.434	1,860	0.58
		English	[60]	0.57	2,560	0.434	1,860	<0.01
		Italians	[61]	0.45	183	0.434	1,860	0.74
		Spaniards	[62]	0.56	118	0.434	1,860	0.01
		Finns	HGVbase	0.46	400	0.434	1,860	0.37
		Finns	HGVbase	0.48	50	0.434	1,860	0.62
		Finns	HGVbase	0.39	42	0.434	1,860	0.68
		Finns	[63]	0.45	400	0.434	1,860	0.60
LEPR	18	English males	HGVbase	0.27	322	0.290	1,839	0.41
		Caucasiansc	[64]	0.28	56	0.290	1,839	0.98
		Swedish males	[65]	0.24	284	0.290	1,839	0.09
	19	Swedish males	[65]	0.15	284	0.182	1,776	0.21
	20	Caucasiansc	[64]	0.38	56	0.353	1,860	0.79
		Caucasiansc	dbSNP	0.4	42	0.353	1,860	0.64
		English males	HGVbase	0.36	322	0.353	1,860	0.77
TGF-ß1	26	Caucasiansc	dbSNP	0.43	31	0.311	1,840	0.21
		Caucasiansc	HGVbase	0.26	42	0.311	1,840	0.59
	27	Caucasiansc	[66]	0.02	304	0.021	1,863	0.98
	28	Danes	[66]	0.027	302	0.008	190	0.33
	29	Caucasianc female DZ twins	[67]	0.25	1,802	0.274	1,852	0.08
TNFR2	31	UK Caucasians	[68]	0.23	192	0.198	1,842	0.39
		Germans	[69]	0.25	94	0.198	1,842	0.33
		Australians	HGVbase	0.26	197	0.198	1,842	0.05
UCP3	32	American Caucasians	[70]	0.21	56	0.265	1,767	0.49
		Danes	[71]	0.27	857	0.265	1,767	0.85
		French	[72]	0.22	894	0.265	1,767	0.01
	33	Caucasiansc	[73]	0.01	165	0.000	295	0.44
	34	American Caucasians	[70]	0.2	58	0.258	1,862	0.38
	35	Caucasiansc	[74]	0.5	501	0.447	1,863	0.04
	36	Caucasiansc	[73]	0.02	165	0.000	295	0.05
VDR	38	American Caucasians	[75]	0.42	49	0.373	1,861	0.61
		English	[76]	0.4	241	0.373	1,861	0.47
		Australians	[77]	0.35	43	0.373	1,861	0.87
	40	Caucasiansc	dbSNP	0.42	107	0.419	1,860	0.94
		English females	[76]	0.44	241	0.419	1,860	0.61
		Australian males	[77]	0.47	39	0.419	1,860	0.66
		Swiss	[78]	0.39	197	0.419	1,860	0.49
		French infants	[79]	0.38	589	0.419	1,860	0.07
	41	Caucasiansc	dbSNP	0.48	31	0.408	1,859	0.51
		French	[80]	0.35	99	0.408	1,859	0.33
		Australians	[81]	0.41	68	0.408	1,859	0.94

For gene abbreviations and SNP designations see Table 1. a Present study. b Frequency of the allele corresponding to a minor allele in the present study. c Caucasians with unidentified ethnic background.

In contrast to the databases, available literature contained data on allele frequencies of the 23 SNPs. We determined significant differences between our and literature data for allele frequencies of four SNPs (SNP16, SNP32, SNP35, and SNP36). We have not found any allele frequency data on the SNPs of the PTHR1 gene (SNP21-SNP25) for Caucasians. Some SNP indicated large discrepancies in allele frequencies between the data from the different databases. For example, the minor allele frequencies of SNP26 for two Caucasian samples reported to dbSNP and HGVbase are 0.43 and 0.26, respectively (Table 3). The differences between the values are not statistically significant due to the small sample sizes (n = 31 and n = 42, respectively), so it is impossible to conclude what actual allele frequencies are in these populations and whether they are indeed significantly different.

African/African-Americans

Table 4 presents our SNP data in comparison with the corresponding data on African and/or African-American samples obtained from literature and the databases. In the public domain, we found in total only 18 SNPs from seven genes, which are shared with our study. The candidate genes for obesity and/or osteoporosis appear to be underrepresented for Africans/African-Americans in the public SNP databases and literature: only 11 shared SNPs were found in the databases and seven in the literature. Among those, one (SNP36) was monomorphic and four had minor allele frequency <0.1. Significant differences in allele frequencies between Caucasians and Africans/African-Americans were found at three SNPs (SNP16, SNP34 and SNP35). Three more SNPs (SNP32, SNP33, and SNP41) manifested nearly significant differences in allele frequencies between these ethnic groups (Table 4). We found no data on the SNPs of three genes (COL1α1, PTHR1, and TNFR2) for Africans/African-Americans in the public domain, which would be shared with our study.

Table 4

Comparison of SNP allele frequencies of the studied Caucasian sample with the data on Africans/African- Americans from databases and literature

Gene	SNP	Reference population	Database or references	Data from a database or reference		Caucasiansa		P
				Frequencyb	n	Frequency	n
APOE	2	African-Americans	dbSNP	0.48	24	0.399	1,839	0.52
	3	African-Americans	dbSNP	0.10	24	0.149	1,858	0.56
	4	African-Americans	dbSNP	0.04	24	0.087	1,783	0.25
ER-α	11	SANGER 12 DNAs of African-American origin	dbSNP	0.08	12	0.257	1,857	0.32
	12	African/African-Americans	dbSNP	0.27	24	0.218	1,859	0.61
	13	TSC panel of African-Americans	dbSNP	0.16	40	0.098	1,863	0.27
	15	African/African-Americans	dbSNP	0.28	23	0.203	1,863	0.42
IL-6	16	African/African-Americans	dbSNP	0.96	24	0.434	1,860	<0.01
LEPR	20	TSC panel of African-Americans	dbSNP	0.45	42	0.353	1,860	0.25
TGF-ß1	26	TSC panel of African-Americans	TSC	0.27	42	0.311	1,840	0.61
		African/African-Americans	dbSNP	0.15	24	0.311	1,840	0.07
UCP3	32	African-American females	[70]	0.15	57	0.265	1,767	0.07
	33	African-Americans	[73]	0.08	18	0.000	295	0.06
	34	African-American females	[70]	0.45	63	0.258	1,862	<0.01
	35	African-American females	[70]	0.247	73	0.447	1,863	<0.01
		African/African-Americans	[74]	0.20	276	0.447	1,863	<0.01
	36	African-Americans	[70]	0.00	18	0.000	295	10.00
VDR	38	African-American males	[77]	0.244	37	0.373	1,861	0.15
	40	African-American males	[77]	0.312	32	0.419	1,860	0.30
		Nigerians	[82]	0.37	93	0.419	1,860	0.41
	41	African/African-Americans	dbSNP	0.35	24	0.408	1,859	0.71
		African-American males	[77]	0.244	39	0.408	1,859	0.06

For gene abbreviations and SNP designations see Table 1. a Present study. b Frequency of the allele corresponding to a minor allele in the present study.

Similar to the case with Caucasians, there are large inconsistencies between some SNP data for Africans/African-Americans from the different databases. For example, the reported minor allele frequencies of SNP26 from the TSC and dbSNP databases are 0.27 and 0.15, respectively (Table 4). The former is not significantly different (P = 0.61) from the corresponding value in Caucasians, while the latter is nearly significantly different (P = 0.07). However, it is hard to determine, which value is true because of the small sample sizes (n = 42 and n = 24, respectively) in the databases.

Asians

As compared to our SNP data, this ethnic group had corresponding data on 22 SNPs (8 genes) in the public domain (Table 5). The data on 13 SNPs of the 6 genes were available from the public databases, and on 10 SNPs of the 5 genes were reported in the literature. Among those SNPs, two (SNP16 and SNP33) were monomorphic and one (SNP36) had minor allele frequency <0.1. Two more SNPs of the VDR gene (SNP40 and SNP41) have a minor allele frequency about this value; however, this frequency varies in different Asian populations. We observed a significant inconsistency between the values of a minor allele frequency for SNP40 in Asians from the database and literature, respectively. In the database, this value (0.41) is about 3–10 times higher than that reported in the literature (0.043–0.13).

Table 5

Comparison of SNP allele frequencies of the studied Caucasian sample with the data on Asian populations from databases and literature

Gene	SNP	Reference population	Database or reference	Data from a database or reference		Caucasiansa		P
				Frequencyb	n	Frequency	n
ER-α	9	Japanese	[83]	0.412	200	0.484	1,861	0.06
		Thai females	[84]	0.337	129	0.484	1,861	<0.01
	10	Korean males	[85]	0.39	219	0.449	1,861	0.11
		Korean females	[86]	0.415	248	0.449	1,861	0.35
		Japanese children	[87]	0.431	102	0.449	1,861	0.80
		Taiwanese	[88]	0.39	246	0.449	1,861	0.09
		Thai females	[89]	0.414	134	0.449	1,861	0.49
	11	SANGER 12 DNAs of Asian origin	dbSNP	0.36	12	0.257	1,857	0.49
	12	Japanese	[83]	0.483	200	0.218	1,859	<0.01
		Japanese females	[90]	0.481	306	0.218	1,859	<0.01
	13	Japanese	JSNP	0.34	744	0.098	1,863	<0.01
		TSC panel of unrelated Asians	TSC	0.27	32	0.098	1,863	0.01
	14	Japanese	JSNP	0.296	752	0.119	1,861	<0.01
	15	Japanese	JSNP	0.18	741	0.203	1,863	0.21
		Japanese	[83]	0.21	200	0.203	1,863	0.88
IL-6	16	Chinese	[31]	0.997	147	0.434	1,860	<0.01
		Chinese males	[31]	0.998	259	0.434	1,860	<0.01
		Japanese	[31]	1.000	388	0.434	1,860	<0.01
LEPR	20	TSC panel of unrelated Asians	dbSNP	0.13	42	0.353	1,860	<0.01
		Japanese	JSNP	0.15	748	0.353	1,860	<0.01
PTHR1	22	Japanese	JSNP	0.43	747	0.398	1,861	0.14
	23	Japanese	JSNP	0.43	748	0.374	1,863	0.01
	24	Japanese	JSNP	0.44	750	0.397	1,855	0.05
	25	Japanese	JSNP	0.43	746	0.380	1,851	0.02
TGF-ß1	26	TSC panel of unrelated Asians	dbSNP	0.47	36	0.311	1,840	0.06
TNFR2	31	Japanese	[91]	0.113	265	0.198	1,842	<0.01
		Thai	[92]	0.127	201	0.198	1,842	0.02
UCP3	33	Asiansc	[73]	0.00	11	0.000	295	10.00
	35	Japanese	JSNP	0.45	750	0.447	1,863	0.92
	36	Asiansc	[73]	0.04	11	0.000	295	0.54
VDR	37	Japanese	JSNP	0.23	750	0.281	1,861	0.01
	38	Japanese	[93]	0.37	195	0.373	1,861	0.99
	40	SANGER 12 DNAs of Asian origin	dbSNP	0.41	12	0.419	1,860	0.77
		Chinese females (Han nationality)	[94]	0.043	162	0.419	1,860	<0.01
		Taiwanese	[95]	0.083	90	0.419	1,860	<0.01
		Japanese	[93]	0.13	195	0.419	1,860	0.02
		Koreans	[96]	0.055	211	0.419	1,860	<0.01
		Thai females	[97]	0.107	84	0.419	1,860	0.01
	41	Chinese (Han nationality)	[98]	0.05	223	0.408	1,859	<0.01
		Japanese females	[99]	0.135	119	0.408	1,859	<0.01
		Koreans	[96]	0.055	120	0.408	1,859	<0.01

For gene abbreviations and SNP designations see Table 1. a Present study. b Frequency of the allele corresponding to a minor allele in the present study. c Asians with unidentified ethnic background.

A comparison of our Caucasian sample and the Asian populations revealed significant differences in allele frequencies of 13 SNPs located in the 6 genes (Table 5). Another SNP representing the TGF-ß1 gene (SNP26) indicated the nearly significant difference (P = 0.06). No such differences were observed only in the UCP3 gene.

Hispanic and Pacific Rim populations

These populations had the poorest data on the SNPs from the databases and literature as compared to the corresponding data from the present study. We found the data on only seven SNPs of the five candidate genes for Hispanics (Table 6) and five SNPs of the four genes for Pacific Rim populations (Table 7). SNP33 and SNP36 for the UCP3 gene were virtually monomorphic in Hispanics, similar to that in our Caucasian sample, while only one significant difference in allele frequencies was observed at SNP16 between these ethnic populations among the seven compared SNPs (Table 6).

Table 6

Comparison of SNP allele frequencies of the studied Caucasian sample with the data on Hispanics from databases and literature

Gene	SNP	Database or references	Data from a database or reference		Caucasiansa		P
			Frequencyb	n	Frequency	n
ER-α	12	dbSNP	0.239	23	0.218	1,859	1.00
	15	dbSNP	0.196	23	0.203	1,863	1.00
IL-6	16	dbSNP	0.795	22	0.434	1,860	<0.01
TGF-ß1	26	dbSNP	0.413	23	0.311	1,840	0.36
UCP3	33	[73]	0.00	27	0.000	295	1.00
	36	[73]	0.02	27	0.000	295	1.00
VDR	41	dbSNP	0.348	23	0.408	1,859	0.83

For gene abbreviations and SNP designations see Table 1. aPresent study. bFrequency of the allele corresponding to a minor allele in the present study.

Table 7

Comparison of SNP allele frequencies of the studied Caucasian sample with the data on Pacific Rim populations from databases

Gene	SNP	Database	Data from a database or reference		Caucasiansa		P
			Frequencyb	n	Frequency	n
ER-α	12	dbSNP	0.478	23	0.218	1,859	0.02
	15	dbSNP	0.326	23	0.203	1,863	0.19
IL-6	16	dbSNP	0.957	23	0.434	1,860	<0.01
TGF-ß1	26	dbSNP	0.458	24	0.311	1,840	0.19
VDR	41	dbSNP	0.13	23	0.408	1,859	<0.01

For gene abbreviations and SNP designations see Table 1. aPresent study. bFrequency of the allele corresponding to a minor allele in the present study

No data on the SNPs of the 10 candidate genes for obesity and/or osteoporosis were found for Pacific Rim populations in the available literature. All the data presented in Table 7 were obtained from the dbSNP database. Three SNPs out of the five differed significantly in the allele frequencies between Pacific Rim populations and the Caucasian sample in our study.

Some new candidate SNPs for genetic studies of obesity and osteoporosis

In this study, we identified a number of SNPs, which may be suitable for genetic studies of obesity and osteoporosis. In particular, the six polymorphic SNPs with minor allele frequency ≥0.1, which were firstly determined in Caucasians (Table 2), may be used in population association studies. A number of SNPs with significant variation in allele frequencies in populations of different ethnicity may be appropriate for studying a genetic basis of between-ethnic differences in the rates of obesity and/or osteoporosis. Examples are SNP34, SNP35, and SNP32 for Caucasians and African-Americans (Table 4) and several SNPs located in the ER-α, LEPR, PTHR1, TNFR2, and VDR genes for Caucasians and Asians (Table 5). Some of the SNPs may also be candidates for such studies, if their allele frequencies reported to the databases are validated by using sufficiently large sample sizes (e.g., SNP26 in Asians).

Discussion

Ethnic variation of SNP allele frequencies of the candidate genes for obesity and/or osteoporosis

We observed significant differences in allele frequencies of some studied SNPs for Caucasians (SNP16, SNP35, SNP32, and SNP36, Table 3) between our results and respective data from the databases and literature. This may be due to various factors, such as population admixture and relatively smaller sample size in the other studies [34]. For example, the sample in the present study consisted of the subjects of various ethnic backgrounds (German, French, Dutch, Swedish, some Portuguese and Italian backgrounds) and the relative proportions of these ethnic groups are not ascertained. In the mixed samples from other studies, these proportions may be different that thus influence the allele frequencies to a larger or lesser extent (e.g., SNP35 and SNP36, Table 3). Likewise, monoethnic samples may have allele frequencies significantly different from those in the mixed samples (e.g., SNP16 in English and Spaniards, SNP32 in French, Table 3). Our results showed that major ethnic groups have significantly different allele distribution at some SNP markers of candidate genes for complex disorders. The largest differences were observed between Asian and Caucasian populations, namely at 13 SNPs out of 22 compared (Table 5). Some genes indicated the significant differences at most SNPs compared (e.g., PTHR1 between Japanese and Caucasians, Table 5). These results are consistent with the previously reported data about high differentiation between Asians and Caucasians at candidate genes for bone mass [35,36].

Another important finding is the observed significant differences between Caucasian and Pacific Rim populations in SNP allele frequencies of three important candidate genes for osteoporosis, ER-α, IL-6, and VDR (Table 7). Interestingly, the respective SNP allele frequencies in Asian and Pacific Rim populations have very similar values (Tables 5 and 7). This may suggest that the Pacific Rim sample from the database has mainly Asian ancestry [37] or, perhaps there are some implications for an evolutionary history of these populations [38]. In any event, much more comprehensive data should be obtained to clarify these issues.

There are abundant data about significant between-ethnic differentiation at loci, which may underlie complex diseases, including obesity and osteoporosis (see, e.g., [36,39-41]). Given that there is well-known different incidence of these disorders in various ethnic groups (e.g., [42,43]), such genetic differentiation may have an important implication for studying a genetic basis of ethnicity-specific definitions of obesity and osteoporosis.

Another important application of the SNP markers with high large between-ethnic and small within-ethnic differences is MALD studies of complex diseases in the admixed populations with known parental ethnicities [33,44]. Our results showed that such SNPs exist in the candidate genes for obesity and osteoporosis. However, the large-scale screening of the candidate genes is necessary to identify such markers with sufficient density.

The results of our study suggest that ethnic heterogeneity of large samples may notably affect the observed allele frequencies. It is supported by the fact that several SNP allele frequencies of some sufficiently large monoethnic samples (e.g., SNP16 in English, SNP32 in French) significantly differ from those determined in our sample of mixed Caucasians (Table 3).

SNP databases: current limitations for studying complex diseases

Our analyses showed that the SNP databases in their current status might have some limitations for studies of complex disorders, especially in different ethnic groups, due to incomplete and/or uneven representation of SNPs and/or candidate genes in these groups. As indicated above, of the ten candidate genes examined here, only four have corresponding but incomplete SNP data in the databases for all the major ethnic groups and may be used in the comparative studies of obesity and/or osteoporosis in the populations of different ethnicity. The SNP data in the databases for the other six genes (60% examined) need substantial updating. To do this, large-scale studies should be performed for the other major ethnic populations. Given that many complex diseases have different rates in different ethnic groups, the extensive volume of the SNP data needs to be updated and validated. This conclusion essentially corresponds to the recently reported results of the SNP databases evaluation regarding their use for whole-genome association studies in humans [21]. Our study provides an example showing the incompleteness of the SNP data in the current databases for studying complex disorders with an ethnic-dependent background.

Another important problem is inconsistency of the data for some SNP markers between either the different databases or databases and literature. It makes selection of a right SNP for a study rather difficult. For example, a SNP is usually considered to be appropriate for association or linkage studies, if its minor allele frequency is ≥0.1 [45,46]. Based on the data of SNP11 from the dbSNP database (Table 3), this SNP is hardly appropriate for the studies in Caucasians, because its minor allele frequency is 0.08. In fact, as determined in the present study, its actual frequency is 0.257 that makes this SNP suitable for the population association and linkage research.

In terms of well-known ethnic differences in incidence of some complex diseases, the discrepancies in the SNP database data may yield wrong conclusion about suitability of particular SNPs for studying genetic basis of these differences. For example, as was mentioned above, SNP26 from the TGF-ß1 gene has the reported minor allele frequency 0.43 (dbSNP) and 0.26 (HGVbase) for Caucasians, 0.27 (TSC) and 0.15 (dbSNP) for African-Americans, and 0.47 (dbSNP) for Asians (Tables 3,4,5). This gene is a candidate for bone mineral density [47] and, thus, may contribute to the ethnic differences in bone mass and osteoporosis. However, from the above data on the SNP allele frequencies, it is impossible to infer whether this particular SNP may be related to these ethnic differences in bone mass, because none of the values significantly differ from any other due to small sample size and, accordingly, limited statistical power.

Large discrepancies in some allele frequency data between the different databases (e.g., SNP4 for Caucasians) are likely due to small sample sizes and respective large sampling errors of the estimates. The sample sizes in the SNP databases are usually smaller than 50. With such a sample size, a relative sampling error of the allele frequency estimates ranges from 0.12 to 0.67 (fig. 1). In this term, the estimates of the SNP allele frequencies in the present study with the sample size of 1,873 are far more reliable than the respective data from the most databases. Estimation of allele frequencies in small samples is notably affected by heterogeneity of the samples. Given that information about an ethnic background of the samples in the SNP databases is usually scarce, the allele frequency estimates in the databases may be significantly biased due to the ethnic admixture.

Figure 1

Relative sampling error of different minor allele frequencies under various sample sizes. The relative sampling error is obtained as , where is the variance of estimated allele frequencies and is computed as [51]. We plotted relative sampling errors of estimated allele frequencies vs. sample sizes.

The small size of the samples in the databases makes it difficult to powerfully test differences in SNP allele frequencies (Δf) between the data from various sources, especially if absolute values of these differences are not large. Fig. 2 illustrates a statistical power of such estimation. Given two samples of size 50 each, the probability to correctly determine Δf = 0.1 is only 27% (fig. 2A). Having a sample size of 1,800 (similar to ours) increases this probability only up to 40% (fig. 2B). Even a 3-fold difference in the absolute values (e.g., SNP11, Table 4) is not statistically significant due to the insufficient sample size in the databases. These examples suggest that some differences in SNP allele frequency data in the public databases may exist due to improper sampling or too small sample sizes [34]. These inconsistencies may be somewhat reduced by pooling the samples. However, this should be done with caution, because pooling data from small samples with a particular genetic background may introduce a large bias to the new values.

Figure 2

Power to detect differences in allele frequencies between two samples at Since the power depends on the sample sizes of the two populations, we fixed one population at a sample size of n, while allowing another population to vary in sample sizes. Δ f is a frequency difference between the sample of given size n and the samples of various sizes. A – n = 50; B – n = 1,800.

Furthermore, in many cases, even after pooling, the sample size may remain too small to gain sufficient reliability of the newly obtained allele frequencies.

The SNP data in literature may somewhat supplement those from the databases. They are sometimes obtained with larger samples and, therefore, have smaller sampling error. However, they are not systematized, which makes their collection and use difficult.

In conclusion, our study demonstrated that, although a large volume of SNP data is available in public databases and literature, a great portion of these data needs comprehensive updating and validating in order to be appropriate for genetic studies of complex disorders, such as obesity and/or osteoporosis. Such large-scale studies of the disease-associated SNPs in various ethnicities may provide important insights into the evolutionary history of human populations as well as in etiology of these diseases.

Methods

Subjects

All the study subjects came from ongoing genetic studies of complex traits that have been approved by the Creighton University Institutional Review Board. All the subjects were Caucasians of western or northern European origin (German, French, Dutch, Swedish), but some had Portuguese and Italian backgrounds. We have recruited 405 nuclear families, each composed of both parents and at least one child. The total sample size was 1,873 subjects, including 840 parents and 1,133 children. All individuals volunteered to participate in the research and signed informed-consent documents before entering the project.

The study SNPs

The candidate SNPs and genes for the present study were selected from publicly available sources based on some of the following criteria: 1) functional relevance and importance for obesity and/or osteoporosis; 2) degree of heterozygosity, i.e., allele frequencies, as reported in literature or databases; 3) position in or around the genes; and 4) their use in previous genetic epidemiology studies. After searching public SNP databases (dbSNP, JSNP, HGVbase, and TSC) and literature, we chose 41 SNPs in 10 genes shown to be associated with obesity and/or osteoporosis. The information about the SNPs is given in Table 1. Among them, SNP27 is an insertion/deletion polymorphism of a cytosine (+/- C) and the others are nucleotide substitutions.

SNP genotyping

DNA was extracted from whole blood using a commercial isolation kit (Gentra Systems, Minneapolis, MN, USA) following the procedure recommended by the manufacturer. The genotyping involved a polymerase chain reaction (PCR) and invader assay (Third Wave Technology, Madison, WI, USA) and was essentially the same for all the SNPs. The PCR mix consisted of 35 ng genomic DNA, 0.2 mM each dNTP, 1 × PCR buffer and 1.5 mM MgCl2, 0.4 µM each primer, and 0.35 units of Taq polymerase (ABI Applied Biosystems, Foster City, CA, USA) in a total 10 µl reaction volume. The PCR primers for all SNPs are given in Table 1. The PCR was performed on PE9700 Thermal Cycler (Perkin Elmer Cetus, Norwalk, CT) using the following profile: 95°C for 5 min, 30 cycles of 94°C for 1 min, 50°C for 1 min, 72°C for 1 min, and then 72°C for 5 min. After the amplification, the product was diluted 1:20 in nuclease-free water. The invader reaction was carried out in 7.5 µl reaction volume containing 3.75 µl diluted PCR product, 1.5 µl probe mix, 1.75 µl Cleavase FRET mix, and 0.5 µl Cleavase enzyme/MgCl2 solution (Third Wave Technology). The reaction mix was overlaid by 15 µl mineral oil and denatured at 95°C for 5 min, and then incubated at 63°C for 20 min on PE9700 Thermal Cycler. After the incubation, the fluorescence intensity for both colors (FAM dye and Red dye) was measured by Cytofluor 4000 (ABI). The data were then loaded to the Invader Analyzer software (Third Wave Technology), and the genotype for every sample was called according to the ratio of the fluorescence intensity of the two dyes.

We initially genotyped all the 41 SNPs in a random sample of 190 to 380 subjects to find SNPs with minor allele frequencies <0.01. Those were then excluded from the further genotyping. Finally, 33 SNPs were genotyped for the whole sample.

Data analysis

PedCheck software [48] was used to verify the accuracy of SNP genotyping in reference to Mendelian inheritance of the alleles within each family. For the 29 polymorphic SNPs the allele frequencies in the nuclear families were estimated by maximum-likelihood method [49] implemented in SOLAR . This method uses all available marker information by accounting the dependence between relatives. Differences in allele frequencies of each SNP among various populations were tested using the χ2 test or Fisher's exact two-tailed test as implemented in Advisor software [50]. SNPs with minor allele frequency less than 0.01 were considered to be monomorphic.

Authors' contributions

VD combined results of the analysis and drafted the manuscript. JRL and DHX performed the literature search, SNP genotyping, and statistical analysis. PYL carried out the statistical analysis. LJZ, HS, YYZ, YJL, SA, and PX performed the SNP genotyping. RRR and HWD conceived of the study, and participated in its design and coordination. All authors read and approved the final manuscript.

Links

Database of Single Nucleotide Polymorphism, dbSNP

Japanese Single Nucleotide Polymorphisms, JSNP

Human Genome Variation Database, HGVbase

The SNP Consortium, TSC