Human genetics. The genetics of Mexico recapitulates Native American substructure and affects biomedical traits

TwinsMX: Exploring the Genetic and Environmental Influences on Health Traits in the Mexican Population

TwinsMX registry is a national research initiative in Mexico that aims to understand the complex interplay between genetics and environment in shaping physical and mental health traits among the country's population. With a multidisciplinary approach, TwinsMX aims to advance our knowledge of the genetic and environmental mechanisms underlying ethnic variations in complex traits and diseases, including behavioral, psychometric, anthropometric, metabolic, cardiovascular and mental disorders. With information gathered from over 2800 twins, this article updates the prevalence of several complex traits; and describes the advances and novel ideas we have implemented such as magnetic resonance imaging. The future expansion of the TwinsMX registry will enhance our comprehension of the intricate interplay between genetics and environment in shaping health and disease in the Mexican population. Overall, this report describes the progress in the building of a solid database that will allow the study of complex traits in the Mexican population, valuable not only for our consortium, but also for the worldwide scientific community, by providing new insights of understudied genetically admixed populations.

Mother's smoking habits affects IL10 methylation but not asthma in Ecuadorian children

There is no evidence evaluating the IL10 epigenetic upregulation among mestizo children in a high-altitude Andean city in Latin America.

OBJECTIVE

To identify polymorphisms and methylation profiles in the IL10 gene associated with asthma in children aged 5 to 11.

METHODS

A case-control study was conducted with asthmatic and non-asthmatic children aged 5 to 11 years in Cuenca-Ecuador. Data on allergic diseases and risk factors were collected through a questionnaire for parents. Atopy was measured by skin prick test (SPT) to relevant aeroallergens. Three IL10 single nucleotide polymorphisms were evaluated in all participants, and methylation analysis was performed in 54 participants. Association between risk factors, allergic diseases and genetic factors were estimated using multivariate logistic regression.

RESULTS

The results of polymorphisms showed no differences between cases and controls when comparing the SNPs rs3024495, rs3024496, rs1800896 allelic and genotypic frequencies. In the methylation analysis, no differences in the IL10 methylation profile were found between cases and controls; however, the multivariate analysis showed an association between the mother's smoking habits and the IL10 methylation profile.

CONCLUSION

Smoking habit could be essential as an environmental exposure factor in regulating gene expression in children with asthma.

STAT4 Gene Variant rs7574865 Is Associated with Rheumatoid Arthritis Activity and Anti-CCP Levels in the Western but Not in the Southern Population of Mexico

Rheumatoid Arthritis (RA) is a multifactorial autoimmune disease. Currently, several genes play an important role in the development of the disease. The objective was to evaluate the association of the STAT4 rs7574865 and rs897200 gene variants with RA susceptibility, DAS28, RF, and anti-CCP in Western and Southern Mexico populations. Genotyping was performed on 476 samples (cases = 240; controls = 236) using the Taqman® system and qPCR probes. Disease activity was assessed using DAS28 and HAQ DI. CRP, ESR, RF, and anti-CCP were determined for clinical assessment. Our study showed there is a statistically significant association with susceptibility to RA for the rs7574865 variant in the Western population for the GT and TT genotypes. The same genotypes also showed a moderate-to-high activity according to DAS28 and positive anti-CCP compared to the control group. This association was not found in the Southern population. This work confirms the association of the rs7574865 variant with RA, as well as a moderate-to-high activity and positive anti-CCP in the Western population but not in the Southern population. No association of the rs897200 variant was found in any of the studied populations.

COMT and SCN9A gene variants do not contribute to chronic low back pain in Mexican-Mestizo patients

BACKGROUND

Chronic low back pain (CLBP) is a complex condition in which genetic factors play a role in its susceptibility. Catechol-O-methyltransferase (COMT) and sodium channel NaV1.7 (SCN9A) genes are implicated in pain perception. The aim is to analyze the association of COMT and SCN9A with CLBP and their interaction, in a Mexican-Mestizo population.

METHODS

A case-control study was conducted. Cases corresponded to adults of both sexes with CLBP. Controls were adults with no CLBP. Variants of SCN9A and COMT were genotyped. Allelic and genotypic frequencies and Hardy-Weinberg equilibrium (HWE) were calculated. Association was tested under codominant, dominant, and recessive models. Multifactor dimensionality reduction was developed to detect epistasis.

RESULTS

Gene variants were in HWE, and there was no association under different inheritance models in the whole sample. In women, in codominant and dominant models, a trend to a high risk was observed for AA of rs4680 of COMT (OR = 1.7 [0.5-5.3] and 1.6 [0.7-3.4]) and for TT of rs4633 (OR = 1.6 [0.7-3.7] and 1.6 [0.7-3.4]). In men, a trend to low risk was observed for AG genotype of rs4680 in the same models (OR = 0.6 [0.2-1.7] and 0.7 [0.3-1.7]), and for TC genotype of rs4633 in the codominant model (OR = 0.6 [0.2-1.7]). In the interaction analysis, a model of the SCN9A and COMT variants showed a CVC of 10/10; however, the TA was 0.4141.

CONCLUSION

COMT and SCN9A variants are not associated with CLBP in the analyzed Mexican-Mestizo population.

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.

Native American ancestry and breast cancer risk in Colombian and Mexican women: ruling out potential confounding through ancestry-informative markers

BACKGROUND

Latin American and Hispanic women are less likely to develop breast cancer (BC) than women of European descent. Observational studies have found an inverse relationship between the individual proportion of Native American ancestry and BC risk. Here, we use ancestry-informative markers to rule out potential confounding of this relationship, estimating the confounder-free effect of Native American ancestry on BC risk.

METHODS AND STUDY POPULATION

We used the informativeness for assignment measure to select robust instrumental variables for the individual proportion of Native American ancestry. We then conducted separate Mendelian randomization (MR) analyses based on 1401 Colombian women, most of them from the central Andean regions of Cundinamarca and Huila, and 1366 Mexican women from Mexico City, Monterrey and Veracruz, supplemented by sensitivity and stratified analyses.

RESULTS

The proportion of Colombian Native American ancestry showed a putatively causal protective effect on BC risk (inverse variance-weighted odds ratio [OR] = 0.974 per 1% increase in ancestry proportion, 95% confidence interval [CI] 0.970-0.978, p = 3.1 × 10-40). The corresponding OR for Mexican Native American ancestry was 0.988 (95% CI 0.987-0.990, p = 1.4 × 10-44). Stratified analyses revealed a stronger association between Native American ancestry and familial BC (Colombian women: OR = 0.958, 95% CI 0.952-0.964; Mexican women: OR = 0.973, 95% CI 0.969-0.978), and stronger protective effects on oestrogen receptor (ER)-positive BC than on ER-negative and triple-negative BC.

