Exploring a Region on Chromosome 8p23.1 Displaying Positive Selection Signals in Brazilian Admixed Populations: Additional Insights Into Predisposition to Obesity and Related Disorders

Abstract

We recently reported a deviation of local ancestry on the chromosome (ch) 8p23.1, which led to positive selection signals in a Brazilian population sample. The deviation suggested that the genetic variability of candidate genes located on ch 8p23.1 may have been evolutionarily advantageous in the early stages of the admixture process. In the present work, we aim to extend the previous work by studying additional Brazilian admixed individuals and examining DNA sequencing data from the ch 8p23.1 candidate region. Thus, we inferred the local ancestry of 125 exomes from individuals born in five towns within the Southeast region of Brazil (São Paulo, Campinas, Barretos, and Ribeirão Preto located in the state of São Paulo and Belo Horizonte, the capital of the state of Minas Gerais), and compared to data from two public Brazilian reference genomic databases, BIPMed and ABraOM, and with information from the 1000 Genomes Project phase 3 and gnomAD databases. Our results revealed that ancestry is similar among individuals born in the five Brazilian towns assessed; however, an increased proportion of sub-Saharan African ancestry was observed in individuals from Belo Horizonte. In addition, individuals from the five towns considered, as well as those from the ABRAOM dataset, had the same overrepresentation of Native-American ancestry on the ch 8p23.1 locus that was previously reported for the BIPMed reference sample. Sequencing analysis of ch 8p23.1 revealed the presence of 442 non-synonymous variants, including frameshift, inframe deletion, start loss, stop gain, stop loss, and splicing site variants, which occurred in 24 genes. Among these genes, 13 were associated with obesity, type II diabetes, lipid levels, and waist circumference (PRAG1, MFHAS1, PPP1R3B, TNKS, MSRA, PRSS55, RP1L1, PINX1, MTMR9, FAM167A, BLK, GATA4, and CTSB). These results strengthen the hypothesis that a set of variants located on ch 8p23.1 that result from positive selection during early admixture events may influence obesity-related disease predisposition in admixed individuals of the Brazilian population. Furthermore, we present evidence that the exploration of local ancestry deviation in admixed individuals may provide information with the potential to be translated into health care improvement.

Introduction

Admixture between different continental populations generates mosaic chromosomes comprised of genomic segments with different ancestry, which is defined as local ancestry (Seldin et al., 2011). As a result, admixed populations may present marked differences in local ancestry patterns (Browning et al., 2016; Deng et al., 2016; Martin et al., 2017; Secolin et al., 2019). These differences may impact disease incidence and genetic risk prediction across populations (Myles et al., 2008; Moonesinghe et al., 2012; Martin et al., 2017). Thus, enhancing our knowledge of the effect of local ancestry is crucial for the development of adequate precision health programs in admixed populations (Aronson and Rehm, 2015; Hindorff et al., 2018).

The Brazilian population was formed via an admixture process comprised mostly of European, sub-Saharan African, and Native-American population ancestry. In terms of global ancestry inference, studies have shown a predominance of European ancestry, followed by sub-Saharan African and Native-American (Kehdy et al., 2015; Lima-Costa et al., 2015; Rodrigues de Moura et al., 2015). Furthermore, a recent study about local ancestry inferences reported that the Native-American component predominated on the chromosome (ch) 8p23.1 due to positive selection (Secolin et al., 2019) (3). Ch 8p23.1 has undergone inversion events stratified across continental populations (Salm et al., 2012), which may influence the recombination landscape (Alves et al., 2014).

Interestingly, the ch 8p23.1 region found to be under positive selection in the Brazilian population has been reported to contain genes previously associated with type 2 diabetes and overweight/obesity in admixed Americans (Dunn et al., 2006; Flores et al., 2016). Indeed, studies taking admixture into account have shown that type 2 diabetes, insulin secretion, body mass index, obesity, and adiposity are the main clinical phenotypes associated with metabolic disorders (Dunn et al., 2006; Hayes et al., 2013; Goetz et al., 2014; Flores et al., 2016; Mehta et al., 2017). Thus, we hypothesize that variants in candidate genes located on ch 8p23.1 could have provided an evolutionary advantage in a restrictive diet environment in the early stages of the Brazilian admixture. However, in the present high caloric diet environment, this genetic variability can result in an increased number of obesity-related traits in admixed Brazilian individuals.

Therefore, our objective was to expand our knowledge of the effects of admixture by describing the genetic variability of ch 8p23.1 from admixed Brazilian exomes compared with global populations. To achieve our aim, we first extended our study to additional admixed exomes from other southeastern Brazil towns. Second, we identified and analyzed sequencing variants identified in the candidate region of ch 8p23.1.

Materials and Methods

Subjects

We evaluated 257 individuals from BIPMed (Rocha et al., 2020), 609 from ABraOM (Naslavsky et al., 2017), and 88 additional exomes from individuals born in the following towns within southeastern Brazil: Barretos (N = 30); Ribeirão Preto (N = 30), located in the state of São Paulo; and Belo Horizonte (N = 28), the capital of the state of Minas Gerais (Supplementary Figure 1). Among individuals included in BIPMed, the birthplace of 193 individuals were included; thus, we were able to extract 21 individuals born in São Paulo city and 37 from Campinas to increase the power of regional comparisons. No information regarding place of birth was obtained from the ABraOM dataset. Permission to use raw, anonymized data from BIPMed and ABraOM public databases and raw, anonymized data associated with the 88 exomes of individuals from Barretos, Ribeirão Preto, and Belo Horizonte was obtained. This study was approved by the University of Campinas’s Research Ethics Committee (UNICAMP, Campinas, São Paulo, Brazil). All methods were performed following relevant guidelines and regulations.

Exome Processing

