Search

Background/objectives: Admixed populations are a resource to study the global genetic architecture of complex phenotypes, which is critical, considering that non-European populations are severely underrepresented in genomic studies. Here, we study the genetic architecture of BMI in children, young adults, and elderly individuals from the admixed population of Brazil. Subjects/methods: Leveraging admixture in Brazilians, whose chromosomes are mosaics of fragments of Native American, European, and African origins, we used genome-wide data to perform admixture mapping/fine-mapping of body mass index (BMI) in three Brazilian population-based cohorts from Northeast (Salvador), Southeast (Bambui), and South (Pelotas). Results: We found significant associations with African-associated alleles in children from Salvador (PALD1 and ZMIZ1 genes), and in young adults from Pelotas (NOD2 and MTUS2 genes). More importantly, in Pelotas, rs114066381, mapped in a potential regulatory region, is significantly associated only in females (p = 2.76e-06). This variant is rare in Europeans but with frequencies of similar to 3% in West Africa and has a strong female-specific effect (95% CI: 2.32-5.65 kg/m(2) per each A allele). We confirmed this sex-specific association and replicated its strong effect for an adjusted fat mass index in the same Pelotas cohort, and for BMI in another Brazilian cohort from Sao Paulo (Southeast Brazil). A meta-analysis confirmed the significant association. Remarkably, we observed that while the frequency of rs114066381-A allele ranges from 0.8 to 2.1% in the studied populations, it attains similar to 9% among women with morbid obesity from Pelotas, Sao Paulo, and Bambui. The effect size of rs114066381 is at least five times higher than the FTO SNPs rs9939609 and rs1558902, already emblematic for their high effects. Conclusions: We identified six candidate SNPs associated with BMI. rs114066381 stands out for its high effect that was replicated and its high frequency in women with morbid obesity. We demonstrate how admixed populations are a source of new relevant phenotype-associated genetic variants.

Genetic and omics analyses frequently require independent observations, which is not guaranteed in real datasets. When relatedness cannot be accounted for, solutions involve removing related individuals (or observations) and, consequently, a reduction of available data. We developed a network-based relatedness-pruning method that minimizes dataset reduction while removing unwanted relationships in a dataset. It uses node degree centrality metric to identify highly connected nodes (or individuals) and implements heuristics that approximate the minimal reduction of a dataset to allow its application to large datasets. NAToRA outperformed two popular methodologies (implemented in software PLINK and KING) by showing the best combination of effective relatedness-pruning, removing all relatives while keeping the largest possible number of individuals in all datasets tested and also, with similar or lesser reduction in genetic diversity. NAToRA is freely available, both as a standalone tool that can be easily incorporated as part of a pipeline, and as a graphical web tool that allows visualization of the relatedness networks. NAToRA also accepts a variety of relationship metrics as input, which facilitates its use. We also present a genealogies simulator software used for different tests performed in the manuscript.

The Covid-19 pandemic has created an unprecedented need for epidemiological monitoring using diverse strategies. We conducted a project combining prevalence, seroprevalence, and genomic surveillance approaches to describe the initial pandemic stages in Betim City, Brazil. We collected 3239 subjects in a population-based age-, sex- and neighbourhood-stratified, household, prospective; cross-sectional study divided into three surveys 21 days apart sampling the same geographical area. In the first survey, overall prevalence (participants positive in serological or molecular tests) reached 0.46% (90% CI 0.12% – 0.80%), followed by 2.69% (90% CI 1.88% – 3.49%) in the second survey and 6.67% (90% CI 5.42% - 7.92%) in the third. The underreporting reached 11, 19.6, and 20.4 times in each survey, respectively. We observed increased odds to test positive in females compared to males (OR 1.88 95% CI 1.25 – 2.82), while the single best predictor for positivity was ageusia/ anosmia (OR 8.12, 95% CI 4.72 – 13.98). Thirty-five SARS-CoV-2 genomes were sequenced, of which 18 were classified as lineage B.1.1.28, while 17 were B.1.1.33. Multiple independent viral introductions were observed. Integration of multiple epidemiological strategies was able to describe Covid-19 dispersion in the city adequately. Presented results have helped local government authorities to guide pandemic management.

