 Methodology article
 Open Access
 Published:
Leveraging local ancestry to detect genegene interactions in genomewide data
BMC Genetics volume 16, Article number: 124 (2015)
Abstract
Background
Although genomewide association studies have successfully identified thousands of variants associated to complex traits, these variants only explain a small amount of the entire heritability of the trait. Genegene interactions have been proposed as a source to explain a significant percentage of the missing heritability. However, detecting genegene interactions has proven to be very difficult due to computational and statistical challenges. The vast number of possible interactions that can be tested induces very stringent multiple hypotheses corrections that limit the power of detection. These issues have been mostly highlighted for the identification of pairwise effects and are even more challenging when addressing higher order interaction effects. In this work we explore the use of local ancestry in recently admixed individuals to find signals of genegene interaction on human traits and diseases.
Results
We introduce statistical methods that leverage the correlation between local ancestry and the hidden unknown causal variants to find distant genegene interactions. We show that the power of this test increases with the number of causal variants per locus and the degree of differentiation of these variants between the ancestral populations. Overall, our simulations confirm that local ancestry can be used to detect genegene interactions, solving the computational bottleneck. When compared to a single nucleotide polymorphism (SNP)based interaction screening of the same sample size, the power of our test was lower on all settings we considered. However, accounting for the dramatic increase in sample size that can be achieve when genotyping only a set of ancestry informative markers instead of the whole genome, we observe substantial gain in power in several scenarios.
Conclusion
Local ancestrybased interaction tests offer a new path to the detection of genegene interaction effects. It would be particularly useful in scenarios where multiple differentiated variants at the interacting loci act in a synergistic manner.
Background
Advances in highthroughput genotyping technologies have enabled largescale studies of genetic variation, from genomewide association studies (GWAS) to inference of population history. The most notable use of highthroughput genotyping has been in GWAS where researchers have reproducibly identified thousands of genetic variants associated with many complex traits and common diseases. Despite the great success in identifying variants that contribute risk to disease, the majority of the genetic component of human traits and diseases remains unexplained. A potential source for this missing heritability is genegene interactions that alter disease risk in a coordinated fashion, for example when several genes are acting synergistically on a trait. Although of potential great interest, robust identification of genegene interactions has largely remained elusive, and despite numerous studies only a few interaction effects have been detected in human data [1–3]. Most genetic association studies of genegene interaction have focused on the joint effect on pairs of single nucleotide polymorphisms (SNPs) and used brute force approaches to evaluate a large number of pairs on homogenous populations (e.g. individuals of European ancestry only), while alternative strategies using heterogeneous populations have been seldom considered [4, 5].
The development of accurate methods for discerning population structure have allowed for studies across different ethnicities including admixed populations (i.e. populations with recent ancestry from more than one continent such as African Americans). In addition to the standard linkage disequilibrium (LD) between nearby markers (used by GWAS to tag hidden causal variants) admixed populations exhibit another form of correlation among variants at a coarser scale due to chromosomal segments of distinct ancestry that is commonly referred to as admixtureLD [6]. This enables admixture mapping to be an effective approach for identifying disease loci that differ in frequency across populations [7–11]. A key component of such studies is the inference of ancestry at each locus in the genome. Several computational and statistical tools, including HAPMIX [12], LAMPLD [13], EILA [14], and LANCCSVs [15] can now be used to reliably call local ancestry. Although local ancestry has been traditionally used in admixture mapping, recent works use analyses of local ancestry to yield novel insights into the dynamics of recombination rate across the genome, to make demographic inferences from genetic data of admixed populations, as well as to understand the genetic basis of complex traits [16–18].
In this work we explore the use of local ancestry in recently admixed populations to find signals of genegene interaction that affect disease risk. We introduce an approach that leverages the correlation between local ancestry and the hidden unknown causal variants to find distant –e.g. on different chromosomes– genegene interactions. Our proposed approach uses multiple linear regression to model the interaction effect between pairs of local ancestry segments. Hence, as opposed to the standard approaches that test all pairs of SNPs assayed in GWAS (e.g. on the order of 10^{12} pairs for a standard GWAS of 2.5 million SNPs), we propose to test interaction only between pairs of local ancestries (on the order of 5×10^{5} pairs for recent admixtures). By performing a much smaller number of statistical tests, our approach solves the computational bottleneck and reduces the multiple testing correction burden. We derive the analytical formulation for our test assuming a single causal variant for each interacting locus and investigate its performance across a wide range of parameter values. Motivated by recent works that show everincreasing evidence for multiple causal variants per locus [19–21], we extend our approach to allow for multiple causal variants at each interacting locus. We find that local ancestry can be used to find genegene interactions, with power increasing with the number of causal variants per locus and the degree of differentiation in the frequency of the causal variants between the ancestral populations. Assuming equal sample size, the test based on pairwise genotyped SNPs appears to be more powerful than the ancestrybased interaction test in most scenarios. However, when accounting for the increase in sample size that can be achieved for a fixed budget when measuring local ancestry only (e.g. based on ancestry informative markers, AIMs), we observed a substantial increase in power under various scenarios.
Results and discussion
Overview of the approach
A standard approach for finding pairwise SNP interactions is to test for nonzero effect size of the product term of the two SNPs considered. The underlying assumption is that the SNPs tested in the model are either the interacting causal variants or correlated to the actual causal variants through LD. Indeed, only a finite number of SNPs are assayed in GWAS (today’s genotyping arrays assay a few million SNPs), with true biologically causal variants likely remaining untyped. While a number of additional SNPs can be imputed on a genomewide scale, the presence of the causal variants in the data can only be assumed for wholegenome sequence data. It is likely the causal variants will only be tagged by the SNPs analyzed. In admixed populations, correlation between SNPs also exists at a coarse level due to the segments of recent ancestry (admixtureLD). Similar to the pairwise SNP interaction screening, we can tag the hidden causal variants using admixtureLD and test for the presence of interaction at hidden causal variants by testing for interaction at the level of local ancestry.
Testing for interaction under a single causal variant per locus assumption
We first considered a scenario where two common SNPs are located on two physically distant segments in the genome, thus independent from each other, and have an interactive effect on a quantitative phenotype, while all other SNPs at the locus harboring these two causals have no effect (Fig. 1a). We derived the performances of three interaction tests, based on full sequence data (S_{S}), genotypic data from a 1 M (1 million) SNPs chip (S_{G}) and local ancestry (S_{L}). Figure 2 presents the sample size required for each of the three interaction tests to achieve a significance level of 5 % with 80 % power after correction for multiple testing. The sample sizes are plotted for a range of correlation levels between the causal variant and the tagging SNP or tagging ancestry segment. More specifically, we refer to ρ _{ GC } for the correlation between the true interaction term and the interaction term derived from the best tag from the genotyping chip, and to ρ _{ LC } for the correlation between the true interaction term and the local ancestry interaction term between the two segments harboring the causals. For simplicity we assumed here that local ancestry is inferred with high accuracy (r ^{2} between true ancestry and inferred ancestry ≥0.99), and therefore does not differ from true local ancestry. This will likely be the case for African Americans [15], but might be too optimistic for other populations such as Latino Americans (see below). Figure 2 shows that for GWASbased test (S _{ G }) to outperform the test based on the true causal variants (S _{ S }), it requires ρ _{ GC } to be above 0.9. For the local ancestry test, ρ _{ LC } has to be above 0.7. Moreover, for ρ _{ LC } above 0.8, the ancestry based interaction test would also outperform the GWASbased interaction test even if the causal variants were genotyped. As expected, the maximum potential gain is achieved when the interaction term is perfectly tagged by either GWAS SNPs or local ancestry. We considered unrealistically large interaction effects in Fig. 2 for illustration purposes. When analyzing 20,000 samples, the smallest interaction effect (as measured by the proportion of variance explained) that can be detected with 80 % power is 0.8, 0.5 and 0.3 % for S _{ S }, S _{ G } and S _{ L }, respectively.
We then estimated the empirical distribution of ρ _{ GC } and ρ _{ LC } using AfricanAmerican individuals simulated using the 1,000 Genomes data (see Methods) [22]. From this simulation we randomly choose 20,000 independent SNPs, and built 10,000,000 hypothetical pairs of interacting SNPs. Bottom panels of Fig. 2 shows the distribution of these two correlation terms when using tagging SNPs from the 1 M Illumina chip and the simulated local ancestry. Despite the large potential increase in power shown in the upper panels of Fig. 2, improvement may actually exists only in very few real situations. For example we observed that the probability of ρ _{ GC } to be above the 0.9 threshold is 0.05. For ρ _{ LC } the “increased power” threshold of 0.7 is achieved only once in 10^{7} times. Hence, even if interaction effects are extremely common in the architecture of complex trait, there is a low probability for the local ancestry–based test to do better than other approaches in the presence of a single causal variant per locus when assuming equal sample size for both tests.
Multiple causal SNPs per segment
Accounting for increasing evidence of multiple causal variants per locus [19–21], we then considered scenarios where genegene interaction effects involved multiple genetic variants per locus. For example when multiple SNPs contribute to gene transcript abundance and the interaction is taking place between the gene products. Such interaction would be challenging to identify using SNP data due to the vast search area among all possible combinations of SNPs. On the other hand, local ancestry offers a more appropriate and natural way to test for such models as it captures a form of an individual’s genetic background at each locus (i.e. genetic variants share the same local ancestry at a given locus in the genome).
To evaluate this assumption we defined a simulation model where multiple SNPs at two independent loci contribute to two latent variables Z _{1} and Z _{2} that have an interaction effect on the outcome (Fig. 1b). The power of the pairwise SNPs test depends on the best tagging SNPs for Z _{1} at locus 1 and for Z _{2} at locus 2. This would be either the strongest causal variants for Z _{1} or Z _{2}, or the best tag of these causals. The power the local ancestrybased test to detect this interaction depends on all parameters influencing \( {\rho}_{L_i{Z}_i} \), the correlation between Z _{ i } and L _{ i }, the latent variable and the local ancestry at locus i, respectively. This includes the number of causal SNPs for Z _{ i }, and the distribution of β _{ i }, the effects of the causal SNPs of Z _{ i }. Assuming the causal SNPs are a random sample of the variants in the segment, power is also bounded by the average difference in minor allele frequency between the two founder populations. Figure 3 presents the empirical distribution of this correlation in a simple scenario, when Z _{ i } depends on 1 to 50 SNPs. Overall, \( {\rho}_{L_i{Z}_i} \), increases with the number of SNPs involved and with increased homogeneity of genetic effect. For example if the β _{ i } = (β _{ i1}, … β _{ iK }) are distributed around the null, the expected value of \( {\rho}_{L_i{Z}_i} \) is null and a local ancestrybased test would have no power. Conversely, if the coded alleles from the causal variants tend to increase the outcome value (while the reference allele has no contribution), \( {\rho}_{L_i{Z}_i} \) can be substantial (e.g. >0.2, Fig. 3b).
We performed a simulation study to compare the performance of the pairwise SNPbased approach (S _{ G }), when using both genotyped and imputed common SNPs (MAF > 1 %), and the local ancestry approach (S _{ L }) while increasing K the number of causal SNPs per locus. For simplicity we assumed the number of causal variants was the same in the two interacting regions, and only considered common variants (minor allele frequency, MAF > 0.10). We explored scenarios where the causal SNPs were either slightly differentiated or highly differentiated between the two ancestral populations. When assuming equal sample size, S _{ L } is underpowered as compared to S _{ G }, despite a dramatic increase in the total number of tests performed (Fig. 4). Hence when GWAS data exists, deriving local ancestry segments would be of limited interest for genegene interaction testing unless the number of causal variants is large (e.g. >10). However, when considering de novo genotyping with a fixed budget, an increase in sample size can be achieved when measuring local ancestry only, S _{ L } can be more powerful than S _{ G }. In particular, assuming a 6 fold decrease in cost, S _{ L } outperform S _{ G } if either the differentiation is moderate or the number of causals is large (>5), or if the causal SNPs are highly differentiated (e.g. correlation between local ancestry and the causal >0.5, Fig. 4c).
We assumed in these simulations that true local ancestry is available. To evaluate the impact of additional noise introduced by ancestry inference, we analyzed the same data but using ancestry inferred using LAMPLD. As shown in Additional file 1: Figure S1, using the inferred ancestry has only minor impact on power. This is partly expected thanks to the high accuracy of inference in African Americans [15]. The accuracy of SNP imputation could also impact the power of the SNPbased test, however the quality of imputation depends on more parameters, varying across allele frequencies and the chips used for genotyping [23]; it would therefore be more difficult to evaluate thoroughly. Instead we applied the SNPbased test using the genotyped SNPs only found on the Illumina Human1MDuo BeadChip (Additional file 1: Figure S1). We observed a substantial decrease in power, highlighting both the importance of using imputed variants and the need for high quality imputation for the SNPbased test, which might be a concern for rare causal variants.
Finally, we evaluated how the relative power of the S _{ L } and S _{ G } tests is impacted when applying a twosteps procedure where SNPs and ancestry segments are first preselected for interaction testing based on their marginal effects [24]. For simplicity we assume that the vast majority of SNPs are not involved in interaction, so that for a pvalue threshold t at step 1, the total number of interaction tests at step 2 can be approximated by \( \left(\begin{array}{c}\hfill n\times t\hfill \\ {}\hfill 2\hfill \end{array}\right) \), where n is the number of predictor (either SNP or local ancestry segment). Additional file 1: Figure S2 shows the results from this strategy when applied to a case similar to Fig. 4, but adding a main effect to each causal SNPs of the same magnitude as the interaction effect, and using t = 0.01 at step 1. In this specific scenario, the 2step approach mostly benefits to the ancestrybased test, which outperform the SNPbased test in many more scenarios, including cases where sample size was the same for the 2 analyses.
Conclusions
We explored the performance of a local ancestrybased interaction test to capture nonlinear effects from two independent loci. The strategy is similar to a standard SNPbased pairwise interaction screening but uses local ancestry segments instead of SNPs. One major underlying motivation for such an approach is that the total number of tests to be conducted is dramatically lower than for a standard pairwise SNP interaction test, reducing both the computational burden and the correction for multiple testing. We demonstrate that such a test would indeed capture interaction effects between two loci as long as the individual effects of the causal variants at each locus do not cancel each other. For existing datasets that only contain local ancestry data derived from AIMs and for de novo genotyping studies looking for the optimal cost/power ratio, our approach (S _{ L }) can be highly relevant as it can outperform the pairwise SNP screening from standard GWAS data (S _{ G }). Conversely, when GWAS genotyping data does exists, in most scenarios we explored, S _{ G } outperforms S _{ L } when the number of causal variants at the locus was small. As the number of causals grows beyond 10, the power of S _{ L } increases but does not substantially exceed S _{ G } unless the differentiation of the causal SNPs between the two populations is very high. Interestingly, as the differentiation increases so does the relative power of S _{ G }, which explains the underperformances of S _{ L }. We found that, as differentiation increases, many genetic variants become good tags for local ancestry, and so S _{ G } benefits from the increase in differentiation as well. Overall, the relative performances of our approach depends on the balance between the gain in power achieved thanks to the decreasing number of test and the decrease in power due to low correlation between local ancestry and the causal variant(s) at the interacting loci.
Furthermore, we used the whole genome sequencebased test (S _{ S }) as a reference to compare the relative performance of the two alternative approaches. While such a test might have higher power than S _{ L } and S _{ G } (Fig. 2), testing all possible pairs of SNPs would requires extremely intensive computational power in practice, and the implementation of such tests, which have been rarely explored to our knowledge, would require substantial software development and hardware structures (e.g. graphics processing units [25]). This confirms that, as of today, GWASbased pairwise interaction tests remain a relevant approach for identifying interactions as compared to whole genome sequencebased approaches.
Regarding power comparison between S _{ L } and S _{ G }, the power of the localancestry based interaction test was derived based on 1,000 local ancestry regions. However, the number of segments depends on the number of ancestral populations and the number of generations since admixture, and will therefore differ across admixed populations. Increase in the total number of segments can impact the correlation between local ancestry and the causal variants within these segments, as well as the total number of tests that have to be performed for an interaction screening. Using the inferred local ancestry had very limited impact on power in our simulation as the accuracy of inference is very high in African Americans. However, for other populations, the impact might be substantial. For example Brown et al. reported squared correlation between true and inferred local ancestry of 0.63, and 0.81 for Mexican and Puerto Rican population when using LANCCSV, which had similar results to other methods [15]. While further analysis might explore such situations, we believe the results described in this study would remain valid. Finally, additional work might also include extensive explorations of scenarios where interaction factors (either single SNPs or single local ancestry segments) are selected based on their marginal association with the phenotype of interest. When applied in our simulation framework we observed a strong improvement of the ancestrybased test over the SNPbased test, however this needs to be confirmed across a broader range of scenarios.
Overall, while our approach shows some limitations when genomewide genotypic data are available and when the number of causal variants per region is small and contains mostly undifferentiated variants, we highlight that genomewide localancestry based interaction screening remains relevant. First, because some datasets only generated local ancestry data through AIMs, and do not have GWAS data. Second, considering budget constraint for de novo genotyping, and assuming a 6 fold decrease in cost for genotyping AIMs as compared to a standard GWAS, substantial additional gain in power can be achieved through localancestry based tests.
Methods
Genetic model
We considered a genetic model similar to the one described in Chatterjee et al. [26], which can be easily adapted to the local ancestry context. It consists in two independent sets of adjacent SNPs from two loci on different genomic regions, which represent in this study two local ancestry segments. Several SNPs within segment i have an indirect association with the outcome of interest through a latent variable Z _{ i }, an unmeasured quantitative biological phenotype partially governed by SNPs within the locus. Interaction effects between the genetic variants on the outcome is introduced through an interaction term between the Z _{1} and Z _{2} variables (i.e. the cumulative effect of the genetic variants within a locus depends on the cumulative effect of the variants in a distant locus). More specifically the outcome Y is defined as follows:
where Z _{1} and Z _{2} are the two latent variables that each depend on K SNPs and θ _{1}, θ _{2}, θ _{12} respectively represent the main effect of Z _{1}, Z _{2}, and the interaction of Z _{1} and Z _{2} on Y; ε is the residual noise and is normally distributed with mean 0. The latent variables Z _{ i } are defined as follows:
where G _{ ik } and β _{ ik } is the standardized genotype of SNP k in locus i and its main effect respectively. The SNP effects β _{ i } = (β _{ i1}, … β _{ iK }) were randomly drawn from either normal or lefttruncated normal depending on the scenario explored. For simplicity, as the main effect of the latent variable has no impact on the interaction test [27], we set θ _{1} and θ _{2} to 0, and set θ _{12} to 1. Except when specified otherwise we scaled the variance of ε so that the proportion of the variance of Y explained by the interaction equals 1 % (to reflect values observed in GWAS for common complex traits).
A single major causal SNP per segment
We first assumed the genetic effect of a segment is driven by a single causal variant, while the effect of other potential SNPs is null or negligible. This is equivalent to assuming K = 1 in equation (2). The effect of the SNPs on Y can then be rewritten:
Assume G _{1} and G _{2} are tagged by L _{1} and L _{2} respectively, where L _{ i } is the local ancestry measured at the segment harboring SNP G _{i}. Note that G _{1} and G _{2} do not necessarily need to be typed to correctly identify local ancestry (local ancestry spans many MB’s in recently admixed populations and can be reliably identified using a small set of variants). For simplicity, we considered only the case of a twoway admixed population, so that local ancestry would be typically coded as an ordinal variable with value corresponding to the number of chromosomes harboring a particular ancestry. Hence, for AfricanAmerican, L _{.} equals 0, 1 or 2. When the population under study is an admixture of more than two ancestries, testing for interaction would be more complex because of additional combinations of ancestries (e.g. for 3 ancestries A, B and C, an individual can have any of the six following ancestries at a given segment: AA, AB, AC, BB, BC and CC). To our knowledge there is no established standard to handle such situations, the simplest solution consists of testing one ancestry versus the rest [17].
We compared the relative performances of the standard test of interaction of β _{ intG } (equation (3)) versus the test of β _{ intL } (equation (4)), the interaction effect observed between L _{1} and L _{2} on Y, which can be obtained from the model:
Note that both equation (3) and (4) are standard 1 ° of freedom tests of interaction effect estimates obtained from multiple linear regression (see our R script example in Additional file 2). The only difference being that the later test uses local ancestry instead of genotyped SNPs. When the causal SNPs are available (e.g. from whole genome sequence data) and have been standardized, the Wald test of β _{ intG } is defined as \( {S}_S={\left({\widehat{\beta}}_{intG}/{\widehat{\sigma}}_{\beta_{intG}}\right)}^2 \). Similarly the test of β _{ intL } is defined as \( {S}_L={\left({\widehat{\beta}}_{intL}/{\widehat{\sigma}}_{\beta_{intL}}\right)}^2 \). Under the null hypothesis of no association both S _{ S } and S _{ L } follow a chisquare distribution with one degree of freedom. Let ρ denote the correlation between two variables. The two scores can be written as:
where ρ _{ LC } ^{2} = ρ ^{2}(G _{1} × G _{2}, L _{1} × L _{2}) is the squaredcorrelation between the true interaction G _{1} × G _{2} term and the local ancestry interaction term L _{1} × L _{2}.
In another scenario, one would test for interaction effects between pairs of SNPs from a standard GWAS chip, which implies that the tested SNPs, say G _{1} ^{*} and G _{2} ^{*} are tagging the two causals, but (likely) at a higher level than the local ancestry. We denote this test S _{G}:
where ρ _{ GC } ^{2} = ρ ^{2}(G _{1} × G _{2}, G _{1} ^{*} × G _{2} ^{*} ) is the squaredcorrelation between the true interaction G _{1} × G _{2} term and the GWAS interaction term G _{1} ^{*} × G _{2} ^{*} . The power of the three tests can be derived as:
where F(χ ^{2}d, S _{.}) is the cumulative probability function of the noncentral chisquare distribution with d degrees of freedom and noncentrality parameter S _{.}; χ _{ d,p,0} ^{2} is the inverse of F under the null, i.e. the quantiles of the noncentral chisquare distribution, and α is the type I error rate. The relative performances of the three strategies can then be evaluated by comparing the sample size N needed to identify the interaction at the 5 % significance level (α) after accounting for n _{ S }, n _{ L }, and n _{ G } the number of tests that has to be performed respectively in the whole genome setting, the local ancestry setting and the GWAS setting, respectively. Therefore, the alpha levels for tests S _{ S } S _{ L } and S _{ G } were set at α/n _{ S }, α/n _{ L }, and α/n _{ G }, respectively.
For the test S _{ S }, we assumed the true causal variants are available as part of a whole genome sequence data of n _{ S } = 20M SNPs, so that the total number of pairwise test equals 2x10^{−14}. We assumed a total of n _{ L } = 1, 000 local ancestry segments for the test S _{ L }, and n _{ G } = 1M SNPs for test S _{G}. We considered values of ρ _{ LC } and ρ _{ GC } in the range [0.3; 1], so that the minimum squaredcorrelation was 0.09. A correlation of 1 corresponding to the highest potential increase in power that can be achieve, since it will be equivalent to testing the interaction with the true causals while dramatically reducing the multiple testing corrections.
Multiple causal SNPs per segment
We used models from equation (1) and (2) and generated outcome data across 1,000 replicates while using simulated genetic and local ancestry data (see further section on simulation). For ease of computation, we considered the number of causal SNPs per segment to be equal between the two interacting loci. The β _{ i } were randomly drawn from a lefttruncated normal (cut at 0) distribution with mean 0 and variance 1. The causal alleles were chosen to be minor in one of the populations so that all effects go in the same direction in one population. We considered scenarios including 1 to 5 causal SNPs per locus, and selected the variants so that correlation between local ancestry and each SNP was either >0.1, >0.2 or >0.5, assuming low, moderate or high differentiation at the causal SNPs, respectively.
We perform both the local ancestry interaction and the pairwise SNP interaction test in a standard linear regression. However, instead of testing all combination of SNP, we first find the best combined tag in each locus (i.e. the single SNP j that maximizes \( {\displaystyle \sum_K}\left({\beta}_{ik}{\rho}_{jk}\right) \) across the K causal variants at locus i, where ρ _{ jk } is the correlation between the tag SNP G _{ j } ^{*} and the causal variant G _{ k }) and then test only the product of that genotype with the best tag genotype at the second locus. Power was defined as the number of replicates for which the pairwise SNP interaction is significant at \( p=0.05/\left(\begin{array}{c}\hfill {n}_G\hfill \\ {}\hfill 2\hfill \end{array}\right) \), where n _{ G } the number of SNPs tested equals 1 M when including genotyped SNPs only and equals 10 M when including common imputed SNPs from 1000 Genomes. For the local ancestrybased test, the threshold was \( p=0.05/\left(\begin{array}{c}\hfill {n}_L\hfill \\ {}\hfill 2\hfill \end{array}\right)=1\times {10}^{7} \), corresponding to the test of n _{ L } = 1, 000 local ancestry segments.
Simulation of admixed populations from the 1000 genome project
Similar to previous work [15], we simulated admixed chromosomes of AfricanAmericans as a random walk over 1000 Genomes haplotypes [22]. CEU and FIN populations were used to represent European haplotypes and the YRI population represented the African haplotypes. We assumed that between any two base pairs there was a 10^{−8} chance of recombination. At a recombination event, the next haplotype was selected with a 20 % chance of being European and an 80 % chance of being African to reflect the estimated admixture proportions in the literature [12, 28]. Haplotypes were sampled with replacement. We simulated 20,000 haplotypes each for chromosomes 19 and 20 in this manner and added them together to form 10,000 unphased genotypes and true local ancestries. Local ancestry was then inferred using LAMPLD [13] with default settings and with the GBR and TSI populations representing Europeans and the LWK representing Africans. LAMPLD was run only using variants found on the Human1MDuo BeadChip.
Availability of supporting data
All supporting data are included as additional files.
Abbreviations
 SNP:

Single nucleotide polymorphisms
 GWAS:

Genomewide association study
 LD:

Linkage disequilibrium
 AIM:

Ancestry informative marker
 MAF:

Minor allele frequency
References
 1.
Aschard H, Lutz S, Maus B, Duell EJ, Fingerlin TE, Chatterjee N, et al. Challenges and opportunities in genomewide environmental interaction (GWEI) studies. Hum Genet. 2012;131(10):1591–613.
 2.
Cordell HJ. Detecting genegene interactions that underlie human diseases. Nat Rev Genet. 2009;10(6):392–404.
 3.
Mackay TF. Epistasis and quantitative traits: using model organisms to study genegene interactions. Nat Rev Genet. 2014;15(1):22–33.
 4.
Aldrich MC, Kumar R, Colangelo LA, Williams LK, Sen S, Kritchevsky SB, et al. Genetic ancestrysmoking interactions and lung function in African Americans: a cohort study. PLoS One. 2012;7(6):e39541.
 5.
Powell R, Davidson D, Divers J, Manichaikul A, Carr JJ, Detrano R, et al. Genetic ancestry and the relationship of cigarette smoking to lung function and per cent emphysema in four race/ethnic groups: a crosssectional study. Thorax. 2013;68(7):634–42.
 6.
Pfaff CL, Parra EJ, Bonilla C, Hiester K, McKeigue PM, Kamboh MI, et al. Population structure in admixed populations: effect of admixture dynamics on the pattern of linkage disequilibrium. Am J Hum Genet. 2001;68(1):198–207.
 7.
