CFTR pathogenic variants spectrum in a cohort of Mexican patients with cystic fibrosis.

Background: Molecular diagnosis of cystic fibrosis (CF) is challenging in Mexico due to the population's high genetic heterogeneity. To date, 46 pathogenic variants (PVs) have been reported, yielding a detection rate of 77%. We updated the spectrum and frequency of PVs responsible for this disease in mexican patients. Methods: We extracted genomic DNA from peripheral blood lymphocytes obtained from 297 CF patients and their parents. First, we analyzed the five most frequent PVs in the Mexican population using PCR-mediated site-directed mutagenesis. In patients with at least one identified allele, CFTR sequencing was performed using next-generation sequencing tools and multiplex ligation-dependent probe amplification. For variants not previously classified as pathogenic, we used a combination of in silico prediction, CFTR modeling, and clinical characteristics to determine a genotype-phenotype correlation. Results: We identified 95 PVs, increasing the detection rate to 87.04%. The most frequent variants were p.(PheF508del) (42.7%), followed by p.(Gly542*) (5.6%), p.(Ser945Leu) (2.9%), p.(Trp1204*) and p.(Ser549Asn) (2.5%), and CFTRdel25-26 and p.(Asn386Ilefs*3) (2.3%). The remaining variants had frequencies of <2.0%, and some were exclusive to one family. We identified 10 novel PVs localized in different exons (frequency range: 0.1-0.8%), all of which produced structural changes, deletions, or duplications in different domains of the protein, resulting in dysfunctional ion flow. The use of different in silico software and American College of Medical Genetics and Genomics (ACMG) and the Association for Molecular Pathology (AMP) criteria allowed us to assume that all of these PVs were pathogenic, causing a severe phenotype. Conclusions: In a highly heterogeneous population, combinations of different tools are needed to identify the variants responsible for CF and enable the establishment of appropriate strategies for CF diagnosis, prevention, and treatment.

Total Antioxidant Capacity in Obese and Non-Obese Subjects and Its Association with Anthropo-Metabolic Markers: Systematic Review and Meta-Analysis.

The total antioxidant capacity (TAC) has been related to the development of and complications associated with chronic diseases, but its importance during obesity is not entirely clear. We conducted a systematic review and meta-analysis to clarify whether there are differences or similarities in the TAC between subjects with obesity (SO) and subjects with normal weight (NW). Following the recommendations of PRISMA and Cochrane, we performed a systematic search in the PubMed, Scopus, Web of Science, Cochrane, and PROSPERO databases, identifying 1607 studies. Among these, 22 studies were included in the final analysis, comprising 3937 subjects (1665 SO and 2272 NW) in whom serum TAC was measured, and from these 19,201 subjects, the correlation of serum TAC with anthropo-metabolic parameters was also estimated. The Newcastle-Ottawa method was used for the evaluation of the risk of bias. Using a random-effect model (REM), TAC was reduced in SO independently of age (SMD, -0.86; 95% CI -1.38 to -0.34; p = 0.0012), whereas malondialdehyde (SMD, 1.50; 95% CI 0.60 to 2.41), oxidative stress index (SMD, 1.0; 95% CI 0.16 to 1.84), and total oxidant status (SMD, 0.80; 0.22 to 1.37) were increased. There were seven significant pooled correlations of TAC with anthropometric and metabolic parameters: weight (r = -0.17), hip circumference (r= -0.11), visceral adipose index (r = 0.29), triglycerides (r = 0.25), aspartate aminotransferase (r = 0.41), alanine aminotransferase (r = 0.38), and uric acid (r = 0.53). Our results confirm a decrease in TAC and an increase in markers of oxidative stress in SO and underpin the importance of these serum biomarkers in obesity.

DNA Methylation Remodeling after Bariatric Surgery Correlates with Clinical Parameters.

