De Novo Genome Assemblies From Two Indigenous Americans from Arizona Identify New Polymorphisms in Non-Reference Sequences.

There is a collective push to diversify human genetic studies by including underrepresented populations. However, analyzing DNA sequence reads involves the initial step of aligning the reads to the GRCh38/hg38 reference genome which is inadequate for non-European ancestries. In this study, using long-read sequencing technology, we constructed de novo genome assemblies from two indigenous Americans from Arizona (IAZ). Each assembly included ∼17 Mb of DNA sequence not present [nonreference sequence (NRS)] in hg38, which consists mostly of repeat elements. Forty NRSs totaling 240 kb were uniquely anchored to the hg38 primary assembly generating a modified hg38-NRS reference genome. DNA sequence alignment and variant calling were then conducted with whole-genome sequencing (WGS) sequencing data from 387 IAZ using both the hg38 and modified hg38-NRS reference maps. Variant calling with the hg38-NRS map identified ∼50,000 single-nucleotide variants present in at least 5% of the WGS samples which were not detected with the hg38 reference map. We also directly assessed the NRSs positioned within genes. Seventeen NRSs anchored to regions including an identical 187 bp NRS found in both de novo assemblies. The NRS is located in HCN2 79 bp downstream of Exon 3 and contains several putative transcriptional regulatory elements. Genotyping of the HCN2-NRS revealed that the insertion is enriched in IAZ (minor allele frequency = 0.45) compared to other reference populations tested. This study shows that inclusion of population-specific NRSs can dramatically change the variant profile in an underrepresented ethnic groups and thereby lead to the discovery of previously missed common variations.

An E115A Missense Variant in CERS2 Is Associated With Increased Sleeping Energy Expenditure and Hepatic Insulin Resistance in American Indians.

Genetic determinants of interindividual differences in energy expenditure (EE) are largely unknown. Sphingolipids, such as ceramides, have been implicated in the regulation of human EE via mitochondrial uncoupling. In this study, we investigated whether genetic variants within enzymes involved in sphingolipid synthesis and degradation affect EE and insulin-related traits in a cohort of American Indians informative for 24-h EE and glucose disposal rates during a hyperinsulinemic-euglycemic clamp. Association analysis of 10,084 genetic variants within 28 genes involved in sphingolipid pathways identified a missense variant (rs267738, A>C, E115A) in exon 4 of CERS2 that was associated with higher sleeping EE (116 kcal/day) and increased rates of endogenous glucose production during basal (5%) and insulin-stimulated (43%) conditions, both indicators of hepatic insulin resistance. The rs267738 variant did not affect ceramide synthesis in HepG2 cells but resulted in a 30% decrease in basal mitochondrial respiration. In conclusion, we provide evidence that the CERS2 rs267738 missense variant may influence hepatic glucose production and postabsorptive sleeping metabolic rate.

Novel signals and polygenic score for height are associated with pubertal growth traits in Southwestern American Indians.

Most genetic variants associated with adult height have been identified through large genome-wide association studies (GWASs) in European-ancestry cohorts. However, it is unclear how these variants influence linear growth during adolescence. This study uses anthropometric and genotypic data from a longitudinal study conducted in an American Indian community in Arizona between 1965-2007. Growth parameters (i.e. height, velocity, and timing of growth spurt) were derived from the Preece-Baines growth model, a parametric growth curve fitted to longitudinal height data, in 787 participants with height measurements spanning the whole period of growth. Heritability estimates suggested that genetic factors could explain 25% to 71% of the variance of pubertal growth traits. We performed a GWAS of growth parameters, testing their associations with 5 077 595 imputed or directly genotyped variants. Six variants associated with height at peak velocity (P < 5 × 10-8, adjusted for sex, birth year and principal components). Implicated genes include NUDT3, previously associated with adult height, and PACSIN1. Two novel variants associated with duration of growth spurt (P < 5 × 10-8) in LOC105375344, an uncharacterized gene with unknown function. We finally examined the association of growth parameters with a polygenic score for height derived from 9557 single nucleotide polymorphisms (SNPs) identified in the GIANT meta-analysis for which genotypic data were available for the American Indian study population. Height polygenic score was correlated with the magnitude and velocity of height growth that occurred before and at the peak of the adolescent growth spurt, indicating overlapping genetic architecture, with no influence on the timing of adolescent growth.

