Association Analysis of Candidate Variants in Admixed Brazilian Patients With Genetic Generalized Epilepsies

Genetic generalized epilepsies (GGEs) include well-established epilepsy syndromes with generalized onset seizures: childhood absence epilepsy, juvenile myoclonic epilepsy (JME), juvenile absence epilepsy (JAE), myoclonic absence epilepsy, epilepsy with eyelid myoclonia (Jeavons syndrome), generalized tonic–clonic seizures, and generalized tonic–clonic seizures alone. Genome-wide association studies (GWASs) and exome sequencing have identified 48 single-nucleotide polymorphisms (SNPs) associated with GGE. However, these studies were mainly based on non-admixed, European, and Asian populations. Thus, it remains unclear whether these results apply to patients of other origins. This study aims to evaluate whether these previous results could be replicated in a cohort of admixed Brazilian patients with GGE. We obtained SNP-array data from 87 patients with GGE, compared with 340 controls from the BIPMed public dataset. We could directly access genotypes of 17 candidate SNPs, available in the SNP array, and the remaining 31 SNPs were imputed using the BEAGLE v5.1 software. We performed an association test by logistic regression analysis, including the first five principal components as covariates. Furthermore, to expand the analysis of the candidate regions, we also interrogated 14,047 SNPs that flank the candidate SNPs (1 Mb). The statistical power was evaluated in terms of odds ratio and minor allele frequency (MAF) by the genpwr package. Differences in SNP frequencies between Brazilian and Europeans, sub-Saharan African, and Native Americans were evaluated by a two-proportion Z-test. We identified nine flanking SNPs, located on eight candidate regions, which presented association signals that passed the Bonferroni correction (rs12726617; rs9428842; rs1915992; rs1464634; rs6459526; rs2510087; rs9551042; rs9888879; and rs8133217; p-values <3.55e^–06). In addition, the two-proportion Z-test indicates that the lack of association of the remaining candidate SNPs could be due to different genomic backgrounds observed in admixed Brazilians. This is the first time that candidate SNPs for GGE are analyzed in an admixed Brazilian population, and we could successfully replicate the association signals in eight candidate regions. In addition, our results provide new insights on how we can account for population structure to improve risk stratification estimation in admixed individuals.

Introduction

Genetic generalized epilepsies (GGEs) are a group of epilepsy syndromes in which the main feature is the recurrence of generalized onset seizures with no known or suspected etiology other than possible genetic predisposition (Berg et al., 2010; Scheffer et al., 2017). GGEs are among the most common types of epilepsy, with an estimated prevalence of 190 per 100,000 individuals (Aaberg et al., 2017). They include well-established syndromes: childhood absence epilepsy, juvenile myoclonic epilepsy (JME), juvenile absence epilepsy (JAE), myoclonic absence epilepsy (a rare form of GGE), epilepsy with eyelid myoclonia (Jeavons syndrome), generalized tonic–clonic seizures, and generalized tonic–clonic seizures alone (Berg et al., 2010). These different GGE syndromes share most genetic susceptibility factors, suggesting an important correlation among the clinical subtypes (International League Against Epilepsy Consortium on Complex Epilepsies (ILAE Consortium on Complex Epilepsies)., 2018). The diagnosis of GGE relies mainly on clinical information and electroencephalographic examination (Scheffer et al., 2017).

Previous genome-wide association studies (GWASs) and exome sequencing analyses have identified 48 single-nucleotide polymorphisms (SNPs) putatively associated with susceptibility to the GGEs (EPICURE Consortium et al., 2012a,b; International League Against Epilepsy Consortium on Complex Epilepsies (ILAE Consortium on Complex Epilepsies), 2014; International League Against Epilepsy Consortium on Complex Epilepsies (ILAE Consortium on Complex Epilepsies)., 2018; Zhang et al., 2014; Wang et al., 2019). These SNPs are located in or near several genes encoding ion channels and synaptic vesicles, making them plausible candidates for the susceptibility to epilepsy (International League Against Epilepsy Consortium on Complex Epilepsies (ILAE Consortium on Complex Epilepsies)., 2018). However, most of these studies evaluated non-admixed populations, including five studies based on Europeans, three based on Asian populations, mainly Chinese, and two based on African populations (EPICURE Consortium et al., 2012a,b; International League Against Epilepsy Consortium on Complex Epilepsies (ILAE Consortium on Complex Epilepsies), 2014; International League Against Epilepsy Consortium on Complex Epilepsies (ILAE Consortium on Complex Epilepsies)., 2018; Zhang et al., 2014; Wang et al., 2019). It is well known that admixed American populations are underrepresented in GWASs, decreasing the accuracy of replicating, predicting, and estimating polygenic risks for complex disorders in these populations (Martin et al., 2017, 2019).

Therefore, this work aims to investigate if a genetic association exists between previously reported candidate SNPs and GGEs in a cohort of admixed Brazilians. To accomplish this goal, we first investigated the population structure of Brazilian patients with GGE. Subsequently, we performed an association study using the 48 previously reported candidate SNPs and their flanking regions.