Fejerman L, Chen GK, Eng C, Huntsman S, Hu D, Williams A, et al. Admixture mapping identifies a locus on 6q25 associated with breast cancer risk in US Latinas. Hum Mol Genet. 2012;21(8):1907–17.
 8.
Cheng CY, Kao WH, Patterson N, Tandon A, Haiman CA, Harris TB, et al. Admixture mapping of 15,280 African Americans identifies obesity susceptibility loci on chromosomes 5 and X. PLoS Genet. 2009;5(5):e1000490.
 9.
Nalls MA, Wilson JG, Patterson NJ, Tandon A, Zmuda JM, Huntsman S, et al. Admixture mapping of white cell count: genetic locus responsible for lower white blood cell count in the Health ABC and Jackson Heart studies. Am J Hum Genet. 2008;82(1):81–7.
 10.
Reich D, Patterson N, Ramesh V, De Jager PL, McDonald GJ, Tandon A, et al. Admixture mapping of an allele affecting interleukin 6 soluble receptor and interleukin 6 levels. Am J Hum Genet. 2007;80(4):716–26.
 11.
Zhu X, Luke A, Cooper RS, Quertermous T, Hanis C, Mosley T, et al. Admixture mapping for hypertension loci with genomescan markers. Nat Genet. 2005;37(2):177–81.
 12.
