Demographic history and biologically relevant genetic variation of Native Mexicans inferred from whole-genome sequencing

Selection scan in Native Americans of Mexico identifies FADS2 rs174616: Evidence of gene-diet interactions affecting lipid levels and Delta-6-desaturase activity

Searching for positive selection signals across genomes has identified functional genetic variants responding to environmental change. In Native Americans of Mexico, we used the fixation index (Fst) and population branch statistic (PBS) to identify SNPs suggesting positive selection. The 103 most differentiated SNPs were tested for associations with metabolic traits, the most significant association was FADS2/rs174616 with body mass index (BMI). This variant lies within a linkage disequilibrium (LD) block independent of previously reported FADS selection signals and has not been clearly associated with metabolic phenotypes. We tested this variant in two independent cohorts with cardiometabolic data. In the Genetics of Atherosclerotic Disease (GEA) cohort, the derived allele (T) was associated with increased BMI, lower LDL-C levels and a decreased risk of subclinical atherosclerosis in women. Significant gene-diet interactions affected lipid, apolipoprotein and adiponectin levels with differences according to sex, involving mainly total and complex dietary carbohydrate%. In the Genotype-related Effects of PUFA trial, the derived allele was associated with lower Δ-6 desaturase activity and erythrocyte membrane dihomo-gamma-linolenic acid (DGLA) levels, and with increased Δ-5 desaturase activity and eicosapentaenoic acid levels. This variant interacted with dietary carbohydrate% affecting Δ-6 desaturase activity. Notably, the relationship of DGLA and other erythrocyte membrane LC-PUFA indices with HOMA-IR differed according to rs174616 genotype, which has implications regarding how these indices should be interpreted. In conclusion, this observational study identified rs174616 as a signal suggesting selection in an independent linkage disequilibrium block, was associated with cardiometabolic and erythrocyte measurements of LC-PUFA in two independent Mexican cohorts and showed significant gene-diet interactions.

Frequencies of BCR::ABL1 Transcripts in Patients with Chronic Myeloid Leukemia: A Meta-Analysis

Chronic myeloid leukemia (CML) is associated with the Philadelphia chromosome and distinct BCR::ABL1 gene transcripts. We assessed the frequencies of these transcripts in Mexico, Latin America, and worldwide. We determined the prevalence of BCR::ABL1 transcripts in CML patients and intercontinental or regional variations using specialized databases and keywords. We analyzed 34 studies from 20 countries, encompassing 5795 patients. Keyword-based searches in specialized databases guided data collection. ANOVA was employed for transcript distribution analysis. The b3a2 transcript was most prevalent globally, followed by b2a2, with e1a2 being the least frequent. Interestingly, Mexico City exhibited a higher incidence of b2a2, while b3a2 predominated in the remaining country. Overall, no significant intercontinental or regional variations were observed. b3a2 was the most common BCR::ABL1 transcript worldwide, with b2a2 following closely; e1a2 was infrequent. Notably, this trend remained consistent in Mexico. Evaluating transcript frequencies holds clinical relevance for CML management. Understanding the frequency of transcript informs personalized CML treatments.

The Mexican Biobank Project promotes genetic discovery, inclusive science and local capacity building

Diversifying genotype-phenotype databases is essential to understanding complex trait and disease etiology across different environments and genetic ancestries. The rise of biobanks across the world is helping reveal the genetic and environmental architecture of multiple disease traits but the diversity they capture remains limited. To help close this gap, the Mexican Biobank (MXB) Project was recently generated, and has already revealed fine-scale genetic ancestries and demographic histories across the country, and their impact on trait-relevant genetic variation. This will help guide future genetic epidemiology and public health efforts, and has also improved polygenic prediction for several traits in Mexican populations compared with using data from other genome-wide association studies, such as the UK Biobank. The MXB illustrates the importance of transnational initiatives and funding calls that prioritize local leadership and capacity building to move towards inclusive genomic science.

Mexican Biobank advances population and medical genomics of diverse ancestries

Latin America continues to be severely underrepresented in genomics research, and fine-scale genetic histories and complex trait architectures remain hidden owing to insufficient data1. To fill this gap, the Mexican Biobank project genotyped 6,057 individuals from 898 rural and urban localities across all 32 states in Mexico at a resolution of 1.8 million genome-wide markers with linked complex trait and disease information creating a valuable nationwide genotype-phenotype database. Here, using ancestry deconvolution and inference of identity-by-descent segments, we inferred ancestral population sizes across Mesoamerican regions over time, unravelling Indigenous, colonial and postcolonial demographic dynamics2-6. We observed variation in runs of homozygosity among genomic regions with different ancestries reflecting distinct demographic histories and, in turn, different distributions of rare deleterious variants. We conducted genome-wide association studies (GWAS) for 22 complex traits and found that several traits are better predicted using the Mexican Biobank GWAS compared to the UK Biobank GWAS7,8. We identified genetic and environmental factors associating with trait variation, such as the length of the genome in runs of homozygosity as a predictor for body mass index, triglycerides, glucose and height. This study provides insights into the genetic histories of individuals in Mexico and dissects their complex trait architectures, both crucial for making precision and preventive medicine initiatives accessible worldwide.

