NAToRA, a relatedness-pruning method to minimize the loss of dataset size in genetic and omics analyses.

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 complex datasets. When compared with two other popular population genetics methodologies (PLINK and KING), NAToRA shows the best combination of removing all relatives while keeping the largest possible number of individuals in all datasets tested and also, with similar effects on the allele frequency spectrum and Principal Component Analysis than PLINK and KING. 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 release a genealogies simulator software used for different tests performed in this study.

Population genetics of PDE4B (phosphodiesterase-4B) in neglected Native Americans: Implications for cancer pharmacogenetics.

PDE4B (phosphodiesterase-4B) has an important role in cancer and in pharmacology of some disorders, such as inflammatory diseases. Remarkably in Native Americans, PDE4B variants are associated with acute lymphoblastic leukemia (ALL) relapse, as this gene modulates sensitivity of glucocorticoids used in ALL chemotherapy. PDE4B allele rs6683977.G, associated with genomic regions of Native American origin in US-Hispanics (admixed among Native Americans, Europeans, and Africans), increases ALL relapse risk, contributing to an association between Native American ancestry and ALL relapse that disappeared with an extra-phase of chemotherapy. This result insinuates that indigenous populations along the Americas may have high frequencies of rs6683977.G, but this has never been corroborated. We studied ancestry and PDE4B diversity in 951 healthy individuals from nine Latin American populations. In non-admixed Native American populations rs6683977.G has frequencies greater than 90%, is in linkage disequilibrium with other ALL relapse associated and regulatory variants in PDE4B-intron-7, conforming haplotypes showing their highest worldwide frequencies in Native Americans (>0.82). Our findings inform the discussion on the pertinence of an extra-phase of chemotherapy in Native American populations, and exemplifies how knowledge generated in US-Hispanics is relevant for their even more neglected and vulnerable Native American ancestors along the American continent.

Seroprevalence, Prevalence, and Genomic Surveillance: Monitoring the Initial Phases of the SARS-CoV-2 Pandemic in Betim, Brazil.

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 neighborhood-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. 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 adequately describe COVID-19 dispersion in the city. Presented results have helped local government authorities to guide pandemic management.

Human-SARS-CoV-2 interactome and human genetic diversity: TMPRSS2-rs2070788, associated with severe influenza, and its population genetics caveats in Native Americans.

For human/SARS-CoV-2 interactome genes ACE2, TMPRSS2 and BSG, there is a 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.

Admixture/fine-mapping in Brazilians reveals a West African associated potential regulatory variant (rs114066381) with a strong female-specific effect on body mass and fat mass indexes.

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 (Bambuí), 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 ~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 women with morbid obesity from Pelotas, São Paulo, and Bambuí. 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.

The genetic structure and adaptation of Andean highlanders and Amazonians are influenced by the interplay between geography and culture.

Western South America was one of the worldwide cradles of civilization. The well-known Inca Empire was the tip of the iceberg of an evolutionary process that started 11,000 to 14,000 years ago. Genetic data from 18 Peruvian populations reveal the following: 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 archeology. 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 migrations during the Inca Empire or the subsequent colonial period. It started at least during the earlier expansion of the Wari Empire (600 to 1,000 years before present). 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, a putative enhancer in HAND2-AS1 (heart and neural crest derivatives expressed 2 antisense RNA1, a noncoding gene related to cardiovascular function) and rs269868-C/Ser1067 in DUOX2 (dual oxidase 2, related to thyroid function and innate immunity) genes and, in the Amazon, the gene encoding for the CD45 protein, essential for antigen recognition by T and B lymphocytes in viral-host interaction.

Origins, Admixture Dynamics, and Homogenization of the African Gene Pool in the Americas.

The Transatlantic Slave Trade transported more than 9 million Africans to the Americas between the early 16th and the mid-19th centuries. We performed a genome-wide analysis using 6,267 individuals from 25 populations to infer how different African groups contributed to North-, South-American, and Caribbean populations, in the context of geographic and geopolitical factors, and compared genetic data with demographic history records of the Transatlantic Slave Trade. We observed that West-Central Africa and Western Africa-associated ancestry clusters are more prevalent in northern latitudes of the Americas, whereas the South/East Africa-associated ancestry cluster is more prevalent in southern latitudes of the Americas. This pattern results from geographic and geopolitical factors leading to population differentiation. However, there is a substantial decrease in the between-population differentiation of the African gene pool within the Americas, when compared with 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 between-population homogenization in the Americas is consistent with the excess of West-Central Africa ancestry (the most prevalent in the Americas) in the United States and Southeast-Brazil, with respect to historical-demography expectations. We also inferred that in most of the Americas, intercontinental admixture intensification occurred between 1750 and 1850, 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, admixture, and the mestizaje process in the Americas.

Genetic structure of pharmacogenetic biomarkers in Brazil inferred from a systematic review and population-based cohorts: a RIBEF/EPIGEN-Brazil initiative.

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.

Population, Epidemiological, and Functional Genetics of Gastric Cancer Candidate Genes in Peruvians with Predominant Amerindian Ancestry.

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.