Materials and Methods

Subjects

We evaluated a total of 87 patients with GGE who were followed up prospectively in the outpatient epilepsy clinic of the University of Campinas (UNICAMP) hospital. All patients had the diagnosis of GGE according to criteria established by the International League Against Epilepsy (ILAE) (Berg et al., 2010; Fisher et al., 2014). Patients were compared with a group of 340 individuals without any neurological disorder from the BIPMed database (Rocha et al., 2020). Both samples are predominantly from the Southeastern region in Brazil. Among the patients with GGE, we found 63 with JME, 10 with JAE, four generalized tonic–clonic seizures alone, two with Jeavons syndrome, one with myoclonic absence epilepsy, one with epilepsy with generalized tonic–clonic seizures, and six patients in whom a specific GGE syndrome could not be determined. All research participants signed an informed consent form previously approved by our Institutional Research Ethics Committee (IRB # 12112913.3.0000.5404).

Single-Nucleotide Polymorphism Quality Control and Population Structure Analysis

We extracted the genotypes for the 48 candidate SNPs (Table 1) from the SNP-array data generated by the Genome-Wide Human SNP Array 6.0 (Affymetrix Inc., Thermo Fisher Scientific, Waltham, MA, United States). These SNP-array data contain 905,171 available SNPs (GRCh37 build). To obtain an unbiased estimation of the population structure of our samples, we processed the SNP-array dataset of the 87 patients with GGE and the 340 BIPMed controls according to previous processing recommendations and pipelines (Anderson et al., 2010; Secolin et al., 2019). First, we removed ambiguous variants (with G/C or A/T alleles) from each dataset. Next, we merged the two datasets into one larger admixed Brazilian dataset (N = 427), maintaining only biallelic SNPs, autosomal SNPs, SNPs without Hardy–Weinberg disequilibrium (p-value <0.000001), and missing data <10%. Then, we estimated the heterozygosity rate for each sample and removed individuals with heterozygosity rates higher or lower than three standard deviations from the mean to avoid individuals with high inbreeding (low heterozygosity rates) or sample contamination (high heterozygosity rates). We also removed pairs of individuals who presented a proportion of identical-by-state (IBS) alleles >0.85, which could indicate duplicated samples, and individuals with genomic relatedness matrix estimations higher than 0.125, which is the expected genomic relatedness for third-degree relatives (Anderson et al., 2010). The merging process, genotyping, and sample filtering were performed using PLINK 1.9 software (Purcell et al., 2007).

TABLE 1

Table 1. Descriptive statistics and logistic regression analysis of the candidate SNPs and the nine flanking SNPs that passed the most stringent Bonferroni correction (p-value = 3.55e^–06; in bold).

Subsequently, we merged the filtered admixed Brazilian sample with the 1000 Genomes Project (1KGP) dataset (The 1000 Genomes Project Consortium et al., 2015), maintaining the SNPs present only in the admixed Brazilian sample. After merging, we removed SNPs with a minor allele frequency (MAF) < 0.01 and SNPs in linkage disequilibrium (LD), using the following parameters: window size = 50 SNPs, shift step = 5 SNPs, and r² = 0.5 (Anderson et al., 2010). We compared our dataset with the 1KGP data by principal component analysis (PCA) using PLINK v1.9 software (Purcell et al., 2007) to evaluate the presence of population-based outliers in the Brazilian samples.

To evaluate whether patients with GGE and BIPMed controls present population stratification, we performed the analysis of molecular variance (AMOVA) (Excoffier et al., 1992) using the poppr.amova R package and the RStudio interface, comparing the genetic distance among the two groups based on a set of 10,000 random SNPs across the genome. The AMOVA partitions the source of genetic variance (σ²) into two components: within-groups and between-groups. The null hypothesis states that the samples were obtained from a global population, with variation due to random sampling in the construction of populations. Thus, we would expect a high heterogeneity within groups (σ² = 100%) and no heterogeneity between groups (σ² = 0%). On the other hand, under the alternative hypothesis, each group was obtained from different populations, and we would expect a low heterogeneity within groups (σ² < 100%) and high heterogeneity between groups (σ² > 0%) (Excoffier et al., 1992). Therefore, to evaluate the significance of σ² components, we generated a Monte Carlo null distribution of 10,000 variance components and tested against the observed variance components by the randtest function in the ade4 R package.

Single-Nucleotide Polymorphism Selection and Imputation

We observed that 31 SNPs were not found in the SNP-array dataset. Therefore, we performed an imputation of all 48 SNPs to obtain the missing SNPs and to evaluate the concordance between the imputed genotypes and the genotypes assessed by the SNP array. Since we analyzed a sample of admixed individuals, we elected to perform the imputation using two approaches. First, we phased and imputed the dataset using SHAPEIT2 v2.r387 (O’Connell et al., 2014) and BEAGLE v5.1 software (Browning et al., 2018) using the default software parameters for phasing and imputation. As a reference for the BEAGLE imputation, we used the 1KGP dataset (GRCh37/hg19 assembly) (The 1000 Genomes Project Consortium et al., 2015). To save on computation time, we imputed only the chromosomes in which the candidate SNPs are located (Table 1). We also evaluated whether the genotypes were successfully imputed by the correlation (in terms of r²) of genotype dosage values between the imputed genotypes and true genotypes used as a reference from the 1KGP provided by the BEAGLE software. For the second imputation approach, we used the TOPMED Imputation Server (Das et al., 2016), with the TOPMed v.R2 on GRCh38 build (Kowalski et al., 2019). The TOPMED server imputation performed the liftover from GRCh37 to GRCh38 and the phasing using the EAGLE v.2.4 algorithm. Finally, imputation was performed by minimac4.

Candidate Single-Nucleotide Polymorphism Association Analysis

After genotype and individual filtering, 360 individuals remained (69 patients with GGE and 291 BIPMed controls), which were used in the association analysis. We estimated the statistical power of our sample by the genpwr package in R (Moore et al., 2020), which analyzes the statistical power under the evaluation between true and test genetic models (Dominant, Additive, Recessive, 2df/unspecified model). In this case, we evaluate the statistical power using a vector of MAFs (from 0.05 to 0.45, by 0.05) and an odds ratios (from 1.5 to 2.0, by 0.1) since not all candidate SNPs presented OR estimations from previous studies. We also set the following parameters for genpwr: model = logistic; N = 360; case/control ratio = 69/291 = 0.237; and alpha = 0.05.

We evaluated candidate SNP association and OR estimation by logistic regression analysis using the PLINK v1.9 software (Purcell et al., 2007), including the first five PCs as covariates. We did not include age, age at seizure onset, and sex since these variables have not been correlated with the GGE phenotype (Berg et al., 2010; Scheffer et al., 2018).

It has been reported that SNPs found to be associated with the phenotype by GWAS in one population may be only nominally associated or non-associated in another population due to difference in LD across populations (Akiyama et al., 2019; Chen et al., 2020; Graff et al., 2021); however, it does not mean that an associated signal in the genomic region cannot be replicated. This is because the SNPs ascertained from GWAS are only tagging variants linked to causal ones. The lack of signals in the replication population could simply be caused by the broken linkage between tagging and causal variants. Therefore, to account for the difference in LD across populations and to investigate the transferability of previous GWAS signals, we used the SNP-array dataset, filtered for population structure and without LD pruning (652,883 SNPs), to interrogate the SNPs flanking the 1 Mb upstream and downstream the candidate SNPs by logistic regression. We assumed a p-value adjusted by the Bonferroni correction to avoid biased results due to the multiple comparisons. In this case, we used two thresholds: the first threshold took into account the 48 SNPs (p-value = 0.05/48 = 0.001), assuming one effective test per region, which is a reasonable assumption and may lead to more informative results. However, this threshold may not be stringent enough. Therefore, we also evaluate the results under a second threshold, considering all the 48 candidate SNPs and the additional flanking SNPs tested, and the results were plotted using the qqman package in R software (Turner, 2014).

Since previous studies of GGE were based on European populations and admixed Brazilians have a large proportion of European ancestry, we decided to evaluate whether the candidate SNP allele frequencies are similar between Brazilian and European populations. We extracted European allele frequencies from the gnomAD database (Karczewski et al., 2020) and performed a two-proportion Z-test using the prop.test function in R. Also, we included African populations from gnomAD in the analysis due to the sub-Saharan African ancestry component present in Brazilian populations. However, since gnomAD does not separate Native American populations in the database, we include the Latin population in the analysis as a proxy.

Results

Population Structure Analysis

The principal components in the PCA plot indicate that both cases and controls clustered together and were spread between Europeans, sub-Saharan Africans, and other admixed American populations (Figure 1). The AMOVA results showed that 99.61% of the genetic variation was observed within groups (patients or controls), and only 0.39% of the genetic variation was observed between groups (Table 2). Because we have one hierarchical level of stratification (patients/controls), the poppr.amova package provided one total φ-statistics = 0.0031, with a p-value = 0.001 (Table 2), indicating evidence of population stratification between patients and controls and the necessity of population structure correction in further association tests.

FIGURE 1

Figure 1. Principal component analysis (PCA) plot of the Brazilian samples and 1000 Genomes Project (1KGP) dataset. The x-axis and y-axis show the first and second principal components (PC1 and PC2) and their respective percentage variability. Each point represents one individual and each color indicates patients with genetic generalized epilepsy (GGE), BIPMed controls, and the continental populations described in 1KGP. as follows: Sub-Saharan Africans (AFR); Europeans (EUR); admixed Americans (AMR); Southwestern Asians (SAS); Eastern Asians (EAS).

TABLE 2

Table 2. AMOVA results.

Single-Nucleotide Polymorphism Selection and Imputation

According to the imputation results from the BEAGLE software (Browning et al., 2018), the correlation between the estimated allele dosage and the true allele dosage from the 1KGP is used as reference (in terms of r²) and presented a minimum value of 95%. In addition, all 17 SNPs genotyped by the SNP array were correctly imputed by the BEAGLE software. However, we observed that the 17 SNPs genotyped by the SNP array presented only 45.7% of matching (on average) with genotypes imputed by the TOPMED server. Thus, we decided to perform further analysis using the imputed genotypes generated by the BEAGLE software.

Candidate Single-Nucleotide Polymorphism Association Analysis

As detailed in Table 1, one candidate SNP (rs1046276) was withdrawn from further association analysis due to the presence of the Hardy–Weinberg disequilibrium (p < 0.000001). According to the analysis performed using the genpwr package (Moore et al., 2020), we observed that the Additive model presented the highest power estimation. We did not observe 80% of statistical power for OR ≤ 1.6 (≥ 0.62 for protection effect) (Figures 2A,B). However, we calculated that our study had 80% power to detect an increased risk in terms of OR ≥ 1.7 (≤0.58 for protection effect) with MAF > 0.25 (Figure 2C), OR ≥ 1.8 (≤0.55 for protection effect) with MAF > 0.2 (Figure 2D), and OR ≥ 1.9 (≤0.52 for protection effect) with MAF > 0.15 (Figures 2E,F).

FIGURE 2

Figure 2. Statistical power estimation. The figure includes six panels (A–F) for the six OR values evaluated. Each panel shows the statistical power (y-axis) estimated by minor allele frequency (MAF) (x-axis) for the combination of the three true genetic models (Dominant, Additive, and Recessive) with the four test models (Dominant, Additive, Recessive, and 2df/unspecified genetic model). Each colored line represents the test model, and each point represents the MAF.

We identified suggestive evidence of a protective effect for the SNP rs1939012^∗A allele (MAF = 0.446; OR = 0.60; 95% CI = 0.40–0.90; nominal p-value = 0.0133) and rs4794333^∗A allele (MAF = 0.356; OR = 0.64; 95% CI = 0.42–0.97; nominal p-value = 0.0361) and an increased risk for rs9788^∗G (MAF = 0.315; OR = 1.57; 95% CI = 1.05–2.34, nominal p-value = 0.0274). However, these results did not pass the corrections for multiple comparisons by Bonferroni (Table 1). Interesting, we found 14,047 flanking SNPs, encompassing 29 candidate regions (Supplementary Data). As shown in Figure 3, under the p-value threshold = 0.001, we observed that the association signals in all candidate regions passed the Bonferroni correction. Adjusting the p-values by Bonferroni for 14,095 tests (p-value = 3.55e^–06), we observed that nine flanking SNPs passed the multiple comparison adjustment: rs12726617 and rs9428842 on chromosome (chr.) 1q43; rs1915992 on chr. 3p14.2; rs1464634 on chr. 3q26.1; rs6459526 on chr. 6p22.3; rs2510087 on chr. 11q22.3; rs9551042 on chr. 13q12.12; rs9888879 on chr. 16p11.2; and rs8133217 on chr. 21q22.3 (Figure 3 and Table 1).

FIGURE 3

Figure 3. Manhattan plot of genetic generalized epilepsy (GGE) candidate regions. The figure shows a plot for each candidate chromosome, including the chromosome position (x-axis) and the –log10(p-value) in the y-axis. The green and black points represent the candidate and the flanking single-nucleotide polymorphisms (SNPs), respectively. The blue line indicates the suggestive association signal based on the p-value adjusted by the Bonferroni correction under the 48 candidate SNPs. The red line indicates the association signal based on the p-value adjusted by Bonferroni under the 48 candidate SNPs plus the 14,047 flanking SNPs.

Since most Brazilian ancestry is derived from European populations (Kehdy et al., 2015; Moura et al., 2015; Secolin et al., 2019), we could hypothesize that effect sizes in terms of OR would present a higher correlation with European effect sizes comparing with Chinese or European/African American samples from previous studies (EPICURE Consortium et al., 2012a,b; International League Against Epilepsy Consortium on Complex Epilepsies (ILAE Consortium on Complex Epilepsies), 2014; International League Against Epilepsy Consortium on Complex Epilepsies (ILAE Consortium on Complex Epilepsies)., 2018; Zhang et al., 2014; Wang et al., 2019). Thus, we show a comparison of the OR estimations of 11 SNPs, which were available from the previous studies, and the OR estimations in our admixed Brazilian samples (Table 3). Remarkably, Chinese and European/African American samples also presented similar OR estimations compared with admixed Brazilians. Two SNPs had different OR estimations for admixed Brazilians compared with European and Chinese samples (rs10496964 and rs11890028).

TABLE 3

Table 3. Odds ratio comparison among studies.

Furthermore, the two-proportion Z-test results showed that 25 candidate SNPs have allele frequencies that were different when comparing admixed Brazilian and the ancestral populations. Among them, 16 SNPs presented differences in allele frequencies comparing admixed Brazilian and European populations. All 25 SNPs presented different allele frequencies when comparing admixed Brazilian and African samples. Remarkably, we also found 15 candidate SNPs with different allele frequencies when comparing admixed Brazilians and the Latin American samples in the gnomAD database (Table 4).

TABLE 4

Table 4. Two-proportion Z-test results comparing Brazilian samples with European, African, and Latin-American samples from gnomAD.

Discussion

The Brazilian population was formed by an admixture of three main ancestry populations: Europeans, sub-Saharan Africans, and Native Americans (Kehdy et al., 2015; Moura et al., 2015; Secolin et al., 2019). In this scenario, it is important to explore whether candidate SNPs previously identified as associated with complex disorders in non-admixed populations also display association signals in the Brazilian admixed population. By doing so, one can better estimate the impact of population structure in estimating polygenic risks, avoiding misinterpretation of risk scores calculated in other populations.

Previous genetic association studies have identified 48 candidate SNPs associated with GGEs (EPICURE Consortium et al., 2012a,b; International League Against Epilepsy Consortium on Complex Epilepsies (ILAE Consortium on Complex Epilepsies), 2014; International League Against Epilepsy Consortium on Complex Epilepsies (ILAE Consortium on Complex Epilepsies)., 2018; Zhang et al., 2014; Wang et al., 2019). These studies were all performed in non-admixed populations, predominantly of European ancestry, raising the question of reproducibility of these results in other populations. Lack of transferability of GWAS results and polygenic risk scores obtained from Europeans and American admixed populations have previously been reported (Martin et al., 2017, 2019), making it important to investigate whether these SNPs are associated with GGEs in our admixed Brazilian sample.

An alternative explanation for the lack of reproducibility among populations relies on the observation that only tagging SNPs are ascertained in GWAS, and the lack of replication in different populations could be due to broken linkage between the tagging SNPs and the causal variants (Akiyama et al., 2019; Chen et al., 2020; Graff et al., 2021). Thus, we searched for SNPs flanking 1 Mb upstream and downstream of the candidate regions to investigate this issue. Indeed, we found 14,047 flanking SNPs, and nine of them presented statistically significant association signals after stringent corrections for multiple comparisons (p-value < 3.55e^–06). These nine SNPs encompass eight candidate regions (Table 2 and Figure 3), which were previously found associated in European samples (1q43; 3p14.2; 3q26.1; 6p22.3; 11q22.3; 13q12.12; 16p.11.2; and 21q22.3), and two of them were also found associated in a mixed sample of European, African, and Asian populations (3q26.1 and 16p.11.2) (EPICURE Consortium et al., 2012a,b; International League Against Epilepsy Consortium on Complex Epilepsies (ILAE Consortium on Complex Epilepsies), 2014; International League Against Epilepsy Consortium on Complex Epilepsies (ILAE Consortium on Complex Epilepsies)., 2018; Wang et al., 2019). Therefore, we may suggest that polygenic risk scores calculated in European populations at these specific loci could indeed be transferable to admixed Brazilian individuals.

However, although all these 29 candidate regions passed the Bonferroni correction based on the 48 candidate SNPs (p-value = 0.001), we understand that this p-value threshold is not stringent. Thus, the lack of association signal cannot be discarded for the 20 remaining candidate regions. Thus, one may still speculate that the lack of reproducibility could be due to the absence of statistical power, population stratification, or the differences in the genomic structure of the admixed sample compared with the previously studied populations.

Although we have identified flanking SNPs in the neighborhood of the candidate regions, which presented 80% of statistical power to detect increased risk or protection allele effect, we acknowledge the limited statistical power provided by the cohort analyzed, with 87 patients with GGE and 340 controls.

Despite the observed high heterogeneity within groups (σ² = 99.61%) and low heterogeneity between patients and controls (σ² = 0.39%), the statistics based on AMOVA results revealed evidence of population stratification between patients with GGE and the BIPMed controls. Thus, we corrected possible spurious association results by taking the first five principal components into account in the logistic regression model (Marchini et al., 2004; Price et al., 2010).

Indeed, the two-proportion Z-test showed that 16 SNPs presented different allele frequencies when comparing admixed Brazilian and European samples, further substantiating the hypothesis of lack of genetic association due to genetic differences when comparing the admixed Brazilians and Europeans.

It is important to note that 31 SNPs were not found in the SNP-array dataset, and we decided to impute them from all populations available in the 1KGP dataset (The 1000 Genomes Project Consortium et al., 2015). Previous studies have demonstrated that imputation accuracy for populations with a high proportion of European ancestry is higher than for populations with African or Native American ancestry (Martin et al., 2017). In addition, the EPIGEN-Brazil Initiative has also imputed admixed Brazilian samples from the 1KGP dataset with high confidence variants (Magalhães et al., 2018). However, the imputation by the TOPMED Consortium has demonstrated improved quality of variant imputation for admixed African and Hispanic/Latin populations compared with the 1KGP dataset (Kowalski et al., 2019). Thus, we also used this approach for comparison. We observed a perfect match between the SNPs genotyped in the SNP-array and their imputed correspondents for the BEAGLE imputation using the 1KGP as reference. By contrast, there was only 45.7% correspondence between the SNPs genotyped and the imputed SNPs using TOPMED. Thus, we can argue that Hispanic/Latin samples included in the TOPMED reference panel (Kowalski et al., 2019) may not represent the genomic structure of admixed Brazilians (Adhikari et al., 2016). This is an important finding and indicates that although allele frequencies of admixed Brazilian populations are different from other populations reported in public databases (Adhikari et al., 2016; Magalhães et al., 2018; Rocha et al., 2020), there is a remarkable accuracy in the SNP imputation for admixed Brazilian individuals based on populations from the 1KGP database, as demonstrated by our results and elsewhere (Magalhães et al., 2018).

In conclusion, we replicated association signals on eight candidate regions previously found in European populations, indicating the possibility of transferability of polygenic risk scores from European studies to admixed Brazilian populations in these specific candidate regions. In addition, we show evidence that differences in the genetic architecture of the population may hinder the replication of association results in admixed Brazilians for the remaining candidate regions, thus supporting the hypothesis of population differences influencing the association results in the present study. Also, we documented the effect of different methods/databases used for genotype imputation in admixed Brazilians. These results could be relevant to improving stratification risk estimation and future precision health applications in admixed Brazilian patients with GGEs and other complex disorders.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ebi.ac.uk/eva/, PRJEB39251; https://www.ebi.ac.uk/ena, PRJEB45235.

Ethics Statement

The studies involving human participants were reviewed and approved by Comitê de Ética em Pesquisa da Universidade Estadual de Campinas. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

FK performed data processing, statistical analysis, and imputation. TA and PA performed data acquisition and SNP array genotyping. MA, CY, and FC performed the clinical analysis of GGE patients. FC and IL-C served as the principal investigators. RS conceptualized the work, created the study design, and served as a principal investigator. All authors reviewed and approved the final version of the manuscript.

Funding

This work was supported by the Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP), Brazil, grant number 2013/07559-3; and Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES, grant number: 001), Brazil. FK was supported by a studentship from Conselho Nacional de Pesquisa (CNPq), Brazil. IL-C is supported by CNPq (grant number 311923/2019-4). RS was supported by FAPESP (grant number 2019/08526-8).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.672304/full#supplementary-material