The current search for host-susceptibility variants for COVID-19 contrasts with the fact that the study of the genetic architecture of Severe Acute Respiratory Syndrome (SARS) has been neglected. For human/SARS-CoV-2 interactome genes ACE2, TMPRSS2 and BSG, there is only one convincing evidence of association in Asians with influenza-induced SARS for TMPRSS2-rs2070788, tag-SNP of the eQTL rs383510. This case illustrates the importance of population genetics and of sequencing data in the design of genetic association studies in different human populations: the high linkage disequilibrium (LD) between rs2070788 and rs383510 is Asian-specific. Leveraging on a combination of genotyping and sequencing data for Native Americans (neglected in genetic studies), we show that while their frequencies of the Asian tag-SNP rs2070788 is, surprisingly, the highest worldwide, it is not in LD with the eQTL rs383510, that therefore, should be directly genotyped in genetic association studies of SARS in populations with Native American ancestry.

Brazil was one of the most affected countries during the SARS-CoV-2 pandemic, with more than 10 million cases and more than 300.000 deaths reported. Together with the increasing rate of cases, the emergence of new variants of concern (VOCs) possessing mutations that affect immunogenicity and transmissibility poses a concerning burden in the healthcare system. In Brazil, three such lineages are currently circulating: P.1, a lineage firstly described in Manaus, P.2, firstly described in Rio de Janeiro, both already reported in other cities in the country, and B.1.1.7, a lineage initially identified in England and recently associated to increasing mortality rates.In this scenario, it urges the necessity of a rapid assessment methodology for identification of variants that could contribute to a better epidemiological surveillance and public healthcare administration decisions in dealing with the pandemic. Here, we describe a PCR genotyping approach using RhAmp technology (IDT Technologies) to correctly assess the circulating emergent variants P.1, P.2 and B.1.1.7 in Brazil. Genotyping is a quicker methodology requiring only a real time PCR machine available in most of laboratories preforming COVID-19 molecular diagnosis.

Brazil currently has one of the fastest-growing severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) epidemics in the world. Because of limited available data, assessments of the impact of nonpharmaceutical interventions (NPIs) on this virus spread remain challenging. Using a mobility-driven transmission model, we show that NPIs reduced the reproduction number from >3 to 1 to 1.6 in São Paulo and Rio de Janeiro. Sequencing of 427 new genomes and analysis of a geographically representative genomic dataset identified >100 international virus introductions in Brazil. We estimate that most (76%) of the Brazilian strains fell in three clades that were introduced from Europe between 22 February and 11 March 2020. During the early epidemic phase, we found that SARS-CoV-2 spread mostly locally and within state borders. After this period, despite sharp decreases in air travel, we estimated multiple exportations from large urban centers that coincided with a 25% increase in average traveled distances in national flights. This study sheds new light on the epidemic transmission and evolutionary trajectories of SARS-CoV-2 lineages in Brazil and provides evidence that current interventions remain insufficient to keep virus transmission under control in this country.

Brazil currently has one of the fastest growing SARS-CoV-2 epidemics in the world. Due to limited available data, assessments of the impact of non-pharmaceutical interventions (NPIs) on virus transmission and epidemic spread remain challenging. We investigate the impact of NPIs in Brazil using epidemiological, mobility and genomic data. Mobility-driven transmission models for São Paulo and Rio de Janeiro cities show that the reproduction number (Rt) reached below 1 following NPIs but slowly increased to values between 1 to 1.3 (1.0–1.6). Genome sequencing of 427 new genomes and analysis of a geographically representative genomic dataset from 21 of the 27 Brazilian states identified >100 international introductions of SARS-CoV-2 in Brazil. We estimate that three clades introduced from Europe emerged between 22 and 27 February 2020, and were already well-established before the implementation of NPIs and travel bans. During this first phase of the epidemic establishment of SARS-CoV-2 in Brazil, we find that the virus spread mostly locally and within-state borders. Despite sharp decreases in national air travel during this period, we detected a 25% increase in the average distance travelled by air passengers during this time period. This coincided with the spread of SARS-CoV-2 from large urban centers to the rest of the country. In conclusion, our results shed light on the role of large and highly connected populated centres in the rapid ignition and establishment of SARS-CoV-2, and provide evidence that current interventions remain insufficient to keep virus transmission under control in Brazil.One Sentence SummaryJoint analysis of genomic, mobility and epidemiological novel data provide unique insight into the spread and transmission of the rapidly evolving epidemic of SARS-CoV-2 in Brazil.