CONCLUSIONS

The present results point to an unconfounded protective effect of Native American ancestry on BC risk in both Colombian and Mexican women which appears to be stronger for familial and ER-positive BC. These findings provide a rationale for personalised prevention programmes that take genetic ancestry into account, as well as for future admixture mapping studies.

Genotyping, sequencing and analysis of 140,000 adults from Mexico City

The Mexico City Prospective Study is a prospective cohort of more than 150,000 adults recruited two decades ago from the urban districts of Coyoacán and Iztapalapa in Mexico City1. Here we generated genotype and exome-sequencing data for all individuals and whole-genome sequencing data for 9,950 selected individuals. We describe high levels of relatedness and substantial heterogeneity in ancestry composition across individuals. Most sequenced individuals had admixed Indigenous American, European and African ancestry, with extensive admixture from Indigenous populations in central, southern and southeastern Mexico. Indigenous Mexican segments of the genome had lower levels of coding variation but an excess of homozygous loss-of-function variants compared with segments of African and European origin. We estimated ancestry-specific allele frequencies at 142 million genomic variants, with an effective sample size of 91,856 for Indigenous Mexican ancestry at exome variants, all available through a public browser. Using whole-genome sequencing, we developed an imputation reference panel that outperforms existing panels at common variants in individuals with high proportions of central, southern and southeastern Indigenous Mexican ancestry. Our work illustrates the value of genetic studies in diverse populations and provides foundational imputation and allele frequency resources for future genetic studies in Mexico and in the United States, where the Hispanic/Latino population is predominantly of Mexican descent.

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.

MEX-PD: A National Network for the Epidemiological & Genetic Research of Parkinson's Disease

Background

Parkinson's Disease (PD) has a complex etiology, involving genetic and environmental factors. Most of our current understanding of the disease comes from studies in populations with mostly European ancestry, representing challenges in generalizing findings to other populations with different genetic, social, and environmental contexts. There are scarce studies focused in Latin American populations. The Mexican population is genetically diverse because its admixture from Native American, European, and African ancestries, coupled with the unique environmental conditions, stressing the relevance of establishing genetic studies in this population. Thus, we have established the Mexican Parkinson's Research Network (MEX-PD), a consortium to research the clinical, genetical, environmental, and neurophysiological bases of the phenotypic diversity in Mexican PD patients.

Objectives

Describing how MEX-PD was established, the methods and instruments and presenting the first results.

Methods

Patients and controls were recruited from medical centers in 20 states of Mexico. Initial recruitment included neurological evaluation, cognitive assessment, and DNA collection.

Results

MEX-PD has registered 302 controls and 262 PD patients with a mean age of diagnosis of 61 years (SD=10.86). There were 19.8% PD patients identified with early onset. Levodopa was the most common pharmacological treatment.

Conclusions

MEX-PD contributes to understand PD nationally. The information gathered here will allow us to understand the prevalence of mental health, neurological symptoms, and cognitive function in the PD Mexican population and how genetical and environmental factors contributes to those outcomes. These will advocate for personalized treatments and improving quality of life in the Mexican population.

Importance of Diversity in Precision Medicine: Generalizability of Genetic Associations Across Ancestry Groups Toward Better Identification of Disease Susceptibility Variants

Genome-wide association studies (GWAS) revolutionized our understanding of common genetic variation and its impact on common human disease and traits. Developed and adopted in the mid-2000s, GWAS led to searchable genotype-phenotype catalogs and genome-wide datasets available for further data mining and analysis for the eventual development of translational applications. The GWAS revolution was swift and specific, including almost exclusively populations of European descent, to the neglect of the majority of the world's genetic diversity. In this narrative review, we recount the GWAS landscape of the early years that established a genotype-phenotype catalog that is now universally understood to be inadequate for a complete understanding of complex human genetics. We then describe approaches taken to augment the genotype-phenotype catalog, including the study populations, collaborative consortia, and study design approaches aimed to generalize and then ultimately discover genome-wide associations in non-European descent populations. The collaborations and data resources established in the efforts to diversify genomic findings undoubtedly provide the foundations of the next chapters of genetic association studies with the advent of budget-friendly whole-genome sequencing.

Prevalence of Parent-Reported Food Allergy in a Mexican Pre-School Population

The magnitude and relevance of food allergies in the preschool population remain unknown in most regions of Mexico and Latin America. Thus, our aim was to estimate the parent-reported prevalence of food allergies in a Mexican preschool population and to characterize their clinical diagnosis and presentation. A cross-sectional survey was conducted in Culiacán City. A validated questionnaire was utilized. A total of 810 parents responded to the questionnaire (valid response rate, 40.7%). The estimated prevalence rates (95% CI) were: "physician-diagnosed Food Allergy (FA), ever" 5.30% (3.86-7.08); "immediate-type FA, ever" 2.96% (1.91-4.38); "immediate-type FA, current" 1.60% (0.86-2.73); and food-dependent anaphylaxis 1.11% (0.51-2.01). The main food allergens were milk (0.49%), strawberries (0.37%), egg, and soy (0.25% each). Atopic diseases and a family history of allergies were significantly associated with immediate-type FA. Among "immediate-type FA, current" cases, 76.9% required emergency room visits, but the prescription of epinephrine autoinjectors was reported in one case only. The food reactions occurred at home (92.35%), in a relative's house (38.5), and at restaurants (23%). Immediate-type FA reactions requiring emergency room visits are not uncommon among the studied population. Actions like proper anaphylaxis management and the prevention of cross-contamination of foods should be encouraged.

Admixture mapping of peripheral artery disease in a Dominican population reveals a putative risk locus on 2q35

Peripheral artery disease (PAD) is a form of atherosclerotic cardiovascular disease, affecting ∼8 million Americans, and is known to have racial and ethnic disparities. PAD has been reported to have a significantly higher prevalence in African Americans (AAs) compared to non-Hispanic European Americans (EAs). Hispanic/Latinos (HLs) have been reported to have lower or similar rates of PAD compared to EAs, despite having a paradoxically high burden of PAD risk factors; however, recent work suggests prevalence may differ between sub-groups. Here, we examined a large cohort of diverse adults in the BioMe biobank in New York City. We observed the prevalence of PAD at 1.7% in EAs vs. 8.5% and 9.4% in AAs and HLs, respectively, and among HL sub-groups, the prevalence was found at 11.4% and 11.5% in Puerto Rican and Dominican populations, respectively. Follow-up analysis that adjusted for common risk factors demonstrated that Dominicans had the highest increased risk for PAD relative to EAs [OR = 3.15 (95% CI 2.33-4.25), p < 6.44 × 10-14]. To investigate whether genetic factors may explain this increased risk, we performed admixture mapping by testing the association between local ancestry and PAD in Dominican BioMe participants (N = 1,813) separately from European, African, and Native American (NAT) continental ancestry tracts. The top association with PAD was an NAT ancestry tract at chromosome 2q35 [OR = 1.96 (SE = 0.16), p < 2.75 × 10-05) with 22.6% vs. 12.9% PAD prevalence in heterozygous NAT tract carriers versus non-carriers, respectively. Fine-mapping at this locus implicated tag SNP rs78529201 located within a long intergenic non-coding RNA (lincRNA) LINC00607, a gene expression regulator of key genes related to thrombosis and extracellular remodeling of endothelial cells, suggesting a putative link of the 2q35 locus to PAD etiology. Efforts to reproduce the signal in other Hispanic cohorts were unsuccessful. In summary, we showed how leveraging health system data helped understand nuances of PAD risk across HL sub-groups and admixture mapping approaches elucidated a putative risk locus in a Dominican population.