Price AL, Tandon A, Patterson N, Barnes KC, Rafaels N, Ruczinski I, et al. Sensitive detection of chromosomal segments of distinct ancestry in admixed populations. PLoS Genet. 2009;5(6):e1000519.
 13.
Baran Y, Pasaniuc B, Sankararaman S, Torgerson DG, Gignoux C, Eng C, et al. Fast and accurate inference of local ancestry in Latino populations. Bioinformatics. 2012;28(10):1359–67.
 14.
Yang JJ, Li J, Buu A, Williams LK. Efficient inference of local ancestry. Bioinformatics. 2013;29(21):2750–6.
 15.
Brown R, Pasaniuc B. Enhanced methods for local ancestry assignment in sequenced admixed individuals. PLoS Comput Biol. 2014;10(4):e1003555.
 16.
Zhang J, Stram DO. The role of local ancestry adjustment in association studies using admixed populations. Genet Epidemiol. 2014;38(6):502–15.
 17.
Seldin MF, Pasaniuc B, Price AL. New approaches to disease mapping in admixed populations. Nat Rev Genet. 2011;12(8):523–8.
 18.
Zaitlen N, Pasaniuc B, Sankararaman S, Bhatia G, Zhang J, Gusev A, et al. Leveraging population admixture to characterize the heritability of complex traits. Nat Genet. 2014;46(12):1356–62.
 19.