Western South America was one of the worldwide cradles of civilization. The well known Inca Empire was the tip of the iceberg of a cultural and biological evolutionary process that started 14-11 thousand years ago. Genetic data from 18 Peruvian populations reveal that: (1) The between-population homogenization of the central-southern Andes and its differentiation with respect to Amazonian populations of similar latitudes do not extend northward. Instead, longitudinal gene flow between the northern coast of Peru, Andes and Amazonia accompanied cultural and socioeconomic interactions revealed by archeological studies. This pattern recapitulates the environmental and cultural differentiation between the fertile north, where altitudes are lower; and the arid south, where the Andes are higher, acting as a genetic barrier between the sharply different environments of the Andes and Amazonia (2). The genetic homogenization between the populations of the arid Andes is not only due to migration during the Inca Empire or the subsequent colonial period. It started at least during the earlier expansion of the pre-Inca Wari Empire (600-1000 YBP) (3) This demographic history allowed for cases of positive natural selection in the high and arid Andes vs. the low Amazon tropical forest: in the Andes, HAND2-AS1 (heart and neural crest derivatives expressed 2 antisense RNA1, related with cardiovascular function) and DUOX2 (dual oxidase 2, related to thyroid function and innate immunity) genes; in the Amazon, the gene encoding for the CD45 protein, essential for antigen recognition by T/B lymphocytes in viral-host interaction, consistent with the host-virus arms race hypothesis.

Admixed populations are a resource to study the global genetic architecture of complex phenotypes, which is critical, considering that non-European populations are severely under-represented in genomic studies. Leveraging admixture in Brazilians, whose chromosomes are mosaics of fragments of Native American, European and African origins, we used genome-wide data to perform admixture mapping/fine-mapping of Body Mass Index (BMI) in three population-based cohorts from Northeast (Salvador), Southeast (Bambuí) and South (Pelotas) of the country. We found significant associations with African-associated alleles in children from Salvador (PALD1 and ZMIZ1 genes), and in young adults from Pelotas (NOD2 and MTUS2 genes). More importantly, in Pelotas, rs114066381, mapped in a potential regulatory region, is significantly associated only in females (p= 2.76 e-06). This variant is very rare in Europeans but with frequencies of ~3% in West Africa, and has a strong female-specific effect (95%CI: 2.32-5.65 kg/m2 per each A allele). We confirmed this sex-specific association and replicated its strong effect for an adjusted fat-mass index in the same Pelotas cohort, and for BMI in another Brazilian cohort from São Paulo (Southeast Brazil). A meta-analysis confirmed the significant association. Remarkably, we observed that while the frequency of rs114066381-A allele ranges from 0.8 to 2.1% in the studied populations, it attains ~9% among morbidly obese women from Pelotas, São Paulo, and Bambuí. The effect size of rs114066381 is at least five-times the effect size of the FTO SNPs rs9939609 and rs1558902, already emblematic for their high effects, and for which we replicated associations in Pelotas. We demonstrate how, after a decade of GWAS mostly performed in European-ancestry populations, non-European and admixed populations are a source of new relevant phenotype-associated genetic variants.

The Transatlantic Slave Trade transported more than 9 million Africans to the Americas between the early 16th and the mid-19th centuries. We performed genome-wide analysis of 6,267 individuals from 22 populations and observed an enrichment in West-African ancestry in northern latitudes of the Americas, whereas South/East African ancestry is more prevalent in southern South-America. This pattern results from distinct geographic and geopolitical factors leading to population differentiation. However, we observed a decrease of 68% in the African gene pool between-population diversity within the Americas when compared to the regions of origin from Africa, underscoring the importance of historical factors favoring admixture between individuals with different African origins in the New World. This is consistent with the excess of West-Central Africa ancestry (the most prevalent in the Americas) in the US and Southeast-Brazil, respect to historical-demography expectations. Also, in most of the Americas, admixture intensification occurred between 1,750 and 1,850, which correlates strongly with the peak of arrivals from Africa. This study contributes with a population genetics perspective to the ongoing social, cultural and political debate regarding ancestry, race, and admixture in the Americas.Significance StatementDifferently from most genetic studies, that have estimated the overall African ancestry in the Americas, we perform a finer geographic analysis and infer how different African groups contributed to North-, South-American and Caribbean populations, in the context of geographic and geopolitical factors. We also perform a formal comparison of information from demographic history records of the Transatlantic Slave Trade with inferences based on genomic diversity of current populations. Our approach reveals the distinct regional African ancestry roots of different populations from North-, South-America and the Caribe and other important aspects of the historical process of mestizaje and its dynamics in the American continent.