Demographic history and genetic structure in pre-Hispanic Central Mexico

Aridoamerica and Mesoamerica are two distinct cultural areas in northern and central Mexico, respectively, that hosted numerous pre-Hispanic civilizations between 2500 BCE and 1521 CE. The division between these regions shifted southward because of severe droughts ~1100 years ago, which allegedly drove a population replacement in central Mexico by Aridoamerican peoples. In this study, we present shotgun genome-wide data from 12 individuals and 27 mitochondrial genomes from eight pre-Hispanic archaeological sites across Mexico, including two at the shifting border of Aridoamerica and Mesoamerica. We find population continuity that spans the climate change episode and a broad preservation of the genetic structure across present-day Mexico for the past 2300 years. Lastly, we identify a contribution to pre-Hispanic populations of northern and central Mexico from two ancient unsampled "ghost" populations.

The role of single nucleotide variant rs3819817 of the Histidine Ammonia-Lyase gene and 25-Hydroxyvitamin D on bone mineral density, adiposity markers, and skin pigmentation, in Mexican population

PURPOSE

Vitamin D (VD) deficiency and osteoporosis have become a global public health problem. A variant in the Histidine Ammonia-Lyase (HAL) gene has been associated with VD levels and bone mineral density (BMD). However, whether this variant has an influence on VD levels and BMD in Mexican adults remain unclear.

METHODS

This cross-sectional analysis included 1,905 adults participating in the Health Worker Cohort Study and 164 indigenous postmenopausal women from the Metabolic Analysis in an Indigenous Sample (MAIS) cohort. The rs3819817 variant was genotyped by TaqMan probe assay. Total 25 hydroxyvitamin D levels were measured by DiaSorin Liaison. BMD at the different sites was assessed through dual-energy X-ray absorptiometry. Linear and logistic regression models were performed to evaluate the associations of interest.

RESULTS

The prevalence of VD deficiency was 41%, showing differences between sexes. Obesity and skin pigmentation were associated with lower levels of VD in males and females. rs3819817-T allele was associated with low levels of 25-hydroxyvitamin D, VD deficiency, and hip and femoral neck BMD values (g/cm²). We found two interactions with VD levels, one between adiposity and rs3819817-T allele (P = 0.017) and another between skin pigmentation and rs3819817-T allele (P = 0.019). In indigenous postmenopausal women, we observed higher VD levels in the southern region compared to the northern region (P < 0.001); however, we did not observe differences by genotype.

CONCLUSION

Our findings confirm that the genetic variant rs3819817 has an essential function in VD levels and BMD and suggests a role in skin pigmentation in the Mexican population.

Aesthetic considerations when treating the Latin American patient: Thriving in diversity international roundtable series

BACKGROUND

The Hispanic/Latin American population is the fastest growing non-Caucasian group in the United States. Within this group, demand for aesthetic procedures is on the rise. High ethnic variability among these patients influences treatment approaches and patient priorities. Understanding these ethnic differences is central to providing optimal care.

AIMS

To discuss similarities and differences in anatomy and treatment preferences of Hispanic/Latin American patients both within the United States and internationally and explore how these differences may influence or inform aesthetic practices.

PATIENTS/METHODS

In support of clinicians who wish to serve a diverse patient population, a 6-part, international roundtable series focused on diversity in aesthetics was conducted from August 24, 2021 to May 16, 2022. In this roundtable, held in Medellin, Columbia, expert clinicians from across Latin America and the United States were invited to contribute and share best practices.

RESULTS

The results of the second roundtable in the series, the Latin American Patient, are described here. A special emphasis is placed on procedures that address the most commonly encountered concerns in these patients.

CONCLUSIONS

Hispanic and Latino patients represent a broad demographic with unique anatomical features, aesthetic preferences, and treatment priorities. Clinicians should consider these differences when treating this patient population.

A review of ancestrality and admixture in Latin America and the caribbean focusing on native American and African descendant populations

Genomics can reveal essential features about the demographic evolution of a population that may not be apparent from historical elements. In recent years, there has been a significant increase in the number of studies applying genomic epidemiological approaches to understand the genetic structure and diversity of human populations in the context of demographic history and for implementing precision medicine. These efforts have traditionally been applied predominantly to populations of European origin. More recently, initiatives in the United States and Africa are including more diverse populations, establishing new horizons for research in human populations with African and/or Native ancestries. Still, even in the most recent projects, the under-representation of genomic data from Latin America and the Caribbean (LAC) is remarkable. In addition, because the region presents the most recent global miscegenation, genomics data from LAC may add relevant information to understand population admixture better. Admixture in LAC started during the colonial period, in the 15th century, with intense miscegenation between European settlers, mainly from Portugal and Spain, with local indigenous and sub-Saharan Africans brought through the slave trade. Since, there are descendants of formerly enslaved and Native American populations in the LAC territory; they are considered vulnerable populations because of their history and current living conditions. In this context, studying LAC Native American and African descendant populations is important for several reasons. First, studying human populations from different origins makes it possible to understand the diversity of the human genome better. Second, it also has an immediate application to these populations, such as empowering communities with the knowledge of their ancestral origins. Furthermore, because knowledge of the population genomic structure is an essential requirement for implementing genomic medicine and precision health practices, population genomics studies may ensure that these communities have access to genomic information for risk assessment, prevention, and the delivery of optimized treatment; thus, helping to reduce inequalities in the Western Hemisphere. Hoping to set the stage for future studies, we review different aspects related to genetic and genomic research in vulnerable populations from LAC countries.

