Genetic determinants of serum bilirubin using inferred native American gene variants in Chilean adolescents.

Gene variants in the UGT1A1 gene are strongly associated with circulating bilirubin levels in several populations, as well as other variants of modest effect across the genome. However, the effects of such variants are unknown regarding the Native American ancestry of the admixed Latino population. Our objective was to assess the Native American genetic determinants of serum bilirubin in Chilean admixed adolescents using the local ancestry deconvolution approach. We measured total serum bilirubin levels in 707 adolescents of the Chilean Growth and Obesity Cohort Study (GOCS) and performed high-density genotyping using the Illumina-MEGA array (>1.7 million genotypes). We constructed a local ancestry reference panel with participants from the 1000 Genomes Project, the Human Genome Diversity Project, and our GOCS cohort. Then, we inferred and isolated haplotype tracts of Native American, European, or African origin to perform genome-wide association studies. In the whole cohort, the rs887829 variant and others near UGT1A1 were the unique signals achieving genome-wide statistical significance (b = 0.30; p = 3.34 × 10-57). After applying deconvolution methods, we found that significance is also maintained in Native American (b = 0.35; p = 3.29 × 10-17) and European (b = 0.28; p = 1.14 × 10-23) ancestry components. The rs887829 variant explained a higher percentage of the variance of bilirubin in the Native American (37.6%) compared to European ancestry (28.4%). In Native American ancestry, carriers of the TT genotype of this variant averaged 4-fold higher bilirubinemia compared to the CC genotype (p = 2.82 × 10-12). We showed for the first time that UGT1A1 variants are the primary determinant of bilirubin levels in Native American ancestry, confirming its pan-ethnic relevance. Our study illustrates the general value of the local ancestry deconvolution approach to assessing isolated ancestry effects in admixed populations.

The first genome-wide association study in the Argentinian and Chilean populations identifies shared genetics with Europeans in Alzheimer's disease.

INTRODUCTION: Genome-wide association studies (GWAS) are fundamental for identifying loci associated with diseases. However, they require replication in other ethnicities. METHODS: We performed GWAS on sporadic Alzheimer's disease (AD) including 539 patients and 854 controls from Argentina and Chile. We combined our results with those from the European Alzheimer and Dementia Biobank (EADB) in a meta-analysis and tested their genetic risk score (GRS) performance in this admixed population. RESULTS: We detected apolipoprotein E ε4 as the single genome-wide significant signal (odds ratio  = 2.93 [2.37-3.63], P = 2.6 × 10-23 ). The meta-analysis with EADB summary statistics revealed four new loci reaching GWAS significance. Functional annotations of these loci implicated endosome/lysosomal function. Finally, the AD-GRS presented a similar performance in these populations, despite the score diminished when the Native American ancestry rose. DISCUSSION: We report the first GWAS on AD in a population from South America. It shows shared genetics modulating AD risk between the European and these admixed populations. HIGHLIGHTS: This is the first genome-wide association study on Alzheimer's disease (AD) in a population sample from Argentina and Chile. Trans-ethnic meta-analysis reveals four new loci involving lysosomal function in AD. This is the first independent replication for TREM2L, IGH-gene-cluster, and ADAM17 loci. A genetic risk score (GRS) developed in Europeans performed well in this population. The higher the Native American ancestry the lower the GRS values.

High Mitochondrial Haplotype Diversity Found in Three Pre-Hispanic Groups from Colombia.

The analysis of mitochondrial DNA (mtDNA) hypervariable region (HVR) sequence data from ancient human remains provides valuable insights into the genetic structure and population dynamics of ancient populations. mtDNA is particularly useful in studying ancient populations, because it is maternally inherited and has a higher mutation rate compared to nuclear DNA. To determine the genetic structure of three Colombian pre-Hispanic populations and compare them with current populations, we determined the haplotypes from human bone remains by sequencing several mitochondrial DNA segments. A wide variety of mitochondrial polymorphisms were obtained from 33 samples. Our results support a high population heterogeneity among pre-Hispanic populations in Colombia.

Bailey-Bloch Congenital Myopathy in Brazilian Patients: A Very Rare Myopathy with Malignant Hyperthermia Susceptibility.