Accuracy of Eye and Hair Color Prediction in Mexican Mestizos from Monterrey City Based on ForenSeqTM DNA Signature Prep

Forensic genomic systems allow simultaneously analyzing identity informative (iiSNPs), ancestry informative (aiSNPs), and phenotype informative (piSNPs) genetic markers. Among these kits, the ForenSeq DNA Signature prep (Verogen) analyzes identity STRs and SNPs as well as 24 piSNPs from the HIrisPlex system to predict the hair and eye color. We report herein these 24 piSNPs in 88 samples from Monterrey City (Northeast, Mexico) based on the ForenSeq DNA Signature prep. Phenotypes were predicted by genotype results with both Universal Analysis Software (UAS) and the web tool of the Erasmus Medical Center (EMC). We observed predominantly brown eyes (96.5%) and black hair (75%) phenotypes, whereas blue eyes, and blond and red hair were not observed. Both UAS and EMC showed high performance in eye color prediction (p ≥ 96.6%), but a lower accuracy was observed for hair color prediction. Overall, UAS hair color predictions showed better performance and robustness than those obtained with the EMC web tool (when hair shade is excluded). Although we employed a threshold (p > 70%), we suggest using the EMC enhanced approach to avoid the exclusion of a high number of samples. Finally, although our results are helpful to employ these genomic tools to predict eye color, caution is suggested for hair color prediction in Latin American (admixed) populations such as those studied herein, principally when no black color is predicted.

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.

Guidelines for genetic ancestry inference created through roundtable discussions

The use of genetic and genomic technology to infer ancestry is commonplace in a variety of contexts, particularly in biomedical research and for direct-to-consumer genetic testing. In 2013 and 2015, two roundtables engaged a diverse group of stakeholders toward the development of guidelines for inferring genetic ancestry in academia and industry. This report shares the stakeholder groups' work and provides an analysis of, commentary on, and views from the groundbreaking and sustained dialogue. We describe the engagement processes and the stakeholder groups' resulting statements and proposed guidelines. The guidelines focus on five key areas: application of genetic ancestry inference, assumptions and confidence/laboratory and statistical methods, terminology and population identifiers, impact on individuals and groups, and communication or translation of genetic ancestry inferences. We delineate the terms and limitations of the guidelines and discuss their critical role in advancing the development and implementation of best practices for inferring genetic ancestry and reporting the results. These efforts should inform both governmental regulation and self-regulation.

Admixture Mapping of Peripheral Artery Disease in a Dominican Population Reveals a Novel Risk Locus on 2q35

Peripheral artery disease (PAD) is a form of atherosclerotic cardiovascular disease, affecting ∼8 million Americans, and is known to have racial and ethnic disparities. PAD has been reported to have significantly higher prevalence in African Americans (AAs) compared to non-Hispanic European Americans (EAs). Hispanic/Latinos (HLs) have been reported to have lower or similar rates of PAD compared to EAs, despite having a paradoxically high burden of PAD risk factors, however recent work suggests prevalence may differ between sub-groups. Here we examined a large cohort of diverse adults in the Bio Me biobank in New York City (NYC). We observed the prevalence of PAD at 1.7% in EAs vs 8.5% and 9.4% in AAs and HLs, respectively; and among HL sub-groups, at 11.4% and 11.5% in Puerto Rican and Dominican populations, respectively. Follow-up analysis that adjusted for common risk factors demonstrated that Dominicans had the highest increased risk for PAD relative to EAs (OR=3.15 (95% CI 2.33-4.25), P <6.44×10 -14 ). To investigate whether genetic factors may explain this increased risk, we performed admixture mapping by testing the association between local ancestry (LA) and PAD in Dominican Bio Me participants (N=1,940) separately for European (EUR), African (AFR) and Native American (NAT) continental ancestry tracts. We identified a NAT ancestry tract at chromosome 2q35 that was significantly associated with PAD (OR=2.05 (95% CI 1.51-2.78), P <4.06×10 -6 ) with 22.5% vs 12.5% PAD prevalence in heterozygous NAT tract carriers versus non-carriers, respectively. Fine-mapping at this locus implicated tag SNP rs78529201 located within a long intergenic non-coding RNA (lincRNA) LINC00607 , a gene expression regulator of key genes related to thrombosis and extracellular remodeling of endothelial cells, suggesting a putative link of the 2q35 locus to PAD etiology. In summary, we showed how leveraging health systems data helped understand nuances of PAD risk across HL sub-groups and admixture mapping approaches elucidated a novel risk locus in a Dominican population.

Single nucleotide variants in microRNA biosynthesis genes in Mexican individuals

Background: MicroRNAs (miRNAs) are important regulators in a variety of biological processes, and their dysregulation is associated with multiple human diseases. Single nucleotide variants (SNVs) in genes involved in the processing of microRNAs may alter miRNA regulation and could present high allele heterogeneity in populations from different ethnic groups. Thus, the aim of this study was to genotype 15 SNVs in eight genes involved in the miRNA processing pathway in Mexican individuals and compare their frequencies across 21 populations from five continental groups.

Methods: Genomic DNA was obtained from 399 healthy Mexican individuals. SNVs in AGO2 (rs2293939 and rs4961280), DGCR8 (rs720012), DICER (rs3742330 and rs13078), DROSHA (rs10719 and rs6877842), GEMIN3 (rs197388 and rs197414), GEMIN4 (rs7813, rs2740349, and rs4968104), TNRC6B (rs9611280), and XP05 (rs11077 and rs34324334) were genotyped using TaqMan probes. The minor allele frequency of each SNV was compared to those reported in the 1,000 Genomes database using chi-squared. Sankey plot was created in the SankeyMATIC package to visualize the frequency range of each variant in the different countries analyzed.