Challenges in selecting admixture models and marker sets to infer genetic ancestry in a Brazilian admixed population

The inference of genetic ancestry plays an increasingly prominent role in clinical, population, and forensic genetics studies. Several genotyping strategies and analytical methodologies have been developed over the last few decades to assign individuals to specific biogeographic regions. However, despite these efforts, ancestry inference in populations with a recent history of admixture, such as those in Brazil, remains a challenge. In admixed populations, proportion and components of genetic ancestry vary on different levels: (i) between populations; (ii) between individuals of the same population, and (iii) throughout the individual's genome. The present study evaluated 1171 admixed Brazilian samples to compare the genetic ancestry inferred by tri-/tetra-hybrid admixture models and evaluated different marker sets from those with small numbers of ancestry informative markers panels (AIMs), to high-density SNPs (HDSNP) and whole-genome-sequence (WGS) data. Analyses revealed greater variation in the correlation coefficient of ancestry components within and between admixed populations, especially for minority ancestral components. We also observed positive correlation between the number of markers in the AIMs panel and HDSNP/WGS. Furthermore, the greater the number of markers, the more accurate the tri-/tetra-hybrid admixture models.

Somatic Mutational Landscape in Mexican Patients: CDH1 Mutations and chr20q13.33 Amplifications Are Associated with Diffuse-Type Gastric Adenocarcinoma

The Hispanic population, compared with other ethnic groups, presents a more aggressive gastric cancer phenotype with higher frequency of diffuse-type gastric adenocarcinoma (GA); this could be related to the mutational landscape of GA in these patients. Using whole-exome sequencing, we sought to present the mutational landscape of GA from 50 Mexican patients who were treated at The Instituto Nacional de Cancerología from 2019 to 2020. We performed a comprehensive statistical analysis to explore the relationship of the genomic variants and clinical data such as tumor histology and presence of signet-ring cell, H. pylori, and EBV. We describe a potentially different mutational landscape between diffuse and intestinal GA in Mexican patients. Patients with intestinal-type GA tended to present a higher frequency of NOTCH1 mutations, copy number gains in cytobands 13.14, 10q23.33, and 12q25.1, and copy number losses in cytobands 7p12, 14q24.2, and 11q13.1; whereas patients with diffuse-type GA tended to present a high frequency of CDH1 mutations and CNV gains in cytobands 20q13.33 and 22q11.21. This is the first description of a mutational landscape of GA in Mexican patients to better understand tumorigenesis in Hispanic patients and lay the groundwork for discovering potential biomarkers and therapeutic targets.

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.

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.

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.

Dietary, Cultural and Pathogens-Related Selective Pressures Shaped Differential Adaptive Evolution Among Native Mexican Populations

Native American genetic ancestry has been remarkably implicated with increased risk of diverse health issues in several Mexican populations, especially in relation to the dramatic changes in environmental, dietary, and cultural settings they have recently undergone. In particular, the effects of these ecological transitions and Westernization of lifestyles have been investigated so far predominantly on Mestizo individuals. Nevertheless, indigenous groups, rather than admixed Mexicans, have plausibly retained the highest proportions of genetic components shaped by natural selection in response to the ancient milieu experienced by Mexican ancestors during their pre-Columbian evolutionary history. These formerly adaptive variants have the potential to represent the genetic determinants of some biological traits that are peculiar to Mexican people, as well as a reservoir of loci with possible biomedical relevance. To test such a hypothesis, we used genome-wide genotype data to infer the unique adaptive evolution of Native Mexican groups selected as reasonable descendants of the main pre-Columbian Mexican civilizations. A combination of haplotype-based and gene-network analyses enabled us to detect genomic signatures ascribable to polygenic adaptive traits plausibly evolved by the main genetic clusters of Mexican indigenous populations to cope with local environmental and/or cultural conditions. Some of these adaptations were found to play a role in modulating the susceptibility/resistance of these groups to certain pathological conditions, thus providing new evidence that diverse selective pressures have contributed to shape the current biological and disease-risk patterns of present-day Native and Mestizo Mexican populations.

Interaction of HLA Class II rs9272219 and TMPO rs17028450 (Arg690Cys) Variants Affects Neuromyelitis Optica Spectrum Disorder Susceptibility in an Admixed Mexican Population