BACKGROUND: Congenital myopathy-13 (CMYP13), also known as Bailey-Bloch congenital myopathy and Native American myopathy (NAM), is a condition caused by biallelic missense pathogenic variants in STAC3, which encodes an important protein necessary for the excitation-relaxation coupling machinery in the muscle. Patients with biallelic pathogenic variants in STAC3 often present with congenital weakness and arthrogryposis, cleft palate, ptosis, myopathic facies, short stature, kyphoscoliosis, and susceptibility to malignant hyperthermia provoked by anesthesia. We present two unrelated cases of Bailey-Bloch congenital myopathy descendants of non-consanguineous parents, which were investigated for delayed psychomotor development and generalized weakness. To the best of our knowledge, these are the first descriptions of CMYP13 in Brazil. In both patients, we found the previously described pathogenic missense variant p.Trp284Ser in homozygosity. CONCLUSION: We seek to highlight the need for screening for CMYP13 in patients expressing the typical phenotype of the disease even in the absence of Lumbee Native American ancestry, and to raise awareness to possible complications like malignant hyperthermia in Bailey-Bloch congenital myopathy.

Distribution of a novel CYP2C haplotype in Native American populations.

The CYP2C19 gene, located in the CYP2C cluster, encodes the major drug metabolism enzyme CYP2C19. This gene is highly polymorphic and no-function (CYP2C19*2 and CYP2C19*3), reduced function (CYP2C19*9) and increased function (CYP2C19*17) star alleles (haplotypes) are commonly used to predict CYP2C19 metabolic phenotypes. CYP2C19*17 and the genotype-predicted rapid (RM) and ultrarapid (UM) CYP2C19 metabolic phenotypes are absent or rare in several Native American populations. However, discordance between genotype-predicted and pharmacokinetically determined CYP2C19 phenotypes in Native American cohorts have been reported. Recently, a haplotype defined by rs2860840T and rs11188059G alleles in the CYP2C cluster has been shown to encode increased rate of metabolism of the CYP2C19 substrate escitalopram, to a similar extent as CYP2C19*17. We investigated the distribution of the CYP2C:TG haplotype and explored its potential impact on CYP2C19 metabolic activity in Native American populations. The study cohorts included individuals from the One Thousand Genomes Project AMR superpopulation (1 KG_AMR), the Human Genome Diversity Project (HGDP), and from indigenous populations living in Brazil (Kaingang and Guarani). The frequency range of the CYP2C:TG haplotype in the study cohorts, 0.469 to 0.598, is considerably higher than in all 1 KG superpopulations (range: 0.014-to 0.340). We suggest that the high frequency of the CYP2C:TG haplotype might contribute to the reported discordance between CYP2C19-predicted and pharmacokinetically verified CYP2C19 metabolic phenotypes in Native American cohorts. However, functional studies involving genotypic correlations with pharmacokinetic parameters are warranted to ascertain the importance of the CYP2C:TG haplotype.

Sequencing errors in Native American mitogenomes: Impact on clade definitions, haplogroup assignation, and beyond.

Available evidence allows the interpretation that some cases of absence of otherwise expected variation, based on phylogenetic expectations in mitogenomes of Native American origin, are due to artificial recombination rather than to homoplasy, while other more complex scenarios involving combination of original Cambridge Reference Sequence mistakes plus incomplete or incorrect scoring of variation are also showed. Several instances of mismatched control and coding regions as well as partially duplicated HV2 are observed in Peruvians, while intra-haplogroup chimaeras of different D1 subhaplogroups are referred to in Mexican Native Americans. A revised definition for haplogroup B2h is proposed, and preventive quality control measures are suggested.

Mucin (MUC) Family Influence on Acute Lymphoblastic Leukemia in Cancer and Non-Cancer Native American Populations from the Brazilian Amazon.

