Volume 5 Supplement 9
Linkage analysis merging replicate phenotypes: an application to three quantitative phenotypes in two African samples
© Hinrichs et al; licensee BioMed Central Ltd. 2011
Published: 29 November 2011
We report two approaches for linkage analysis of data consisting of replicate phenotypes. The first approach is specifically designed for the unusual (in human data) replicate structure of the Genetic Analysis Workshop 17 pedigree data. The second approach consists of a standard linkage analysis that, although not specifically tailored to data consisting of replicate genotypes, was envisioned as providing a sounding board against which our novel approach could be assessed. Both approaches are applied to the analysis of three quantitative phenotypes (Q1, Q2, and Q4) in two sets of African families. All analyses were carried out blind to the generating model (i.e., the “answers”). Using both methods, we found numerous significant linkage signals for Q1, although population colocalization was absent for most of these signals. The linkage analysis of Q2 and Q4 failed to reveal any strong linkage signals.
The simulated mini-exome data for Genetic Analysis Workshop 17 were derived from the pilot3 study of the 1000 Genomes Project (http://www.1000genomes.org). The data consist of information on 3,205 genes containing 24,487 single-nucleotide polymorphisms (SNPs) on 697 unrelated individuals drawn from 7 population samples. The family data consist of 8 extended pedigrees (varying in size from 73 to 128 individuals), also with a total of 697 individuals. The pedigrees were seeded with 202 founders, who were described as having been drawn at random from the 697 unrelated individuals. Two hundred replicates were provided for both the unrelated individuals data set and the family data set. A noteworthy feature of these data is that the SNP genotypes from replicate to replicate are identical. In other words, for any given individual over all replicates, each person is represented as though he or she is one of 200 monozygotic twins.
Because of the unusual nature of the data sets, we elected to explore an innovative approach that takes into account the data’s replicate structure. As a complement to this approach, we also carried out an old-fashioned linkage analysis that does not take into account the fact that the families are not genetically independent from one replicate to another.
Because our analyses were blind to the generating model, we began our investigation by obtaining some descriptive statistics. We started with the genotypes of the N = 697 unrelated individuals and quickly discovered that a large proportion (38.5%) of the 24,487 SNPs have a frequency of 1/(2N); that is, only one copy of the variant exists in the entire set of unrelated individuals, and fully 87.2% have a minor allele frequency less than 5.0%. Under usual circumstances, this collection of SNPs would have such a low polymorphic information content as to render ordinary linkage analysis powerless to identify genes involved in the quantitative phenotypes of interest. Fortunately for the family data set, the data simulators provided fully informative identity-by-descent (IBD) matrices for each pair of individuals at each gene. These IBD scores were given as 0, 0.5, and 1, denoting the sharing of 0, 1, or both genes identical by descent. Both methodological approaches reported here used these IBD matrices exclusively for the linkage analyses.
We note that in the process of simulation, recombination could occur only once per chromosome per meiosis and that this recombination could not occur within a gene. This allows us to refer to IBD at a particular gene. In real data, we would not have perfect information, although with extensive genotyping and sequencing we would be able to approach this. Recombination within a gene is also possible, although it is less likely than in the larger intergenic regions. Furthermore, a multipoint analysis is unnecessary because the single point information is completely informative.
Before performing the linkage analyses, we began by making sure that the data had the properties we expected. We believe that this is an important but often neglected step for both real and simulated data. While attempting to identify the population origin of the 202 founders of the 8 pedigrees, we discovered to our surprise that 4,239 SNPs had been altered, although these alterations did not prevent the unambiguous identity and population origin of the founders. We were even more surprised to discover that the founders were not a random sample from the unrelated individuals data set. Instead, all but three of the founders were drawn from the two African samples. Moreover, three of the eight pedigrees were established by approximately equal numbers of founders from the Yoruba and the Luhya population samples. The remaining five pedigrees were primarily founded by one ethnic group or the other. Under ordinary circumstances when using genetic markers for a linkage analysis, it is more cumbersome to analyze families composed of individuals drawn from various ethnic groups because differences in allele frequencies need to be taken into account (see Suarez et al.  for an example of combining families of different ancestries in a linkage analysis). Using IBD matrices somewhat obviates this problem, although the possibility looms large that different ancestral populations may be segregating different causative genes.
Accordingly, we decided to concentrate our attention on two groups of kindreds: (A) pedigree 2 (deleting founder 2005 [a Tuscan] and all of her descendants), pedigree 3, and pedigree 4 (deleting founder 4021 [a Yoruba] and all of her descendants); and (B) pedigrees 5 and 6. Individuals in kindred group A are all sampled from the Luhya population, a Bantu ethnic group living primarily in Kenya. It is unclear whether the Luhya sample from the 1000 Genomes Project was randomly sampled from all of the 16 tribes that make up the Luhya ethnic group or from a smaller subset. Individuals in kindred group B are all Yoruba, a large heterogeneous group of about 30 million people from West Africa, predominantly Nigeria.
as the quantitative phenotype, where and , the maximum-likelihood estimators (MLEs) of the parameters α and β, are computed from the data and G represents the sum of the genetic contributions to the phenotype.
In this way, each unique genetic observation is counted only once in the analysis, but full information from the repeated measures is used.
For the standard linkage analysis (which does not take into account the repeated measures feature of these data), we used SOLAR , and the p-values were combined over multiple replicates using Fisher’s method . Smoking, Sex, and Age were used as covariates, although the nonsignificant covariates were not in the final model. We chose not to include Population among the covariates because our analyses were performed separately for the Luhya pedigrees and the Yoruba pedigrees.
True positives and false positives from Fisher’s method and GLMM analyses
A striking result of our analyses is the lack of colocalization of linkage signals between the two African groups. Moreover, the maximum signal is often a plateau that occurs over a region of a number of contiguous genes. This appears to be an artifact of the simulation, which allowed only a single recombination per chromosome for each meiosis. Consequently, the position of the recombinations is constant for all of the replicates.
Figures 2 and 3 indicate that the linkage signals obtained with Fisher’s method tend to be higher than those obtained using the GLMM, which, we speculate, is due to the fact that Fisher’s method treats the identical genotypic structure of the 50 replicates as independent observations.
Discussion and conclusions
The most striking feature revealed by our linkage analysis of the quantitative trait Q1 is that, with one possible exception, none of the five signals seen in the Luhya sample colocalize with any of the four signals found in the Yoruba sample. The exception is the signal seen on chromosome 8 at CDCA2, although only the Yoruba signal is highly significant. Had we analyzed each of the three Luhya pedigrees and each of the two Yoruba pedigrees separately, we would have found family-specific signals, all of which would have been stronger than those reported in Table 1 because of our pooling of pedigrees from the two ethnic groups. Clearly, for private polymorphisms or in cases of locus heterogeneity (see, e.g., Morton  for the first published example of locus heterogeneity documented by a linkage analysis), linkage analysis of individual families often is preferable to combining all families from a single ethnic group, as we originally did. Indeed, Culverhouse et al. , analyzing the unrelated individuals data set, also found that pooling populations could obscure a strong signal even when Population was included as a covariate.
Our Eigenstrat analysis of just the two African samples from the unrelated individuals file (Figure 1) suggests much greater genetic similarity among the Luhya sample than among the Yoruba sample. One possible reason for the apparent heterogeneity of the Yoruba sample is their ancestors’ involvement in the slave trade. Indeed, the Yoruba were one of the largest African groups enslaved and taken by Europeans to the New World, especially Cuba.
It would be naïve to expect any single method to pinpoint all the genetic factors that were simulated, especially for private polymorphisms (i.e., polymorphisms found only in a single family or ethnic group). For instance, Hill et al.  recently were able to localize a linkage signal for pleuropulmonary blastoma on the basis of a few small families. However, it took a different method (sequencing) to identify the lesions in the responsible gene (DICER1).
The linkage results for Q2 and Q4 (not shown in the figures) did not produce any compelling evidence for the involvement of a well-demarked chromosomal region for these phenotypes. Fisher’s method identified a modest true positive for Q2 on chromosome 8 at PLAT and a modest false positive for Q4 on chromosome 20 at LOC643659. Similarly, the GLMM identified modest true positives for Q2 at SREBF1 on chromosome 17 and at PLAT and a modest false positive for Q4 on chromosome 5 at RNF145.
Other analyses by our group  indicate that the four covariates (Sex, Age, Smoking status, and Population) are highly significant for Q4 and explain a large proportion of the variance of Q4. (Over all 200 unrelated replicates, the mean R2 is 0.787 with a range of 0.759–0.811.) We concluded that Q4 was primarily environmentally determined but that the generating model did include unlinked heritable factors.
Although we correctly located most of the simulated signals for Q1, we also found a number of false-positive signals (Table 1). Two of these (the regions on chromosomes 7 and 9) attained p-values of approximately 10−4, which normally would not be considered strong evidence of linkage. The other two signals (on chromosomes 8 and 17), however, would be considered persuasive evidence of linkage if the families were independent.
Finally, we note that the present situation does not reflect any situation observed in human genetics. However, these analyses are appropriate for plant and animal genetics. In particular, homozygous inbred strains and their F1 hybrids offer investigators the ability to observe the characteristics of hundreds of genetically identical organisms in a variety of contexts. Furthermore, the analyses of monozygotic twin data and longitudinal phenotype data are also amenable to analysis with generalized linear models.
The Genetic Analysis Workshops are supported by National Institutes of Health (NIH) grant R01 GM031575 from the National Institute of General Medical Sciences. Additional support was obtained from the Urological Research Foundation and from NIH grants K01 AA015572, K25 GM069590, R03 DA023166, and IRG-58-010-50 from the American Cancer Society.
This article has been published as part of BMC Proceedings Volume 5 Supplement 9, 2011: Genetic Analysis Workshop 17. The full contents of the supplement are available online at http://www.biomedcentral.com/1753-6561/5?issue=S9.
- Suarez BK, Duan J, Sanders AR, Hinrichs AL, Jin CH, Hou C, Buccola NG, Hale N, Weilbaecher AN, Nertney DA, et al: Genomewide linkage scan of 409 European-ancestry and African American families with schizophrenia: suggestive evidence of linkage at 8p23.3–p21.2 and 11p13.1–q14.1 in the combined sample. Am J Hum Genet. 2006, 78: 315-333. 10.1086/500272.PubMed CentralView ArticlePubMedGoogle Scholar
- SAS Institute Inc: SAS 9.2. Edited by: Cary, NC. 2008, SAS Institute IncGoogle Scholar
- Almasy L, Blangero J: Multipoint quantitative trait linkage analysis in general pedigrees. Am J Hum Genet. 1998, 62: 1198-1211. 10.1086/301844.PubMed CentralView ArticlePubMedGoogle Scholar
- Province MA: The significance of not finding a gene. Am J Hum Genet. 2001, 69: 660-663. 10.1086/323316.PubMed CentralView ArticlePubMedGoogle Scholar
- Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D: Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006, 38: 904-909. 10.1038/ng1847.View ArticlePubMedGoogle Scholar
- Entine J: Taboo: why black athletes dominate sports and why we are afraid to talk about it. 2000, New York, Public AffairsGoogle Scholar
- Morton NE: The detection and estimation of linkage between the genes for elliptocytosis and the Rh blood type. Am J Hum Genet. 1956, 8: 80-96.PubMed CentralPubMedGoogle Scholar
- Culverhouse R, Hinrichs AL, Suarez BK: Stratify or adjust? Dealing with multiple populations when evaluating rare variants. BMC Proc. 2011, 5 (suppl 9): S101-10.1186/1753-6561-5-S9-S101.PubMed CentralView ArticlePubMedGoogle Scholar
- Hill DA, Ivanovich J, Priest JR, Gurnett CA, Dehner LP, Desruisseau D, Jarzembowski JA, Wikenheiser-Brokamp KA, Suarez BK, Whelan AJ, et al: DICER1 mutations in familial pleuropulmonary blastoma. Science. 2009, 325: 965-10.1126/science.1174334.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.