We present allele frequencies involving 39 pharmacogenetic biomarkers studied in Brazil, and their distribution on self-reported race/color categories that: (1) involve a mix of perceptions about ancestry, morphological traits, and cultural/identity issues, being social constructs pervasively used in Brazilian society and medical studies; (2) are associated with disparities in access to health services, as well as in their representation in genetic studies, and (3), as we report here, explain a larger portion of the variance of pharmaco-allele frequencies than geography. We integrated a systematic review of studies on healthy volunteers (years 1968–2017) and the analysis of allele frequencies on three population-based cohorts from northeast, southeast, and south, the most populated regions of Brazil. Cross-validation of results from these both approaches suggest that, despite methodological heterogeneity of the 120 studies conducted on 51,747 healthy volunteers, allele frequencies estimates from systematic review are reliable. We report differences in allele frequencies between color categories that persist despite the homogenizing effect of >500 years of admixture. Among clinically relevant variants: CYP2C9*2 (null), CYP3A5*3 (defective), SLCO1B1-rs4149056(C), and VKORC1-rs9923231(A) are more frequent in Whites than in Blacks. Brazilian Native Americans show lower frequencies of CYP2C9*2, CYP2C19*17 (increased activity), and higher of SLCO1B1-rs4149056(C) than other Brazilian populations. We present the most current and informative database of pharmaco-allele frequencies in Brazilian healthy volunteers.

A iniciativa EPIGEN-Brasil é um projeto estratégico do Ministério da Saúde que constitui a maior iniciativa latino-americana para o estudo da diversidade genômica humana, integrando duas importantes tradições científicas brasileiras: Saúde Coletiva e Genética de Populações. No contexto dos resultados dessa iniciativa, neste artigo discutimos a importância dos estudos de genética de populações para o entendimento do processo de miscigenação no Brasil e para a Saúde Pública. Mostramos como a alta resolução genômica do projeto EPIGEN-Brasil e metodologias computacionais avançadas possibilitaram novas inferências sobre as origens e a dinâmica da miscigenação e como a história social brasileira influenciou a estrutura genética da população.

While South Americans are underrepresented in human genomic diversity studies, Brazil has been a classical model for population genetics studies on admixture. We present the results of the EPIGEN Brazil Initiative, the most comprehensive up-to-date genomic analysis of any Latin-American population. A population-based genome-wide analysis of 6,487 individuals was performed in the context of worldwide genomic diversity to elucidate how ancestry, kinship, and inbreeding interact in three populations with different histories from the Northeast (African ancestry: 50%), Southeast, and South (both with European ancestry >70%) of Brazil. We showed that ancestry-positive assortative mating permeated Brazilian history. We traced European ancestry in the Southeast/South to a wider European/Middle Eastern region with respect to the Northeast, where ancestry seems restricted to Iberia. By developing an approximate Bayesian computation framework, we infer more recent European immigration to the Southeast/South than to the Northeast. Also, the observed low Native-American ancestry (6–8%) was mostly introduced in different regions of Brazil soon after the European Conquest. We broadened our understanding of the African diaspora, the major destination of which was Brazil, by revealing that Brazilians display two within-Africa ancestry components: one associated with non-Bantu/western Africans (more evident in the Northeast and African Americans) and one associated with Bantu/eastern Africans (more present in the Southeast/South). Furthermore, the whole-genome analysis of 30 individuals (42-fold deep coverage) shows that continental admixture rather than local post-Columbian history is the main and complex determinant of the individual amount of deleterious genotypes.

Background. Gastric adenocarcinoma is associated with chronic infection by Helicobacter pylori and with the host inflammatory response triggered by it, with substantial inter-person variation in the immune response profile due to host genetic factors. Aim. To investigate the diversity of the proinflammatory genes IL8, its receptors and PTGS2 in Amerindians; to test whether candidate SNPs in these genes are associated with gastric cancer in an admixed population with high Amerindian ancestry from Lima, Peru; and to assess whether an IL8RB promoter-derived haplotype affects gene expression. Methods. We performed a Sanger-resequencing population survey, a candidate-gene association study (220 cases, 288 controls) and meta-analyses. We also performed an in vitro validation by a reporter gene assay of IL8RB promoter. Results. The diversity of the promoter of studied genes in Native Americans is similar to Europeans. Although an association between candidate SNPs and gastric cancer was not found in Peruvians, trend in our data is consistent with meta-analyses results that suggest PTGS2-rs689466-A is associated with H. pylori-associated gastric cancer in East Asia. IL8RB promoter-derived haplotype (rs3890158-A/rs4674258-T), common in Peruvians, was up-regulated by TNF-α unlike the ancestral haplotype (rs3890158-G/rs4674258-C). Bioinformatics analysis suggests that this effect stemmed from creation of a binding site for the FOXO3 transcription factor by rs3890158G>A. Conclusions. Our updated meta-analysis reinforces the role of PTGS2-rs689466-A in gastric cancer in Asians, although more studies that control for ancestry are necessary to clarify its role in Latin Americans. Finally, we suggest that IL8RB-rs3890158G>A is a cis-regulatory SNP.