Results: In Mexican individuals, all 15 SNVs were found in Hardy-Weinberg equilibrium, with frequencies ranging from 0.04 to 0.45. The SNVs rs4961280, rs2740349, rs34324334, and rs720012 in Mexican individuals had the highest minor allele frequencies worldwide, whereas the minor allele frequencies of rs197388, rs10719, rs197414, and rs1107 were among the lowest in Mexican individuals. The variants had high allele heterogeneity among the sub-continental populations, ranging from monomorphic, as was the case for rs9611280 and rs34324334 in African groups, to &gt;0.50, which was the case for variants rs11077 and rs10719 in most of the populations. Importantly, the variants rs197388, rs720012, and rs197414 had FST values &gt; 0.18, indicating a directional selective process. Finally, the SNVs rs13078 and rs10719 significantly correlated with both latitude and longitude.

Conclusion: These data indicate the presence of high allelic heterogeneity in the worldwide distribution of the frequency of SNVs located in components of the miRNA processing pathway, which could modify the genetic susceptibility associated with human diseases in populations with different ancestry.

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.

Pathogenic variant profile in DNA damage response genes correlates with metastatic breast cancer progression-free survival in a Mexican-mestizo population

Introduction

Metastatic breast cancer causes the most breast cancer-related deaths around the world, especially in countries where breast cancer is detected late into its development. Genetic testing for cancer susceptibility started with the BRCA 1 and 2 genes. Still, recent research has shown that variations in other members of the DNA damage response (DDR) are also associated with elevated cancer risk, opening new opportunities for enhanced genetic testing strategies.

Methods

We sequenced BRCA1/2 and twelve other DDR genes from a Mexican-mestizo population of 40 metastatic breast cancer patients through semiconductor sequencing.

Results

Overall, we found 22 variants -9 of them reported for the first time- and a strikingly high proportion of variations in ARID1A. The presence of at least one variant in the ARID1A, BRCA1, BRCA2, or FANCA genes was associated with worse progression-free survival and overall survival in our patient cohort.

Discussion

Our results reflected the unique characteristics of the Mexican-mestizo population as the proportion of variants we found differed from that of other global populations. Based on these findings, we suggest routine screening for variants in ARID1A along with BRCA1/2 in breast cancer patients from the Mexican-mestizo population.

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.

Human leukocyte antigen Class II alleles associated with acral lentiginous melanoma in Mexican Mestizo patients: A case-control study

Background Melanoma is an aggressive cutaneous cancer. Acral lentiginous melanoma is a melanoma subtype arising on palms, soles, and nail-units. The incidence, prevalence and prognosis differ among populations. The link between expression of major histocompatibility complex Class II alleles and melanoma progression is known. However, available studies report variable results regarding the association of melanoma with specific HLA Class II loci. Aims The aim of the study was to determine HLA Class II allele frequencies in acral lentiginous melanoma patients and healthy Mexican Mestizo individuals. Methods Eighteen patients with acral lentiginous melanoma and 99 healthy controls were recruited. HLA Class II typing was performed based on the sequence-specific oligonucleotide method. Results Three alleles were associated with increased susceptibility to develop acral lentiginous melanoma, namely: HLA-DRB1*13:01; pC = 0.02, odds ratio = 6.1, IC95% = 1.4-25.5, HLA-DQA1*01:03; pC = 0.001, odds ratio = 9.3, IC95% = 2.7-31.3 and HLA-DQB1*02:02; pC = 0.01, odds ratio = 3.7, IC95% = 1.4-10.3. Limitations The small sample size was a major limitation, although it included all acral lentiginous melanoma patients seen at the dermatology department of Dr. Manuel Gea González General Hospital during the study period. Conclusion HLA-DRB1*13:01, HLA-DQB1*02:02 and HLA-DQA*01:03 alleles are associated with increased susceptibility to develop acral lentiginous melanoma in Mexican Mestizo patients.

Recommendations on the use and reporting of race, ethnicity, and ancestry in genetic research: Experiences from the NHLBI TOPMed program

How race, ethnicity, and ancestry are used in genomic research has wide-ranging implications for how research is translated into clinical care and incorporated into public understanding. Correlation between race and genetic ancestry contributes to unresolved complexity for the scientific community, as illustrated by heterogeneous definitions and applications of these variables. Here, we offer commentary and recommendations on the use of race, ethnicity, and ancestry across the arc of genetic research, including data harmonization, analysis, and reporting. While informed by our experiences as researchers affiliated with the NHLBI Trans-Omics for Precision Medicine (TOPMed) program, these recommendations are applicable to basic and translational genomic research in diverse populations with genome-wide data. Moving forward, considerable collaborative effort will be required to ensure that race, ethnicity, and ancestry are described and used appropriately to generate scientific knowledge that yields broad and equitable benefit.

The immunogenetic impact of European colonization in the Americas

The introduction of pathogens originating from Eurasia into the Americas during early European contact has been associated with high mortality rates among Indigenous peoples, likely contributing to their historical and precipitous population decline. However, the biological impacts of imported infectious diseases and resulting epidemics, especially in terms of pathogenic effects on the Indigenous immunity, remain poorly understood and highly contentious to this day. Here, we examine multidisciplinary evidence underpinning colonization-related immune genetic change, providing contextualization from anthropological studies, paleomicrobiological evidence of contrasting host-pathogen coevolutionary histories, and the timings of disease emergence. We further summarize current studies examining genetic signals reflecting post-contact Indigenous population bottlenecks, admixture with European and other populations, and the putative effects of natural selection, with a focus on ancient DNA studies and immunity-related findings. Considering current genetic evidence, together with a population genetics theoretical approach, we show that post-contact Indigenous immune adaptation, possibly influenced by selection exerted by introduced pathogens, is highly complex and likely to be affected by multifactorial causes. Disentangling putative adaptive signals from those of genetic drift thus remains a significant challenge, highlighting the need for the implementation of population genetic approaches that model the short time spans and complex demographic histories under consideration. This review adds to current understandings of post-contact immunity evolution in Indigenous peoples of America, with important implications for bettering our understanding of human adaptation in the face of emerging infectious diseases.

Genetic Diversity of Mycobacterium tuberculosis Isolates From an Amerindian Population in Chiapas, México