References

The 1000 Genomes Project Consortium, Auton, A., Brooks, L. D., Durbin, R. M., Garrison, E. P., Kang, H. M., et al. (2015). A global reference for human genetic variation. Nature 526, 68–74. doi: 10.1038/nature15393

PubMed Abstract | CrossRef Full Text | Google Scholar

Aaberg, K. M., Surén, P., Søraas, C. L., Bakken, I. J., Lossius, M. I., Stoltenberg, C., et al. (2017). Seizures, syndromes, and etiologies in childhood epilepsy: the International League Against Epilepsy 1981, 1989, and 2017 classifications used in a population-based cohort. Epilepsia 58, 1880–1891. doi: 10.1111/epi.13913

PubMed Abstract | CrossRef Full Text | Google Scholar

Adhikari, K., Mendoza-Revilla, J., Chacón-Duque, J. C., Fuentes-Guajardo, M., Ruiz-Linares, A., Chacón-Duque, J. C., et al. (2016). Admixture in Latin America. Curr. Opin. Genet. Dev. 41, 106–114. doi: 10.1016/j.gde.2016.09.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Akiyama, M., Ishigaki, K., Sakaue, S., Momozawa, Y., Horikoshi, M., Hirata, M., et al. (2019). Characterizing rare and low-frequency height-associated variants in the Japanese population. Nat. Commun. 10:4393. doi: 10.1038/s41467-019-12276-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Anderson, C., Pettersson, F., Clarke, G., Cardon, L., Morris, A., and Zondervan, K. (2010). Data quality control in genetic case-control association studies. Nat. Protoc. 5, 1564–1573. doi: 10.1038/nprot.2010.116

