Low Frequency Variants, Collapsed Based on Biological Knowledge, Uncover Complexity of Population Stratification in 1000 Genomes Project Data
[['In the field of human genetics research , there has been increasing interest in the role of low frequency variation in complex human disease ( defined in this text as variants with a minor allele frequency between 0 . 5%–3% ) .', 'This is in many ways a response to changing technology , but more importantly a response to the inability to completely explain heritability in common complex diseases and recognition of the true multifactorial mechanisms of genetic inheritance [1] .', 'Since low frequency variants are likely essential in understanding the etiology of common , complex traits , it is critical to elucidate the genetic architecture and population substructure of low frequency variants for future work in this field .', 'Factors such as rapid population growth and weak purifying selection have allowed ancestral populations to accumulate an excess of low frequency variants across the genome .', 'This affects genomic analyses in two ways: proportion of deleterious versus neutral variation expected in low frequency variants and population stratification .', 'It has been suggested that slightly deleterious single nucleotide variants ( SNVs ) subjected to weak purifying selection are major players in common disease susceptibility [2] , [3] .', 'For example , Nelson et al . found that in 202 drug target genes , 2/3 of the low frequency variants were nonsynonymous mutations .', 'This is a much higher ratio than found for common variants , and reflect the expected proportion given random mutation and degenerate coding .', 'This ratio also suggests low frequency variants are only weakly filtered by selection [2] , [4] .', 'In addition , low frequency variants represent a considerable proportion of the genome due to recent explosive population growth [3] .', 'Gorlov estimates up to 60% of SNVs in the genome are variants with an allele frequency <5% [5] .', 'Since the allele frequency distribution is skewed towards more low frequency variants , a higher number of low frequency deleterious variants are expected .', 'Subsequently , low frequency variants appear to be enriched for functional variation , including protein coding changes and altered function [6] .', 'Further , low frequency variants exhibit extreme population stratification .', 'Demonstrating the magnitude of low frequency population stratification between two populations , Tennessen et al . identified more than 500 , 000 SNVs using 15 , 585 protein-coding genes from 2 , 440 individuals .', 'Of these SNVs , 86% had a MAF<0 . 5% and 82% were population specific between European Americans and African Americans [3] .', 'Low frequency allele sharing between populations on the same continent can be between 70% and 80% .', 'In contrast , low frequency allele sharing between populations on different continents can be lower than 30% and variants are often unique to a single population .', 'This extreme geographic stratification can lead to higher false positives and difficulty in replicating associations across genetic studies when not considered as part of the experimental design for low frequency SNV analyses [6] .', 'To study the “landscape” of low frequency variant stratification across populations , we grouped low frequency variants across pertinent genome-wide biological features in a series of pairwise population comparisons across multiple ancestries .', 'We define the boundaries of grouping by features , which consist of genomic regions ( one or many ) that belong to a genomic category , for example , a gene or a set of genes in a pathway .', 'Methods that aggregate variants have been shown to be much more powerful than single-variant association testing for low frequency variants [7]–[10] , and thus are reliable to detect population stratification .', 'Our collapsing method , BioBin , provides the opportunity to cast a broader net and uncover stratification across meaningful elements such as genes , pathways , and evolutionary conserved regions by aggregating low frequency variants based on expert biological knowledge .', 'Herein we have applied BioBin to individuals from 1000 Genomes Project Phase I data; we defined “cases” and “controls” randomly between exhaustive pairwise population comparisons .', 'Our goal was to identify features across the genome with differences in low frequency burden between populations; specifically , to look for aggregate differences in low frequency variation between populations , not to detect individual population-specific variants .', 'We show that BioBin is effective in identifying differences in low frequency variant burden centered on biological criteria and highlights the considerable differences in low frequency variants across ancestry groups .', 'These results further emphasize the critical importance of considering low frequency population substructure in future rare and low frequency variant analyses .'], ['We applied BioBin to whole-genome population data using the 1000 Genomes Project Phase I data .', 'The populations , sample sizes , and total number of loci , variants , low frequency variants , and private variants are listed in Table 1 .', 'Although the Iberian population ( IBS ) is listed in Table 1 , this population was not used in the analyses presented in this paper .', 'There was not a sufficient sample size to meet our low frequency criteria ( N\u200a=\u200a14 , MAF cutoff\u200a=\u200a0 . 03 ) .', 'In addition to the differences in overall magnitude of variation seen in Table 1 between these population groups , there were also differences in the distribution of this variation .', 'In Figure 1 , we present an allele frequency density distribution plot of all autosomal chromosomes for all 13 populations .', 'African descent populations have the highest density of low frequency variation .', 'Others have found a similar trend genome-wide [11] .', 'In general , the African ancestral populations not only have more variants overall than other ancestral groups ( see Table 1 ) , these populations also have a higher distribution of low frequency variants than other ancestral groups ( see Figure 1 ) .', 'Although low coverage next generation sequence data is prone to errors , we found no evidence that sequence technology led to differential bias in a way that could explain the trends found in this paper ( Text S1 , Table S1 , Figures S1 , S2 ) .', 'We investigated sample-relatedness with respect to common and low frequency variants using both identity-by-descent ( IBD ) and identity-by-state ( IBS ) estimations , and in each analysis , we found evidence of increased relatedness in ASW ( African ancestry , USA ) , CHB ( Han Chinese Beijing , China ) , CHS ( Han Chinese Shanghai , China ) , CLM ( Medellin , Columbia ) , GBR ( England and Scotland ) , JPT ( Japan ) , LWK ( Luhya , Kenya ) , and MXL ( Mexican Ancestry , California ) .', 'We performed iterative IBD calculations in plink to eliminate related individuals from continental groups .', 'Seventy-five individuals of 1080 total individuals were parsimoniously removed to achieve a pi_hat<\u200a=\u200a0 . 3 in each continental population .', 'The remaining 1 , 005 individuals were used for the binning analyses presented in this paper .', 'An alternate allele sharing method described by Abecasis et al . uses IBS rather than IBD to review allele sharing [12] , [13] .', 'In the case of low frequency or rare variants , IBS approximates IBD .', 'Figure 2 shows within population IBS for all 13 populations for variants with a MAF<3% , where each point represents a pairwise IBS calculation within the same population ( i . e . YRI-YRI but not YRI-CEU ) .', 'In Figure 2A , the pairs with average IBS calculations that fall outside of the cluster are cryptically related individuals with increased allele sharing .', 'Figure 2B shows the IBS calculations after removing 75 individuals with cryptic relatedness .', 'Complete details of these and additional sample-relatedness analyses are available in Text S2 , Figure S3 , Figure S4 , and Figure S5 .', 'Knowledge of population substructure in low frequency variants is critical for genomic studies .', 'We applied BioBin to test for low frequency ( MAF≤0 . 03 ) variant burden differences between 13 populations from the 1000 Genomes Project across different genomic features: genes ( intronic and exonic variants , filtered nonsynonymous and predicted damaging variants ) , intergenic regions , ORegAnno annotated regulatory regions , pathways , pathway-exons , evolutionary conserved regions , and regions considered to be under natural selection .', 'Results are shown in Figure 3 , Figure 4 , Figure 5 , Figure 6 , Figure S6 , and Figure S7 .', 'In each matrix plot , we have indicated the proportion of significant bins ( after Bonferroni correction ) out of the total number of bins generated between two populations .', 'The color intensity represents the proportion of total bins that were significant [0 , 1] .', 'Overall , there are large differences across populations with regard to low frequency variant burden and the distribution of low frequency variants is not random across the genome .', 'The magnitude of stratification corresponds to the mutational landscape of the region .', 'We chose NCBI Entrez to provide the boundaries for gene regions and created a custom role file of intron and exon boundaries using data provided from UCSC Genome Browser [14] .', 'In Figure 3 , the top matrix corresponds to bins created using gene-exon boundaries , the middle matrix corresponds to bins created using gene-intron boundaries , and the bottom matrix corresponds to bins created using regions between genes ( intergenic ) .', 'The values and color intensity within each block represent the proportion of significant bins after Bonferroni correction out of the total number of bins generated between two populations .', 'The coding regions show a trend of a lower proportion of significant bins with low frequency variant burden differences than either the intron or intergenic bins .', 'For example , in the CEU ( Northern/Western European Ancestry , USA ) −YRI ( Yoruba African ) comparison , approximately 44% of the gene-exon bins had significant differences in low frequency variant burden .', 'In contrast , the noncoding region bins , gene-introns and intergenic bins had 66% and 70% of bins with significant differences in low frequency variant burden .', 'The coding regions appear to be under more constraint across populations than noncoding regions .', 'Comparing only the noncoding regions , introns tend to have slightly fewer variation differences than intergenic bins , most likely because introns are by default nearest neighbors to the selective pressures on coding regions .', 'We filtered the gene-exon bins using annotations from the Variant Effect Predictor Software ( VEP ) [15] .', 'We created gene bins with only nonsynonymous variants and a second analysis using only predicted damaging variants annotated by SIFT or PolyPhen2 [15]–[17] .', 'The results in Figure 4 indicate that these potentially functional and significant changes are even more conserved between populations than coding regions ( Figure 3A ) .', 'We used ORegAnno ( Open Regulatory Annotation database ) to define regulatory region boundaries for the bin analysis .', 'The top matrix of Figure 5 shows the 78 population comparisons for the ORegAnno regulatory feature analysis .', 'In comparison to Figure 3 , the annotated regulatory regions have fewer significant bins .', 'For example , in gene-exon analysis shown in Figure 3 , approximately 44% of the ASW-CHB gene-exon bins contained significant differences in low frequency burden .', 'However , in Figure 5 , only 28% of the ASW-CHB annotated regulatory bins contained significant differences in low frequency burden .', 'This trend is consistent across the matrix of population comparisons; regulatory regions have fewer significant bins than the coding or noncoding features of the same population comparison .', 'Several biological pathway and group sources from LOKI ( the Library of Knowledge Integration , which is described in detail in the methods ) were used to generate low frequency variant bins; including , PFAM , KEGG , NetPath , PharmGKB , MINT , GO , dbSNP , Entrez , and Reactome .', 'The Figure 5B shows the 78 population comparisons for the pathway group feature analysis .', 'Of all of the feature analyses , pathway bins consistently show the highest proportion of significant differences in low frequency variant burden between populations .', 'There are several potential explanations .', 'First , since pathway bins are generally much larger than the other feature types , it is possible that large bins increase the false positive rate .', 'Second , the same genes and regions can recur in multiple pathways .', 'If the region has significant differences in low frequency variant burden , then each pathway or group containing that region will have a higher chance of having significant differences in low frequency variant burden .', 'Following this logic , a pathway containing many genes has a higher chance of having at least one gene with extreme low frequency variant stratification .', 'To compare only coding regions within a pathway , we filtered the pathway analysis to include only variants within exons .', 'The proportions are reduced ( shown in Figure S6 ) but still higher than the gene-exon proportions shown in Figure 3A .', 'PhastCons output downloaded from UCSC Genome Browser was used to derive evolutionary conserved feature boundaries for primates , mammals , and more than 40 species of vertebrates .', 'Figure S7 shows the 78 population comparisons for the ECR feature analysis .', 'Of all of the feature analyses , ECR bins had the smallest proportion of significant bins .', 'More ancestrally similar populations tended to have negligible low frequency burden differences in these conserved segments .', 'For example , approximately 7% of the ECR region bins ( vertebrate alignment ) were significantly different between FIN ( Finnish ) and JPT ( Japanese ) individuals .', 'However , the significant number of bins between the two ancestrally similar GBR ( British ) and CEU individuals was less than 1% .', 'To investigate regions of natural selection , we created a feature list using regions recently identified/confirmed by Grossman et al . with the Composite Multiple Signals algorithm on the 1000 Genomes Project data [18] .', 'In addition , a publication by Barreiro et al . provided a list of specific genes with the strongest signatures of positive selection; i . e . genes that contained at least one nonsynonymous or 5′ UTR mutation with an FST value greater than 0 . 65 [19] .', 'After lifting positions to build 37 , there were only 368 regions from the regions identified by Grossman et al . The results are shown in Figure 6 .', 'The top plot corresponds to regions identified in African ancestry , the middle plot corresponds to regions identified in populations of Asian ancestry , and the bottom plot corresponds to regions identified specifically in populations of European ancestry .', 'The trends in these three matrix plots are distinctly different from the trends shown in Figures 3–5 .', 'The blocks of comparisons within a continental group ( shown in black boxes on each matrix plot ) still have very little color , which means that the low frequency variant burden between populations within a continental group is very similar .', 'The main difference is the gain of intensity outside of the continental groups .', 'For example , in Figure 6B ( regions identified in Asian populations ) , the European continental group and Spanish continental group mostly have proportions over 60% when compared to populations of Asian descent .', 'In the same plot , the populations in the African group have proportions over 85% when compared to populations in the Asian group .', 'In general , we found regions considered to be under natural selection unlikely to have significant differences in low frequency burden between ancestrally similar populations , and very likely to have significant differences in regions considered to be under natural selection between ancestrally distant populations ( see Figure 6 ) .', 'Additional analyses were performed using regions identified in other publications and can be found in Text S3 , Table S2 , and Table S3 .', 'Although low frequency variants are commonly assumed as independent ( in low linkage disequilibrium ( LD ) with other variants ) , there are rare haplotypes within related individuals and populations [20] .', 'In Figure 7 , three pairwise population comparisons are shown .', 'We investigated the top 10 ranked bins from the CEU-CHB ( A ) , CHB-YRI ( B ) , and CEU-YRI ( C ) coding and noncoding analyses for presence of LD ( r2>0 . 3 ) between two variants in the same bin .', 'Figure 7 shows bins predominately filled with white-space indicating low to no pairwise LD between variants in those bins .', 'In the top ten bins from these three analyses , rare haplotypes do not appear to be driving the significant differences seen in low frequency variant burden .', 'Since the proportion of significant bins in the feature analyses is considerably higher for pathway bins than any other feature , we wanted to investigate the correlation between pathway p-value and bin size .', 'We chose to assess the correlation between significance and several characteristics of the pathways using the pathway feature CEU-YRI population comparison .', 'Figure S8 and Figure S9 show the correlations between six untransformed and transformed variables ( with outliers removed ) , where each pairwise correlation is significant ( p-value<0 . 05 ) .', 'A bin was considered an outlier if the number of loci in the bin was more than 2 . 5 standard deviations from the mean transformed loci value .', 'The most interesting correlations were the nonlinear correlations between the loci/variants/genomic coverage and p-values .', 'Figure S9B is a higher magnification of the highlighted correlation in Figure S9A , specifically; we plotted the correlation between −log10 p-values and log10 variants .', 'The lowess smoothing function is shown in red , and the function appears to change slope twice .', 'From x\u200a=\u200a1 to x\u200a=\u200a3 , the slope is increasing with increasing number of variants .', 'From x\u200a=\u200a3 to x\u200a=\u200a4 , the slope is near 0 .', 'From x\u200a=\u200a4 to x\u200a=\u200a5 , the slope is increasing with increasing numbers of variants .', 'When the log10-transformed value of the number of variants is less than 3 or greater than 4 , there appears to be a positive correlation between the number of variants in a bin and increasing significance of that bin .', 'However , the data is not uniformly distributed and is sparse in those same areas .', 'Therefore , the trends in the tails are most likely unreliable .', 'We created boxplots describing certain characteristics from each data source .', 'Figure S10 shows that specific sources ( i . e . KEGG ) consistently have larger bin characteristics ( number of loci , number of genes , coverage ( kb ) , etc .', ') and also have much more significant bin p-values ( Figure S10B ) .', 'It appears that certain sources might inherently have more significant groups by nature of the information that source provides , or because of the size of groups found in the source .', 'From the matrix plots shown above , there is undoubtedly a functional component that influences the evolution of low frequency variants .', 'However , from the correlations in the pathway analyses , it is also clear that larger bins can contain more stratification and thus more likely to have significance differences in low frequency burden .', 'In more traditional case/control analyses , large bins are less likely to be significant because increasing binsize generally means more noise to mitigate the signal .', 'However , in this study , when diverse populations are compared , larger bin sizes have more opportunity to capture population stratification .', 'Figure S11 shows the relative number of loci across tested features and varying interregion parameters .', 'Boxplots in Figure S11A represent each feature tested in the population comparison .', 'The small inset figure shows the magnitude of difference between the numbers of loci in pathways ( peak ) versus other feature types .', 'The main plot in Figure S11A shows the same information , but is limited to 2000 loci .', 'In general , ECRs/exonic regions/nonsynonymous gene variants/ORegAnno annotated regions/predicted deleterious gene variants/UTRs are very small bins .', 'Pathway bins have a broader distribution , but in general are much larger .', 'For comparison , we varied the size of intergenic regions ( only noncoding regions ) between 10 kb and 200 kb , those results are shown in Figure S11B .', 'We also split the entire genome ( including coding and noncoding regions ) by various windows between 10 kb–200 kb .', 'Figure S11C represents a genome “average , ” and both Figure S11B and S11C can be used as comparison for feature tests .', 'Figure S11B and S11C show increasing bin size as windows increase , the proportion of significant bins increase as window size increases as well ( see Figure S12 and S13 ) .', 'For example , Figure S13 shows matrix plots from whole-genome “average” analyses ( A–G correspond to 10 kb , 25 kb , 50 kb , 75 kb , 100 kb , 150 kb , and 200 kb respectively ) .', 'According to Figure S11A and S11C , exon bins from the original feature analysis are roughly comparable in size to 10 kb bins from the whole-genome “average” analysis .', 'In the gene analysis results , approximately 43 . 9% of bins are significant after Bonferroni correction between CEU-YRI .', 'Comparatively , the genome average between CEU-YRI for 10 kb bins is 57 . 64% .', 'This supports the idea that coding regions are presumably more functional and perhaps more conserved than other regions in the genome of comparative size .', 'According to Figure S11A and S11C , pathway bins from the original feature analysis are roughly comparable in size to 150 kb bins from the whole-genome “average” analysis .', 'In the pathway analysis results , approximately 81 . 28% of bins are significant after Bonferroni correction between CEU-YRI .', 'Comparatively , the genome average between CEU-YRI for 150 kb bins is 86 . 09% .', 'The gap between pathway bins and “average” genome stratification given similar size is much smaller for pathways than it is for exons .', 'This particular pathway analysis includes introns ( which typically have more variation than coding regions and larger bins are expected to collect more stratification .', 'However , there are still fewer significant bins than expected on average .'], ['Since the reference genome is predominantly of European ancestry [21]–[23] , populations with non-European ancestry generally have more variation with respect to the reference genome than those of European ancestry ( see Table 1 ) .', 'Therefore , to interpret the results of this study , one might conclude that non-European populations have higher rates of sequencing error than European descent populations .', 'However , in the most recent 1000 Genomes Project publication , the authors report an accuracy of individual genotype calls at heterozygous sites more than 99% for common SNVs and 95% for SNVs at a frequency of 0 . 5% .', 'Furthermore , the authors found that variation in genotype accuracy was most related to sequencing depth and technical issues than population-level characteristics [11] .', 'Therefore , neither the sequencing error nor the predominantly European reference genome adequately explains the trends seen in the genomic feature exploration ( see Text S1 , Table S1 , Figures S1 , S2 ) .', 'Both sequence generation ( technology and/or site ) and population identity strongly contributes to underlying stratification in next-generation sequence data .', 'After removing individuals with cryptic relatedness , 4 out of 13 Phase I populations were sequenced entirely using a single sequence technology ( CHB , CHS , JPT , and TSI ) .', 'The other 9 populations had between 3–18 individuals or ∼5%–57% of the population sequenced on technologies other than Illumina ( ABI SOLID or LS454 ) .', 'Note: all three of the Asian populations ( after removing individuals with cryptic relatedness ) were sequenced only with Illumina technologies .', 'In our IBD analysis using variants with a minor allele frequency of 5% or greater and linkage disequilibrium r2<\u200a=\u200a0 . 2 , we identified and eliminated 75 individuals of various population backgrounds .', 'In addition to the previously documented relatedness in 1000 Genomes Project [http://www . 1000genomes . org/phase1-analysis-results-directory] , we also found additional cryptic relatedness seen in other work [24] , [25] .', 'The differences are likely because we used continental groups ( not a single population or the entire 1080 individuals ) to identify cryptically related individuals and in our analysis that could include variants with fixed opposite frequencies and are overall common .', 'This is infrequent in populations of the same continental group , but could be stratification introduced by different sequencing technologies .', 'The major goal of this study was to investigate population stratification across multiple biological features .', 'We created matrix plots to illustrate the proportion of significant bins in each comparison ( shown in Figure 3 , Figure 4 , Figure 5 , Figure 6 , Figure S6 , and Figure S7 ) .', 'Our results show an interesting trend between functional regions of the genome and variant tolerance .', 'Mutations appear to be less tolerated in functional regions .', 'Similarly , ECRs , which are known to be conserved among species , are also the features least likely to have variation burden differences between two populations .', 'There is some debate about selection and functional significance in these conserved regions , it is unknown what factors have the largest effect on mutation rates [26] , but it is possible that consistently low mutation rates in these features have generated conserved regions throughout evolution [27] .', 'There are two potential explanations:', '1 ) additional level of repair of DNA damage in transcriptional active regions by transcription coupled repair ( TCR ) ,', '2 ) approximately 3% of the genome is subject to negative selection , however it is estimated that functionally dense regions contain up to 20% of the sites under selection [26] , [27] .', 'A number of the top results in each comparison have an interesting context , particularly in light of natural selection .', 'Perhaps one of the most notable is SLC24A5 ( Ensembl ID:ENSG00000188467 ) , which is one of the top ten results in 19 out of 78 populations comparisons in the gene feature analysis .', 'European specific selective sweeps estimated in the last 20 , 000 years suggest that SLC24A5 is key in skin pigmentation and Zebrafish with “golden” mutations exhibit melanosomal changes [28]–[30] .', 'The presence of selection in particular populations due to environmental factors such as distance to the equator has led to the evolution and expansion of low frequency variants in some populations but not others .', 'A second notable top result is DARC ( Ensembl ID:ENSG00000213088 ) , which encodes the Duffy antigen .', 'The DARC gene bin was in the top ten results in 14 out of 78 population comparisons in the gene feature analysis .', 'It has long been known that populations of African descent have increased diversity due to natural selection at this location , which prevents Plasmodium vivax infection .', 'The top result from the regulatory region analysis was a region on chromosome 20 ( chr20:45395536–45396346 ) which was in the top ten bins in 24 out of 78 populations comparisons in the ORegAnno feature analysis .', 'This region also overlaps ENCODE transcription factor binding sites in multiple cell lines , including: CTCF , POLR2A , NFYA , E2F1 , FOS , and more .', 'It was also annotated as an insulator in multiple cell lines in ENCODE Chromatin State Segmentation analyses using Hidden Markov Models [14] , [31] .', 'One last example , chr15 . 968 , contains variants in the genome location chr15:48400199–48412256 .', 'This bin is one of the top ten bins in 17 out of 78 population comparisons in the intergenic analysis .', 'The region covered by the chr15 . 968 bin is less than 1 kb upstream of SCL24A5 on chromosome 15 and overlaps with several transcription factor-binding sites ( including CTCF ) , regions thought to be weak enhancers , and regions thought to be insulators .', 'According to Grossman et al . , there are defined regions under natural selection before and after this region ( chr15:45145764–45258860 and chr15:48539026–48633153 ) , and all are very likely to participate in the transcriptional regulation of SLC24A5 [18] .', 'The natural selection features require knowledge of three things for interpretation:', '1 ) population A ,', '2 ) population B , and', '3 ) the population where this signature was identified .', 'When all three of these are within the same ancestral or continental group , we expect very few differences in low frequency burden .', 'However , if population A is the same or similar to the population possessing the selection signature and population B is different , we expect significant differences in low frequency burden between population A and population B . In our results , we found that the vast majority of regions considered to be under natural selection had significant differences in low frequency burden between disparate ancestral populations , which support the theory of selection in these regions .', 'In general , size of bins can influence the number of stratified variants contained and thus the significance of that bin .', 'It is important to prove that this is because larger bins have a greater opportunity to “collect” variants that are stratified and not because of inflated type I error .', 'We have tested type I error rates in bins between approximately 40 variants to over 100 , 000 variants , which covers all analyses presented in this paper , and found no correlation between bin size and Type I error rate ( unpublished data ) .', 'However , it should also be noted that while larger bins have more chances to collect stratified variants , there is also a larger capacity to collect neutral variants that contribute noise and decrease the signal .', 'Using CEU-YRI pathway burden analysis , we reviewed the correlation between pathway size and significance .', 'The number of genes in pathways ranged from 1 to over 700 genes , with the average around 5 genes per group .', 'Correlations for this data are shown in Figure S9 .', 'Not surprisingly , there were very linear and positive correlations between number of loci , number of variants , and genomic coverage .', 'However , each of these had a nonlinear and somewhat complex relationship with the log-transformed p-value .', 'This is highlighted in Figure S9B , which shows the relationship between the −log10 transformed p-value and the log10 transformed number of variants in the bin .', 'The trend indicates that p-values are positively correlated ( become more significant ) with numbers of variants in a bin when the numbers of variants are relatively small or very large .', 'Two reasons could explain this correlation:', '1 ) the false-positive rate is influenced by bin size ( number of variants per bin ) , and', '2 ) true signals from gene bins with burden differences perpetuate higher numbers of significant pathway bins .', 'After extensive simulation testing ( unpublished data ) and recent publications in the literature , we believe the later is true [32] , [33] .', 'A single or small number of child bins ( gene bins in this example ) , can drive parent bins ( pathways in this example ) to be significant even if no other child bin contains stratification .', 'The comparison in Figure S10 between group sources available in LOKI suggests KEGG , NetPATH , PharmGKB , and Reactome have consistently larger bins ( higher number of loci , variants , and coverage ) .', 'On average , these same four sources also tend to have bins with smaller p-values .', 'Therefore , larger pathways are more likely to contain a gene with extreme low frequency variant stratification .', 'Population stratification is incredibly important in genomic analyses , particularly when low frequency variants are being studied .', 'Expected stratification and potential bi
['Introduction', 'Results', 'Discussion', 'Methods']
['Analyses investigating low frequency variants have the potential for explaining additional genetic heritability of many complex human traits .', 'However , the natural frequencies of rare variation between human populations strongly confound genetic analyses .', 'We have applied a novel collapsing method to identify biological features with low frequency variant burden differences in thirteen populations sequenced by the 1000 Genomes Project .', 'Our flexible collapsing tool utilizes expert biological knowledge from multiple publicly available database sources to direct feature selection .', 'Variants were collapsed according to genetically driven features , such as evolutionary conserved regions , regulatory regions genes , and pathways .', 'We have conducted an extensive comparison of low frequency variant burden differences ( MAF<0 . 03 ) between populations from 1000 Genomes Project Phase I data .', 'We found that on average 26 . 87% of gene bins , 35 . 47% of intergenic bins , 42 . 85% of pathway bins , 14 . 86% of ORegAnno regulatory bins , and 5 . 97% of evolutionary conserved regions show statistically significant differences in low frequency variant burden across populations from the 1000 Genomes Project .', 'The proportion of bins with significant differences in low frequency burden depends on the ancestral similarity of the two populations compared and types of features tested .', 'Even closely related populations had notable differences in low frequency burden , but fewer differences than populations from different continents .', 'Furthermore , conserved or functionally relevant regions had fewer significant differences in low frequency burden than regions under less evolutionary constraint .', 'This degree of low frequency variant differentiation across diverse populations and feature elements highlights the critical importance of considering population stratification in the new era of DNA sequencing and low frequency variant genomic analyses .']
['Low frequency variants are likely to play an important role in uncovering complex trait heritability; however , they are often continent or population specific .', 'This specificity complicates genetic analyses investigating low frequency variants for two reasons: low frequency variant signals in an association test are often difficult to generalize beyond a single population or continental group , and there is an increase in false positive results in association analyses due to underlying population stratification .', 'In order to reveal the magnitude of low frequency population stratification , we performed pairwise population comparisons using the 1000 Genomes Project Phase I data to investigate differences in low frequency variant burden across multiple biological features .', 'We found that low frequency variant confounding is much more prevalent than one might expect , even within continental groups .', 'The proportion of significant differences in low frequency variant burden was also dependent on the region of interest; for example , annotated regulatory regions showed fewer low frequency burden differences between populations than intergenic regions .', 'Knowledge of population structure and the genomic landscape in a region of interest are important factors in determining the extent of confounding due to population stratification in a low frequency genomic analysis .']
['genome-wide association studies', 'genome sequencing', 'genome analysis tools', 'biology', 'genomics', 'computational biology']