The altered functions of adipose tissue are one of the main issues in obesity. Bariatric surgery is associated with improvement of obesity associated comorbidities. Here DNA methylation remodeling in adipose tissue after bariatric surgery is examined. After six months postoperative, DNA methylation shows changes in 1155 CpG sites, 66 of these sites correlate with body mass index. Some sites also show correlation with LDL-C, HDL-C, total cholesterol, and triglycerides. CpG sites are located in genes that have not previously been linked to obesity or metabolic diseases. GNAS complex locus is one of those that presented CpG site with the greatest changes after surgery, and the most significant correlation with BMI and lipid profiles. These results show that epigenetic regulation may be involved in the alteration of adipose tissue functions in obesity.

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 >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 > 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.

Frequent copy number variants in a cohort of Mexican-Mestizo individuals.

BACKGROUND: The human genome presents variation at distinct levels, copy number variants (CNVs) are DNA segments of variable lengths that range from several base pairs to megabases and are present at a variable number of copies in human genomes. Common CNVs have no apparent influence on the phenotype; however, some rare CNVs have been associated with phenotypic traits, depending on their size and gene content. CNVs are detected by microarrays of different densities and are generally visualized, and their frequencies analysed using the HapMap as default reference population. Nevertheless, this default reference is inadequate when the samples analysed are from people from Mexico, since population with a Hispanic genetic background are minimally represented. In this work, we describe the variation in the frequencies of four common CNVs in Mexican-Mestizo individuals. RESULTS: In a cohort of 147 unrelated Mexican-Mestizo individuals, we found that the common CNVs 2p11.2 (99.6%), 8p11.22 (54.5%), 14q32.33 (100%), and 15q11.2 (71.1%) appeared with unexpectedly high frequencies when contrasted with the HapMap reference (ChAS). Yet, while when comparing to an ethnically related reference population, these differences were significantly reduced or even disappeared. CONCLUSION: The findings in this work contribute to (1) a better description of the CNVs characteristics of the Mexican Mestizo population and enhance the knowledge of genome variation in different ethnic groups. (2) emphasize the importance of contrasting CNVs identified in studied individuals against a reference group that-as best as possible-share the same ethnicity while keeping this relevant information in mind when conducting CNV studies at the population or clinical level.

Brain radiotoxicity-related 15CAcBRT gene expression signature predicts survival prognosis of glioblastoma patients.

BACKGROUND: Glioblastoma is the most common and devastating primary brain cancer. Radiotherapy is standard of care; however, it is associated with brain radiation toxicity (BRT). This study used a multi-omics approach to determine whether BRT-related genes (RGs) harbor survival prognostic value and whether their encoded proteins represent novel therapeutic targets for glioblastoma. METHODS: RGs were identified through analysis of single-nucleotide variants associated with BRT (R-SNVs). Functional relationships between RGs were established using Protein-Protein Interaction networks. The influence of RGs and their functional groups on glioblastoma prognosis was evaluated using clinical samples from the Glioblastoma Bio-Discovery Portal database and validated using the Chinese Glioma Genome Atlas dataset. The identification of clusters of radiotoxic and putative pathogenic variants in proteins encoded by RGs was achieved by computational 3D structural analysis. RESULTS: We identified the BRT-related 15CAcBRT molecular signature with prognostic value in glioblastoma, by analysis of the COMT and APOE protein functional groups. Its external validation confirmed clinical relevance independent of age, MGMT promoter methylation status, and IDH mutation status. Interestingly, the genes IL6, APOE, and MAOB documented significant gene expression levels alteration, useful for drug repositioning. Biological networks associated with 15CAcBRT signature involved pathways relevant to cancer and neurodegenerative diseases. Analysis of 3D clusters of radiotoxic and putative pathogenic variants in proteins coded by RGs unveiled potential novel therapeutic targets in neuro-oncology. CONCLUSIONS: 15CAcBRT is a BRT-related molecular signature with prognostic significance for glioblastoma patients and represents a hub for drug repositioning and development of novel therapies.

Unraveling Signatures of Local Adaptation among Indigenous Groups from Mexico.

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

DNA methylation and gene expression analysis in adipose tissue to identify new loci associated with T2D development in obesity.