The mucin (MUC) family includes several genes aberrantly expressed in multiple carcinomas and mediates diverse pathways essentials for oncogenesis, in both solid and hematological malignancies. Acute Lymphoblastic Leukemia (ALL) can have its course influenced by genetic variants, and it seems more frequent in the Amerindian population, which has been understudied. Therefore, the present work aimed to investigate the MUC family exome in Amerindian individuals from the Brazilian Amazon, in a sample containing healthy Native Americans (NAMs) and indigenous subjects with ALL, comparing the frequency of polymorphisms between these two groups. The population was composed of 64 Amerindians from the Brazilian Amazon, from 12 different isolated tribes, five of whom were diagnosed with ALL. We analyzed 16 genes from the MUC family and found a total of 1858 variants. We compared the frequency of each variant in the ALL vs. NAM group, which led to 77 variants with a significant difference and, among these, we excluded those with a low impact, resulting in 63 variants, which were distributed in nine genes, concentrated especially in MUC 19 (n = 30) and MUC 3A (n = 18). Finally, 11 new variants were found in the NAM population. This is the first work with a sample of native Americans with cancer, a population which is susceptible to ALL, but remains understudied. The MUC family seems to have an influence on the development of ALL in the Amerindian population and especially MUC19 and MUC3A are shown as possible hotspots. In addition, the 11 new variants found point to the need to have their clinical impact analyzed.

Unraveling Signatures of Local Adaptation among Indigenous Groups from Mexico.

Few studies have addressed how selective pressures have shaped the genetic structure of the current Native American populations, and they have mostly limited their inferences to admixed Latin American populations. Here, we searched for local adaptation signals, based on integrated haplotype scores and population branch statistics, in 325 Mexican Indigenous individuals with at least 99% Native American ancestry from five previously defined geographical regions. Although each region exhibited its own local adaptation profile, only PPARG and AJAP1, both negative regulators of the Wnt/ß catenin signaling pathway, showed significant adaptation signals in all the tested regions. Several signals were found, mainly in the genes related to the metabolic processes and immune response. A pathway enrichment analysis revealed the overrepresentation of selected genes related to several biological phenotypes/conditions, such as the immune response and metabolic pathways, in agreement with previous studies, suggesting that immunological and metabolic pressures are major drivers of human adaptation. Genes related to the gut microbiome measurements were overrepresented in all the regions, highlighting the importance of studying how humans have coevolved with the microbial communities that colonize them. Our results provide a further explanation of the human evolutionary history in response to environmental pressures in this region.

Molecular Epidemiology in Amerindians of the Brazilian Amazon Reveals New Genetic Variants in DNA Repair Genes.

Native American populations from the Brazilian Amazon have a low genetic diversity and a different genetic profile when compared to people from other continents. Despite this, few studies have been conducted in this group, and there is no description of their genetic data in the various currently existent international databases. The characterization of the genomic profile of a population not only has an impact in studies of population genetics, but also helps to advance diagnostic and therapeutic response studies, leading to the optimization of clinical applicability. Genetic variations in DNA repair genes have been associated with the modulation of susceptibility to various pathologies, as well as in their prognosis and therapy. This is the first study to investigate DNA repair genes in Amerindians from the Brazilian Amazon region. We investigated 13 important DNA repair genes in the exome of 63 Native Americans, comparing our results with those found in 5 continental populations, whose data are available in the Genome Aggregation Database. Our results showed that 57 variants already described in literature were differentially distributed in the Amerindian populations in relation to the continental populations, 7 of which have significant clinical relevance. In addition, 9 new variants were described, suggesting that they are unique to these populations. Our study reinforces the understanding that the Amazonian Native American population presents a unique genetic profile, and our findings may collaborate with the creation of public policies that optimize the quality of life of these groups as well as the Brazilian population, which presents a high degree of interethnic mixing with Amerindian groups.

HLA molecular study of patients in a public kidney transplant program in Guatemala.

Guatemala is a country located in Central America, and while it is one of the most populated countries in the region, the genetic diversity of the population has been poorly analyzed. Currently, there are no analyses of the distribution of human leukocyte antigen (HLA) system alleles in mixed ancestry (i.e., ladino) populations in Guatemala. The HLA system exhibits the most extensive polymorphism in the human genome and has been extensively analyzed in a large number of studies related to disease association, transplantation, and population genetics (with particular importance in the understanding of diversity in the human population). Here, we present HLA typing data from 127 samples of unrelated individuals from the kidney transplant program of the San Juan de Dios General Hospital (Guatemala City) using a PCR-SSOP-based (PCR-sequence specific oligonucleotide probes) typing method. We found 16 haplotypes that accounted for 39.76 % of the total haplotype diversity, of which thirteen have been reported previously in Native American populations and three have been reported in European populations. The analyses showed no deviations from Hardy-Weinberg equilibrium, and admixture estimates calculated with k = 3 ancestral components showed that Native American was the most represented component, followed by the European component. The African component was less prominent in the Guatemala mixed ancestry sample in comparison to samples from other countries in Central America. The HLA-based admixture results for Central America showed a continuum in the distribution of Native American, European and African ancestries throughout the region, which is consistent with the complex demographic history of the region.