Neuromyelitis Optica Spectrum Disorder (NMOSD) is a demyelinating autoimmune disease of the central nervous system, more prevalent in individuals of non-European ancestry. Few studies have analyzed genetic risk factors in NMOSD, and HLA class II gene variation has been associated NMOSD risk in various populations including Mexicans. Thymopoietin (TMPO) has not been tested as a candidate gene for NMOSD or other autoimmune disease, however, experimental evidence suggests this gene may be involved in negative selection of autoreactive T cells and autoimmunity. We thus investigated whether the missense TMPO variant rs17028450 (Arg630Cys, frequent in Latin America) is associated with NMOSD, and whether this variant shows an interaction with HLA-class II rs9272219, previously associated with NMOSD risk. A total of 119 Mexican NMOSD patients, 1208 controls and 357 Native Mexican individuals were included. The HLA rs9272219 “T” risk allele frequency ranged from 21 to 68%, while the rs17028450 “T” minor allele frequency was as high as 18% in Native Mexican groups. Both rs9272219 and rs17028450 were significantly associated with NMOSD risk under additive models (OR = 2.48; p = 8 × 10^–10 and OR = 1.59; p = 0.0075, respectively), and a significant interaction between both variants was identified with logistic regression models (p = 0.048). Individuals bearing both risk alleles had an estimated 3.9-fold increased risk of NMOSD. To our knowledge, this is the first study reporting an association of TMPO gene variation with an autoimmune disorder and the interaction of specific susceptibility gene variants, that may contribute to the genetic architecture of NMOSD in admixed Latin American populations.

Allele-dosage genetic polymorphisms of cannabinoid receptor 1 predict attention, but not working memory performance in humans

Attention and working memory (WM) are under high genetic regulation. Single nucleotide polymorphisms (SNPs) of the CNR1 gene, that encode for CB1R, have previously been shown to be related with individual differences in attentional control and WM. However, it remains unclear whether there is an allele-dosage or a dominant contribution of polymorphisms of CNR1 affecting attention and WM performance. This study evaluated the associations between attention and WM performance and three SNPs of CNR1: rs1406977, rs2180619, and rs1049353, previously associated with both processes. Healthy volunteers (n = 127) were asked to perform the Attention Network Task (ANT) to evaluate their overall attention and alerting, orienting, and executive systems, and the n-back task for evaluating their WM. All subjects were genotyped using qPCR with TaqMan assays; and dominant and additive models were assessed using the risk alleles of each SNP as the predictor variable. Results showed an individual association of the three SNPs with attention performance, but the composite genotype by the three alleles had the greatest contribution. Moreover, the additive-dosage model showed that for each G-allele added to the genotypic configuration, there was an increase in the percentage of correct responses respect to carriers who have no risk alleles in their genotypic configuration. The number of risk alleles in the genotypic configurations did not predict efficiency in any of the attention systems, nor in WM performance. Our model showed a contribution of three single nucleotide polymorphisms of the CNR1 gene to explain 9% of the variance of attention in an additive manner.

Whole genome variation in 27 Mexican indigenous populations, demographic and biomedical insights

There has been limited study of Native American whole genome diversity to date, which impairs effective implementation of personalized medicine and a detailed description of its demographic history. Here we report high coverage whole genome sequencing of 76 unrelated individuals, from 27 indigenous groups across Mexico, with more than 97% average Native American ancestry. On average, each individual has 3.26 million Single Nucleotide Variants and short indels, that together comprise a catalog of 9,737,152 variants, 44,118 of which are novel. We report 497 common Single Nucleotide Variants (with allele frequency > 5%) mapped to drug responses and 316,577 in enhancer or promoter elements; interestingly we found some of these enhancer variants in PPARG, a nuclear receptor involved in highly prevalent health problems in Mexican population, such as obesity, diabetes, and insulin resistance. By detecting signals of positive selection we report 24 enriched key pathways under selection, most of them related to immune mechanisms. No missense variants in ACE2, the receptor responsible for the entry of the SARS CoV-2 virus, were found in any individual. Population genomics and phylogenetic analyses demonstrated stratification in a Northern-Central-Southern axis, with major substructure in the Central region. The Seri, a northern group with the most genetic divergence in our study, showed a distinctive genomic context with the most novel variants, and the most population specific genotypes. Genome-wide analysis showed that the average haplotype blocks are longer in Native Mexicans than in other world populations. With this dataset we describe previously undetected population level variation in Native Mexicans, helping to reduce the gap in genomic data representation of such groups.

Population History and Gene Divergence in Native Mexicans Inferred from 76 Human Exomes

Native American genetic variation remains underrepresented in most catalogs of human genome sequencing data. Previous genotyping efforts have revealed that Mexico’s Indigenous population is highly differentiated and substructured, thus potentially harboring higher proportions of private genetic variants of functional and biomedical relevance. Here we have targeted the coding fraction of the genome and characterized its full site frequency spectrum by sequencing 76 exomes from five Indigenous populations across Mexico. Using diffusion approximations, we modeled the demographic history of Indigenous populations from Mexico with northern and southern ethnic groups splitting 7.2 KYA and subsequently diverging locally 6.5 and 5.7 KYA, respectively. Selection scans for positive selection revealed BCL2L13 and KBTBD8 genes as potential candidates for adaptive evolution in Rarámuris and Triquis, respectively. BCL2L13 is highly expressed in skeletal muscle and could be related to physical endurance, a well-known phenotype of the northern Mexico Rarámuri. The KBTBD8 gene has been associated with idiopathic short stature and we found it to be highly differentiated in Triqui, a southern Indigenous group from Oaxaca whose height is extremely low compared to other Native populations.