BACKGROUND: Obesity is accompanied by excess adipose fat storage, which may lead to adipose dysfunction, insulin resistance, and type 2 diabetes (T2D). Currently, the tendency to develop T2D in obesity cannot be explained by genetic variation alone-epigenetic mechanisms, such as DNA methylation, might be involved. Here, we aimed to identify changes in DNA methylation and gene expression in visceral adipose tissue (VAT) that might underlie T2D susceptibility in patients with obesity. METHODS: We investigated DNA methylation and gene expression in VAT biopsies from 19 women with obesity, without (OND = 9) or with T2D (OD = 10). Differences in genome-scale methylation (differentially methylated CpGs [DMCs], false discovery rate < 0.05; and differentially methylated regions [DMRs], p value < 0.05) and gene expression (DEGs, p value <0.05) between groups were assessed. We searched for overlap between altered methylation and expression and the impact of altered DNA methylation on gene expression, using bootstrap Pearson correlation. The relationship of altered DNA methylation to T2D-related traits was also tested. RESULTS: We identified 11 120 DMCs and 96 DMRs distributed across all chromosomes, with the greatest density of epigenomic alterations at the MHC locus. These alterations were found in newly and previously T2D-related genes. Several of these findings were supported by validation and extended multi-ethnic analyses. Of 252 DEGs in the OD group, 68 genes contained DMCs (n = 88), of which 24 demonstrated a significant relationship between gene expression and methylation (p values <0.05). Of these, 16, including ATP11A, LPL and EHD2 also showed a significant correlation with fasting glucose and HbA1c levels. CONCLUSIONS: Our results revealed novel candidate genes related to T2D pathogenesis in obesity. These genes show perturbations in DNA methylation and expression profiles in patients with obesity and diabetes. Methylation profiles were able to discriminate OND from OD individuals; DNA methylation is thus a potential biomarker.

Kazak mitochondrial genomes provide insights into the human population history of Central Eurasia.

As a historical nomadic group in Central Asia, Kazaks have mainly inhabited the steppe zone from the Altay Mountains in the East to the Caspian Sea in the West. Fine scale characterization of the genetic profile and population structure of Kazaks would be invaluable for understanding their population history and modeling prehistoric human expansions across the Eurasian steppes. With this mind, we characterized the maternal lineages of 200 Kazaks from Jetisuu at mitochondrial genome level. Our results reveal that Jetisuu Kazaks have unique mtDNA haplotypes including those belonging to the basal branches of both West Eurasian (R0, H, HV) and East Eurasian (A, B, C, D) lineages. The great diversity observed in their maternal lineages may reflect pivotal geographic location of Kazaks in Eurasia and implies a complex history for this population. Comparative analyses of mitochondrial genomes of human populations in Central Eurasia reveal a common maternal genetic ancestry for Turko-Mongolian speakers and their expansion being responsible for the presence of East Eurasian maternal lineages in Central Eurasia. Our analyses further indicate maternal genetic affinity between the Sherpas from the Tibetan Plateau with the Turko-Mongolian speakers.

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.

Ancestry-dependent genetic structure of the Xq28 risk haplotype in the Mexican population and its association with childhood-onset systemic lupus erythematosus.