PubMed Abstract | CrossRef Full Text | Google Scholar

Berg, A. T., Berkovic, S. F., Brodie, M. J., Buchhalter, J., Cross, J. H., Van Emde Boas, W., et al. (2010). Revised terminology and concepts for organization of seizures and epilepsies: report of the ILAE Commission on Classification and Terminology, 2005-2009. Epilepsia 51, 676–685. doi: 10.1111/j.1528-1167.2010.02522.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Browning, B. L., Zhou, Y., and Browning, S. R. (2018). A One-Penny Imputed Genome from Next-Generation Reference Panels. Am. J. Hum. Genet. 103, 338–348. doi: 10.1016/j.ajhg.2018.07.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, M. H., Raffield, L. M., Mousas, A., Sakaue, S., Huffman, J. E., Moscati, A., et al. (2020). Trans-ethnic and Ancestry-Specific Blood-Cell Genetics in 746,667 Individuals from 5 Global Populations. Cell 182, 1198–1213.e14. doi: 10.1016/j.cell.2020.06.045

PubMed Abstract | CrossRef Full Text | Google Scholar

Das, S., Forer, L., Schönherr, S., Sidore, C., Locke, A. E., Kwong, A., et al. (2016). Next-generation genotype imputation service and methods. Nat. Genet. 48, 1284–1287. doi: 10.1038/ng.3656