Differences in MTHFR and LRRK2 variant's association with sporadic Parkinson's disease in Mexican Mestizos correlated to Native American ancestry

Parkinson's disease (PD), a common neurodegenerative disorder, has a complex etiology where environmental and genetic factors intervene. While a number of genes and variants have been identified in recent decades as causative or protective agents of this condition, a limited number of studies have been conducted in mixed populations, such as Mexican Mestizos. The historical convergence of two founding groups and three ethnicities, and the increasing north-to-south gradient of Native American ancestry in Mexico resulted in a subpopulation structure with considerable genetic diversity. In this work, we investigate the influence of 21 known susceptibility variants for PD. Our case-control study, with a cohort of 311 Mexican Mestizo subjects, found a significant risk association for the variant rs1491942 in LRRK2. However, when stratification by ancestry was performed, a risk effect for MTHFR rs1801133 was observed only in the group with the highest percentage of European ancestry, and the PD risk effect for LRRK2 rs1491942 was significant in subjects with a higher ratio of Native American ancestry. Meta-analyses of these SNP revealed the effect of LRRK2 rs1491942 to be even more significant than previously described in populations of European descent. Although corroboration is necessary, our findings suggest that polymorphism rs1491942 may be useful as a risk marker of PD in Mexican Mestizos with greater Native American ancestry. The absence of associations with the remaining known risk factors is, in itself, a relevant finding and invites further research into the shared risk factors' role in the pathophysiological mechanisms of this neurodegenerative disorder.

Phenotypic and Genotypic Antibiotic Resistance Patterns in Helicobacter pylori Strains From Ethnically Diverse Population in México

Helicobacter pylori strains carry a range of mutations in genes that confer antimicrobial resistance and restrict the available options to treat the infection. Latin America is a region that conserve a large number of indigenous communities relatively isolated that practice a traditional medicine without consumption of drugs. We hypothesized that rates of antibiotic resistance are lower in these communities. Recent progress in whole-genome sequencing has allowed the study of drug susceptibility by searching for the known mutations associated with antibiotic resistance. The aim of this work was to study trends of antibiotic resistance over a 20-year period in Mexican H. pylori strains and to compare susceptibility between strains from Mexican mestizos and from indigenous population; we also aimed to learn the prevalence of mutational patterns in genes gyrA, gyrB, rdxA, frxA, rpsU, omp11, dppA, and 23S rRNA and its association with phenotypic tests. Resistance to clarithromycin, metronidazole, amoxicillin and levofloxacin was determined in167 H. pylori isolates by E-test, and the occurrence of mutational patterns in specific genes was determined by whole genome sequencing (WGS). The trend of resistance over 20 years in mestizo isolates showed significant resistant increase for clarithromycin and levofloxacin to frequencies that banned its clinical use. Resistance in H. pylori isolates of native communities was lower for all antibiotics tested. Phenotypic resistance showed good to moderate correlation with genotypic tests. Genetic methods for characterizing antibiotic resistance require further validation in each population.

Allele Frequency of ACE2 Intron Variants and Its Association with Blood Pressure

Angiotensin-converting enzyme 2 (ACE2) is known as the counter-regulator of the renin-angiotensin system, it cleaves angiotensin II to render Ag 1-7, a potent vasodilator with multiple roles in cardiovascular protection. A few studies have pinpointed ACE2 polymorphisms and their relationship with heart function and hypertension in a sex-dependent manner. These observations still lack replication mostly for admixed populations. This study aimed to report minor allele frequencies of four ACE2 intron variants, rs2285666, rs2048683, rs2106809, and rs4240157, derived from previous research using the GSA, v1.0, microarray in 1231 hypertensive and nonhypertensive patients. Logistic and multiple linear regression models were developed to identify potential associations with hypertension status and systolic and diastolic blood pressure (SBP and DBP). Allele frequency differences were identified for ACE2 rs2048683 and rs4240157 in populations with European ancestry and people of the Americas. Regression analyses identified a significant association of ACE2 rs2048683 and rs4240157 with SBP/DBP in males or females. Our observations confirm sex differences in ACE2 genetic associations with SBP and DBP and contribute to the collection of genetic variation in ACE2 for admixed populations.

Native American ancestry significantly contributes to neuromyelitis optica susceptibility in the admixed Mexican population