Objective: Here we aimed to investigate the association of the Xq28 risk haplotype (H1) with susceptibility to childhood-onset systemic lupus erythematosus (SLE), and to compare its frequency and genetic structure in the Mexican population with those in other continental populations. Methods: We genotyped 15 single-nucleotide variants (SNVs) that form the H1 haplotype, using TaqMan real-time PCR. The association analysis [case-control and transmission disequilibrium test (TDT)] included 376 cases and 400 adult controls, all of whom were mestizos (MEZ). To identify risk alleles in Mexican Indigenous individuals, SNVs were imputed from whole-exome sequencing data of 1,074 individuals. The allelic frequencies determined in MEZ and Indigenous individuals were compared with those of the continental populations from the 1,000 Genomes database phase 3. Linkage disequilibrium (LD) analysis of risk alleles was performed on all populations. Interleukin-1 receptor associated kinase 1 (IRAK1) and methyl CpG binding protein 2 (MECP2) mRNA levels were determined using real-time PCR. Results: Case-control analysis revealed genetic association with childhood-onset SLE for all 15 SNVs (OR = 1.49-1.75; p = 0.0095 to 1.81 × 10-4) and for the Xq28 risk haplotype (OR = 1.97, p = 4 × 10-6). Comparing with individuals of European ancestry (0.14-0.16), the frequencies of the risk alleles were significantly higher in the MEZ individuals (0.55-0.68) and even higher in Indigenous individuals (0.57-0.83). LD analysis indicated a differential haplotype structure within the Indigenous groups, which was inherited to the MEZ population as a result of genetic admixture. Individuals homozygous for the Xq28 risk haplotype exhibited decreased levels of both MECP2A and B transcripts. Conclusion: We found that the H1 risk haplotype differs in its conformation in the Mexican population. This difference could be attributed to positive selection within the Indigenous population, with its inheritance now having an autoimmune health impact in both the Mexican Indigenous and MEZ populations.

Replication of Integrative Data Analysis for Adipose Tissue Dysfunction, Low-Grade Inflammation, Postprandial Responses and OMICs Signatures in Symptom-Free Adults.

We previously reported preliminary characterization of adipose tissue (AT) dysfunction through the adiponectin/leptin ratio (ALR) and fasting/postprandial (F/P) gene expression in subcutaneous (SQ) adipose tissue (AT) biopsies obtained from participants in the GEMM study, a precision medicine research project. Here we present integrative data replication of previous findings from an increased number of GEMM symptom-free (SF) adults (N = 124) to improve characterization of early biomarkers for cardiovascular (CV)/immunometabolic risk in SF adults with AT dysfunction. We achieved this goal by taking advantage of the rich set of GEMM F/P 5 h time course data and three tissue samples collected at the same time and frequency on each adult participant (F/P blood, biopsies of SQAT and skeletal muscle (SKM)). We classified them with the presence/absence of AT dysfunction: low (<1) or high (>1) ALR. We also examined the presence of metabolically healthy (MH)/unhealthy (MUH) individuals through low-grade chronic subclinical inflammation (high sensitivity C-reactive protein (hsCRP)), whole body insulin sensitivity (Matsuda Index) and Metabolic Syndrome criteria in people with/without AT dysfunction. Molecular data directly measured from three tissues in a subset of participants allowed fine-scale multi-OMIC profiling of individual postprandial responses (RNA-seq in SKM and SQAT, miRNA from plasma exosomes and shotgun lipidomics in blood). Dynamic postprandial immunometabolic molecular endophenotypes were obtained to move towards a personalized, patient-defined medicine. This study offers an example of integrative translational research, which applies bench-to-bedside research to clinical medicine. Our F/P study design has the potential to characterize CV/immunometabolic early risk detection in support of precision medicine and discovery in SF individuals.

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.

Alterations of DNA methylation during adipogenesis differentiation of mesenchymal stem cells isolated from adipose tissue of patients with obesity is associated with type 2 diabetes.

Adipogenesis regulation is crucial for mature adipocyte function. In obesity, a major driver of type 2 diabetes (T2D), this process is disrupted and remains poorly characterized. Here we identified altered DNA methylation profiles in diabetic obese patients, during three adipocytes differentiation stages. We isolated mesenchymal cells from visceral adipose tissue of obese patients with and without T2D to analyse DNA methylation profiles at 0, 3, and 18 days of ex vivo differentiation and documented their impact on gene expression. Methylation and gene expression were analysed with EPIC and Clarion S arrays, respectively. Patients with T2D had epigenetic alterations in all the analysed stages, and these were mainly observed in genes important in adipogenesis, insulin resistance, cell death programming, and immune effector processes. Importantly, at 3 days, we found six-fold more methylated CpG alterations than in the other stages. This is the first study to document epigenetic markers that persist through all three adipogenesis stages and their impact on gene expression, which could be a cellular metabolic memory involved in T2D. Our data provided evidence that, throughout the adipogenesis process, alterations occur in methylation that might impact mature adipocyte function, cause tissue malfunction, and potentially, lead to the development of T2D.

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

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