PubMed Abstract | CrossRef Full Text | Google Scholar

EPICURE Consortium, EMINet Consortium, Steffens, M., Leu, C., Ruppert, A. K., Zara, F., et al. (2012a). Genome-wide association analysis of genetic generalized epilepsies implicates susceptibility loci at 1q43, 2p16.1, 2q22.3 and 17q21.32. Hum. Mol. Genet. 21, 5359–5372. doi: 10.1093/hmg/dds373

PubMed Abstract | CrossRef Full Text | Google Scholar

EPICURE Consortium, Leu, C., de Kovel, C., Zara, F., Striano, S., Pezzella, M., et al. (2012b). Genome-wide linkage meta-analysis identifies susceptibility loci at 2q34 and 13q31.3 for genetic generalized epilepsies. Epilepsia 53, 308–318. doi: 10.1111/j.1528-1167.2011.03379.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Excoffier, L., Smouse, P. E., and Quattro, J. M. (1992). Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics 131, 479–491. doi: 10.5962/bhl.title.86657

PubMed Abstract | CrossRef Full Text | Google Scholar

Fisher, R. S., Acevedo, C., Arzimanoglou, A., Bogacz, A., Cross, J. H., Elger, C. E., et al. (2014). ILAE Official Report: a practical clinical definition of epilepsy. Epilepsia 55, 475–82. doi: 10.1111/epi.12550