Neuromyelitis Optica (NMO) is an autoimmune disease with a higher prevalence in non-European populations. Because the Mexican population resulted from the admixture between mainly Native American and European populations, we used genome-wide microarray, HLA high-resolution typing and AQP4 gene sequencing data to analyze genetic ancestry and to seek genetic variants conferring NMO susceptibility in admixed Mexican patients. A total of 164 Mexican NMO patients and 1,208 controls were included. On average, NMO patients had a higher proportion of Native American ancestry than controls (68.1% vs 58.6%; p = 5 × 10^-6). GWAS identified a HLA region associated with NMO, led by rs9272219 (OR = 2.48, P = 8 × 10^-10). Class II HLA alleles HLA-DQB1*03:01, -DRB1*08:02, -DRB1*16:02, -DRB1*14:06 and -DQB1*04:02 showed the most significant associations with NMO risk. Local ancestry estimates suggest that all the NMO-associated alleles within the HLA region are of Native American origin. No novel or missense variants in the AQP4 gene were found in Mexican patients with NMO or multiple sclerosis. To our knowledge, this is the first study supporting the notion that Native American ancestry significantly contributes to NMO susceptibility in an admixed population, and is consistent with differences in NMO epidemiology in Mexico and Latin America.

The immunogenetic diversity of the HLA system in Mexico correlates with underlying population genetic structure