Epistasis Between HLA-DRB1*16:02:01 and SLC16A11 T-C-G-T-T Reduces Odds for Type 2 Diabetes in Southwest American Indians.

We sought to identify genetic/immunologic contributors of type 2 diabetes (T2D) in an indigenous American community by genotyping all study participants for both high-resolution HLA-DRB1 alleles and SLC16A11 to test their risk and/or protection for T2D. These genes were selected based on independent reports that HLA-DRB1*16:02:01 is protective for T2D and that SLC16A11 associates with T2D in individuals with BMI <35 kg/m2. Here, we test the interaction of the two loci with a more complete data set and perform a BMI sensitivity test. We defined the risk protection haplotype of SLC16A11, T-C-G-T-T, as allele 2 of a diallelic genetic model with three genotypes, SLC16A11*11, *12, and *22, where allele 1 is the wild type. Both earlier findings were confirmed. Together in the same logistic model with BMI ≥35 kg/m2, DRB1*16:02:01 remains protective (odds ratio [OR] 0.73), while SLC16A11 switches from risk to protection (OR 0.57 [*22] and 0.78 [*12]); an added interaction term was statistically significant (OR 0.49 [*12]). Bootstrapped (b = 10,000) statistical power of interaction, 0.4801, yielded a mean OR of 0.43. Sensitivity analysis demonstrated that the interaction is significant in the BMI range of 30-41 kg/m2. To investigate the epistasis, we used the primary function of the HLA-DRB1 molecule, peptide binding and presentation, to search the entire array of 15-mer peptides for both the wild-type and ancient human SLC16A11 molecules for a pattern of strong binding that was associated with risk and protection for T2D. Applying computer binding algorithms suggested that the core peptide at SLC16A11 D127G, FSAFASGLL, might be key for moderating risk for T2D with potential implications for type 1 diabetes.

Genetic factors associated with suicidal behaviors and alcohol use disorders in an American Indian population.

American Indians (AI) demonstrate the highest rates of both suicidal behaviors (SB) and alcohol use disorders (AUD) among all ethnic groups in the US. Rates of suicide and AUD vary substantially between tribal groups and across different geographical regions, underscoring a need to delineate more specific risk and resilience factors. Using data from over 740 AI living within eight contiguous reservations, we assessed genetic risk factors for SB by investigating: (1) possible genetic overlap with AUD, and (2) impacts of rare and low-frequency genomic variants. Suicidal behaviors included lifetime history of suicidal thoughts and acts, including verified suicide deaths, scored using a ranking variable for the SB phenotype (range 0-4). We identified five loci significantly associated with SB and AUD, two of which are intergenic and three intronic on genes AACSP1, ANK1, and FBXO11. Nonsynonymous rare and low-frequency mutations in four genes including SERPINF1 (PEDF), ZNF30, CD34, and SLC5A9, and non-intronic rare and low-frequency mutations in genes OPRD1, HSD17B3 and one lincRNA were significantly associated with SB. One identified pathway related to hypoxia-inducible factor (HIF) regulation, whose 83 nonsynonymous rare and low-frequency variants on 10 genes were significantly linked to SB as well. Four additional genes, and two pathways related to vasopressin-regulated water metabolism and cellular hexose transport, also were strongly associated with SB. This study represents the first investigation of genetic factors for SB in an American Indian population that has high risk for suicide. Our study suggests that bivariate association analysis between comorbid disorders can increase statistical power; and rare and low-frequency variant analysis in a high-risk population enabled by whole-genome sequencing has the potential to identify novel genetic factors. Although such findings may be population specific, rare functional mutations relating to PEDF and HIF regulation align with past reports and suggest a biological mechanism for suicide risk and a potential therapeutic target for intervention.