This is the first report of the genetic diversity of the Mycobacterium tuberculosis complex isolates found in a Mexican-Amerindian setting. In this study, we analyzed isolates collected from the Highlands region of Chiapas, Mexico, by using spoligotyping and whole-genome sequencing analyses. Seventy-three M. tuberculosis isolates were analyzed initially by spoligotyping; no new spoligotypes were identified. Nineteen percent of the isolates were identified as SIT53 (T1) (n = 14), followed by SIT42 (14%, n = 10, LAM9) and SIT119 (11%; n = 8, X1). SIT53, SIT42, and orphan isolates (16.4%, n = 12) constituted about 50% of the isolates studied and were subjected to whole-genome sequencing (WGS) analysis. Most SIT53 (10/12) isolates belonged to the Euro-American sub-lineage 4.8. Most SIT42 isolates (4/7) as .well as most orphan isolates (5/8) belonged to the lineage 4.3.3 LAM group. By comparing the single-nucleotide polymorphism (SNP) patterns of the SIT53 isolates, we found one clone (&lt;7 SNPs) and four clustered isolates (&lt;15 SNPs). In isolates from the SIT42 and orphan groups, we did not find any clones or clusters. This work demonstrates the success of sub-lineage 4.8 to predominate in Mexico and confirms the dominion of sub-lineage 4.3.3 in Central and South America.

Variability in Susceptibility to Type I Interferon Response and Subgenomic RNA Accumulation Between Clinical Isolates of Dengue and Zika Virus From Oaxaca Mexico Correlate With Replication Efficiency in Human Cells and Disease Severity

Dengue and Zika viruses cocirculate annually in endemic areas of Mexico, causing outbreaks of different magnitude and severity every year, suggesting a continuous selection of Flavivirus variants with variable phenotypes of transmissibility and virulence. To evaluate if Flavivirus variants with different phenotypes cocirculate during outbreaks, we isolated dengue and Zika viruses from blood samples of febrile patients from Oaxaca City during the 2016 and 2019 epidemic years. We compared their replication kinetics in human cells, susceptibility to type I interferon antiviral response, and the accumulation of subgenomic RNA on infected cells. We observed correlations between type I interferon susceptibility and subgenomic RNA accumulation, with high hematocrit percentage and thrombocytopenia. Our results suggest that Flaviviruses that cocirculate in Oaxaca, Mexico, have variable sensitivity to the antiviral activity of type I interferons, and this phenotypic trait correlates with the severity of the disease.

Hispanic ethnicity and the rs4880 variant in SOD2 are associated with elevated liver enzymes and bilirubin levels in children receiving asparaginase-containing chemotherapy for acute lymphoblastic leukemia

Asparaginase is an integral component of acute lymphoblastic leukemia (ALL)³ treatment. Hepatotoxicity related to asparaginase is one of the most common treatment-related toxicities in ALL therapy. Hispanic children are at higher risk of developing ALL, and toxicities from ALL therapy. The rs4880 variant in the superoxide dismutase 2 (SOD2)⁴ gene, a critical mitochondrial enzyme that protects cells against oxidative stress, was found to be associated with increased incidence of asparaginase-related hepatotoxicity in adult cohort of largely White non-Hispanics patients with ALL. The risk genotype (rs4880-CC) is more frequent among adult Hispanic patients with ALL. To assess the prevalence of hepatotoxicity and risk genotype among pediatric patients with ALL, particularly of Hispanic ethnicity, we conducted a prospective study of 143 pediatric patients with ALL (62.2% Hispanic). Bilirubin and hepatic transaminase levels were collected at different times during multiagent therapy including asparaginase treatment. Germline DNA blood samples were genotyped for the SOD2 rs4880. We found that the frequency of hepatotoxicity and the rs4880-CC risk genotype are higher in Hispanic patients than non-Hispanic. Patients with the CC genotype exhibit higher bilirubin and hepatic transaminase levels compared with patients with the TT and CT genotypes. In a multivariate Cox analysis, Hispanic ethnicity was identified as a strong predictor of hepatotoxicity (hazard ratio [HR] = 1.9, 95% confidence interval [95% CI] 1.0-3.5, p = 0.05). Altogether, these findings demonstrate that hepatotoxicity is highly prevalent among Hispanic pediatric patients with ALL, and those with rs4880-CC genotype.

Exome Sequencing Data Analysis and a Case-Control Study in Mexican Population Reveals Lipid Trait Associations of New and Known Genetic Variants in Dyslipidemia-Associated Loci

Background: Plasma lipid levels are a major risk factor for cardiovascular diseases. Although international efforts have identified a group of loci associated with the risk of dyslipidemia, Latin American populations have been underrepresented in these studies.

Objective: To know the genetic variation occurring in lipid-related loci in the Mexican population and its association with dyslipidemia.

Methods: We searched for single-nucleotide variants in 177 lipid candidate genes using previously published exome sequencing data from 2838 Mexican individuals belonging to three different cohorts. With the extracted variants, we performed a case-control study. Logistic regression and quantitative trait analyses were implemented in PLINK software. We used an LD pruning using a 50-kb sliding window size, a 5-kb window step size and a r2 threshold of 0.1.

Results: Among the 34251 biallelic variants identified in our sample population, 33% showed low frequency. For case-control study, we selected 2521 variants based on a minor allele frequency ≥1% in all datasets. We found 19 variants in 9 genes significantly associated with at least one lipid trait, with the most significant associations found in the APOA1/C3/A4/A5-ZPR1-BUD13 gene cluster on chromosome 11. Notably, all 11 variants associated with hypertriglyceridemia were within this cluster; whereas variants associated with hypercholesterolemia were located at chromosome 2 and 19, and for low high density lipoprotein cholesterol were in chromosomes 9, 11, and 19. No significant associated variants were found for low density lipoprotein. We found several novel variants associated with different lipemic traits: rs3825041 in BUD13 with hypertriglyceridemia, rs7252453 in CILP2 with decreased risk to hypercholesterolemia and rs11076176 in CETP with increased risk to low high density lipoprotein cholesterol.

Conclusions: We identified novel variants in lipid-regulation candidate genes in the Mexican population, an underrepresented population in genomic studies, demonstrating the necessity of more genomic studies on multi-ethnic populations to gain a deeper understanding of the genetic structure of the lipemic traits.

The Tibetan-Yi region is both a corridor and a barrier for human gene flow

The Tibetan-Yi Corridor (TYC) region between Tibet and the rest of east Asia has served as a crossroads for human migrations for thousands of years. The lack of whole-genome sequencing data specific to the TYC populations has hindered the understanding of the fundamental patterns of migration and divergence between humans in east Asia and southeast Asia. Here, we provide 248 individual whole genomes from the 16 TYC and 3 outgroup populations to elucidate historical relationships. We find that the Tibetan plateau forms an important barrier to gene flow, with a more Tibetan-like ancestry in northern populations and a southern east Asian-related ancestry in south populations. An isolated population, Achang, shows a prolonged isolation and genetic drift compared to other TYC populations. We also note that previous claims regarding the history and structure of TYC populations inferred by linguistics are incompatible with the genetic evidence.