We studied HLA class I (HLA-A, -B) and class II (HLA-DRB1, -DQB1) allele groups and alleles by PCR-SSP based typing in a total of 15,318 mixed ancestry Mexicans from all the states of the country divided into 78 sample sets, providing information regarding allelic and haplotypic frequencies and their linkage disequilibrium, as well as admixture estimates and genetic substructure. We identified the presence of 4268 unique HLA extended haplotypes across Mexico and find that the ten most frequent (HF > 1%) HLA haplotypes with significant linkage disequilibrium (Δ'≥0.1) in Mexico (accounting for 20% of the haplotypic diversity of the country) are of primarily Native American ancestry (A*02~B*39~DRB1*04~DQB1*03:02, A*02~B*35~DRB1*08~DQB1*04, A*68~B*39~DRB1*04~DQB1*03:02, A*02~B*35~DRB1*04~DQB1*03:02, A*24~B*39~DRB1*14~DQB1*03:01, A*24~B*35~DRB1*04~DQB1*03:02, A*24~B*39~DRB1*04~DQB1*03:02, A*02~B*40:02~DRB1*04~DQB1*03:02, A*68~B*35~DRB1*04~DQB1*03:02, A*02~B*15:01~DRB1*04~DQB1*03:02). Admixture estimates obtained by a maximum likelihood method using HLA-A/-B/-DRB1 as genetic estimators revealed that the main genetic components in Mexico as a whole are Native American (ranging from 37.8% in the northern part of the country to 81.5% in the southeastern region) and European (ranging from 11.5% in the southeast to 62.6% in northern Mexico). African admixture ranged from 0.0 to 12.7% not following any specific pattern. We were able to detect three major immunogenetic clusters correlating with genetic diversity and differential admixture within Mexico: North, Central and Southeast, which is in accordance with previous reports using genome-wide data. Our findings provide insights into the population immunogenetic substructure of the whole country and add to the knowledge of mixed ancestry Latin American population genetics, important for disease association studies, detection of demographic signatures on population variation and improved allocation of public health resources.

The Identification of Admixture Patterns Could Refine Pharmacogenetic Counseling: Analysis of a Population-Based Sample in Mexico

Pharmacogenetic analysis has generated translational data that could be applied to guide treatments according to individual genetic variations. However, pharmacogenetic counseling in some mestizo (admixed) populations may require tailoring to different patterns of admixture. The identification and clustering of individuals with related admixture patterns in such populations could help to refine the practice of pharmacogenetic counseling. This study identifies related groups in a highly admixed population-based sample from Mexico, and analyzes the differential distribution of actionable pharmacogenetic variants. A subsample of 1728 individuals from the Mexican Genomic Database for Addiction Research (MxGDAR/Encodat) was analyzed. Genotyping was performed with the commercial PsychArray BeadChip, genome-wide ancestry was estimated using EIGENSOFT, and model-based clustering was applied to defined admixture groups. Actionable pharmacogenetic variants were identified with a query to the Pharmacogenomics Knowledge Base (PharmGKB) database, and functional prediction using the Variant Effect Predictor (VEP). Allele frequencies were compared with chi-square tests and differentiation was estimated by F_ST. Seven admixture groups were identified in Mexico. Some, like Group 1, Group 4, and Group 5, were found exclusively in certain geographic areas. More than 90% of the individuals, in some groups (Group 1, Group 4 and Group 5) were found in the Central-East and Southeast region of the country. MTRR p.I49M, ABCG2 p.Q141K, CHRNA5 p.D398N, SLCO2B1 rs2851069 show a low degree of differentiation between admixture groups. ANKK1 p.G318R and p.H90R, had the lowest allele frequency of Group 1. The reduction in these alleles reduces the risk of toxicity from anticancer and antihypercholesterolemic drugs. Our analysis identified different admixture patterns and described how they could be used to refine the practice of pharmacogenetic counseling for this admixed population.

Variation in Actionable Pharmacogenetic Markers in Natives and Mestizos From Mexico

The identification and characterization of pharmacogenetic variants in Latin American populations is still an ongoing endeavor. Here, we investigated SNVs on genes listed by the Pharmacogenomics Knowledge Base in 1284 Mestizos and 94 Natives from Mexico. Five institutional cohorts with NGS data were retrieved from different research projects at INMEGEN, sequencing files were filtered for 55 pharmacogenes present in all cohorts to identify novel and known variation. Bioinformatic tools VEP, PROVEAN, and FATHMM were used to assess, in silico, the functional impact of this variation. Next, we focused on 17 genes with actionable variants that have been clinically implemented. Allele frequencies were compared with major continental groups and differences discussed in the scope of a pharmacogenomic impact. We observed a wide genetic variability for known and novel SNVs, the largest variation was on UGT1A > ACE > COMT > ABCB1 and the lowest on APOE and NAT2. Although with allele frequencies around 1%, novel variation was observed in 16 of 17 PGKB genes. In Natives we identified 59 variants and 58 in Mestizos. Several genes did not show novel variation, on CYP2B6, CYP2D6, and CYP3A4 in Natives; and APOE, UGT1A, and VKORC1 in Mestizos. Similarities in allele frequency, comparing major continental groups for VIP pharmacogenes, hint towards a comparable PGx for drugs metabolized by UGT1A1, DPYD, ABCB1, CBR3, COMT, and TPMT; in contrast to variants on CYP3A5 and CYP2B6 for which significant MAF differences were identified. Our observations offer some discernment into the extent of pharmacogenetic variation registered up-to-date in Mexicans and contribute to quantitatively dissect actionable pharmacogenetic variants in Natives and Mestizos.

HLA concordance between hematopoietic stem cell transplantation patients and umbilical cord blood units: Implications for cord blood banking in admixed populations

Umbilical cord blood stem cell transplantation is an important choice for treating a variety of hematopoietic, neoplastic, and genetic disorders. The optimal size for a cord blood bank to provide matching units for 80% of patients requiring a stem cell transplantation procedure depends on the particular characteristics of each population. In this study, we analyzed the immunogenetic diversity of a sample set of Mexican patients suffering from blood, hematopoietic, and immunological diseases, to assess the best strategy for cord blood banking. For achieving that, we analyzed HLA-A, HLA-B, HLA-DRB1, and HLA-DQB1 genotype and allele frequencies of both units from the bioarchive of the Umbilical Cord Blood Bank from La Raza and patients requiring a stem cell transplant and compared these variables with data from the same geographic and genetic context. We were able to detect significant differences for at least half of the alleles were observed for HLA class I and class II genes between units and patients. Five Native American haplotypes had lower frequencies in patients sample than in the cord blood units. Genetic admixture estimations for both groups showed a higher contribution of Native American component in the cord blood units. Differences in ancestral components in the Umbilical Cord Blood Bank from La Raza and six virtual banks modeled from a pool of Mexican mixed ancestry individuals show that genetic background is important in cord blood collection. In conclusion, increasing diversity over quantity of new cord blood units will improve the cost effectiveness of cord blood banking and health policies regarding hematopoietic stem cell transplantation in admixed populations such as those present in Latin American countries.

Microbiome-MX 2018: microbiota and microbiome opportunities in Mexico, a megadiverse country

A weekly conference series paired with lectures entitled "Microbiome-MX: exploring the Microbiota and Microbiome Research in Mexico" was organized to provide a multidisciplinary overview of the most recent research done in Mexico using high-throughput sequencing. Scientists and postgraduate students from several disciplines such as microbiology, bioinformatics, virology, immunology, nutrition, and medical genomics gathered to discuss state of the art in each of their respective subjects of expertise, as well as advances, applications and new opportunities on microbiota/microbiome research. In particular, high-throughput sequencing is a crucial tool to understand the challenges of a megadiverse developing country as Mexico, and moreover to know the scientific capital and capabilities available for collaboration. The conference series addressed three main topics important for Mexico: i) the complex role of microbiota in health and prevalent diseases such as obesity, diabetes, inflammatory bowel disease, tuberculosis, HIV, autoimmune diseases and gastric cancer; ii) the use of local, traditional and prehispanic products as pre/probiotics to modulate the microbiota and improve human health; and iii) the impact of the microbiota in shaping the biodiversity of economically important terrestrial and marine ecosystems. Herein, we summarize the contributions that Mexican microbiota/microbiome research is making to the global trends, describing the highlights of the conferences and lectures, rather than a review of the state-of-the-art of this research. This meeting report also presents the efforts of a multidisciplinary group of scientist to encourage collaborations and bringing this research field closer for younger generations.

Whole Genome Sequence, Variant Discovery and Annotation in Mapuche-Huilliche Native South Americans