Pharmacogenomic Profile of Amazonian Amerindians.

Given the role of pharmacogenomics in the large variability observed in drug efficacy/safety, an assessment about the pharmacogenomic profile of patients prior to drug prescription or dose adjustment is paramount to improve adherence to treatment and prevent adverse drug reaction events. A population commonly underrepresented in pharmacogenomic studies is the Native American populations, which have a unique genetic profile due to a long process of geographic isolation and other genetic and evolutionary processes. Here, we describe the pharmacogenetic variability of Native American populations regarding 160 pharmacogenes involved in absorption, distribution, metabolism, and excretion processes and biological pathways of different therapies. Data were obtained through complete exome sequencing of individuals from 12 different Amerindian groups of the Brazilian Amazon. The study reports a total of 3311 variants; of this, 167 are exclusive to Amerindian populations, and 1183 are located in coding regions. Among these new variants, we found non-synonymous coding variants in the DPYD and the IFNL4 genes and variants with high allelic frequencies in intronic regions of the MTHFR, TYMS, GSTT1, and CYP2D6 genes. Additionally, 332 variants with either high or moderate (disruptive or non-disruptive impact in protein effectiveness, respectively) significance were found with a minimum of 1% frequency in the Amazonian Amerindian population. The data reported here serve as scientific basis for future design of specific treatment protocols for Amazonian Amerindian populations as well as for populations admixed with them, such as the Northern Brazilian population.

Characterization of PCLO Gene in Amazonian Native American Populations.

Genetic variations in PCLO have been associated with different pathologies in global literature, but there are no data regarding this gene in Native American populations. The Amazonian Native American populations have lower genetic diversity and are more different from other continental groups. We investigated 18 genetic variants in the PCLO gene in Amazonian indigenous and compared our results with the ones found in global populations, which were publicly available in the 1000 Genomes Project, gnmAD and ABraOM databases. The results demonstrated that the variants of the PCLO, especially rs17156844, rs550369696, rs61741659 and rs2877, have a significantly higher frequency in Amerindian populations in comparison with other continental populations. These data outline the singular genetic profile of the Native American population from the Brazilian Amazon region.

Genetic Distribution of Five Spinocerebellar Ataxia Microsatellite Loci in Mexican Native American Populations and Its Impact on Contemporary Mestizo Populations.

Spinocerebellar ataxias (SCAs) conform a heterogeneous group of neurodegenerative disorders with autosomal dominant inheritance. Five of the most frequent SCAs are caused by a CAG repeat expansion in the exons of specific genes. The SCAs incidence and the distribution of polymorphic CAG alleles vary among populations and ethnicities. Thus, characterization of the genetic architecture of ethnically diverse populations, which have undergone recent admixture and demographic events, could facilitate the identification of genetic risk factors. Owing to the great ethnic diversity of the Mexican population, this study aimed to analyze the allele frequencies of five SCA microsatellite loci (SCA1, SCA2, SCA3, SCA6, and SCA7) in eleven Mexican Native American (MNA) populations. Data from the literature were used to compare the allelic distribution of SCA loci with worldwide populations. The SCA loci allelic frequencies evidenced a certain genetic homogeneity in the MNA populations, except for Mayans, who exhibited distinctive genetic profiles. Neither pathological nor large normal alleles were found in MNA populations, except for the SCA2 pre-mutated allele in the Zapotec population. Collectively, our findings demonstrated the contribution of the MNA ancestry in shaping the genetic structure of contemporary Mexican Mestizo populations. Our results also suggest that Native American ancestry has no impact on the origin of SCAs in the Mexican population. Instead, the acquisition of pathological SCA alleles could be associated with European migration.

Uniparental Lineages from the Oldest Indigenous Population of Ecuador: The Tsachilas.