PubMed Abstract | CrossRef Full Text | Google Scholar

Graff, M., Justice, A. E., Young, K. L., Marouli, E., Zhang, X., Fine, R. S., et al. (2021). Discovery and fine-mapping of height loci via high-density imputation of GWASs in individuals of African ancestry. Am. J. Hum. Genet. 108, 564–582. doi: 10.1016/j.ajhg.2021.02.011

PubMed Abstract | CrossRef Full Text | Google Scholar

International League Against Epilepsy Consortium on Complex Epilepsies (ILAE Consortium on Complex Epilepsies) (2014). Genetic determinants of common epilepsies: a meta-analysis of genome-wide association studies. Lancet Neurol. 13, 893–903. doi: 10.1016/S1474-4422(14)70171-1

CrossRef Full Text | Google Scholar

International League Against Epilepsy Consortium on Complex Epilepsies (ILAE Consortium on Complex Epilepsies). (2018). Genome-wide mega-analysis identifies 16 loci and highlights diverse biological mechanisms in the common epilepsies. Nat. Commun. 9:5269. doi: 10.1038/s41467-018-07524-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Karczewski, K. J., Francioli, L. C., Tiao, G., Cummings, B. B., Alföldi, J., Wang, Q., et al. (2020). The mutational constraint spectrum quantified from variation in 141,456 humans. Nature 581, 434–443. doi: 10.1038/s41586-020-2308-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Kehdy, F. S. G., Gouveia, M. H., Machado, M., Magalhães, W. C. S., Horimoto, A. R., Horta, B. L., et al. (2015). Origin and dynamics of admixture in Brazilians and its effect on the pattern of deleterious mutations. Proc. Natl. Acad. Sci.U. S. A. 112, 8696–8701. doi: 10.1073/pnas.1504447112