Meyer KB, O’Reilly M, Michailidou K, Carlebur S, Edwards SL, French JD, et al. Finescale mapping of the FGFR2 breast cancer risk locus: putative functional variants differentially bind FOXA1 and E2F1. Am J Hum Genet. 2013;93(6):1046–60.
 20.
Trynka G, Hunt KA, Bockett NA, Romanos J, Mistry V, Szperl A, et al. Dense genotyping identifies and localizes multiple common and rare variant association signals in celiac disease. Nat Genet. 2011;43(12):1193–201.
 21.
Gusev A, Bhatia G, Zaitlen N, Vilhjalmsson BJ, Diogo D, Stahl EA, et al. Quantifying missing heritability at known GWAS loci. PLoS Genet. 2013;9(12):e1003993.
 22.
Genomes Project C, Abecasis GR, Auton A, Brooks LD, DePristo MA, Durbin RM, et al. An integrated map of genetic variation from 1,092 human genomes. Nature. 2012;491(7422):56–65.
 23.
International HapMap C, Altshuler DM, Gibbs RA, Peltonen L, Altshuler DM, Gibbs RA, et al. Integrating common and rare genetic variation in diverse human populations. Nature. 2010;467(7311):52–8.
 24.
Kooperberg C, Leblanc M. Increasing the power of identifying gene x gene interactions in genomewide association studies. Genet Epidemiol. 2008;32(3):255–63.
 25.