Exome data were stored in variant call format (VCF) files created using the GRCh37 assembly. We used PLINK 1.9 (Purcell et al., 2007) software to convert VCF to PLINK files, variant and individual filtering, and data merging (Anderson et al., 2010). First, we removed ambiguous variants (with G/C or A/T alleles) from VCFs associated with each town, BIPMed, and ABraOM. Next, we merged all Brazilian VCFs, maintaining only biallelic variants, autosomal variants, variants in Hardy–Weinberg equilibrium (Anderson et al., 2010), and removal of missing data (> 10%). These filters were used only to analyze population structure and local ancestry and were removed in the analysis to identify variants in the candidate region at ch 8p23.1.

We evaluated the heterozygosity rate of each individual to search for inbreeding (low heterozygosity rate) or sample contamination (high heterozygosity rate) (Anderson et al., 2010), and individuals with a heterozygosity rate higher or lower than three standard deviations (SDs) from the mean were removed. We also removed individuals with genomic relatedness matrix estimations higher than 0.125, which is the expected genomic relatedness of third-degree relatives (Anderson et al., 2010). The genomic relatedness matrix estimation used a greedy algorithm implemented using the PLINK 1.9 software to maximize the sample size.

After genotype and individual filtering, a total of 893 exomes and 661,617 variants remained in the Brazilian datasets analyzed. We merged this dataset with the 1000 Genome project data phase 3 (1KGP) (1000 Genomes Project Consortium et al., 2015) and removed SNPs with a minor allele frequency (MAF) < 0.01. As a result, 225,997 variants for local ancestry inference were used. We also removed variants in linkage disequilibrium (LD) from the MAF-filtered dataset (parameters: window size = 50 SNPs; shift step = 5 SNPs; and r2 = 0.5) (Anderson et al., 2010), which left 127,172 SNPs for an investigation of population structure.

Population Structure

To evaluate whether our Brazilian sample (BRS) presents a geographical substructure based on birthplace, we performed the analysis of molecular variance (AMOVA) using the poppr.amova package in R software (Excoffier et al., 1992), which compares the genetic distance among birthplace/town groups based on a set of 10,000 random SNPs across the genome. In addition, we compared the BRS data classified by birthplace to the 1KGP dataset via principal component analysis (PCA) using PLINK v1.9 software to evaluate the presence of population-based outliers in the BRS dataset.

Local Ancestry Inference and Positive Selection Test

We phased SNPs without LD pruning using the SHAPEIT2 v2.r387 software with default parameters (O’Connell et al., 2014). After phasing, we converted the output data from SHAPEIT2 to input files required by RFMix v.1.5.4 software (Maples et al., 2013) using a pipeline previously reported¹ (Martin et al., 2017).

Previous studies showed that using Peruvian individuals from the 1KGP with a high degree of Native-American ancestry as a Native-American reference produced the same result as using Native-American indigenous individuals (Secolin et al., 2019). Therefore, we inferred the local ancestry of 23 Peruvian individuals who possessed a > 0.95 proportion of Native-American ancestry (NAT) (Secolin et al., 2019), 23 random Europeans (EUR), and 23 random sub-Saharan Africans from the 1KGP (AFR). The size sample of ancestry references was selected based on the 23 NAT to avoid biases due to unbalanced reference panel sizes of ancestry references, according to the RFMix v.1.5.4. Manual (Maples et al., 2013). We ran RFMix in PopPhased mode with a minimum window size of 0.2 cM, using one EM iteration and node size 5. The reference panel was maintained after the initial inference step, and forward-backward probabilities were saved. We analyzed the proportion of EUR, AFR, and NAT ancestry in the BRS dataset for each variant across the genome using in-house-developed R scripts (Secolin et al., 2019), and results were plotted using the man package in R software (Turner, 2014). In order to evaluate the presence of ch 8p32.1 inversions, we performed an inversion inference using the invClust package in R software (Cáceres and González, 2015), as performed in our previous work (Secolin et al., 2019). Since we have individuals that overlap the previous paper, we decided to compare the inversion inference between the SNP array data (Secolin et al., 2019) and the exome data from the same individuals to evaluate whether the inversion analysis generated a perfect match.