Two novel variants in DYRK1B causative of AOMS3: expanding the clinical spectrum.

BACKGROUND: We investigated pathogenic DYRK1B variants causative of abdominal obesity-metabolic syndrome 3 (AOMS3) in a group of patients originally diagnosed with type 2 diabetes. All DYRK1B exons were analyzed in a sample of 509 unrelated adults with type 2 diabetes and 459 controls, all belonging to the DMS1 SIGMA-cohort (ExAC). We performed in silico analysis on missense variants using Variant Effect Predictor software. To evaluate co-segregation, predicted pathogenic variants were genotyped in other family members. We performed molecular dynamics analysis for the co-segregating variants. RESULTS: After filtering, Mendelian genotypes were confirmed in two probands bearing two novel variants, p.Arg252His and p.Lys68Gln. Both variants co-segregated with the AOMS3 phenotype in classic dominant autosomal inheritance with full penetrance. In silico analysis revealed impairment of the DYRK1B protein function by both variants. For the first time, we describe age-dependent variable expressivity of this entity, with central obesity and insulin resistance apparent in childhood; morbid obesity, severe hypertriglyceridemia, and labile type 2 diabetes appearing before 40 years of age; and hypertension emerging in the fifth decade of life. We also report the two youngest individuals suffering from AOMS3. CONCLUSIONS: Monogenic forms of metabolic diseases could be misdiagnosed and should be suspected in families with several affected members and early-onset metabolic phenotypes that are difficult to control. Early diagnostic strategies and medical interventions, even before symptoms or complications appear, could be useful.

Reconstruction of ancient microbial genomes from the human gut.

Loss of gut microbial diversity1-6 in industrial populations is associated with chronic diseases7, underscoring the importance of studying our ancestral gut microbiome. However, relatively little is known about the composition of pre-industrial gut microbiomes. Here we performed a large-scale de novo assembly of microbial genomes from palaeofaeces. From eight authenticated human palaeofaeces samples (1,000-2,000 years old) with well-preserved DNA from southwestern USA and Mexico, we reconstructed 498 medium- and high-quality microbial genomes. Among the 181 genomes with the strongest evidence of being ancient and of human gut origin, 39% represent previously undescribed species-level genome bins. Tip dating suggests an approximate diversification timeline for the key human symbiont Methanobrevibacter smithii. In comparison to 789 present-day human gut microbiome samples from eight countries, the palaeofaeces samples are more similar to non-industrialized than industrialized human gut microbiomes. Functional profiling of the palaeofaeces samples reveals a markedly lower abundance of antibiotic-resistance and mucin-degrading genes, as well as enrichment of mobile genetic elements relative to industrial gut microbiomes. This study facilitates the discovery and characterization of previously undescribed gut microorganisms from ancient microbiomes and the investigation of the evolutionary history of the human gut microbiota through genome reconstruction from palaeofaeces.

Towards precision medicine: defining and characterizing adipose tissue dysfunction to identify early immunometabolic risk in symptom-free adults from the GEMM family study.

Interactions between macrophages and adipocytes are early molecular factors influencing adipose tissue (AT) dysfunction, resulting in high leptin, low adiponectin circulating levels and low-grade metaflammation, leading to insulin resistance (IR) with increased cardiovascular risk. We report the characterization of AT dysfunction through measurements of the adiponectin/leptin ratio (ALR), the adipo-insulin resistance index (Adipo-IRi), fasting/postprandial (F/P) immunometabolic phenotyping and direct F/P differential gene expression in AT biopsies obtained from symptom-free adults from the GEMM family study. AT dysfunction was evaluated through associations of the ALR with F/P insulin-glucose axis, lipid-lipoprotein metabolism, and inflammatory markers. A relevant pattern of negative associations between decreased ALR and markers of systemic low-grade metaflammation, HOMA, and postprandial cardiovascular risk hyperinsulinemic, triglyceride and GLP-1 curves was found. We also analysed their plasma non-coding microRNAs and shotgun lipidomics profiles finding trends that may reflect a pattern of adipose tissue dysfunction in the fed and fasted state. Direct gene differential expression data showed initial patterns of AT molecular signatures of key immunometabolic genes involved in AT expansion, angiogenic remodelling and immune cell migration. These data reinforce the central, early role of AT dysfunction at the molecular and systemic level in the pathogenesis of IR and immunometabolic disorders.