Global expression profiling of CD10 + /CD19 + pre-B lymphoblasts from Hispanic B-ALL patients correlates with comparative TARGET database analysis

Mexico City has one of the highest incidences of acute lymphoblastic leukemia (ALL) globally, with patients showing low survival, and high relapse rates. To gain more insight into the molecular features of B-ALL in Mexican children, we isolated CD10 + /CD19 + precursor B lymphoblasts from four bone marrow and nine peripheral blood samples of B-ALL patients using a fluorescence-activated cell sorting protocol. The global gene expression profile (BM vs PB) revealed 136 differentially expressed genes; 62 were upregulated (45.6%) and 74 were downregulated (54.4%). Pearson's correlation coefficient was calculated to determine the similarity between pre-B lymphoblast populations. We selected 26 highly significant genes and validated 21 by RT-qPCR (CNN3, STON2, CALN1, RUNX2, GADD45A, CDC45, CDC20, PLK1, AIDA, HCK, LY86, GPR65, PIK3CG, LILRB2, IL7R, TCL1A, DOCK1, HIST1H3G, PTPN14, CD72, and NT5E). The gene set enrichment analysis of the total expression matrix and the ingenuity pathway analysis of the 136 differentially expressed genes showed that the cell cycle was altered in the bone marrow with four overexpressed genes (PLK1, CDC20, CDC45, and GADD45A) and a low expression of IL7R and PIK3CG, which are involved in B cell differentiation. A comparative bioinformatics analysis of 15 bone marrow and 10 peripheral blood samples from Hispanic B-ALL patients collected by the TARGET program, corroborated the genes observed, except for PIK3CG. We conclude the Mexican and the Hispanic B-ALL patients studied present common driver alterations and histotype-specific mutations that could facilitate risk stratification and diagnostic accuracy and serve as potential therapeutic targets.

The genetic legacy of the Manila galleon trade in Mexico

The population of Mexico has a considerable genetic substructure due to both its pre-Columbian diversity and due to genetic admixture from post-Columbian trans-oceanic migrations. The latter primarily originated in Europe and Africa, but also, to a lesser extent, in Asia. We analyze previously understudied genetic connections between Asia and Mexico to infer the timing and source of this genetic ancestry in Mexico. We identify the predominant origin within Southeast Asia—specifically western Indonesian and non-Negrito Filipino sources—and we date its arrival in Mexico to approximately 13 generations ago (1620 CE). This points to a genetic legacy from the seventeenth century Manila galleon trade between the colonial Spanish Philippines and the Pacific port of Acapulco. Indeed, within Mexico we observe the highest level of this trans-Pacific ancestry in Acapulco, located in the state of Guerrero. This colonial Spanish trade route from East Asia to Europe was centred on Mexico and appears in historical records, but its legacy has been largely ignored. Identities and stories were suppressed due to slavery, assimilation of the immigrants as ‘Indios’ and incomplete historical records. Here we characterize this understudied Mexican ancestry.

This article is part of the theme issue ‘Celebrating 50 years since Lewontin's apportionment of human diversity’.

Human genetic admixture through the lens of population genomics

Over the past 50 years, geneticists have made great strides in understanding how our species' evolutionary history gave rise to current patterns of human genetic diversity classically summarized by Lewontin in his 1972 paper, ‘The Apportionment of Human Diversity’. One evolutionary process that requires special attention in both population genetics and statistical genetics is admixture: gene flow between two or more previously separated source populations to form a new admixed population. The admixture process introduces ancestry-based structure into patterns of genetic variation within and between populations, which in turn influences the inference of demographic histories, identification of genetic targets of selection and prediction of complex traits. In this review, we outline some challenges for admixture population genetics, including limitations of applying methods designed for populations without recent admixture to the study of admixed populations. We highlight recent studies and methodological advances that aim to overcome such challenges, leveraging genomic signatures of admixture that occurred in the past tens of generations to gain insights into human history, natural selection and complex trait architecture.

This article is part of the theme issue ‘Celebrating 50 years since Lewontin's apportionment of human diversity’.

Ancestral diversity improves discovery and fine-mapping of genetic loci for anthropometric traits-The Hispanic/Latino Anthropometry Consortium

Hispanic/Latinos have been underrepresented in genome-wide association studies (GWAS) for anthropometric traits despite their notable anthropometric variability, ancestry proportions, and high burden of growth stunting and overweight/obesity. To address this knowledge gap, we analyzed densely imputed genetic data in a sample of Hispanic/Latino adults to identify and fine-map genetic variants associated with body mass index (BMI), height, and BMI-adjusted waist-to-hip ratio (WHRadjBMI). We conducted a GWAS of 18 studies/consortia as part of the Hispanic/Latino Anthropometry (HISLA) Consortium (stage 1, n = 59,771) and generalized our findings in 9 additional studies (stage 2, n = 10,538). We conducted a trans-ancestral GWAS with summary statistics from HISLA stage 1 and existing consortia of European and African ancestries. In our HISLA stage 1 + 2 analyses, we discovered one BMI locus, as well as two BMI signals and another height signal each within established anthropometric loci. In our trans-ancestral meta-analysis, we discovered three BMI loci, one height locus, and one WHRadjBMI locus. We also identified 3 secondary signals for BMI, 28 for height, and 2 for WHRadjBMI in established loci. We show that 336 known BMI, 1,177 known height, and 143 known WHRadjBMI (combined) SNPs demonstrated suggestive transferability (nominal significance and effect estimate directional consistency) in Hispanic/Latino adults. Of these, 36 BMI, 124 height, and 11 WHRadjBMI SNPs were significant after trait-specific Bonferroni correction. Trans-ancestral meta-analysis of the three ancestries showed a small-to-moderate impact of uncorrected population stratification on the resulting effect size estimates. Our findings demonstrate that future studies may also benefit from leveraging diverse ancestries and differences in linkage disequilibrium patterns to discover novel loci and additional signals with less residual population stratification.

Disentangling Signatures of Selection Before and After European Colonization in Latin Americans

Throughout human evolutionary history, large-scale migrations have led to intermixing (i.e., admixture) between previously separated human groups. Although classical and recent work have shown that studying admixture can yield novel historical insights, the extent to which this process contributed to adaptation remains underexplored. Here, we introduce a novel statistical model, specific to admixed populations, that identifies loci under selection while determining whether the selection likely occurred post-admixture or prior to admixture in one of the ancestral source populations. Through extensive simulations, we show that this method is able to detect selection, even in recently formed admixed populations, and to accurately differentiate between selection occurring in the ancestral or admixed population. We apply this method to genome-wide SNP data of ∼4,000 individuals in five admixed Latin American cohorts from Brazil, Chile, Colombia, Mexico, and Peru. Our approach replicates previous reports of selection in the human leukocyte antigen region that are consistent with selection post-admixture. We also report novel signals of selection in genomic regions spanning 47 genes, reinforcing many of these signals with an alternative, commonly used local-ancestry-inference approach. These signals include several genes involved in immunity, which may reflect responses to endemic pathogens of the Americas and to the challenge of infectious disease brought by European contact. In addition, some of the strongest signals inferred to be under selection in the Native American ancestral groups of modern Latin Americans overlap with genes implicated in energy metabolism phenotypes, plausibly reflecting adaptations to novel dietary sources available in the Americas.