We tested our exome sample for positive selection by the same approach used previously (Patin et al., 2017). Briefly, this approach combines the results of five neutrality statistics (intrapopulation absolute integrated haplotype scores (| iHS|, | ΔiHH|) (Voight et al., 2006; Sabeti et al., 2007), interpopulation integrated haplotype score (| ΔiHH_derived|) (Grossman et al., 2010), interpopulation extended haplotype homozygosity (XP-EHH) (Sabeti et al., 2007), and population branch statistics (PBS) (Yi et al., 2010) based on Hudson’s Fst (Bhatia et al., 2013) in a single Fisher combined score (FCS) (Deschamps et al., 2016). The variants with values of FCS higher than 99% of the SNP FCS values across the genome (i.e., the 1% highest FCS values) were defined as outliers. Then, we split the genome into 100-variants blocks. Finally, we estimated the proportion of outliers within each block. If a block presents a proportion of outliers higher than the 99.5th percentile (the highest 0.5%), it was defined as a region under positive selection.

Analysis of Chromosome 8p23.1

We extracted the ch 8p23.1 region (8092025–11859740 bp) (Secolin et al., 2019) from the VCF file of each sample using vcftools (Danecek et al., 2011). Variant consequences from each gene region were annotated using the ANNOVAR software (version 2019Oct24) (Wang et al., 2010) with the following flags: -other info (to include our sample AF); -one transcript; -buildver hg19; -remove; -protocol refGene, gnomad211_exome, ALL.sites.2015_08,EUR.sites.2015_ 08,AFR.sites.2015_08,AMR.sites.2015_08,EAS.sites.2015_08,SAS. sites.2015_08,dbnsfp35a; -operation g,f,f,f,f,f,f,f,f; and -nastring.

We included the allele frequency (AF) information from African/African-American (AFR/AFA), Latino/admixed American (LAT/AMR), East Asian (EAS), non-Finish European (NFE), and South Asian (SAS) populations from the gnomAD exome dataset (Karczewski et al., 2020); sub-Saharan African (AFR), Europeans (EUR), admixed Americans (AMR), East Asians (EAS), and South Asians (SAS) from 1KGP, which are publicly available in ANNOVAR resource data. In addition, we annotated variants that were not identified via ANNOVAR using the Variant Effect Prediction (VEP) algorithm (McLaren et al., 2016) with the following parameters: a buffer_size 500; –canonical; –distance 5000; –regulatory; –species homo_sapiens; –symbol.

To predict the impact of non-synonymous variants identified on protein function, we analyzed the information provided by the use of 12 algorithms, which included PolyPhen2 (Adzhubei et al., 2013), Sort Intolerant from Tolerant (SIFT) (Sim et al., 2012), MutationTaster (Schwarz et al., 2010), PROVEAN (Choi et al., 2012), Combined Annotation Dependent Depletion (CADD) (Rentzsch et al., 2019), MutPred². Functional Analysis through Hidden Markov Models (FATHMM) (Shihab et al., 2015), PhD-SNPg (Capriotti and Fariselli, 2017), Condel (González-Pérez and López-Bigas, 2011), PANTHER (Mi et al., 2013), Align Grantham Variation/Grantham Difference score (GVGD) (Tavtigian et al., 2006), and SNPs&GO (Calabrese et al., 2009). For Align-GVGD, variants graded higher than C35 were classified as deleterious. For MutPred2, variants with a score higher than 0.5 were considered pathogenic. For all other algorithms, we used default classifications.

Associated trait information for genes located on ch 8p23.1 was accessed from the NHGRI-EBI GWAS Catalog (Buniello et al., 2019) on October 30, 2020, and results were complemented by a search of the PubMed^® database.

Results

Population Structure

AMOVA results obtained using 10,000 random SNPs showed that 99.21% of observed variation occurred within groups and 0.79% occurred among groups (total φ-statistics = 0.0079; p = 0.0001), indicating the absence of a population substructure, which is consistent with the lack of clusters observed in the PCA of the BRS sample that was based on birthplace (Figure 1). In addition, the global PCA of the 1KGP dataset (Supplementary Figure 2) showed that our sample consisted of a mixture of European, sub-Saharan African, and Native-American/East Asian individuals, similar to other admixed American populations. However, the population was distributed mainly between Europeans and sub-Saharans rather than Native-Americans, consistent with previous studies (Ruiz-Linares et al., 2014; Kehdy et al., 2015; Secolin et al., 2019; Rocha et al., 2020).

FIGURE 1

Local Ancestry Inference

The proportion of the BRS sample, which included individuals born in the different towns and two public datasets, had an average local ancestry proportion for its EUR component of 74.6% (SD = 1.4%). The proportion of the sample that comprised the AFR component was 16.0% (SD = 1.1%), and the NAT component was 9.4% (SD = 1.1%) (Figure 2A and Supplementary Figure 3). We observed differences in EUR and AFR ancestry proportions among towns, with the Belo Horizonte population containing the lowest EUR component (mean = 66.8%; SD = 6.0%), and the highest AFR component (mean = 26.9%; SD = 5.6%). São Paulo, in contrast, had the greatest proportion of EUR ancestry (mean = 87.8%; SD = 5.7%) and the lowest AFR proportion (mean = 6.2%; SD = 4.1%). The NAT component of the sample remained constant among individuals from the different towns and the two public Brazilian databases, ranging from a mean of 4.6% (SD = 2.8%) in Barretos to 8.2% (SD = 1.3%) in the BIPMed sample (Figure 2A). Our assessment revealed a decreased EUR component on ch 8p23.1, and an elevated NAT in individuals born in Campinas, São Paulo, Barretos, and Belo Horizonte, as well as in those included in the ABraOM dataset (Figure 2B and Supplementary Figures 3–9).

FIGURE 2

We inferred the inversion events on ch 8p23.1 from the exome data by the same approach used in our previous report (Secolin et al., 2019), with the invClust package in R (Cáceres and González, 2015). The results showed that 48.9% of the inferred inversions in the exome data matched the results previously obtained with the SNP-array dataset of the BIPMed sample used in our previous work (Secolin et al., 2019; Supplementary Table 1).

We tested for positive selection in the exome dataset using the Fisher combined scores (FCS). FCS, which includes PBS tests in the calculation, this is the same approach used in our previous work (Secolin et al., 2019). However, the results did not recapture the same positive selection signal on ch 8p23.1 previously observed (Secolin et al., 2019) (Supplementary Figures 3–10).

Analysis of Chromosome 8p23.1

We found 17,536 variants within ch 8p23.1. We focused on the following variants with the potential to impact gene function: 414 non-synonymous variants, ten frameshifts, eight inframe deletions, one start loss, five stop gains, one stop loss, and five splicing sites. The variants affected 24 genes and two open reading frames (Supplementary Data Sheet 1). Among these variants, 355 were also found in gnomAD and/or 1KGP databases, and 44 such variants were determined to be common with an alternative allele frequency (AAF) > 0.01 in the BRS dataset but rare (AAF < 0.01) in gnomAD and 1KGP (Table 1). Also, we identified nine common variants (AAF > 0.01) among the 89 variants exclusive to the Brazilian population. The AF comparison of these 89 variants separated by Brazilian cities and datasets showed that the RP1L1 gene in the ABraOM database contained the largest number of exclusive Brazilian variants (Supplementary Figure 11).

We observed that 374 of the 414 non-synonymous variants, in genes located at ch 8p23.1, were classified as deleterious via in silico prediction of at least one algorithm (Supplementary Data Sheet 2), and 19 of these were predicted to be pathogenic with an 80% concordance among the different algorithms; these were present in five genes PRSS55, RP1L1, SOX7, GATA4, and CTSB. As shown in Figure 3 and Supplementary Table 2, we found 167 non-synonymous variants predicted to be benign by at least one algorithm; these were present in 16 different genes. Also, there were 309 variants predicted to be deleterious and 50 variants predicted to be benign when considering less than 20% concordance among the algorithms.

FIGURE 3

Interestingly, 140 of the variants predicted to be benign by at least one algorithm (140/167) were found in 13 genes, which were previously associated with metabolic phenotypes such as type 1 diabetes mellitus (T1DM), type 2 diabetes mellitus (T2DM), obesity, insulin resistance, body mass index (BMI), body fat distribution, waist circumference, and diet measurement (MFHAS1, ERI1, TNKS, PRSS55, RP1L1, PINX1, XKR6, FAM167A, BLK, GATA4, and CTSB genes), Table 2. Also, 45 of these variants were located in the following eight genes with an AAF > 0.01, considering the BRSs and gnomAD/1KGP databases: MFHAS1, ERI1, PRSS55, RP1L1, PINX1, FAM167A, GATA4, and CTSB (Supplementary Data Sheet 2).

Discussion

The sequencing analysis of ch 8p23.1 performed in the current study revealed the presence of 442 non-synonymous variants, including frameshift, inframe deletion, start loss, stop gain, stop loss, and splicing site variants, which occurred in 24 genes and two open reading frames. Among the genes, 13 were associated with obesity, type II diabetes, lipid levels, and waist circumference (PRAG1, MFHAS1, PPP1R3B, TNKS, MSRA, PRSS55, RP1L1, PINX1, MTMR9, FAM167A, BLK, GATA4, and CTSB).

The inversion event on ch 8p23.1 generated a large haplotype, which was able to be traced through continental populations globally (Salm et al., 2012). It presented us with an opportunity to investigate how admixture events and evolutionary processes have affected variability within non-inverted and inverted haplotypes in admixed populations. Previously, two independent studies reported local ancestry deviation on ch 8p23.1 in admixed American populations, likely due to inversion events (Guan, 2014; Secolin et al., 2019). Furthermore, our own work using SNP-array data demonstrated that the proportion of non-inverted haplotypes inherited from Native-Americans is higher than those inherited from Europeans in admixed Brazilian individuals (Secolin et al., 2019). Here, we replicated these findings using exome datasets in populations originating from an extended geographic region in the southeastern region of Brazil. Besides, since we evaluated individuals with unknown information regarding the presence of obesity-related disorders, our study is not biased toward a specific phenotype, and it is suitable to assess the genetic variability of the candidate region on ch 8p23.1.

We observed that the results from the inversion inference on ch 8p23.1 obtained in the present work, using the exome data, did not completely match that resulted from the analysis using the SNP-array dataset (Secolin et al., 2019). However, it is noteworthy that the invClust package was developed to be used with SNP-array data, and to our knowledge, there is no reference to its use with exome datasets. Thus, it is possible that inferences of chromosomal inversions using exome data may not be accurate with the invClust package. As pointed out in our previous work (Secolin et al., 2019), it is not likely that inversion bias would influence the high NAT proportions observed in the sample. However, we agree that this is a limitation of our current work. Further analysis, in which inverted and non-inverted genotypes are unequivocally identified, would help evaluate the distribution of the inflation in NAT ancestry in inverted and non-inverted genotypes.

There is evidence that the deviation towards Native-American ancestry on ch 8p23.1 could be due to positive selection events after the Brazilian admixture (Secolin et al., 2019). Indeed, previous studies suggested that Native-American ancestry was admixed early in the European colonization in Brazil (approximately 18 to 16 generations ago) and was followed by the posterior depletion of NAT (Kehdy et al., 2015). This early admixture could have catalyzed positive selection events among the first admixed Brazilian individuals. Although environmental causes that drove this positive selection remain unknown, studies had identified variants associated with type 2 diabetes mellitus, insulin secretion, body mass index, obesity, and adiposity, when admixture was considered (Dunn et al., 2006; Hayes et al., 2013; Goetz et al., 2014; Flores et al., 2016; Mehta et al., 2017). Therefore, the large number of genes located on ch 8p23.1 related to diet and metabolic traits suggest that positive selection may have occurred due to the restrictive diet environment and severe famine periods in early admixed Brazilian individuals (Davis, 2001).

In the present work, our results did not recapture the same positive selection signal detected previously (Secolin et al., 2019) (Supplementary Figures 3–10). However, since FCS has only been used with whole-genome sequencing (Deschamps et al., 2016; Patin et al., 2017) and SNP-array datasets (Secolin et al., 2019), we believe that the decrease in genetic variability present in exome data may render FCS less suitable for this type of analysis.

Increasing fat and glucose storage could increase body fat, glucose storage, and obesity-related diseases in individuals who eat a fat and glucose-rich diet today; findings consistent with previous association studies (Pulit et al., 2019; Rask-Andersen et al., 2019; Schlauch et al., 2020). Indeed, we identified 89 variants with the potential to impact gene function that were found exclusively in the admixed Brazilian sample. Unfortunately, we cannot define the correct phase for the allele variants and the ancestry block by the RFMix algorithm. However, we know that the 89 variants are in the region, presenting 60.69% of EUR ancestry proportion, followed by 15.47% of AFR and 23.83% of NAT ancestry proportions. The comparison of these proportions with the average EUR ancestry proportion among Brazilian genomes genomes (74.6%), AFR (16.0%), and NAT (9.4%) suggests that these variants present exclusively in the Brazilian samples could, most likely, be the main contributors to the signals of selection identified, and are possibly influencing obesity-related phenotypes. Therefore, we consider the region of ch 8p23.1 a hotspot for genetic variants that predispose individuals to obesity disorders. It may be useful, as a first strategy, to concentrate efforts on studying effects of non-synonymous variants identified within the 13 candidate genes of the region, PRAG1, MFHAS1, PPP1R3B, TNKS, MSRA, PRSS55, RP1L1, PINX1, MTMR9, FAM167A, BLK, GATA4, and CTSB. It may also be useful to expand genetic studies to include patients with obesity-related phenotypes and studying the expression levels of candidate genes in relevant tissue may also give additional clues regarding their roles in disease-related phenotypes.

In Table 2, we present the list of 17 genes that have been linked to diet patterns in large association studies and are located in the candidate region on ch 8p23.1. Seven of these large studies included Hispanic, Native-American, and Caribbean populations (Comuzzie et al., 2012; Divers et al., 2017; Pulit et al., 2019; Rask-Andersen et al., 2019; Niarchou et al., 2020; Schlauch et al., 2020; Vujkovic et al., 2020), and three contained association signals in the ERI1 gene (Pulit et al., 2019; Rask-Andersen et al., 2019; Schlauch et al., 2020), which is located within the region and was determined to possess the greatest degree of positive selection in admixed Brazilians in our previous work (Secolin et al., 2019).

Finally, it is also important to study the non-synonymous variants identified in the candidate genes on ch 8p23 and predicted to be benign. These variants were identified in 11 of the candidate genes listed in Table 2. Currently, we cannot exclude the possibility that even though these variants are not predicted to affect protein function individually, they may contribute to a polygenic phenotype.

Furthermore, when considering a polygenic phenotype, one aspect that we should take into account is the presence of epistatic interactions. Thus, we could argue that an increase in the frequency of the genes of NAT ancestry on ch 8p23.1 could be due to the breakup of negative epistatic interactions among genes on other regions from NAT genomes and the genes on ch 8p23.1, which are currently coupled with AFR and EUR ancestry tracts, and could lead to increased fitness. We count the number of AFR-NAT-AFR, EUR-NAT-EUR, and NAT-NAT-NAT haplotypes, including ch 8p23.1 and adjacent regions (approximately 3.7Mb upstream and downstream ch 8p23.1). However, we did not observe an overrepresentation of AFR-NAT-AFR (n = 32) or EUR-NAT-EUR (n = 98) ancestry haplotypes compared to NAT-NAT-NAT haplotypes (n = 118). Therefore, our results suggest that interactions among EUR or AFR ancestry genes in adjacent regions on ch 8p23.1 with the NAT ancestry core seems not to be enough to boost negative or positive selection in our sample. However, gene interactions can occur among genes on ch 8p23.1 and genes on other regions of the genome, and further studies should be performed to clarify this issue.

Conclusion

We successfully replicated previous results that identified local ancestry deviation on ch 8p23.1, which seems to have occurred in populations from the southeastern region of Brazil, including the states of São Paulo and Minas Gerais. Thus, the candidate region on ch 8p23.1 emerges as a hotspot for obesity-related genes in admixed Brazilians, which should be further explored. In particular, the information presented here could be used in the future to support risk stratification and implement personalized public health policies and preventive medical treatments.

References

• 1

  1000 Genomes Project Consortium, AutonA.BrooksL. D.DurbinR. M.GarrisonE. P.KangH. M.et al (2015). A global reference for human genetic variation.Nature52668–74. 10.1038/nature15393

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  AdzhubeiI.JordanD. M.SunyaevS. R. (2013). Predicting functional effect of human missense mutations using PolyPhen-2.Curr. Protoc. Hum. Genet.767.20.1–7.20.41. 10.1002/0471142905.hg0720s76

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 3

  AlvesJ. M.ChikhiL.AmorimA.LopesA. M. (2014). The 8p23 inversion polymorphism determines local recombination heterogeneity across human populations.Genome Biol. Evol.6921–930. 10.1093/gbe/evu064

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 4

  AndersonC.PetterssonF.ClarkeG.CardonL.MorrisA.ZondervanK. (2010). Data quality control in genetic case-control association studies.Nat. Protoc.51564–1573. 10.1038/nprot.2010.116

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 5

  AndreasenC. H.MogensenM. S.Borch-JohnsenK.SandbaekA.LauritzenT.AlmindK.et al (2009). Studies of CTNNBL1 and FDFT1 variants and measures of obesity: analyses of quantitative traits and case-control studies in 18,014 Danes.BMC Med. Genet.10:17. 10.1186/1471-2350-10-17

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 6

  AronsonS. J.RehmH. L. (2015). Building the foundation for genomics in precision medicine.Nature526336–342. 10.1038/nature15816

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  BhatiaG.PattersonN.SankararamanS.PriceA. L. (2013). Estimating and interpreting FST: the impact of rare variants.Genome Res.231514–1521. 10.1101/gr.154831.113

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 8

  BorowiecM.LiewC. W.ThompsonR.BoonyasrisawatW.HuJ.MlynarskiW. M.et al (2009). Mutations at the BLK locus linked to maturity onset diabetes of the young and beta-cell dysfunction.Proc. Natl. Acad. Sci. U.S.A.10614460–14465. 10.1073/pnas.0906474106

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 9

  BrowningS. R.GrindeK.PlantingaA.GogartenS. M.StilpA. M.KaplanR. C.et al (2016). Local ancestry inference in a large US-Based Hispanic/Latino study: hispanic community health study/study of latinos (HCHS/SOL).G3 (Bethesda)61525–1534. 10.1534/g3.116.028779

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  BunielloA.MacArthurJ. A. L.CerezoM.HarrisL. W.HayhurstJ.MalangoneC.et al (2019). The NHGRI-EBI GWAS catalog of published genome-wide association studies, targeted arrays and summary statistics 2019.Nucleic Acids Res.47D1005–D1012. 10.1093/nar/gky1120

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 11

  CáceresA.GonzálezJ. R. (2015). Following the footprints of polymorphic inversions on SNP data: from detection to association tests.Nucleic Acids Res.43:e53. 10.1093/nar/gkv073

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 12

  CalabreseR.CapriottiE.FariselliP.MartelliP. L.CasadioR. (2009). Functional annotations improve the predictive score of human disease-related mutations in proteins.Hum. Mutat.301237–1244. 10.1002/humu.21047

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  CapriottiE.FariselliP. (2017). PhD-SNPg: a webserver and lightweight tool for scoring single nucleotide variants.Nucleic Acids Res.45W247–W252. 10.1093/nar/gkx369

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  ChalasaniN.GuoX.LoombaR.GoodarziM. O.HarituniansT.KwonS.et al (2010). Genome-wide association study identifies variants associated with histologic features of nonalcoholic Fatty liver disease.Gastroenterology139:e1-e6. 10.1053/j.gastro.2010.07.057

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  ChoiY.SimsG. E.MurphyS.MillerJ. R.ChanA. P. (2012). Predicting the functional effect of amino acid substitutions and indels.PLoS One7:e46688. 10.1371/journal.pone.0046688

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  ComuzzieA. G.ColeS. A.LastonS. L.VorugantiV. S.HaackK.GibbsR. A.et al (2012). Novel genetic loci identified for the pathophysiology of childhood obesity in the Hispanic population.PLoS One7:e51954. 10.1371/journal.pone.0051954

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  DanecekP.AutonA.AbecasisG.AlbersC. A.BanksE.DePristoM. A.et al (2011). The variant call format and VCFtools.Bioinformatics272156–2158. 10.1093/bioinformatics/btr330

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 18

  DavisM. (2001). Late Victorian Holocausts: El Niño Famines and the Making of the Third World.London: Verso.

  • Google Scholar
  • View reference in article

• 19

  DengL.Ruiz-linaresA.XuS.WangS. (2016). Ancestry variation and footprints of natural selection along the genome in Latin American populations.Sci. Rep.61–7. 10.1038/srep21766

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  DeschampsM.LavalG.FagnyM.ItanY.AbelL.CasanovaJ. L.et al (2016). Genomic signatures of selective pressures and introgression from archaic hominins at human innate immunity genes.Am. J. Hum. Genet.985–21. 10.1016/j.ajhg.2015.11.014

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 21

  DiversJ.PalmerN. D.LangefeldC. D.BrownW. M.LuL.HicksP. J.et al (2017). Genome-wide association study of coronary artery calcified atherosclerotic plaque in African Americans with type 2 diabetes.BMC Genet.18:105. 10.1186/s12863-017-0572-9

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  DunnJ. S.MlynarskiW. M.PezzolesiM. G.BorowiecM.PowersC.KrolewskiA. S.et al (2006). Examination of PPP1R3B as a candidate gene for the type 2 diabetes and MODY loci on chromosome 8p23.Ann. Hum. Genet.70587–593. 10.1111/j.1469-1809.2005.00248.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  ExcoffierL.SmouseP. E.QuattroJ. M. (1992). Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data.Genetics131479–491. 10.5962/bhl.title.86657

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 24

  FloresY. N.Velázquez-CruzR.RamírezP.BañuelosM.ZhangZ. F.YeeH. F.et al (2016). Association between PNPLA3 (rs738409), LYPLAL1 (rs12137855), PPP1R3B (rs4240624), GCKR (rs780094), and elevated transaminase levels in overweight/obese Mexican adults.Mol. Biol. Rep.431359–1369. 10.1007/s11033-016-4058-z

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  GoetzL. H.Uribe-BruceL.QuarlessD.LibigerO.SchorkN. J. (2014). Admixture and clinical phenotypic variation.Hum. Hered.7773–86. 10.1159/000362233

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 26

  González-PérezA.López-BigasN. (2011). Improving the assessment of the outcome of nonsynonymous SNVs with a consensus deleteriousness score.Condel. Am. J. Hum. Genet.88440–449. 10.1016/j.ajhg.2011.03.004

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 27

  GrossmanS. R.ShylakhterI.KarlssonE. K.ByrneE. H.MoralesS.FriedenG.et al (2010). A composite of multiple signals distinguishes causal variants in regions of positive selection.Science327883–886. 10.1126/science.1183863

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 28

  GuanY. (2014). Detecting structure of haplotypes and local ancestry.Genetics196625–642. 10.1534/genetics.113.160697

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 29

  HayesM. G.UrbanekM.HivertM. F.ArmstrongL. L.MorrisonJ.GuoC.et al (2013). Identification of HKDC1 and BACE2 as genes influencing glycemic traits during pregnancy through genome-wide association studies.Diabetes623282–3291. 10.2337/db12-1692

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 30

  HindorffL. A.BonhamV. L.BrodyL. C.GinozaM. E. C.HutterC. M.ManolioT. A.et al (2018). Prioritizing diversity in human genomics research.Nat. Rev. Genet.19175–185. 10.1038/nrg.2017.89

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 31

  KarczewskiK. J.FrancioliL. C.TiaoG.CummingsB. B.AlföldiJ.WangQ.et al (2020). The mutational constraint spectrum quantified from variation in 141,456 humans.Nature581434–443. 10.1038/s41586-020-2308-7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 32

  KehdyF. S. G.GouveiaM. H.MachadoM.MagalhãesW. C. S.HorimotoA. R.HortaB. 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.1128696–8701. 10.1073/pnas.1504447112

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 33

  Lima-CostaM. F.RodriguesL. C.BarretoM. L.GouveiaM.HortaB. L.MambriniJ.et al (2015). Genomic ancestry and ethnoracial self-classification based on 5,871 community-dwelling Brazilians (The Epigen Initiative).Sci. Rep.5:9812. 10.1038/srep09812

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 34

  MahajanA.TaliunD.ThurnerM.RobertsonN. R.TorresJ. M.RaynerN. W.et al (2018). Fine-mapping type 2 diabetes loci to single-variant resolution using high-density imputation and islet-specific epigenome maps.Nat. Genet.501505–1513. 10.1038/s41588-018-0241-6

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 35

  MaplesB. K.GravelS.KennyE. E.BustamanteC. D. (2013). RFMix: a discriminative modeling approach for rapid and robust local-ancestry inference.Am. J. Hum. Genet.93278–288. 10.1016/j.ajhg.2013.06.020

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 36

  MartinA. R.GignouxC. R.WaltersR. K.WojcikG. L.NealeB. M.GravelS.et al (2017). Human demographic history impacts genetic risk prediction across diverse populations.Am. J. Hum. Genet.100635–649. 10.1016/j.ajhg.2017.03.004

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 37

  McLarenW.GilL.HuntS. E.RiatH. S.RitchieG. R. S.ThormannA.et al (2016). The ensembl variant effect predictor.Genome Biol.17:122. 10.1186/s13059-016-0974-4

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 38

  MehtaM. B.ShewaleS. V.SequeiraR. N.MillarJ. S.HandN. J.RaderD. J. (2017). Hepatic protein phosphatase 1 regulatory subunit 3B (Ppp1r3b) promotes hepatic glycogen synthesis and thereby regulates fasting energy homeostasis.J. Biol. Chem.29210444–10454. 10.1074/jbc.M116.766329

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 39

  MiH.MuruganujanA.ThomasP. D. (2013). PANTHER in 2013: modeling the evolution of gene function, and other gene attributes, in the context of phylogenetic trees.Nucleic Acids Res.41D377–D386. 10.1093/nar/gks1118

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 40

  MoonesingheR.IoannidisJ. P. A.FlandersW. D.YangQ.TrumanB. I.KhouryM. J. (2012). Estimating the contribution of genetic variants to difference in incidence of disease between population groups.Eur. J. Hum. Genet.20831–836. 10.1038/ejhg.2012.15

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 41

  MylesS.DavisonD.BarrettJ.StonekingM.TimpsonN. (2008). Worldwide population differentiation at disease-associated SNPs.BMC Med. Genomics1:22. 10.1186/1755-8794-1-22

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 42

  NaslavskyM. S.YamamotoG. L.de AlmeidaT. F.EzquinaS. A. M.SunagaD. Y.PhoN.et al (2017). Exomic variants of an elderly cohort of Brazilians in the ABraOM database.Hum. Mutat.38751–763. 10.1002/humu.23220

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 43

  NiarchouM.ByrneE. M.TrzaskowskiM.SidorenkoJ.KemperK. E.McGrathJ. J.et al (2020). Genome-wide association study of dietary intake in the UK biobank study and its associations with schizophrenia and other traits.Transl. Psychiatry10:51. 10.1038/s41398-020-0688-y

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 44

  O’ConnellJ.GurdasaniD.DelaneauO.PirastuN.UliviS.CoccaM.et al (2014). A general approach for haplotype phasing across the full spectrum of relatedness.PLoS Genet.10:e1004234. 10.1371/journal.pgen.1004234

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 45

  PatinE.LopezM.GrollemundR.VerduP.HarmantC.QuachH.et al (2017). Dispersals and genetic adaptation of Bantu-speaking populations in Africa and North America.Science356543–546. 10.1126/science.aal1988

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 46

  PeltolaP.PihlajamäkiJ.KoutnikovaH.RuotsalainenE.SalmenniemiU.VauhkonenI.et al (2006). Visceral obesity is associated with high levels of serum squalene.Obesity (Silver Spring)141155–1163. 10.1038/oby.2006.132

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 47

  PulitS. L.StonemanC.MorrisA. P.WoodA. R.GlastonburyC. A.TyrrellJ.et al (2019). Meta-analysis of genome-wide association studies for body fat distribution in 694 649 individuals of European ancestry.Hum. Mol. Genet.28166–174. 10.1093/hmg/ddy327

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 48

  PurcellS.NealeB.Todd-BrownK.ThomasL.FerreiraM. A. R.BenderD.et al (2007). PLINK: a tool set for whole-genome association and population-based linkage analyses.Am. J. Hum. Genet.81559–575. 10.1086/519795

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 49

  Rask-AndersenM.KarlssonT.EkW. E.JohanssonÅ (2019). Genome-wide association study of body fat distribution identifies adiposity loci and sex-specific genetic effects.Nat. Commun.10:339. 10.1038/s41467-018-08000-4

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 50

  RentzschP.WittenD.CooperG. M.ShendureJ.KircherM. (2019). CADD: predicting the deleteriousness of variants throughout the human genome.Nucleic Acids Res.47D886–D894. 10.1093/nar/gky1016

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 51

  RochaC. S.SecolinR.RodriguesM. R.CarvalhoB. S.Lopes-CendesI. (2020). The Brazilian Initiative on Precision Medicine (BIPMed): fostering genomic data-sharing of underrepresented populations.NPJ Genomic Med.5:42. 10.1038/s41525-020-00149-6

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 52

  Rodrigues de MouraR.CoelhoA. V. C.de Queiroz BalbinoV.CrovellaS.BrandãoL. A. C. (2015). Meta-analysis of Brazilian genetic admixture and comparison with other Latin America countries.Am. J. Hum. Biol.27674–680. 10.1002/ajhb.22714

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 53

  Ruiz-LinaresA.AdhikariK.Acuña-AlonzoV.Quinto-SanchezM.JaramilloC.AriasW.et al (2014). Admixture in Latin America: geographic structure, phenotypic diversity and self-perception of ancestry based on 7,342 individuals.PLoS Genet.10:e1004572. 10.1371/journal.pgen.1004572

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 54

  SabetiP. C.VarillyP.FryB.LohmuellerJ.HostetterE.CotsapasC.et al (2007). Genome-wide detection and characterization of positive selection in human populations.Nature449913–918. 10.1038/nature06250

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 55

  SalmM. P. A.HorswellS. D.HutchisonC. E.SpeedyH. E.YangX.LiangL.et al (2012). The origin, global distribution, and functional impact of the human 8p23 inversion polymorphism.Genome Res.221144–1153. 10.1101/gr.126037.111

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 56

  SartoriD. J.WilburC. J.LongS. Y.RankinM. M.LiC.BradfieldJ. P.et al (2014). GATA factors promote ER integrity and β-cell survival and contribute to type 1 diabetes risk.Mol. Endocrinol.2828–39. 10.1210/me.2013-1265

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 57

  ScheragA.DinaC.HinneyA.VatinV.ScheragS.VogelC. I.et al (2010). Two new Loci for body-weight regulation identified in a joint analysis of genome-wide association studies for early-onset extreme obesity in French and german study groups.PLoS Genet.6:e1000916. 10.1371/journal.pgen.1000916

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 58

  SchlauchK. A.ReadR. W.LombardiV. C.ElhananG.MetcalfW. J.SlonimA. D.et al (2020). A comprehensive genome-wide and phenome-wide examination of BMI and obesity in a Northern Nevadan Cohort.G310645–664. 10.1534/g3.119.400910

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 59

  SchwarzJ. M.RödelspergerC.SchuelkeM.SeelowD. (2010). MutationTaster evaluates disease-causing potential of sequence alterations.Nat. Methods7575–576. 10.1038/nmeth0810-575

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 60

  SecolinR.Mas-SandovalA.AraunaL. R.TorresF. R.de AraujoT. K.SantosM. L.et al (2019). Distribution of local ancestry and evidence of adaptation in admixed populations.Sci. Rep.9:13900. 10.1038/s41598-019-50362-2

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 61

  SeldinM. F.PasaniucB.PriceA. L. (2011). New approaches to disease mapping in admixed populations.Nat. Rev. Genet.12523–528. 10.1038/nrg3002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 62

  Shaw-SmithC.De FrancoE.Lango AllenH.BatlleM.FlanaganS. E.BorowiecM.et al (2014). GATA4 mutations are a cause of neonatal and childhood-onset diabetes.Diabetes632888–2894. 10.2337/db14-0061

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 63

  ShihabH. A.RogersM. F.GoughJ.MortM.CooperD. N.DayI. N. M.et al (2015). An integrative approach to predicting the functional effects of non-coding and coding sequence variation.Bioinformatics311536–1543. 10.1093/bioinformatics/btv009

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 64

  SimN. L.KumarP.HuJ.HenikoffS.SchneiderG.NgP. C. (2012). SIFT web server: predicting effects of amino acid substitutions on proteins.Nucleic Acids Res.40W452–W457. 10.1093/nar/gks539

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 65

  TavtigianS. V.DeffenbaughA. M.YinL.JudkinsT.SchollT.SamollowP. B.et al (2006). Comprehensive statistical study of 452 BRCA1 missense substitutions with classification of eight recurrent substitutions as neutral.J. Med. Genet.43295–305. 10.1136/jmg.2005.033878

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 66

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

  • CrossRef
  • Google Scholar

• 67

  VoightB. F.KudaravalliS.WenX.PritchardJ. K. (2006). A map of recent positive selection in the human genome.PLoS Biol.4:e72. 10.1371/journal.pbio.0040072

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 68

  VujkovicM.KeatonJ. M.LynchJ. A.MillerD. R.ZhouJ.TcheandjieuC.et al (2020). Discovery of 318 new risk loci for type 2 diabetes and related vascular outcomes among 1.4 million participants in a multi-ancestry meta-analysis.Nat. Genet.52680–691. 10.1038/s41588-020-0637-y

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 69

  WangH.ZhangF.ZengJ.WuY.KemperK. E.XueA.et al (2019). Genotype-by-environment interactions inferred from genetic effects on phenotypic variability in the UK Biobank.Sci. Adv.5:eaaw3538. 10.1126/sciadv.aaw3538

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 70

  WangK.LiM.HakonarsonH. (2010). ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data.Nucleic Acids Res.38:e164. 10.1093/nar/gkq603

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 71

  WillerC. J.SchmidtE. M.SenguptaS.PelosoG. M.GustafssonS.KanoniS.et al (2013). Discovery and refinement of loci associated with lipid levels.Nat. Genet.451274–1283. 10.1038/ng.2797

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 72

  XueA.WuY.ZhuZ.ZhangF.KemperK. E.ZhengZ.et al (2018). Genome-wide association analyses identify 143 risk variants and putative regulatory mechanisms for type 2 diabetes.Nat. Commun.9:2941. 10.1038/s41467-018-04951-w

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 73

  YiX.LiangY.Huerta-SanchezE.JinX.CuoZ. X. P.PoolJ. E.et al (2010). Sequencing of 50 human exomes reveals adaptation to high altitude.Science32975–78. 10.1126/science.1190371

  • Pubmed Abstract
  • CrossRef
  • Google Scholar