Whole human genome sequencing initiatives help us understand population history and the basis of genetic diseases. Current data mostly focuses on Old World populations, and the information of the genomic structure of Native Americans, especially those from the Southern Cone is scant. Here we present annotation and variant discovery from high-quality complete genome sequences of a cohort of 11 Mapuche-Huilliche individuals (HUI) from Southern Chile. We found approximately 3.1 × 106 single nucleotide variants (SNVs) per individual and identified 403,383 (6.9%) of novel SNVs events. Analyses of large-scale genomic events detected 680 copy number variants (CNVs) and 4,514 structural variants (SVs), including 398 and 1,910 novel events, respectively. Global ancestry composition of HUI genomes revealed that the cohort represents a sample from a marginally admixed population from the Southern Cone, whose main genetic component derives from Native American ancestors. Additionally, we found that HUI genomes contain variants in genes associated with 5 of the 6 leading causes of noncommunicable diseases in Chile, which may have an impact on the risk of prevalent diseases in Chilean and Amerindian populations. Our data represents a useful resource that can contribute to population-based studies and for the design of early diagnostics or prevention tools for Native and admixed Latin American populations.

Latin Americans show wide-spread Converso ancestry and imprint of local Native ancestry on physical appearance

Historical records and genetic analyses indicate that Latin Americans trace their ancestry mainly to the intermixing (admixture) of Native Americans, Europeans and Sub-Saharan Africans. Using novel haplotype-based methods, here we infer sub-continental ancestry in over 6,500 Latin Americans and evaluate the impact of regional ancestry variation on physical appearance. We find that Native American ancestry components in Latin Americans correspond geographically to the present-day genetic structure of Native groups, and that sources of non-Native ancestry, and admixture timings, match documented migratory flows. We also detect South/East Mediterranean ancestry across Latin America, probably stemming mostly from the clandestine colonial migration of Christian converts of non-European origin (Conversos). Furthermore, we find that ancestry related to highland (Central Andean) versus lowland (Mapuche) Natives is associated with variation in facial features, particularly nose morphology, and detect significant differences in allele frequencies between these groups at loci previously associated with nose morphology in this sample.

Genetic diversity in populations across Latin America: implications for population and medical genetic studies

Hispanic/Latino (H/L) populations, although linked by culture and aspects of shared history, reflect the complexity of history and migration influencing the Americas. The original settlement by indigenous Americans, followed by postcolonial admixture from multiple continents, has yielded localized genetic patterns. In addition, numerous H/L populations appear to have signatures of pre-colonization and post-colonization bottlenecks, indicating that tens of millions of H/Ls may harbor signatures of founder effects today. Based on both population and medical genetic findings we highlight the extreme differentiation across the Americas, providing evidence for why H/Ls should not be considered a single population in modern human genetics. We highlight the need for additional sampling of understudied H/L groups, and ramifications of these findings for genomic medicine in one-tenth of the world's population.

Exploring Cuba's population structure and demographic history using genome-wide data

Cuba is the most populated country in the Caribbean and has a rich and heterogeneous genetic heritage. Here, we take advantage of dense genomic data from 860 Cuban individuals to reconstruct the genetic structure and ancestral origins of this population. We found distinct admixture patterns between and within the Cuban provinces. Eastern provinces have higher African and Native American ancestry contributions (average 26% and 10%, respectively) than the rest of the Cuban provinces (average 17% and 5%, respectively). Furthermore, in the Eastern Cuban region, we identified more intense sex-specific admixture patterns, strongly biased towards European male and African/Native American female ancestries. Our subcontinental ancestry analyses in Cuba highlight the Iberian population as the best proxy European source population, South American and Mesoamerican populations as the closest Native American ancestral component, and populations from West Central and Central Africa as the best proxy sources of the African ancestral component. Finally, we found complex admixture processes involving two migration pulses from both Native American and African sources. Most of the inferred Native American admixture events happened early during the Cuban colonial period, whereas the African admixture took place during the slave trade and more recently as a probable result of large-scale migrations from Haiti.

Blood Groups Distribution and Gene Diversity of the ABO and Rh (D) Loci in the Mexican Population

Objective

To determine the frequency and distribution of ABO and Rh (D) antigens and, additionally, investigate gene diversity and the structure of Mexican populations.

Materials and Methods

Blood groups were tested in 271,164 subjects from 2014 to 2016. The ABO blood group was determined by agglutination using the antibodies anti-A, Anti-B, and Anti-D for the Rh factor, respectively.

Results

The overall distribution of ABO and Rh (D) groups in the population studied was as follows: O: 61.82%; A: 27.44%; B: 8.93%; and AB: 1.81%. For the Rh group, 95.58% of people were Rh (D), and 4.42% were Rh (d). Different distributions of blood groups across regions were found; additionally, genetic analysis revealed that the I^O and I^D allele showed an increasing trend from the north to the center, while the I^A and I^d allele tended to increase from the center to the north. Also, we found more gene diversity in both loci in the north compared with the center, suggesting population structure in Mexico.

Conclusion

This work could help health institutions to identify where they can obtain blood products necessary for medical interventions. Moreover, this piece of information contributes to the knowledge of the genetic structure of the Mexican populations which could have significant implications in different fields of biomedicine.