- Open Access
Selecting cases from nuclear families for case-control association analysis
BMC Genetics volume 6, Article number: S105 (2005)
We examine the efficiency of a number of schemes to select cases from nuclear families for case-control association analysis using the Genetic Analysis Workshop 14 simulated dataset. We show that with this simulated dataset comparing all affected siblings with unrelated controls is considerably more powerful than all of the other approaches considered. We find that the test statistic is increased by almost 3-fold compared to the next best sampling schemes of selecting all affected sibs only from families with affected parents (AFaff), one affected sib with most evidence of allele-sharing from each family (SF), and all affected sibs from families with evidence for linkage (AFL). We consider accounting for biological relatedness of samples in the association analysis to maintain the correct type I error. We also discuss the relative efficiencies of increasing the ratio of unrelated cases to controls, methods to confirm associations and issues to consider when applying our conclusions to other complex disease datasets.
Case-control association studies are regaining popularity in the challenge to identify markers conferring susceptibility to complex diseases. A sample of affected cases is compared to a sample of suitable controls to test for association between allelic variants and disease status. In the recent past, family-based association designs were advocated to protect against spurious associations arising from population substructure. However, such designs are 2- to 5-fold less efficient than using unrelated controls . Furthermore, methods such as genomic control and structured association have since been developed to detect and account for population stratification. These methods rely on the premise that stratification would lead to differences in allele frequencies between two or more populations and that these differences could be detected by analyzing anonymous markers [2–4].
Further improvements in power can be obtained by including sibships with multiple affected sibs that are readily available from prior linkage studies . Most of this gain is generally attributable to an increased allele frequency difference between related cases and unrelated controls. When the number of affected relatives increases, the expected allele frequency of the high-risk allele increases in the cases but remains the same in the unrelated controls. In contrast, the frequency of the high-risk allele also increases in the control group when related controls are used.
Where genotyping more than one sibling from a family is cost prohibitive, it may be useful to select the most informative sib for association analysis. A recent study used allele sharing to select the most informative sib from sibships of various sizes and found that choosing the sib showing the greatest allele sharing from each sibship increased the efficiency of case-control associations under a variety of genetic models .
When using related subjects in case-control studies the correlations among relatives must be accounted for in the statistical analysis to avoid an increase in type I error. A number of tests have been proposed that take account of the sampling of biologically related subjects in the variance of test statistic. Risch and Teng  propose a transmission disequilibrium test- (TDT) like statistic for sibling data; Slager and Schaid  advocate an adjusted trend test that allows cousin data to be used as well as sibling data; Bourgain et al.  suggest a quasi-likelihood trend test, particularly when cases are selected from complex inbred pedigrees.
Here we examine the efficiency of a number of strategies for selecting cases from nuclear families with multiple affected subjects and comparing with unrelated controls to identify a known Kofendrerd Personality Disorder (KPD) disease susceptibility marker on a region of chromosome 5. We examine the efficiencies of a number of case selection strategies including those proposed by Fingerlin et al.  and Risch and Teng . The test statistic at the disease locus for each selection scheme is compared with the maximum test statistic we observed, and the number of other associated markers identified is also considered. We discuss the impact of over-sampling controls relative to cases and present approaches for confirming putative associations.
The Genetic Analysis Workshop 14 (GAW14) simulated dataset was used for this analysis. A region of chromosome 5 was known to us to contain a susceptibility locus for KPD and was chosen for investigation. The actual disease locus was originally blinded from those doing the association analysis. However, it became clear from the analysis which marker was the true association and hence results are reported with reference to the known answer. Data packages 206–210 containing 100 markers were used. The Aipotu family dataset (001) with KPD affection status was initially used for the analysis. Unrelated control populations of various sizes (50, 100, 200, 400, and 1,000) were created by randomly combining control replication sets (replicates 001–020). The control data sets for each scheme contained approximately the same number of unrelated controls as affected cases, with different randomly selected controls used in each scheme.
Case selection schemes
Seven case selection strategies were compared. We selected one affected sibling at random from each family (RF), all affected siblings from each family (AF), and all affected sibs from families with one or both parents affected (AFaff). We also selected sibs on the basis of linkage and allele sharing information. Families with evidence for linkage were defined as those with a multipoint linkage NPL score >0 in the chromosome 5 region (centered at D05S0172) (calculated using GENEHUNTER ). We employed schemes using one sibling selected at random from families with evidence for linkage (RFL) and all siblings from those families (AFL). The selection of siblings with the most evidence for allele-sharing was achieved by using the IBD probabilities for each family . The schemes considered used one affected sib from each family with the most evidence for allele sharing (SF) and one affected sib with the most evidence for allele sharing in families with evidence for linkage (SFL). The schemes are summarized in Table 1.
Comparison of selection strategies
Because some of the selection schemes include multiple siblings per family, the test statistic proposed by Risch and Teng  was used to test for single-nucleotide polymorphism (SNP) marker associations in the region, with the adaptation of using the average number of siblings per family. The test statistic is of the form , where and are the estimated allele frequencies of the SNP marker in the cases and controls, and is the estimated standard deviation of the difference.
The test statistic is defined as
The average number of affected sibs per family is denoted r, n is the number of families, and u represents the number of unrelated controls per family. and represent the number of marker alleles carried by case i and control i,, respectively. The test statistic is distributed as a chi-squared with 1 degree of freedom.
The test statistic was calculated for all of the markers in the region using SAS ® (SAS version 8.2, SAS Institute Inc., Cary, NC, USA) and markers yielding p < 0.05 were considered further. To compare the efficiencies of the selection strategies the ratio of the test statistic at the disease locus for the given selection scheme to the maximum test statistic obtained (using all affected sibs and 5 times the number of controls, test statistic = 19.47) was used. The number of other markers in the region showing evidence of association (p < 0.05) to KPD but not identified as the causal variant was also noted. Linkage disequilibrium (LD) across the region was visualized using the program HELIXTREE . The program TRANSMIT  was used to calculate the TDT test for the disease susceptibility marker.
Results and discussion
Identification of disease locus
The efficiencies of the seven selection schemes using equal ratios of cases to controls are shown in Table 1. The most efficient scheme was sampling all affected individuals (n = 233) (AF) and 5 times the number of controls, which gave a test statistic of 19.47. However, using an equal number of controls resulted in a mere 7% reduction in efficiency. The AF scheme was considerably more powerful than the any of the other selection strategies, with a greater than 2.5-fold increase in the test statistic than the next best approaches; AFaff (n = 118), SF (n = 100), and AFL (n = 110), with efficiencies of 33%, 29%, and 27%, respectively. The remaining approaches were particularly lacking in power with efficiencies in the range of 22% to 9%. The number of other positives (p < 0.05) identified was approximately the same across all of the selection schemes.
The higher efficiency of the AF scheme relative to the other schemes may be attributable to the larger sample size alone. However, others have shown that the gain in efficiency of such a design is due to an increase in the disease allele frequency in the case group rather than sample size . In this study, the average number of affected siblings per family was 2.33 (range 2 to 7). It would be interesting to explore the impact of sampling a greater number of siblings per family but we were unable to do this due to the limitations of the simulated dataset. Table 1 shows the risk allele frequency for the case sampling schemes; the AF and AFaff schemes result in the greatest enrichment for the risk allele followed by AFL, SF, SFL, RFL, and RF. This ranking corresponds well with ranking on the difference in the risk allele between cases and controls (differences were observed in the risk allele frequency for each of the control groups as controls were selected randomly for comparing with cases from each of the selection schemes).
Although using all available cases from families may be preferable, from a genotyping perspective, this may not be feasible. When restricted to using only one affected sibling from each family we show that selecting the sib with the greatest allele sharing in a family results in a greater efficiency than randomly selecting an affected sib from each family, as observed by Fingerlin et al. .
Although these findings hold true for the specific genetic and phenotypic models used in the simulated data set, it is not clear how robust these findings are across a range of genetic models (e.g., allele frequency, dominance, epistasis, penetrance, LD between marker and causal variant(s)) and phenotypic traits (e.g., continuous and categorical). However, it is beyond the scope of this article to investigate the impact of these factors on the selection strategies considered here.
Confirmation of the disease locus
The risk locus for KPD on chromosome 5, B05T4136, was identified in the Aipotu population. In addition to the disease locus a number of other markers were found to be statistically associated with KPD. We attempted to strengthen confidence in the association by using a number of methods described below.
We examined the pair-wise marker LD across the region (Figure 1) using D' and r and found it to be low. This was not surprising given an average marker density of 1 per 0.3 cM. The known susceptibility marker was not found to be in strong LD with any of the markers genotyped in the region.
We next considered the potential existence of population stratification. We used a family-based TDT test of association to control for population stratification. The disease locus was found not to be associated, although a trend of increased transmissions of the risk allele was observed. Although the TDT is robust to population stratification, it is somewhat less powerful than case-control sampling . The reduction in power is not only due to the smaller difference in allele frequency between cases and controls, as discussed above, but also because of the lower sample size. The reduction in sample size is attributable to only heterozygous parents being informative in the TDT, resulting in at least a third greater sampling and genotyping being needed compared to case-control sampling.
In addition to attempting to replicate the association in the same population, we also considered replicating in the available Karangar population. This marker was found not to be associated in the Karangar population. There are a number of reasons for not being able to replicate across populations. There could be important differences in allele frequency or LD structure across populations, resulting in the risk allele exhibiting a different pattern of association with marker alleles and haplotypes in the different populations . Hidden population stratification can further complicate this situation by producing spurious association or changing the pattern of a true association [4, 13–16]. Locus and allelic heterogeneity are also possible explanations.
In the GAW14 simulated dataset we have shown that comparing all available cases from nuclear families with unrelated controls in an association study is considerably more powerful than any of the case selection schemes considered. However, strategies using all affected siblings with affected parents or cases with strong allele-sharing result in comparable enrichment of the risk allele but with fewer cases being selected. Although our results are similar to those published by other investigators, we suggest a degree of caution when generalizing all of these findings. Further investigation into the robustness of the results over a range of genetic and phenotypic models is required.
All affected sibs selected from each family
All affected sibs selected from families with affected parents
All affected sibs from families with evidence for linkage
Genetic Analysis Workshop 14
Kofendrerd Personality Disorder
One affected sib selected randomly from each family
One affected sib selected randomly from families with evidence for linkage
Transmission disequilibrium test
One affected sib with most evidence of allele-sharing from each family
One affected sib with most evidence of allele-sharing from families with evidence for linkage
Risch N, Teng J: The relative power of family-based and case-control designs for linkage disequilibrium studies of complex human diseases. I. DNA pooling. Genome Res. 1998, 8: 1273-1288.
Devlin B, Roeder K: Genomic control for association studies. Biometrics. 1999, 55: 997-1004. 10.1111/j.0006-341X.1999.00997.x.
Satten GA, Flanders WD, Yang Q: Accounting for unmeasured population substructure in case-control studies of genetic association using a novel latent-class model. Am J Hum Genet. 2001, 68: 466-477. 10.1086/318195.
Pritchard JK, Stephens M, Rosenberg NA, Donnelly P: Association mapping in structured populations. Am J Hum Genet. 2000, 67: 170-181. 10.1086/302959.
Fingerlin TE, Boehnke M, Abecasis GR: Increasing the power and efficiency of disease-marker case-control association studies through use of allele-sharing information. Am J Hum Genet. 2004, 74: 432-443. 10.1086/381652.
Slager SL, Schaid DJ: Evaluation of candidate genes in case-control studies: a statistical method to account for related subjects. Am J Hum Genet. 2001, 68: 1457-1462. 10.1086/320608.
Bourgain C, Hoffjan S, Nicolae R, Newman D, Steiner L, Walker K, Reynolds R, Ober C, McPeek MS: Novel case-control test in a founder population identifies P-selectin as an atopy-susceptibility locus. Am J Hum Genet. 2003, 73: 612-626. 10.1086/378208.
Kruglyak L, Daly MJ, ReeveDaly MP, Lander ES: Parametric and nonparametric linkage analysis: a unified multipoint approach. Am J Hum Genet. 1996, 58: 1347-1363.
Cardon LR, Bell JI: Association study designs for complex diseases. Nat Rev Genet. 2001, 2: 91-99. 10.1038/35052543.
Neale BM, Sham PC: The future of association studies: gene-based analysis and replication. Am J Hum Genet. 2004, 75: 353-362. 10.1086/423901.
Morton NE, Collins A: Tests and estimates of allelic association in complex inheritance. Proc Natl Acad Sci USA. 1998, 95: 11389-11393. 10.1073/pnas.95.19.11389.
Thomas DC, Witte JS: Population stratification: a problem for case-control studies of candidate-gene associations?. Cancer Epidemiol Biomarkers Prev. 2002, 11: 505-512.
Stumpf MP, Goldstein DB: Demography, recombination hotspot intensity, and the block structure of linkage disequilibrium. Curr Biol. 2003, 13: 1-8. 10.1016/S0960-9822(02)01404-5.
Freedman ML, Reich D, Penney KL, McDonald GJ, Mignault AA, Patterson N, Gabriel SB, Topol EJ, Smoller JW, Pato CN, Pato MT, Petryshen TL, Kolonel LN, Lander ES, Sklar P, Henderson B, Hirschhorn JN, Altshuler D: Assessing the impact of population stratification on genetic association studies. Nat Genet. 2004, 36: 388-393. 10.1038/ng1333.
About this article
Cite this article
Moore, R.M., Pinel, T., Zhao, J.H. et al. Selecting cases from nuclear families for case-control association analysis. BMC Genet 6, S105 (2005) doi:10.1186/1471-2156-6-S1-S105
- Risk Allele
- Population Stratification
- Simulated Dataset
- Transmission Disequilibrium Test
- Allele Sharing