Together with Cayapas, the Tsachilas constitute the oldest population in the country of Ecuador and, according to some historians, they are the last descendants of the ancient Yumbos. Several anthropological issues underlie the interest towards this peculiar population: the uncertainty of their origin, their belonging to the Barbacoan linguistic family, which is still at the center of an intense linguistic debate, and the relations of their Yumbo ancestors with the Inca invaders who occupied their ancient territory. Our contribution to the knowledge of their complex past was the reconstruction of their genetic maternal and paternal inheritance through the sequencing of 70 entire mitochondrial genomes and the characterization of the non-recombinant region of the Y chromosome in 26 males. For both markers, we built comprehensive datasets of various populations from the surrounding geographical area, northwestern South America, NW, with a known linguistic affiliation, and we could then compare our sample against the overall variability to infer relationships with other Barbacoan people and with other NW natives. We found contrasting patterns of genetic diversity for the two markers, but generally, our results indicated a possible common origin between the Tsachilas, the Chachi, and other Ecuadorian and Colombian Barbacoans and are suggestive of an interesting ancient linkage to the Inca invaders in Yumbo country.

Genetic characterization of a collection of Tsantsas from Ecuadorian museums.

Tsantsas are shrunken human heads originally made for ceremonial purposes by Amazonian indigenous groups of the Shuar and Achuar family, previously called Jivaroan tribes. A significant demand of these objects during the first half of the 20th century led to the manufacture of counterfeit shrunken heads for commercial purposes. For museums where these collections are held, as well as for the indigenous groups who claim their ownership, it is important to identify the origin and authenticity of these tsantsas. We hypothesized that a collection of 14 tsantsas from 3 different museum collections in Ecuador are human and aimed to characterize their sex and potential origin. We amplified the amelogenin gene and performed a high resolution melting analysis to determine their human origin and characterize their sex. We also analyzed a fragment (16209-16402) from the HVR-1 region to identify the mtDNA haplogroups present in the tsantsa collection. Our exploratory results show that all the tsantsas are human and that the collection is comprised of 13 males and 1 female. A total of seven mtDNA haplogroups were found among the tsantsa collection using the mtDNA EMPOP database. These results show a predominance of the Amerindian mtDNA haplogroups B, C and D. Additional principal component analysis, genetic distance tree and haplotype network analyses suggest a relationship between the tsantsa specimens and Native American groups.

Study of 17 X-STRs in Native American and Mestizo populations of Central America for forensic and population purposes.

In the present work, an extensive analysis of the X-chromosomal pool of Native American and Mestizo groups of Central America (Guatemala, El Salvador, Nicaragua, and Panama) has been carried out. Allele and haplotype frequency databases, as well as other forensic parameters for these populations, are presented. The admixture analysis supports the tri-hybrid composition in terms of ancestry in the Mestizo populations, with a predominant Native American contribution (54-69%), followed by European (19-28%) and African contributions (12-19%). Pairwise FST genetic distances highlight the genetic proximity between the northernmost Central American populations, especially among admixed populations. The unique and complex nature of this area, where populations from different origins intercrossed, as well as the informativity of X-STR data, highpoint the great interest of this genetic study. Furthermore, the X-chromosome databases for Central American populations here provided will be not only useful for forensic and population purposes not only in the target countries but also in the host countries.

Maternal ancestry and hematological cancer risk: case-control study in an Argentinean population.

Aim: We investigated the role of maternal ancestry in neoplastic hematological malignancies (HMs) risk in a population from Central Argentina. Materials & methods: We analyzed 125 cases with HMs and 310 controls from a public hospital, and a set of 202 colorectal, breast, lung, and hematologic cancer patients from a private hospital. Results: A decreased risk for HMs was associated with the Native American haplogroup B2 (odds ratio = 0.49; 95% CI: 0.25-0.92; p = 0.02). The sub-Saharan African parahaplogroup L was associated with higher susceptibility for disease (odds ratio = 3.10; 95% CI: 1.04-9.31; p = 0.043). Although the mean ancestral proportions in the total studied population was as published (61.7% Native American, 34.6% European and 3.7% African), an unequal distribution was observed between hospitals. Conclusion: We confirmed the tri-hybrid nature of the Argentinean population, with proportions varying within the country. Our finding supports the notion that associated haplogroup is population and cancer specific.

Ancestral genetic legacy of the extant population of Argentina as predicted by autosomal and X-chromosomal DIPs.