Development of birthweight and length for gestational age and sex references in Yucatan, Mexico

OBJECTIVE

To develop sex- and gestational age specific reference percentiles and curves for birth weight and length for Yucatec neonates using data from birth registers of infants born during 2015-2019.

MATERIAL AND METHODS

Observational, descriptive, epidemiologic study in a 5-year period including every registered birth in the state of Yucatan, Mexico using birth registries. A total of 158 432 live, physically healthy singletons (76 442 females and 81 990 males) between 25 and 42 weeks of gestation were included in the analysis. We used the LMS method to construct smoothed reference centiles (3rd, 10th, 25th, 50th, 75th, 95th, and 97th) and curves for males and females separately.

RESULTS

Mean maternal age was 26 (SD = 6.22) years. Fifty-two percent of births occurred by vaginal delivery, 37% were firstborn and similar proportions were second (33%) and third or more (30%) born. 5.5% of newborns included in the references corresponds to neonates born before 37 weeks of gestation (5.9% boys and 5.1% girls). In both sexes, the percentage of infants with a birthweight less than 2500 g was 6.7%. The birthweight at the 50th percentile for males and females at 40 weeks of gestation in this cohort was 3256 and 3167 g, respectively, and the corresponding values for birth length were 50.23 and 49.84 cm (mean differences between sexes: 89 g and 0.40 cm, respectively).

CONCLUSION

The reference percentile and curves developed in this study are useful for research purposes and can help health practitioners to assess the biological status of infants born in Yucatán.

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.

A Bibliometric Analysis of Mexican Bioinformatics: A Portrait of Actors, Structure, and Dynamics

Bioinformatics is a very important informatics tool for health and biological sciences, focusing on biological data management. The objective of this work was to perform a bibliometric analysis regarding the development of Mexican bioinformatics. An exhaustive revision of the literature associated with Mexican bioinformatics in a period of 25-years was performed. Bibliometric tools, such as performance analysis and science mapping were included in the analysis. We identified the main actors as well as the structure and dynamics of Mexican bioinformatics. Some of the main findings were as follows: the thematic structure in the field is defined by the research lines of outstanding authors; the outstanding collaborations of Mexican institutions with foreign countries and institutions are influenced by the geographic proximity and binational agreements, as well as philanthropic and academic programs that promote collaborations, and there is an inclination for health issues promoted by public health financing and philanthropic organizations. It is identified that publications had an explosion since 2012, we consider that this growth may be influenced by the democratization of data, derived from the mass sequencing of biological molecules stored in public databases.

Is there still evolution in the human population?

It is often claimed that humanity has stopped evolving because modern medicine erased all selection on survival. Even if that would be true, and it is not, there would be other mechanisms of evolution which could still led to changes in allelic frequencies. Here I show, by applying basic evolutionary genetics knowledge, that we expect humanity to evolve. The results from genome sequencing projects have repeatedly affirmed that there are still recent signs of selection in our genomes. I give some examples of such adaptation. Then I briefly discuss what our evolutionary future has in store for us.

Dyslipidemia, Obesity, and Ethnicity in Mexican Children

The aim of this study was to assess lipid disorders in children from five ethnic groups, both urban and indigenous, from northern and central Mexico. We measured the lipid profile to determine the ability of the body mass index (BMI) to discriminate an abnormally high lipid level using receiving operating characteristics (ROC). We analyzed the association and interaction of obesity and ethnicity with lipid disorders using generalized linear models in 977 children. The highest prevalence of lipid disorders (high TG, high TC, high LDL, high APOB, and dyslipidemia) was found in central Mexico-Mexico City and urban northern Mexico. The BMI performed better at predicting low HDL in Seris, a northern indigenous group (0.95, CI: 0.69–0.85), and Mexico City (0.75, CI: 0.69–0.82), and high LDL in Puebla (central Mexico, 0.80, CI: 0.69–0.85). Obesity significantly (p < 0.05) increases lipid disorders by around two times (OR~2) for almost all lipid markers. Obesity and ethnic interaction increase the lipid disorders by more than five times for different lipid markers and ethnic groups (high total cholesterol OR = 5.31; low HDL OR = 5.11, and dyslipidemia OR = 5.68). Lipid disorders are not restricted to children with high BMIs, but obesity exacerbates these. The emerging lipid disorder risk depends on the ethnic group.

Neuropsychiatric Genetics of Psychosis in the Mexican Population: A Genome-Wide Association Study Protocol for Schizophrenia, Schizoaffective, and Bipolar Disorder Patients and Controls

No large-scale genome-wide association studies (GWASs) of psychosis have been conducted in Mexico or Latin America to date. Schizophrenia and bipolar disorder in particular have been found to be highly heritable and genetically influenced. However, understanding of the biological basis of psychosis in Latin American populations is limited as previous genomic studies have almost exclusively relied on participants of Northern European ancestry. With the goal of expanding knowledge on the genomic basis of psychotic disorders within the Mexican population, the National Institute of Psychiatry Ramón de la Fuente Muñiz (INPRFM), the Harvard T.H. Chan School of Public Health, and the Broad Institute's Stanley Center for Psychiatric Research launched the Neuropsychiatric Genetics Research of Psychosis in Mexican Populations (NeuroMex) project to collect and analyze case-control psychosis samples from 5 states across Mexico. This article describes the planned sample collection and GWAS protocol for the NeuroMex study. The 4-year study will span from April 2018 to 2022 and aims to recruit 9,208 participants: 4,604 cases and 4,604 controls. Study sites across Mexico were selected to ensure collected samples capture the genomic diversity within the Mexican population. Blood samples and phenotypic data will be collected during the participant interview process and will contribute to the development of a local biobank in Mexico. DNA extraction will be done locally and genetic analysis will take place at the Broad Institute in Cambridge, MA. We will collect extensive phenotypic information using several clinical scales. All study materials including phenotypic instruments utilized are openly available in Spanish and English. The described study represents a long-term collaboration of a number of institutions from across Mexico and the Boston area, including clinical psychiatrists, clinical researchers, computational biologists, and managers at the 3 collaborating institutions. The development of relevant data management, quality assurance, and analysis plans are the primary considerations in this protocol article. Extensive management and analysis processes were developed for both the phenotypic and genetic data collected. Capacity building, partnerships, and training between and among the collaborating institutions are intrinsic components to this study and its long-term success.