PubMed Abstract | CrossRef Full Text | Google Scholar

Kowalski, M. H., Qian, H., Hou, Z., Rosen, J. D., Tapia, A. L., Shan, Y., et al. (2019). Use of >100,000 NHLBI Trans-Omics for Precision Medicine (TOPMed) Consortium whole genome sequences improves imputation quality and detection of rare variant associations in admixed African and Hispanic/Latino populations. PLoS Genet. 15:e1008500. doi: 10.1371/journal.pgen.1008500

PubMed Abstract | CrossRef Full Text | Google Scholar

Magalhães, W. C. S., Araujo, N. M., Leal, T. P., Araujo, G. S., Viriato, P. J. S., Kehdy, F. S., et al. (2018). EPIGEN-Brazil Initiative resources: a Latin American imputation panel and the Scientific Workflow. Genome Res. 28, 1090–1095. doi: 10.1101/gr.225458.117

PubMed Abstract | CrossRef Full Text | Google Scholar

Marchini, J., Cardon, L. R., Phillips, M. S., and Donnelly, P. (2004). The effects of human population structure on large genetic association studies. Nat. Genet. 36, 512–517. doi: 10.1038/ng1337

PubMed Abstract | CrossRef Full Text | Google Scholar

Martin, A. R., Gignoux, C. R., Walters, R. K., Wojcik, G. L., Neale, B. M., Gravel, S., et al. (2017). Human Demographic History Impacts Genetic Risk Prediction across Diverse Populations. Am. J. Hum. Genet. 100, 635–649. doi: 10.1016/j.ajhg.2017.03.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Martin, A. R., Kanai, M., Kamatani, Y., Okada, Y., Neale, B. M., and Daly, M. J. (2019). Clinical use of current polygenic risk scores may exacerbate health disparities. Nat. Genet. 51, 584–591. doi: 10.1038/s41588-019-0379-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Moore, C. M., Jacobson, S. A., and Fingerlin, T. E. (2020). Power and Sample Size Calculations for Genetic Association Studies in the Presence of Genetic Model Misspecification. Hum. Hered. 84, 1–16. doi: 10.1159/000508558