Hemani G, Theocharidis A, Wei W, Haley C. EpiGPU: exhaustive pairwise epistasis scans parallelized on consumer level graphics cards. Bioinformatics. 2011;27(11):1462–5.
 26.
Chatterjee N, Kalaylioglu Z, Moslehi R, Peters U, Wacholder S. Powerful multilocus tests of genetic association in the presence of genegene and geneenvironment interactions. Am J Hum Genet. 2006;79(6):1002–16.
 27.
Aschard H. A Perspective on Interaction Tests in Genetic Association Studies. bioRxiv 2015. doi:10.1101/019661".
 28.
Pasaniuc B, Zaitlen N, Lettre G, Chen GK, Tandon A, Kao WH, et al. Enhanced statistical tests for GWAS in admixed populations: assessment using African Americans from CARe and a Breast Cancer Consortium. PLoS Genet. 2011;7(4):e1001371.
Acknowledgements
We acknowledge funding and support from the National Institute of Health, R03HG006720 (HA), T32HG002536 (RB), R01GM053275 (BP,RB), and R01HG006399 (BP).
Author information
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
HA and BP conceived the study. HA, BP, and AG participated in its design and coordination. HA, AG and RB conducted all analyses. HA, BP, AG and RB were involved in the drafting of the manuscript. All authors read and approved the final manuscript.
Additional files
Additional file 1: Figures S1 and S2.
Power comparison between SNPbased and local ancestrybased interaction tests when using inferred ancestry and 1 M genotyped SNPs. Figure S2. Power comparison between SNPbased and local ancestrybased interaction tests when using a two steps approach. (PDF 555 kb)
Additional file 2:
R script and corresponding input files to perform 1 ° of freedom interaction test between ancestry segments. (ZIP 40 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Aschard, H., Gusev, A., Brown, R. et al. Leveraging local ancestry to detect genegene interactions in genomewide data. BMC Genet 16, 124 (2015). https://doi.org/10.1186/s128630150283z
Received:
Accepted:
Published:
Keywords
 Genegene interaction
 GWAS
 Local ancestry
 Statistical genetics