Prevalence of protective haplotypes of the SLCO1B1 gene for statin transport in Mexican populations

Aim: To evaluate the genetic distribution of the rs4149056 and rs2306283 variants in the SLCO1B1 gene in Mexican Mestizo (admixed) and Native American groups. Materials & methods: We recruited 360 volunteers who were qPCR-genotyped with TaqMan probes. Results: Allele and genotype frequencies are reported. Among the expected rs4149056-rs2306283 haplotypes, T-A (42.35-58.47%) was the most prevalent which relates to the normal activity of the OATP1B1 transporter. This was followed by the T-G haplotype associated with further statin transport and cholesterol reduction (32.49-43.76%). Conclusion: Based on these SLCO1B1 gene variants, we confirmed that a minimum fraction of the Mexican study populations would be at risk from decreasing simvastatin transport and the development of statin-induced myopathy.

The population genetics characteristics of a 90 locus panel of microhaplotypes

Single-nucleotide polymorphisms (SNPs) and small genomic regions with multiple SNPs (microhaplotypes, MHs) are rapidly emerging as novel forensic investigative tools to assist in individual identification, kinship analyses, ancestry inference, and deconvolution of DNA mixtures. Here, we analyzed information for 90 microhaplotype loci in 4009 individuals from 79 world populations in 6 major biogeographic regions. The study included multiplex microhaplotype sequencing (mMHseq) data analyzed for 524 individuals from 16 populations and genotype data for 3485 individuals from 63 populations curated from public repositories. Analyses of the 79 populations revealed excellent characteristics for this 90-plex MH panel for various forensic applications achieving an overall average effective number of allele values (A_e) of 4.55 (range 1.04–19.27) for individualization and mixture deconvolution. Population-specific random match probabilities ranged from a low of 10^–115 to a maximum of 10^–66. Mean informativeness (I_n) for ancestry inference was 0.355 (range 0.117–0.883). 65 novel SNPs were detected in 39 of the MHs using mMHseq. Of the 3018 different microhaplotype alleles identified, 1337 occurred at frequencies > 5% in at least one of the populations studied. The 90-plex MH panel enables effective differentiation of population groupings for major biogeographic regions as well as delineation of distinct subgroupings within regions. Open-source, web-based software is available to support validation of this technology for forensic case work analysis and to tailor MH analysis for specific geographical regions.

Genetic ancestry and ethnic identity in Ecuador

We investigated the ancestral origins of four Ecuadorian ethnic groups-Afro-Ecuadorian, Mestizo, Montubio, and the Indigenous Tsáchila-in an effort to gain insight on the relationship between ancestry, culture, and the formation of ethnic identities in Latin America. The observed patterns of genetic ancestry are largely concordant with ethnic identities and historical records of conquest and colonization in Ecuador. Nevertheless, a number of exceptional findings highlight the complex relationship between genetic ancestry and ethnicity in Ecuador. Afro-Ecuadorians show far less African ancestry, and the highest levels of Native American ancestry, seen for any Afro-descendant population in the Americas. Mestizos in Ecuador show high levels of Native American ancestry, with substantially less European ancestry, despite the relatively low Indigenous population in the country. The recently recognized Montubio ethnic group is highly admixed, with substantial contributions from all three continental ancestries. The Tsáchila show two distinct ancestry subgroups, with most individuals showing almost exclusively Native American ancestry and a smaller group showing a Mestizo characteristic pattern. Considered together with historical data and sociological studies, our results indicate the extent to which ancestry and culture interact, often in unexpected ways, to shape ethnic identity in Ecuador.

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.

Relevance of NR1I2 variants on carbamazepine therapy in Mexican Mestizos with epilepsy at a tertiary-care hospital

Aim: We evaluated the potential influence of genetic (CYP3A5, EPHX1, NR1I2, HNF4A, ABCC2, RALBP1, SCN1A, SCN2A and GABRA1) and nongenetic factors on carbamazepine (CBZ) response, adverse drug reactions and CBZ plasma concentrations in 126 Mexican Mestizos (MM) with epilepsy. Subjects & methods: Patients were genotyped for 27 variants using TaqMan^® assays. Results: CBZ response was associated with NR1I2 variants and lamotrigine cotreatment. CBZ-induced adverse drug reactions were related to antiepileptic polytherapy and SCN1A rs2298771/rs3812718 haplotype. CBZ plasma concentrations were influenced by NR1I2-rs2276707 and -rs3814058, and by phenytoin cotreatment. CBZ daily dose was also influenced by NR1I2-rs3814055 and EPHX1-rs1051740. Conclusion: Interindividual variability in CBZ treatment was partly explained by NR1I2, EPHX1 and SCN1A variants, as well as antiepileptic cotreatment in MM with epilepsy.

Association between genetic variations at 8q24 and prostate cancer risk in Mexican Men

BACKGROUND

Variants of 8q24 locus have been associated with prostate cancer (PCa) susceptibility. This study aims to analyze the genetic basis of PCa susceptibility in Mexican men by analyzing SNPs in the 8q24 locus for the first time.

METHODS

A case-control study was performed in 875 men recruited from the Mexican Social Security Institute, 326 patients with PCa, and 549 non-PCa patients (88 with benign prostatic hyperplasia BPH and 461 healthy controls). The 8q24 locus SNPs: rs16901979, rs16983267, rs1447295, and rs7837328 were genotyped by allelic discrimination assays using TaqMan probes. Statistical analysis was performed using Epi Info statistical 7.0 and SNPstats softwares.

RESULTS

All genotype frequencies were in Hardy-Weinberg Equilibrium. No differences were observed in genotype distribution between PCa and non-PCa patients for rs6983267. Under different inheritance models, the rs16901979, rs1447295, and rs7837328 SNPs were associated with PCa (OR = 2.8, 1.8, and 1.72, respectively; Pc < 0.001) when comparing PCa patients against controls. This association remains between PCa and BPH patients under different models (OR = 8.5, 2.2, and 1.9, respectively; Pc < 0.001). There were no significant differences in allele and genotype distribution among BPH patients and controls. The combined effect of the alleles CGAA for the SNPs rs16901979, rs6983267, rs1447295, and rs7837328 showed significant differences between PCa patients and controls (OR = 2.9, 95% CI = 1.48-5.83, Pc = 0.008). Four 8q24 variants were not associated with D'Amico score, age at diagnosis, and bone metastases.

CONCLUSIONS

Our study provides the first confirmation that variants rs16901979, rs1447295, and 7837328 at 8q24 locus are associated with PCa susceptibility in Mexican men.