Blood DNA methylation and liver cancer in American Indians: evidence from the Strong Heart Study.

PURPOSE: Liver cancer incidence among American Indians/Alaska Natives has risen over the past 20 years. Peripheral blood DNA methylation may be associated with liver cancer and could be used as a biomarker for cancer risk. We evaluated the association of blood DNA methylation with risk of liver cancer. METHODS: We conducted a prospective cohort study in 2324 American Indians, between age 45 and 75 years, from Arizona, Oklahoma, North Dakota and South Dakota who participated in the Strong Heart Study between 1989 and 1991. Liver cancer deaths (n = 21) were ascertained using death certificates obtained through 2017. The mean follow-up duration (SD) for non-cases was 25.1 (5.6) years and for cases, 11.0 (8.8) years. DNA methylation was assessed from blood samples collected at baseline using MethylationEPIC BeadChip 850 K arrays. We used Cox regression models adjusted for age, sex, center, body mass index, low-density lipoprotein cholesterol, smoking, alcohol consumption, and immune cell proportions to examine the associations. RESULTS: We identified 9 CpG sites associated with liver cancer. cg16057201 annotated to MRFAP1) was hypermethylated among cases vs. non-cases (hazard ratio (HR) for one standard deviation increase in methylation was 1.25 (95% CI 1.14, 1.37). The other eight CpGs were hypomethylated and the corresponding HRs (95% CI) ranged from 0.58 (0.44, 0.75) for cg04967787 (annotated to PPRC1) to 0.77 (0.67, 0.88) for cg08550308. We also assessed 7 differentially methylated CpG sites associated with liver cancer in previous studies. The adjusted HR for cg15079934 (annotated to LPS1) was 1.93 (95% CI 1.10, 3.39). CONCLUSIONS: Blood DNA methylation may be associated with liver cancer mortality and may be altered during the development of liver cancer.

Identification and characterization of the long non-coding RNA NFIA-AS2 as a novel locus for body mass index in American Indians.

BACKGROUND: Genome-wide association studies have shown that body mass index (BMI), an estimate of obesity, is highly polygenic. Individual variants typically have small effect sizes, making it challenging to identify unique loci in under-represented ethnic groups which lack statistical power due to their small sample size. Yet obesity is a major health disparity and is particularly prevalent in southwestern American Indians. Here, we identify and characterize a new locus for BMI that was detected by analyzing moderate associations with BMI obtained in a population-based sample of southwestern American Indians together with the well-powered GIANT dataset. METHODS: Genotypes for 10.5 million variants were tested for association with BMI in 5870 American Indians and 2600 variants that showed an association P < 10-3 in the American Indian sample were combined in a meta-analysis with the BMI data reported in GIANT (N = 240,608). The newly identified gene, NFIA-AS2 was functionally characterized, and the impact of its lead associated variant rs1777538 was studied both in-silico and in-vitro. RESULTS: Rs1777538 (T/C; C allele frequency = 0.16 in American Indians and 0.04 in GIANT, meta-analysis P = 5.0 × 10-7) exhibited a large effect in American Indians (1 kg/m2 decrease in BMI per copy of C allele). NFIA-AS2 was found to be a nuclear localized long non-coding RNA expressed in tissues pertinent to human obesity. Analysis of this variant in human brown preadipocytes showed that NFIA-AS2 transcripts carrying the C allele had increased RNA degradation compared to the T allele transcripts (half-lives = 9 h, 13 h respectively). During brown adipogenesis, NFIA-AS2 featured a stage-specific regulation of nearby gene expression where rs1777538 demonstrated an allelic difference in regulation in the mature adipocytes (the strongest difference was observed for L1TD1, P = 0.007). CONCLUSION: Our findings support a role for NFIA-AS2 in regulating pathways that impact BMI.

Spatiotemporal fluctuations of population structure in the Americas revealed by a meta-analysis of the first decade of archaeogenomes.

Objectives: Since 2010, genome-wide data from hundreds of ancient Native Americans have contributed to the understanding of Americas' prehistory. However, these samples have never been studied as a single dataset, and distinct relationships among themselves and with present-day populations may have never come to light. Here, we reassess genomic diversity and population structure of 223 ancient Native Americans published between 2010 and 2019. Materials and Methods: The genomic data from ancient Americas was merged with a worldwide reference panel of 278 present-day genomes from the Simons Genome Diversity Project and then analyzed through ADMIXTURE, D-statistics, PCA, t-SNE, and UMAP. Results: We find largely similar population structures in ancient and present-day Americas. However, the population structure of contemporary Native Americans, traced here to at least 10,000 years before present, is noticeably less diverse than their ancient counterparts, a possible outcome of the European contact. Additionally, in the past there were greater levels of population structure in North than in South America, except for ancient Brazil, which harbors comparatively high degrees of structure. Moreover, we find a component of genetic ancestry in the ancient dataset that is closely related to that of present-day Oceanic populations but does not correspond to the previously reported Australasian signal. Lastly, we report an expansion of the Ancient Beringian ancestry, previously reported for only one sample. Discussion: Overall, our findings support a complex scenario for the settlement of the Americas, accommodating the occurrence of founder effects and the emergence of ancestral mixing events at the regional level.

Functional characterization of a novel p.Ser76Thr variant in IGFBP4 that associates with body mass index in American Indians.

Insulin-like growth factor binding protein 4 (IGFBP4) is involved in adipogenesis, and IGFBP4 null mice have decreased body fat through decreased PPAR-γ expression. In the current study, we assessed whether variation in the IGFBP4 coding region influences body mass index (BMI) in American Indians who are disproportionately affected by obesity. Whole exome sequence data from a population-based sample of 6779 American Indians with longitudinal measures of BMI were used to identify variation in IGFBP4 that associated with BMI. A novel variant that predicts a p.Ser76Thr in IGFBP4 (Thr-allele frequency = 0.02) was identified which associated with the maximum BMI measured during adulthood (BMI 39.8 kg/m2 for Thr-allele homozygotes combined with heterozygotes vs. 36.2 kg/m2 for Ser-allele homozygotes, ß = 6.7% per Thr-allele, p = 8.0 × 10-5, adjusted for age, sex, birth-year and the first five genetic principal components) and the maximum age- and sex-adjusted BMI z-score measured during childhood/adolescence (z-score 0.70 SD for Thr-allele heterozygotes vs. 0.32 SD for Ser-allele homozygotes, ß = 0.37 SD per Thr-allele, p = 8.8 × 10-6). In vitro functional studies showed that IGFBP4 with the Thr-allele (BMI-increasing) had a 55% decrease (p = 0.0007) in FOXO-induced transcriptional activity, reflecting increased activation of the PI3K/AKT pathway mediated through increased IGF signaling. Over-expression and knock-down of IGFBP4 in OP9 cells during differentiation showed that IGFBP4 upregulates adipogenesis through PPARγ, CEBPα, AGPAT2 and SREBP1 expression. We propose that this American Indian specific variant in IGFBP4 affects obesity via an increase of IGF signaling.

The genomic landscape of Mexican Indigenous populations brings insights into the peopling of the Americas.

The genetic makeup of Indigenous populations inhabiting Mexico has been strongly influenced by geography and demographic history. Here, we perform a genome-wide analysis of 716 newly genotyped individuals from 60 of the 68 recognized ethnic groups in Mexico. We show that the genetic structure of these populations is strongly influenced by geography, and our demographic reconstructions suggest a decline in the population size of all tested populations in the last 15-30 generations. We find evidence that Aridoamerican and Mesoamerican populations diverged roughly 4-9.9 ka, around the time when sedentary farming started in Mesoamerica. Comparisons with ancient genomes indicate that the Upward Sun River 1 (USR1) individual is an outgroup to Mexican/South American Indigenous populations, whereas Anzick-1 was more closely related to Mesoamerican/South American populations than to those from Aridoamerica, showing an even more complex history of divergence than recognized so far.

Population diversity of three variants of the SLC47A2 gene (MATE2-K transporter) in Mexican Mestizos and Native Americans.

BACKGROUND: MATE2-K is an efflux transporter protein of organic cation expressed mainly in the kidney and encoded by the SLC47A2 gene. Different variants of this gene have shown an impact on the pharmacokinetics of various drugs, including metformin, which represents one of the most widely used drugs in treating type 2 diabetes. The SLC47A2 gene variants have been scarcely studied in Mexican populations, especially in Native American groups. For this reason, we analyzed the distribution of the variants rs12943590, rs35263947, and rs9900497 within the SLC47A2 gene in 173 Native Americans (Tarahumara, Huichol, Maya, Puerépecha) and 182 Mestizos (admixed) individuals from Mexico. METHODS AND RESULTS: Genotypes were determined through TaqMan probes (qPCR). The Hardy-Weinberg agreement was confirmed for all three SLC47A2 gene variants in all the Mexican populations analyzed. When worldwide populations were included for comparison purposes, for alleles and genotypes a relative interpopulation homogeneity was observed for rs35263947 (T allele; range 23.3-51.1%) and rs9900497 (T allele; range 18.6-40.9%). Conversely, heterogeneity was evident for rs12943590 (A allele, range 22.1-59.1%), where the most differentiated population was the Huichol, with high frequencies of the risk genotype associated with decreased response to metformin treatment (A/A = 40.9%). CONCLUSIONS: Although the SLC47A2 gene variants allow predicting favorable response to the metformin treatment in Mexican populations, the probable high frequency of ineffectiveness should be discarded in Huichols.

Association of rs9939609-FTO with metabolic syndrome components among women from Mayan communities of Chiapas, Mexico.

BACKGROUND: Metabolic syndrome (MetS) is a complex cluster of risk factors, considered as a polygenic and multifactorial entity. The objective of this study was to determine the association of rs9939609-FTO polymorphism and MetS components in adult women of Mayan communities of Chiapas. METHODS: In a cross-sectional study, sociodemographic, anthropometric, clinical, and biochemical data were obtained from 291 adult women from three regions of Chiapas, Mexico. The prevalence of MetS and the allele and genotype frequencies of the rs9939609-FTO were estimated. Multivariate logistic regression models were used to assess the association of the single nucleotide polymorphism (SNP) with each of the MetS components. RESULTS: The MetS prevalence was 60%. We found a statistically significant association between rs9939609-FTO and hyperglycemia in the dominant model (OR 2.6; 95% CI 1.3-5.3; p = 0.007). CONCLUSIONS: Women from Mayan communities of Chiapas presented a high prevalence of MetS and a relevant association of the FTO variant with hyperglycemia. This is the first study carried out in these Mayan indigenous communities from Chiapas.

The L125F MATE1 variant enriched in populations of Amerindian origin is associated with increased plasma levels of metformin and lactate.

Genetic factors that affect variability in metformin response have been poorly studied in the Latin American population, despite its being the initial drug therapy for type 2 diabetes, one of the most prevalent diseases in that region. Metformin pharmacokinetics is carried out by members of the membrane transporters superfamily (SLCs), being the multidrug and toxin extrusion protein 1 (MATE1), one of the most studied. Some genetic variants in MATE1 have been associated with reduced in vitro metformin transport. They include rs77474263 p.[L125F], a variant present at a frequency of 13.8% in Latin Americans, but rare worldwide (less than 1%). Using exome sequence data and TaqMan genotyping, we revealed that the Mexican population has the highest frequency of this variant: 16% in Mestizos and 27% in Amerindians, suggesting a possible Amerindian origin. To elucidate the metformin pharmacogenetics, a children cohort was genotyped, allowing us to describe, for the first time, a MATE1 rs77474263 TT homozygous individual. An additive effect of the L125F variant was observed on blood metformin accumulation, revealing the highest metformin and lactate serum levels in the TT homozygote, and intermediate metformin values in the heterozygotes. Moreover, a molecular dynamics analysis suggested that the genetic variant effect on metformin efflux could be due to a decreased protein permeability. We conclude that pharmacogenetics could be useful in enhancing metformin pharmacovigilance in populations having a high frequency of the risk genotype, especially considering that these populations also have a higher susceptibility to the diseases for which metformin is the first-choice drug.

Impact of Amerind ancestry and FADS genetic variation on omega-3 deficiency and cardiometabolic traits in Hispanic populations.

Long chain polyunsaturated fatty acids (LC-PUFAs) have critical signaling roles that regulate dyslipidemia and inflammation. Genetic variation in the FADS gene cluster accounts for a large portion of interindividual differences in circulating and tissue levels of LC-PUFAs, with the genotypes most strongly predictive of low LC-PUFA levels at strikingly higher frequencies in Amerind ancestry populations. In this study, we examined relationships between genetic ancestry and FADS variation in 1102 Hispanic American participants from the Multi-Ethnic Study of Atherosclerosis. We demonstrate strong negative associations between Amerind genetic ancestry and LC-PUFA levels. The FADS rs174537 single nucleotide polymorphism (SNP) accounted for much of the AI ancestry effect on LC-PUFAs, especially for low levels of n-3 LC-PUFAs. Rs174537 was also strongly associated with several metabolic, inflammatory and anthropomorphic traits including circulating triglycerides (TGs) and E-selectin in MESA Hispanics. Our study demonstrates that Amerind ancestry provides a useful and readily available tool to identify individuals most likely to have FADS-related n-3 LC-PUFA deficiencies and associated cardiovascular risk.

LCAT deficiency: a systematic review with the clinical and genetic description of Mexican kindred.

BACKGROUND: LCAT (lecithin-cholesterol acyltransferase) deficiency is characterized by two distinct phenotypes, familial LCAT deficiency (FLD) and Fish Eye disease (FED). This is the first systematic review evaluating the ethnic distribution of LCAT deficiency, with particular emphasis on Latin America and the discussion of three Mexican-Mestizo probands. METHODS: A systematic review was conducted following the PRISMA (Preferred Reporting Items for Systematic review and Meta-Analysis) Statement in Pubmed and SciELO. Articles which described subjects with LCAT deficiency syndromes and an assessment of the ethnic group to which the subject pertained, were included. RESULTS: The systematic review revealed 215 cases (154 FLD, 41 FED and 20 unclassified) pertaining to 33 ethnic/racial groups. There was no association between genetic alteration and ethnicity. The mean age of diagnosis was 42 ± 16.5 years, with fish eye disease identified later than familial LCAT deficiency (55 ± 13.8 vs. 41 ± 14.7 years respectively). The prevalence of premature coronary heart disease was significantly greater in FED vs. FLD. In Latin America, 48 cases of LCAT deficiency have been published from six countries (Argentina (1 unclassified), Brazil (38 FLD), Chile (1 FLD), Columbia (1 FLD), Ecuador (1 FLD) and Mexico (4 FLD, 1 FED and 1 unclassified). Of the Mexican probands, one showed a novel LCAT mutation. CONCLUSIONS: The systematic review shows that LCAT deficiency syndromes are clinically and genetically heterogeneous. No association was confirmed between ethnicity and LCAT mutation. There was a significantly greater risk of premature coronary artery disease in fish eye disease compared to familial LCAT deficiency. In FLD, the emphasis should be in preventing both cardiovascular disease and the progression of renal disease, while in FED, cardiovascular risk management should be the priority. The LCAT mutations discussed in this article are the only ones reported in the Mexican- Amerindian population.

Autosomal recessive diseases among the Athabaskans of the southwestern United States: anthropological, medical, and scientific aspects.

The peopling of the Americas by Native Americans occurred in 4 waves of which the last was Nadene language speakers of whom Athabaskans are the largest group. As the Europeans were entering the Southwestern states of the USA, Athabaskan hunting-gathering tribes were migrating South from Canada along the Rocky Mountains and undergoing potential bottlenecks reflected in autosomal recessive diseases shared by Apaches and Navajos. About 300 years ago, the Navajo developing a sedentary culture learned from Pueblo Indians while the Apache remained hunter-gathers. Although most of the tribe was rounded up and forced to relocate to Bosque Redondo, the adult breeding population was large enough to prevent a genetic bottleneck. However, some Navajo underwent further population bottlenecks while hiding from the brutal US Army action (under Kit Carson's guidance). This led to an increased frequency of other autosomal recessive diseases. Recent advances in population genetics, pathophysiology of the diseases, and social/ethical issues concerning their study are reviewed.

Applicability of ancestral genotyping in pharmacogenomic research with hormonal contraception.

To compare etonogestrel pharmacokinetic and pharmacodynamic outcomes by both self-reported race/ethnicity and genetically determined ancestry among contraceptive implant users. We conducted a secondary analysis of our parent pharmacogenomic study of 350 implant users. We genotyped these reproductive-aged (18-45 years) women for 88 ancestry-informative single nucleotide polymorphisms. We then assigned each participant a proportion value for African (AFR), European (EUR), and Indigenous American (AMR) ancestry based on reference population data. We correlated genetic ancestry with self-reported race/ethnicity and utilized genetic ancestry proportion values as variables for previously performed association analyses with serum etonogestrel concentrations and progestin-related side effects (e.g., bothersome bleeding and subjective weight gain). We successfully estimated genetically determined ancestry for 332 participants. EUR, AFR, and AMR ancestry were each highly correlated with self-reported White/non-Hispanic race (r = 0.64, p = 4.14 × 10-40 ), Black/African American race (r = 0.88, p = 1.36 × 10-107 ), and Hispanic/Latina ethnicity (r = 0.68, p = 4.03 × 10-47 ), respectively. Neither genetically determined ancestry nor self-reported race/ethnicity were significantly associated with serum etonogestrel concentrations. AFR ancestry and self-reported Black race had similar associations with reporting monthly periods (odds ratio [OR] 2.18, p = 0.09 vs. OR 2.22, p = 0.02) and having received treatment for bothersome bleeding (OR 5.19, p = 0.005 vs. OR 4.73, p = 2.0 × 10-4 ). In multivariable logistic regression for subjective weight gain, AMR ancestry dropped out of the model in preference for self-reported Hispanic/Latina ethnicity. We found no new associations between genetically determined ancestry and contraceptive implant pharmacodynamics/pharmacokinetics. Self-reported race/ethnicity were strong surrogates for genetically determined ancestry among this population of contraceptive implant users. Our data suggest that self-reported race/ethnicity, capturing societal and cultural aspects, remain important to the investigation of progestin-related side effects.

Genetic diversity, structure, and admixture in Mayans from Guatemala and Mexico based on 15 short tandem repeats.

OBJECTIVE: To analyze the genetic origin, relationships, structure, and admixture in Mayan Native American groups from Guatemala and Mexico based on 15 autosomal short tandem repeats (STRs) loci commonly used in human identification (HID). METHODS: We genotyped 513 unrelated Mayan samples from Guatemala based on 15 STR loci (AmpFlSTR® Identifiler kit). Moreover, we included 4408 genotypes previously reported, as following: Mayas from Guatemala and Mexico (n = 1666) and from Latin American, European, and African (n = 2742) populations. Forensic parameters, genetic distances, admixture, and population structure were assessed. RESULTS: Forensic parameters of the 15 STRs in different Mayan groups from Guatemala were reported. Low (Fst = 0.78%; p = 0.000) and non-significant differentiation (Fst = 1.8%; p = 0.108) were observed in Mayas from Guatemala and Mexico, respectively. The relative homogeneity observed among Mayan groups supported theories of extensive pre-Columbian gene flow and trade throughout the Mayan Empire. The distribution of the three Native American ancestries among these Mayan groups did not support the presumable Guatemalan origin of Tojolabal and Lacandon people (South, Mexico). The nonsignificant differentiation between Ladinos and Mayas suggests a relative panmixia in Guatemala. Mestizos from southeastern Mexico and Guatemala constitute a core of Native American ancestry in Latin America related to the Mayan Empire in Central America. CONCLUSIONS: The higher European admixture and homogeneity in Mexican Mayas of the Yucatan Peninsula suggest more intensive post-Columbian gene flow in this region than in Guatemalan Mayas.

Characterization of CYP3A pharmacogenetic variation in American Indian and Alaska Native communities, targeting CYP3A4*1G allele function.

The frequencies of genetic variants in the CYP3A4 and CYP3A5 genes differ greatly across global populations, leading to profound differences in the metabolic activity of these enzymes and resulting drug metabolism rates, with important consequences for therapeutic safety and efficacy. Yet, the impact of genetic variants on enzyme activity are incompletely described, particularly in American Indian and Alaska Native (AIAN) populations. To characterize genetic variation in CYP3A4 and CYP3A5 and its effect on enzyme activity, we partnered with AIAN people living in two regions of Alaska: Yup'ik Alaska Native people living in the Yukon-Kuskokwim Delta region of rural southwest Alaska and AIAN people receiving care at the Southcentral Foundation in Anchorage, Alaska. We identified low frequencies of novel and known variation in CYP3A4 and CYP3A5, including low frequencies of the CYP3A4*1G and CYP3A5*1 variants, and linkage disequilibrium patterns that differed from those we previously identified in an American Indian population in western Montana. We also identified increased activity of the CYP3A4*1G allele in vitro and in vivo. We demonstrated that the CYP3A4*1G allele confers increased protein content in human lymphoblastoid cells and both increased protein content and increased activity in human liver microsomes. We confirmed enhanced CYP3A4-mediated 4ß-vitamin D hydroxylation activity in Yup'ik people with the CYP3A4*1G allele. AIAN people in Alaska and Montana who carry the CYP3A4*1G allele-coupled with low frequency of the functional CYP3A5*1 variant-may metabolize CYP3A substrates more rapidly than people with the reference CYP3A4 allele.

"Distribution of paternal lineages in Mestizo populations throughout Mexico: an in silico study based on Y-STR haplotypes".

The Mexican-Mestizo population arose following European contact with the Americas due to the admixture of principally Spaniards, Native Americans, and Africans around 500 years ago. Because the paternal lineage distribution of the Mexican population has been poorly investigated, this study inferred the haplogroups of ten populations based on 1859 haplotypes (Y-STR data) using two haplogroup predictor programs. In the Mexican population sample, we found predominantly European ancestry (50.1%), followed by Native American (32.5%), Eurasian (13.4%), African (2.1%), East African-South Eurasian (1.3%), and Asian (0.6%) ancestries. In general, our results support a contrary north-to-south gradient throughout the Mexican territory of European and Native-American ancestries, respectively. Moreover, the presence of West-European R1b and Sub-Saharan African E1b1a haplogroups agrees with historical and genetic data of gene flow during the European conquest. This study represents the effort to analyze these paternal lineages on a large scale by taking advantage of Y-STR haplotype data to determine the distribution and ancestry proportions in this country.