Genetic diversity of HLA system in six populations from Mexico City Metropolitan Area, Mexico: Mexico City North, Mexico City South, Mexico City East, Mexico City West, Mexico City Center and rural Mexico City

Heterogeneity of Genetic Admixture Determines SLE Susceptibility in Mexican

Systemic Lupus Erythematosus (SLE) is an autoimmune inflammatory disorder for which Major Histocompatibility Complex (MHC) genes are well identified as risk factors. SLE patients present different clinical phenotypes, which are partly explained by admixture patterns variation among Mexicans. Population genetic has insight into the high genetic variability of Mexicans, mainly described through HLA gene studies with anthropological and biomedical importance. A prospective, case-control study was performed. In this study, we recruited 146 SLE patients, and 234 healthy individuals were included as a control group; both groups were admixed Mexicans from Mexico City. The HLA typing methods were based on Next Generation Sequencing and Sequence-Based Typing (SBT). The data analysis was performed with population genetic programs and statistical packages. The admixture estimations based on HLA-B and -DRB1 revealed that SLE patients have a higher Southwestern European ancestry proportion (48 ± 8%) than healthy individuals (30 ± 7%). In contrast, Mexican Native American components are diminished in SLE patients (44 ± 1%) and augmented in Healthy individuals (63 ± 4%). HLA alleles and haplotypes' frequency analysis found variants previously described in SLE patients from Mexico City. Moreover, a conserved extended haplotype that confers risk to develop SLE was found, the HLA-A^∗29:02∼C^∗16:01∼B^∗44:03∼DRB1^∗07:01∼DQB1^∗02:02, pC = 0.02, OR = 1.41. Consistent with the admixture estimations, the origin of all risk alleles and haplotypes found in this study are European, while the protection alleles are Mexican Native American. The analysis of genetic distances supported that the SLE patient group is closer to the Southwestern European parental populace and farthest from Mexican Native Americans than healthy individuals. Heterogeneity of genetic admixture determines SLE susceptibility and protection in Mexicans. HLA sequencing is helpful to determine susceptibility alleles and haplotypes restricted to some populations.

A bioinformatic prediction of antigen presentation from SARS-CoV-2 spike protein revealed a theoretical correlation of HLA-DRB1*01 with COVID-19 fatality in Mexican population: An ecological approach

SARS‐CoV‐2 infection is causing a pandemic disease that is reflected in challenging public health problems worldwide. Human leukocyte antigen (HLA)‐based epitope prediction and its association with disease outcomes provide an important base for treatment design. A bioinformatic prediction of T cell epitopes and their restricted HLA Class I and II alleles was performed to obtain immunogenic epitopes and HLA alleles from the spike protein of the severe acute respiratory syndrome coronavirus 2 virus. Also, a correlation with the predicted fatality rate of hospitalized patients in 28 states of Mexico was done. Here, we describe a set of 10 highly immunogenic epitopes, together with different HLA alleles that can efficiently present these epitopes to T cells. Most of these epitopes are located within the S1 subunit of the spike protein, suggesting that this area is highly immunogenic. A statistical negative correlation was found between the frequency of HLA‐DRB1*01 and the fatality rate in hospitalized patients in Mexico.

First HLA association study for COVID‐19 in Mexico

An epitope prediction for HLA Class I and II provided a list of highly immunogenic epitopes from the S protein of SARS‐CoV2 with potential use for vaccine development.

A multi‐level approach revealed a correlation of HLA‐DRB1*01 frequency with fatality in Mexican hospitalized patients at ecological level.

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.

Genome-Wide Association Study of Body Mass Index and Body Fat in Mexican-Mestizo Children

BACKGROUND

Childhood obesity is a major health problem in Mexico. Obesity prevalence estimated by body mass index (BMI) is almost half than that estimated by percent body fat (%BF) in the Childhood Obesity pediatric cohort (COIPIS).

OBJECTIVE

We performed a genome-wide association study (GWAS) of BMI and %BF in 828 children from the COIPIS to identify markers of predisposition to high values for both phenotypes used for obesity classification.

METHODS

For the GWAS we used the LAT Axiom 1, Affymetrix and 2.5 million single loci from the 1000 Genomes Phase 3 imputation panel. We used a linear model, adjusted by age, sex, and Amerindian ancestry assuming an additive inheritance model.

RESULTS

Genome-wide significance (p ≤ 5.0 × 10-8) and 80% of statistical power was reached for associations of two loci in two genes (CERS3 and CYP2E1) to BMI. Also, 11 loci in six genes (ANKS1B, ARNTL2, KCNS3, LMNB1, SRGAP3, TRPC7) reached genome-wide significance for associations to %BF, though not 80% of statistical power.

DISCUSSION

None of the SNPs were previously reported as being associated to BMI or %BF. In addition, different loci were found for BMI and %BF. These results highlight the importance of gaining deeper understanding of genetic markers of predisposition to high values for the phenotypes used for obesity diagnosis.