Aiming to determine their ancestry diagnostic potential, we selected two sets of nuclear deletion/insertion polymorphisms (DIPs), including 30 located on autosomal chromosomes and 33 on the X chromosome. We analysed over 200 unrelated Argentinean individuals living in urban areas of Argentina. As in most American countries, the extant Argentinean population is the result of tricontinental genetic admixture. The peopling process within the continent was characterised by mating bias involving Native American and enslaved African females and European males. Differential results were detected between autosomal DIPs and X-DIPs. The former showed that the European component was the largest (77.8%), followed by the Native American (17.9%) and African (4.2%) components, in good agreement with the previously published results. In contrast, X-DIPs showed that the European genetic contribution was also predominant but much smaller (52.9%) and considerably larger Native American and African contributions (39.6% and 7.5%, respectively). Genetic analysis revealed continental genetic contributions whose associated phenotypic traits have been mostly lost. The observed differences between the estimated continental genetic contribution proportions based on autosomal DIPs and X-DIPs reflect the effects of autosome and X-chromosome transmission behaviour and their different recombination patterns. This work shows the ability of the tested DIP panels to infer ancestry and confirm mating bias. To the best of our knowledge, this is the first study focusing on ancestry-informative autosomal DIP and X-DIP comparisons performed in a sample representing the entire Argentinean population.

Pain Prevalence, Chronicity and Impact Within Subpopulations Based on Both Hispanic Ancestry and Race: United States, 2010-2017.

We provide national surveillance estimates of pain chronicity, severity and impact in adult subpopulations defined by both Hispanic Ancestry and Race. Data are from 144,434 adults who completed validated questionnaires in the 2010-2017 National Health Interview Survey asking about pain status within the last 3 (N = 84,664) or 6 months (N = 59,770). Multivariable logistic regression was used to assess the relationship between pain and ethnicity/race. Compared to White Puerto Rican participants, White participants with Central/South American and Mexican ancestry had reduced odds of reporting Category 3-4 pain and High-Impact Chronic Pain (HICP), while those of Cuban ancestry had reduced odds of only HICP - eg, White participants with Mexican ancestry had 32% lower odds of having Category 3-4 pain and 50% lower odds of having HICP. While no differences were seen between White Puerto Rican and White Non-Hispanic participants for Category 3-4 pain, White Non-Hispanics had 40% lower odds of reporting HICP. Asian Non-Hispanic and Black Non-Hispanic participants had significantly lower odds of reporting Category 3-4 pain and HICP compared to White Puerto Rican participants, eg, Black Non-Hispanic participants had 26% lower odds off having Category 3-4 pain and 42% lower odds of having HICP. Perspective: By examining pain status in discrete demographic groups based on Hispanic Ancestry and Race, this report further documents substantial difference in health status among underserved populations and provides a baseline for continuing surveillance research on pain, with the eventual goal of eliminating disparities in pain assessment and treatment.

Imputation Performance in Latin American Populations: Improving Rare Variants Representation With the Inclusion of Native American Genomes.

Current Genome-Wide Association Studies (GWAS) rely on genotype imputation to increase statistical power, improve fine-mapping of association signals, and facilitate meta-analyses. Due to the complex demographic history of Latin America and the lack of balanced representation of Native American genomes in current imputation panels, the discovery of locally relevant disease variants is likely to be missed, limiting the scope and impact of biomedical research in these populations. Therefore, the necessity of better diversity representation in genomic databases is a scientific imperative. Here, we expand the 1,000 Genomes reference panel (1KGP) with 134 Native American genomes (1KGP + NAT) to assess imputation performance in Latin American individuals of mixed ancestry. Our panel increased the number of SNPs above the GWAS quality threshold, thus improving statistical power for association studies in the region. It also increased imputation accuracy, particularly in low-frequency variants segregating in Native American ancestry tracts. The improvement is subtle but consistent across countries and proportional to the number of genomes added from local source populations. To project the potential improvement with a higher number of reference genomes, we performed simulations and found that at least 3,000 Native American genomes are needed to equal the imputation performance of variants in European ancestry tracts. This reflects the concerning imbalance of diversity in current references and highlights the contribution of our work to reducing it while complementing efforts to improve global equity in genomic research.