PubMed Abstract | CrossRef Full Text | Google Scholar

Moura, R. R., de, Coelho, A. V. C., de Queiroz Balbino, V., Crovella, S., Brandão, L. A. C., et al. (2015). Meta-analysis of Brazilian genetic admixture and comparison with other Latin America countries. Am. J. Hum. Biol. 27, 674–680. doi: 10.1002/ajhb.22714

PubMed Abstract | CrossRef Full Text | Google Scholar

O’Connell, J., Gurdasani, D., Delaneau, O., Pirastu, N., Ulivi, S., Cocca, M., et al. (2014). A General Approach for Haplotype Phasing across the Full Spectrum of Relatedness. PLoS Genet. 10:e1004234. doi: 10.1371/journal.pgen.1004234

PubMed Abstract | CrossRef Full Text | Google Scholar

Price, A. L., Zaitlen, N. A., Reich, D., and Patterson, N. (2010). New approaches to population stratification in genome-wide association studies. Nat. Rev. Genet. 11, 459–463. doi: 10.1038/nrg2813

PubMed Abstract | CrossRef Full Text | Google Scholar

Purcell, S., Neale, B., Todd-Brown, K., Thomas, L., Ferreira, M. A. R., Bender, D., et al. (2007). PLINK: a Tool Set for Whole-Genome Association and Population-Based Linkage Analyses. Am. J. Hum. Genet. 81, 559–575. doi: 10.1086/519795

PubMed Abstract | CrossRef Full Text | Google Scholar

Rocha, C. S., Secolin, R., Rodrigues, M. R., Carvalho, B. S., and Lopes-Cendes, I. (2020). The Brazilian Initiative on Precision Medicine (BIPMed): fostering genomic data-sharing of underrepresented populations. NPJ Genom. Med. 5:42. doi: 10.1038/s41525-020-00149-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Scheffer, I. E., Berkovic, S., Capovilla, G., Connolly, M. B., French, J., Guilhoto, L., et al. (2017). ILAE classification of the epilepsies: position paper of the ILAE Commission for Classification and Terminology. Epilepsia 58, 512–521. doi: 10.1111/epi.13709

PubMed Abstract | CrossRef Full Text | Google Scholar

Scheffer, I. E., Berkovic, S., Capovilla, G., Connolly, M. B., French, J., Guilhoto, L., et al. (2018). ILAE classification of the epilepsies: position paper of the ILAE Commission for Classification and Terminology. Z. Epileptol. 31, 296–306. doi: 10.1007/s10309-018-0218-6

CrossRef Full Text | Google Scholar

Secolin, R., Mas-Sandoval, A., Arauna, L. R., Torres, F. R., Araujo, T. K., Santos, M. L., et al. (2019). Distribution of local ancestry and evidence of adaptation in admixed populations. Sci. Rep. 9:13900. doi: 10.1038/s41598-019-50362-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Turner, S. (2014). qqman: an R package for visualizing GWAS results using Q-Q and manhattan plots. BioRxiv [Preprint]. doi: 10.1101/005165

CrossRef Full Text | Google Scholar

Wang, M., Greenberg, D. A., and Stewart, W. C. L. (2019). Replication, reanalysis, and gene expression: ME2 and genetic generalized epilepsy. Epilepsia 60, 539–546. doi: 10.1111/epi.14654

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Qu, J., Mao, C. X., Wang, Z. B., Mao, X. Y., Zhou, B. T., et al. (2014). Novel Susceptibility Loci were Found in Chinese Genetic Generalized Epileptic Patients by Genome-wide Association Study. CNS Neurosci. Ther. 20, 1008–1010. doi: 10.1111/cns.12328

PubMed Abstract | CrossRef Full Text | Google Scholar