Metabolic syndrome in indigenous communities in Mexico: a descriptive and cross-sectional study.

BACKGROUND: An Amerindian genetic background could play an important role in susceptibility to metabolic diseases, which have alarmingly increased in recent decades. Mexico has one of the highest prevalences of metabolic disease worldwide. The purpose of this study was to determine the prevalence of metabolic syndrome and its components in a population with high Amerindian ancestry. METHODS: We performed a descriptive, quantitative, and analytical cross-sectional study of 2596 adult indigenous volunteers from 60 different ethnic groups. Metabolic syndrome and its components were evaluated using the American Heart Association/National Heart, Lung, and Blood Institute Scientific Statement criteria. RESULTS: The overall prevalence of metabolic syndrome in the indigenous Mexican population was 50.3%. Although females had a higher prevalence than males (55.6% vs. 38.2%), the males presented with combinations of metabolic syndrome components that confer a higher risk of cardiovascular disease. The most frequent metabolic syndrome component in both genders was low HDL-cholesterol levels (75.8%). Central obesity was the second most frequent component in females (61%), though it had a low prevalence in males (16.5%). The overall prevalence of elevated blood pressure was 42.7% and was higher in males than females (48.8 vs. 40%). We found no gender differences in the overall prevalence of elevated triglycerides (56.7%) or fasting glucose (27.9%). CONCLUSIONS: We documented that individuals with Amerindian ancestry have a high prevalence of metabolic syndrome. Health policies are needed to control the development of metabolic disorders in a population with high genetic risk.

Genetic variability of five ADRB2 polymorphisms among Mexican Amerindian ethnicities and the Mestizo population.

The Mexican population is characterized by high and particular admixture, and the picture of variants associated with disease remains unclear. Here we investigated the distribution of single nucleotide polymorphisms (SNPs) in the Mexican population. We focused on two non-synonymous and three synonymous SNPs in the beta-2 adrenergic receptor gene (ADRB2), which plays key roles in energy balance regulation. These SNPs were genotyped in 2,011 Mexican Amerindians (MAs) belonging to 62 ethnic groups and in 1,980 geographically matched Mexican Mestizos (MEZs). The frequency distribution of all five ADRB2 variants significantly differed between MAs, MEZs, and other continental populations (CPs) from the 1000 Genomes database. Allele frequencies of the three synonymous SNPs rs1042717A, rs1042718A, and rs1042719C were significantly higher in Mexican individuals, particularly among MAs, compared to in the other analyzed populations (P<0.05). The non-synonymous ADRB2 Glu27 allele (rs1042714G), which is associated with several common conditions, showed the lowest frequency in MAs (0.03) compared to other populations worldwide. Among MEZs, this allele showed a frequency of 0.15, intermediate between that in MAs and in Iberians (0.43). Moreover, Glu27 was the only SNP exhibiting a geographic gradient within the MEZ population (from 0.22 to 0.11), reflecting admixed mestizo ancestry across the country. Population differentiation analysis demonstrated that Glu27 had the highest FST value in MAs compared with Europeans (CEU) (0.71), and the lowest between MAs and Japanese (JPT) (0.01), even lower than that observed between MAs and MEZs (0.08). This analysis demonstrated the genetic diversity among Amerindian ethnicities, with the most extreme FST value (0.34) found between the Nahuatls from Morelos and the Seris. This is the first study of ADRB2 genetic variants among MA ethnicities. Our findings add to our understanding of the genetic contribution to variability in disease susceptibility in admixed populations.