Search

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

Abstract The Nav1.7 voltage-gated sodium channel plays a key role in nociception. Three functional variants in the SCN9A gene (encoding M932L, V991L, and D1908G in Nav1.7), have recently been identified as stemming from Neanderthal introgression and to associate with pain symptomatology in UK BioBank data. In 1000 genomes data, these variants are absent in Europeans but common in Latin Americans. Analysing high-density genotype data from 7594 Latin Americans, we characterized Neanderthal introgression in SCN9A. We find that tracts of introgression occur on a Native American genomic background, have an average length of ~123 kb and overlap the M932L, V991L, and D1908G coding positions. Furthermore, we measured experimentally six pain thresholds in 1623 healthy Colombians. We found that Neanderthal ancestry in SCN9A is significantly associated with a lower mechanical pain threshold after sensitization with mustard oil and evidence of additivity of effects across Nav1.7 variants. Our findings support the reported association of Neanderthal Nav1.7 variants with clinical pain, define a specific sensory modality affected by archaic introgression in SCN9A and are consistent with independent effects of the Neanderthal variants on Nav1.7 function.

Abstract We report a genome-wide association study of facial features in >6000 Latin Americans based on automatic landmarking of 2D portraits and testing for association with inter-landmark distances. We detected significant associations (P-value

Postmenopausal women are at an increased risk of developing metabolic syndrome (MetS) due to hormonal changes and lifestyle factors. Gut microbiota (GM) have been linked to the development of MetS, and they are influenced by dietary habits. However, the interactions between dietary patterns (DP) and the GM of postmenopausal women, as well as their influence on MetS, still need to be understood. The present study evaluated the DP and microbiota composition of postmenopausal Mexican women with MetS and those in a control group. Diet was assessed using a food frequency questionnaire, and the GM were profiled using 16S rRNA gene sequencing. Greater adherence to a “healthy” DP was significantly associated with lower values of MetS risk factors. GM diversity was diminished in women with MetS, and it was negatively influenced by an “unhealthy” DP. Moreover, a higher intake of fats and proteins, as well as lower amounts of carbohydrates, showed a reduction in some of the short-chain fatty acid-producing genera in women with MetS, as well as increases in some harmful bacteria. Furthermore, Roseburia abundance was positively associated with dietary fat and waist circumference, which may explain 7.5% of the relationship between this macronutrient and MetS risk factors. These findings suggest that GM and diet interactions are important in the development of MetS in postmenopausal Mexican women.

Abstract Background Elevations of circulating branched-chain amino acids (BCAA) are observed in humans with obesity and metabolic comorbidities, such as insulin resistance. Although it has been described that microbial metabolism contributes to the circulating pool of these amino acids, studies are still scarce, particularly in pediatric populations. Thus, we aimed to explore whether in early adolescents, gut microbiome was associated to circulating BCAA and in this way to insulin resistance. Methods Shotgun sequencing was performed in DNA from fecal samples of 23 early adolescents (10–12 years old) and amino acid targeted metabolomics analysis was performed by LC–MS/MS in serum samples. By using the HUMAnN2 algorithm we explored microbiome functional profiles to identify whether bacterial metabolism contributed to serum BCAA levels and insulin resistance markers. Results We identified that abundance of genes encoding bacterial BCAA inward transporters were negatively correlated with circulating BCAA and HOMA-IR (P

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.

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

Abstract Background In the last decade, increasing evidence has shown that changes in human gut microbiota are associated with diseases, such as obesity. The excreted/secreted proteins (secretome) of the gut microbiota affect the microbial composition, altering its colonization and persistence. Furthermore, it influences microbiota-host interactions by triggering inflammatory reactions and modulating the host's immune response. The metatranscriptome is essential to elucidate which genes are expressed under diseases. In this regard, little is known about the expressed secretome in the microbiome. Here, we use a metatranscriptomic approach to delineate the secretome of the gut microbiome of Mexican children with normal weight (NW) obesity (O) and obesity with metabolic syndrome (OMS). Additionally, we performed the 16S rRNA profiling of the gut microbiota. Results Out of the 115,712 metatranscriptome genes that codified for proteins, 30,024 (26%) were predicted to be secreted, constituting the Secrebiome of the gut microbiome. The 16S profiling confirmed an increased abundance in Firmicutes and decreased in Bacteroidetes in the obesity groups, and a significantly higher richness and diversity than the normal weight group. We found novel biomarkers for obesity with metabolic syndrome such as increased Coriobacteraceae, Collinsela, and Collinsella aerofaciens; Erysipelotrichaceae, Catenibacterium and Catenibacterium sp., and decreased Parabacteroides distasonis, which correlated with clinical and anthropometric parameters associated to obesity and metabolic syndrome. Related to the Secrebiome, 16 genes, homologous to F. prausniitzi, were overexpressed for the obese and 15 genes homologous to Bacteroides, were overexpressed in the obesity with metabolic syndrome. Furthermore, a significant enrichment of CAZy enzymes was found in the Secrebiome. Additionally, significant differences in the antigenic density of the Secrebiome were found between normal weight and obesity groups. Conclusions These findings show, for the first time, the role of the Secrebiome in the functional human-microbiota interaction. Our results highlight the importance of metatranscriptomics to provide novel information about the gut microbiome’s functions that could help us understand the impact of the Secrebiome on the homeostasis of its human host. Furthermore, the metatranscriptome and 16S profiling confirmed the importance of treating obesity and obesity with metabolic syndrome as separate conditions to better understand the interplay between microbiome and disease.

Background: Gut microbiota is a complex organized collection of microorganisms that confers multiple metabolic advantages to the host. The reduced diversity and proportion of specific gut microbial species have been associated with obesity and metabolic disorders. Multidimensional interventions, including modifications in dietary and physical activity habits, are associated with favorable changes in microbiota composition. This pilot study aimed to evaluate changes in the gut microbiota composition of Mexican children with obesity before and after a 6-week multidimensional intervention. Methods: Blood and stool samples were collected, and anthropometric measurements were obtained from six children with obesity before and after the intervention. The intervention consisted of modeling a hypo energetic diet and giving nutritional and physical activation recommendations. DNA from stool samples was used to characterize gut microbial composition by sequencing the 16S rRNA gene. Results: The decrease in waist circumference was associated with increased Odoribacter relative abundance. However, gut microbiota composition and diversity remained unchanged. Conclusions: Although no modifications in the body mass index, body fat, composition, or diversity of the gut microbiota were observed with the intervention, it was possible to associate the reduction in waist circumference with the presence of Odoribacter after a multidimensional intervention in Mexican children with obesity.

The main roles of adipose tissue include triglycerides storage and adipokine secretion, which regulate energy balance and inflammation status. In obesity, adipocyte dysfunction leads to proinflammatory cytokine production and insulin resistance. Bariatric surgery is the most effective treatment for obesity, the gold-standard technique being Roux-en-Y gastric bypass (RYGB). Since metabolic improvements after RYGB are clear, a better understanding of adipose tissue molecular modifications could be derived from this study. Thus, the aim of this systematic review was to find differentially expressed genes in subcutaneous adipose tissue of lean, obese and post-RYGB (distinct timepoints). To address this objective, publications from 2015–2022 reporting gene expression (candidate genes or transcriptomic approach) of subcutaneous adipose tissue from lean and obese individuals before and after RGYB were searched in PubMed, Elsevier, and Springer Link. Excluded publications were reviews, studies analyzing serum, other types of tissues, or bariatric procedures. A risk-of-bias summary was created for each paper using Robvis, to finally include 17 studies. Differentially expressed genes in post-RYGB vs. obese and lean vs. obese were obtained and the intersection among these groups was used for analysis and gene classification by metabolic pathway. Results showed that the lean state as well as the post-RYGB is similar in terms of increased expression of insulin-sensitizing molecules, inducing lipogenesis over lipolysis and downregulating leukocyte activation, cytokine production and other factors that promote inflammation. Thus, massive weight loss and metabolic improvements after RYGB are accompanied by gene expression modifications reverting the “adipocyte dysfunction” phenomenon observed in obesity conditions.

Metabolic syndrome (MetS) is a group of several metabolic conditions predisposing to chronic diseases. Individuals diagnosed with MetS are physiologically heterogeneous, with significant sex-specific differences. Therefore, we aimed to investigate the potential sex-specific serum modifications of amino acids and acylcarnitines (ACs) and their relationship with MetS in the Mexican population. This study included 602 participants from the Health Workers Cohort Study. Forty serum metabolites were analyzed using a targeted metabolomics approach. Multivariate regression models were used to test associations of clinical and biochemical parameters with metabolomic profiles. Our findings showed a serum amino acid signature (citrulline and glycine) and medium-chain ACs (AC14:1, AC10, and AC18:10H) associated with MetS. Glycine and AC10 were specific metabolites representative of discrimination according to sex-dependent MetS. In addition, we found that glycine and short-chain ACs (AC2, AC3, and AC8:1) are associated with age-dependent MetS. We also reported a significant correlation between body fat and metabolites associated with sex-age-dependent MetS. In conclusion, the metabolic profile varies by MetS status, and these differences are sex-age-dependent in the Mexican population.

Summary: Changes in the human gut microbiome are associated with obesity and metabolic syndrome, but the role of the gut virome in both diseases remains largely unknown. We characterized the gut dsDNA virome of 28 school-aged children with healthy normal-weight (NW, n = 10), obesity (O, n = 10), and obesity with metabolic syndrome (OMS, n = 8), using metagenomic sequencing of virus-like particles (VLPs) from fecal samples. The virome classification confirmed the bacteriophages' dominance, mainly composed of Caudovirales. Notably, phage richness and diversity of individuals with O and OMS tended to increase, while the VLP abundance remained the same among all groups. Of the 4,611 phage contigs composing the phageome, 48 contigs were highly prevalent in ≥80% of individuals, suggesting high inter-individual phage diversity. The abundance of several contigs correlated with gut bacterial taxa; and with anthropometric and biochemical parameters altered in O and OMS. To our knowledge, this gut phageome represents one of the largest datasets and suggests disease-specific phage alterations.

Pigmentation variation in humans is influenced by complex genetic architecture in different populations. Here, the authors perform a genome-wide association analysis involving > 6,000 Latin Americans for pigmentation of skin and eyes, and identify known and novel genetic associations.

Gut microbiota has been suggested to modulate circulating lipids. However, the relationship between the gut microbiota and atherogenic dyslipidemia (AD), defined as the presence of both low HDL-C and hypertriglyceridemia, is not fully understood. Moreover, because obesity is among the main causes of secondary AD, it is important to analyze the effect of gut microbiota composition on lipid profiles after a weight loss intervention. We compared the microbial diversity and taxonomic composition in patients with AD (n = 41) and controls (n = 38) and sought correlations of genera abundance with serum lipid levels in 20 patients after weight loss induced by Roux-en-Y gastric bypass (RYGB) surgery. Gut microbiota composition was profiled using next-generation sequencing of 16S rRNA. Gut microbiota diversity was significantly lower in atherogenic dyslipidemia. Moreover, relative abundance of two genera with LDA score >3.5 (Megasphaera and LPS-producing Escherichia-Shigella), was significantly higher in AD subjects, while the abundance of four short chain fatty acids (SCFA) producing-genera (Christensenellaceae R-7, Ruminococcaceae UCG-014; Akkermansia and [Eubacterium] eligens group) was significantly higher in controls. Notably, [Eubacterium] eligens group abundance was also significantly associated with higher HDL-C levels in RYGB patients one year after surgery. Although dietary polyunsaturated fatty acid/saturated fatty acid (PUFA/SFA) ratio and PUFA intake were higher in controls than in AD subjects, of the four genera differentiated in cases and controls, only Akkermansia abundance showed a positive and significant correlation with PUFA/SFA ratio. Our results suggest that SCFA-producing bacteria promote a healthy lipid homeostasis, while the presence of LPS-producing bacteria such Escherichia-Shigella may contribute to the development of atherogenic dyslipidemia.

Latin Americans trace their ancestry to the admixture of Native Americans, Europeans and Sub-Saharan Africans. Here, the authors develop a novel haplotype-based approach and analyse over 6,500 Latin Americans to infer the geographically-detailed genetic structure of this population.

Background. The gene for patatin-like phospholipase domain containing 3 (PNPLA3) is associated with non-alcoholic fatty liver disease (NAFLD) development. We previously found that Mexican indigenous population had the highest frequency reported of the PNPLA3 148M risk allele. Further, we observed a relationship between M148M genotype with elevated ALT levels in individuals with normal weight, overweight and obese. We sought to investigate whether PNPLA3 polymorphism is associated with NAFLD development in Mexicans.Material and methods. We enrolled 189 Mexican patients with NAFLD and 201 healthy controls. Anthropometric, metabolic, and biochemical variables were measured, and rs738409 (Ile148Met substitution) polymorphism was genotyped by sequencing.Results. Logistic regression analysis, using a recessive model, suggested that PNPLA3 polymorphism in Mexican population is significantly associated (OR = 1.711, 95% CI: 1.014-2.886; P = 0.044) with NAFLD.Conclusions: The PNPLA3 gene is associated with NAFLD in Mexican population. More studies are required to explain the high prevalence of PNPLA3 polymorphism in Mexican-Americans, Mexican-Indians, and Mexican-Mestizos.

Abstract Objectives crAssphage is a newly found phage described as the most abundant virus in the human gut microbiome. The majority of the crAssphage proteins are unknown in sequences databases, and its pathogenicity and epidemiology in humans are yet unclear. Hence, being one of the most abundant phages in the human gut microbiome more investigation at the genomic level is necessary to improve our understanding, especially in the Latin American population. Data description In this article, we provide the whole genome of a crAssphage isolated from the human gut microbiome of the Mexican population, which was named Mexican-crAssphage. The genome consists of 96,283 bp, G+C content of 29.24% and 87 coding sequences. Notably, we did not find any transfer RNA genes in the genome sequence. We also sequenced viral-like enriched particles from 28 fecal samples, and we detected the presence of the Mexican-crAssphage genome in 8 samples (28.5%). To our knowledge, our data is the first whole genome report of the crAssphage isolated from the Latin American Population and provides valuable information for the experimental characterization of the most abundant human gut bacteriophage. The whole genome shotgun project of the Mexican-crAssphage is available at DDBJ/ENA/GenBank under the GenBank MK069403.

Dietary fiber (DF) is a major substrate for the gut microbiota that contributes to metabolic health. Recent studies have shown that diet–metabolic phenotype effect might be related to individual gut microbial profiles or enterotypes. Thus, the aim of this study was to examine whether microbial enterotypes modify the association between DF intake and metabolic traits. This cross-sectional study included 204 children (6–12 years old) and 75 adults (18–60 years old). Habitual DF intake was estimated with a Food Frequency Questionnaire and biochemical, clinical and anthropometric data were obtained. Gut microbiota was assessed through 16S sequencing and participants were stratified by enterotypes. Correlations adjusting for age and sex were performed to test the associations between dietary fiber components intake and metabolic traits. In children and adults from the Prevotella enterotype, a nominal negative correlation of hemicellulose intake with insulin and HOMA-IR levels was observed (p < 0.05), while in individuals of the other enterotypes, these associations were not observed. Interestingly, the latter effect was not related to the fecal short-chain-fatty acids profile. Our results contribute to understanding the enterotype influence on the diet–phenotype interaction, which ultimate could provide evidence for their use as potential biomarkers for future precision nutrition strategies.

Recent evidence shows that obesity correlates negatively with bone mass. However, traditional anthropometric measures such as body mass index could not discriminate visceral adipose tissue from subcutaneous adipose tissue. The visceral adiposity index (VAI) is a reliable sex-specified indicator of visceral adipose distribution and function. Thus, we aimed to identify metabolomic profiles associated with VAI and low bone mineral density (BMD). A total of 602 individuals from the Health Workers Cohort Study were included. Forty serum metabolites were measured using the targeted metabolomics approach, and multivariate regression models were used to test associations of metabolomic profiles with anthropometric, clinical, and biochemical parameters. The analysis showed a serum amino acid signature composed of glycine, leucine, arginine, valine, and acylcarnitines associated with high VAI and low BMD. In addition, we found a sex-dependent VAI in pathways related to primary bile acid biosynthesis, branched-chain amino acids, and the biosynthesis of pantothenate and coenzyme A (CoA). In conclusion, a metabolic profile differs by VAI and BMD status, and these changes are gender-dependent.

People of Mexico have diverse historical and genetic background. Here, Romero-Hidalgo and colleagues sequence whole genomes of Native Americans of Mexico, and show demographic history and genetic variation shared among subgroups of Native Americans.

Abstract Childhood obesity is associated with a number of metabolic abnormalities leading to increased cardiovascular risk. Metabolites can be useful as early biomarkers and new targets to promote early intervention beginning in school age. Thus, we aimed to identify metabolomic profiles associated with obesity and obesity-related metabolic traits. We used data from the Obesity Research Study for Mexican children (ORSMEC) in Mexico City and included a case control (n = 1120), cross-sectional (n = 554) and a longitudinal study (n = 301) of 6–12-year-old children. Forty-two metabolites were measured using electrospray MS/MS and multivariate regression models were used to test associations of metabolomic profiles with anthropometric, clinical and biochemical parameters. Principal component analysis showed a serum amino acid signature composed of arginine, leucine/isoleucine, phenylalanine, tyrosine, valine and proline significantly associated with obesity (OR = 1.57; 95%CI 1.45–1.69, P = 3.84 × 10−31) and serum triglycerides (TG) (β = 0.067, P = 4.5 × 10−21). These associations were validated in the cross-sectional study (P

Abstract Background The aim of this study was to explore whether interactions between FTO rs9939609 and ABCA1 rs9282541 affect BMI and waist circumference (WC), and could explain previously reported population differences in FTO-obesity and FTO-BMI associations in the Mexican and European populations. Methods A total of 3938 adults and 636 school-aged children from Central Mexico were genotyped for both polymorphisms. Subcutaneous and visceral adipose tissue biopsies from 22 class III obesity patients were analyzed for FTO and ABCA1 mRNA expression. Generalized linear models were used to test for associations and gene-gene interactions affecting BMI, WC and FTO expression. Results FTO and ABCA1 risk alleles were not individually associated with higher BMI or WC. However, in the absence of the ABCA1 risk allele, the FTO risk variant was significantly associated with higher BMI (P = 0.043) and marginally associated with higher WC (P = 0.067), as reported in Europeans. The gene-gene interaction affecting BMI and WC was statistically significant only in adults. FTO mRNA expression in subcutaneous abdominal adipose tissue according to ABCA1 genotype was consistent with these findings. Conclusions This is the first report showing evidence of FTO and ABCA1 gene variant interactions affecting BMI, which may explain previously reported population differences. Further studies are needed to confirm this interaction.

The human face represents a combined set of highly heritable phenotypes, but knowledge on its genetic architecture remains limited, despite the relevance for various fields. A series of genome-wide association studies on 78 facial shape phenotypes quantified from 3-dimensional facial images of 10,115 Europeans identified 24 genetic loci reaching study-wide suggestive association (p < 5 × 10−8), among which 17 were previously unreported. A follow-up multi-ethnic study in additional 7917 individuals confirmed 10 loci including six unreported ones (padjusted < 2.1 × 10−3). A global map of derived polygenic face scores assembled facial features in major continental groups consistent with anthropological knowledge. Analyses of epigenomic datasets from cranial neural crest cells revealed abundant cis-regulatory activities at the face-associated genetic loci. Luciferase reporter assays in neural crest progenitor cells highlighted enhancer activities of several face-associated DNA variants. These results substantially advance our understanding of the genetic basis underlying human facial variation and provide candidates for future in-vivo functional studies.

Humans show great diversity in facial appearance and this variation is highly heritable. Here, Andres Ruiz-Linares and colleagues examined facial features in admixed Latin Americans and identify genome-wide associations for 14 facial traits, including four gene loci (RUNX2, GLI3, DCHS2 and PAX1) influencing nose morphology.

By examining Latin American individuals of mixed European, Native American and African ancestry, Adhikari et al. identify novel loci influencing various features of facial and scalp hair. The study also provides experimental evidence that one of the implicated genes (PRSS53) is expressed in the hair follicle and that the top associated variant alters processing of this enzyme.

The effect of microbiota composition and its health on bone tissue is a novel field for research. However, their associations with bone mineral density (BMD) have not been established in postmenopausal women. The present study investigates the relation of diet, the microbiota composition, and the serum metabolic profile in postmenopausal women with normal-BMD or with low-BMD. Ninety-two Mexican postmenopausal women were classified into normal-BMD (n = 34) and low-BMD (n = 58). The V4 hypervariable region was sequenced using the Miseq platform. Serum vitamin D was determined by chemiluminescence immunoassay. Serum concentrations of acyl-carnitines and amino acids were determined by electrospray tandem mass spectrometry. Diet was assessed by a food frequency questionnaire. The low-BMD group had fewer observed species, higher abundance of γ-Proteobacteria, lower consumption of lycopene, and lower concentrations of leucine, valine, and tyrosine compared with the normal-BMD group. These amino acids correlated positively with the abundance of Bacteroides. Lycopene consumption positively correlated with Oscillospira and negatively correlated with Pantoea genus abundance. Finally, the intestinal microbiota of women with vitamin D deficiency was related to Erysipelotrichaceae and Veillonellaceae abundance compared to the vitamin D non-deficient group. Associations mediated by the gut microbiota between diet and circulating metabolites with low-BMD were identified.

Hair plays an important role in primates and is clearly subject to adaptive selection. While humans have lost most facial hair, eyebrows are a notable exception. Eyebrow thickness is heritable and widely believed to be subject to sexual selection. Nevertheless, few genomic studies have explored its genetic basis. Here, we performed a genome-wide scan for eyebrow thickness in 2961 Han Chinese. We identified two new loci of genome-wide significance, at 3q26.33 near SOX2 (rs1345417: P = 6.51×10(-10)) and at 5q13.2 near FOXD1 (rs12651896: P = 1.73×10(-8)). We further replicated our findings in the Uyghurs, a population from China characterized by East Asian-European admixture (N = 721), the CANDELA cohort from five Latin American countries (N = 2301), and the Rotterdam Study cohort of Dutch Europeans (N = 4411). A meta-analysis combining the full GWAS results from the three cohorts of full or partial Asian descent (Han Chinese, Uyghur and Latin Americans, N = 5983) highlighted a third signal of genome-wide significance at 2q12.3 (rs1866188: P = 5.81×10(-11)) near EDAR. We performed fine-mapping and prioritized four variants for further experimental verification. CRISPR/Cas9-mediated gene editing provided evidence that rs1345417 and rs12651896 affect the transcriptional activity of the nearby SOX2 and FOXD1 genes, which are both involved in hair development. Finally, suitable statistical analyses revealed that none of the associated variants showed clear signals of selection in any of the populations tested. Contrary to popular speculation, we found no evidence that eyebrow thickness is subject to strong selective pressure.

Background In spite of the emergence of RNA sequencing (RNA-seq), microarrays remain in widespread use for gene expression analysis in the clinic. There are over 767,000 RNA microarrays from human samples in public repositories, which are an invaluable resource for biomedical research and personalized medicine. The absolute gene expression analysis allows the transcriptome profiling of all expressed genes under a specific biological condition without the need of a reference sample. However, the background fluorescence represents a challenge to determine the absolute gene expression in microarrays. Given that the Y chromosome is absent in female subjects, we used it as a new approach for absolute gene expression analysis in which the fluorescence of the Y chromosome genes of female subjects was used as the background fluorescence for all the probes in the microarray. This fluorescence was used to establish an absolute gene expression threshold, allowing the differentiation between expressed and non-expressed genes in microarrays. Methods We extracted the RNA from 16 children leukocyte samples (nine males and seven females, ages 6–10 years). An Affymetrix Gene Chip Human Gene 1.0 ST Array was carried out for each sample and the fluorescence of 124 genes of the Y chromosome was used to calculate the absolute gene expression threshold. After that, several expressed and non-expressed genes according to our absolute gene expression threshold were compared against the expression obtained using real-time quantitative polymerase chain reaction (RT-qPCR). Results From the 124 genes of the Y chromosome, three genes (DDX3Y, TXLNG2P and EIF1AY) that displayed significant differences between sexes were used to calculate the absolute gene expression threshold. Using this threshold, we selected 13 expressed and non-expressed genes and confirmed their expression level by RT-qPCR. Then, we selected the top 5% most expressed genes and found that several KEGG pathways were significantly enriched. Interestingly, these pathways were related to the typical functions of leukocytes cells, such as antigen processing and presentation and natural killer cell mediated cytotoxicity. We also applied this method to obtain the absolute gene expression threshold in already published microarray data of liver cells, where the top 5% expressed genes showed an enrichment of typical KEGG pathways for liver cells. Our results suggest that the three selected genes of the Y chromosome can be used to calculate an absolute gene expression threshold, allowing a transcriptome profiling of microarray data without the need of an additional reference experiment. Discussion Our approach based on the establishment of a threshold for absolute gene expression analysis will allow a new way to analyze thousands of microarrays from public databases. This allows the study of different human diseases without the need of having additional samples for relative expression experiments.

En este trabajo se evalúa la relación entre la asimetría fluctuante facial (AFF) y los tratamientos hormonales, cirugías maxilofaciales, ortodoncia, traumatismos y malformaciones. En el marco del proyecto CANDELA, se tomaron cinco fotografías faciales de 3162 voluntarios entre los 18 y 85 años. Por fotogrametría se colocaron 34 landmarks o puntos en 3D y mediante el método Procrustes ANOVA se obtuvieron valores individuales de asimetría fluctuante facial. Se realizó una prueba de ANOVA de una vía y la prueba de Welch y Levene para conocer las diferencias entre media y varianza de los valores de asimetría facial y las variables respuesta. También, se caracterizó la variación morfológica del componente asimétrico de la forma facial mediante técnicas multivariadas sobre los grupos que resultaran diferentes significativamente. Las mujeres que reportaron haber recibido algún tipo de tratamiento hormonal mostraron mayores valores de asimetría fluctuante facial respecto al grupo sin tratamiento. Esta asociación se mantuvo una vez removido el efecto de la ancestría genética y sin interactuar con el resto de variables incluidas en el análisis. Los cambios morfológicos asociados a este factor se concentran en el mentón, maxilar labio inferior, región perifrontal, región nasal y orejas. Algunos trabajos anteriores dieron cuenta de la posible relación entre la asimetría facial y los niveles de hormonas, pero no hay estudios que sustenten la relación causal o directa entre la asociación aquí planteada. El presente trabajo es una evidencia más de la asociación entre el consumo de hormonas y modificaciones de caracteres faciales en para poblaciones urbanas mestizas latinoamericanas.

Background and aims. Secreted frizzled-related protein 5 (SFRP5) was recently described as a new adipokine protective for hepatic steatosis and other obesity-related complications in the mouse model. To date, SFRP5 expression in non-alcoholic fatty liver disease (NAFLD) has not been fully assessed in humans. We measured circulating SFRP5 levels and its expression in liver and adipose tissue, and evaluated its association with NAFLD in morbidly obese women.Material and methods. Fifty-four morbidly obese women undergoing bariatric surgery were included in the study. Liver biopsies were used for histology and hepatic triglyceride content quantification. Circulating SFRP5 levels were measured through enzyme-linked immunoabsorbent assay, and SFRP5 expression was performed in hepatic and adipose tissue (subcutaneous and visceral).Results. Although circulating SFRP5 levels showed a tendency to decrease with NAFLD progression, no significant differences were observed among non-alcoholic steatosis, steatohepatitis, and control subjects. Hepatic SFRP5 expression showed a negative correlation with hepatic triglyceride content (r = -0.349, P = 0.016 for mRNA and r = -0.291, P = 0.040 for SRFP5 protein) and ALT serum levels (r = -0.437, P = 0.001 for SRFP5 protein). In addition, hepatic SFRP5 protein levels were significantly lower in NASH than in control subjects (P = 0.006). Conclusion. This is the first study reporting an association of hepatic SFRP5 expression with NAFLD in humans.

The advances in experimental methods and the development of high performance bioinformatic tools have substantially improved our understanding of microbial communities associated with human niches. Many studies have documented that changes in microbial abundance and composition of the human microbiome is associated with human health and diseased state. The majority of research on human microbiome is typically focused in the analysis of one level of biological information, i.e., metagenomics or metatranscriptomics. In this review, we describe some of the different experimental and bioinformatic strategies applied to analyze the 16S rRNA gene profiling and shotgun sequencing data of the human microbiome. We also discuss how some of the recent insights in the combination of metagenomics, metatranscriptomics and viromics can provide more detailed description on the interactions between microorganisms and viruses in oral and gut microbiomes. Recent studies on viromics have begun to gain importance due to the potential involvement of viruses in microbial dysbiosis. In addition, metatranscriptomic combined with metagenomic analysis have shown that a substantial fraction of microbial transcripts can be differentially regulated relative to their microbial genomic abundances. Thus, understanding the molecular interactions in the microbiome using the combination of metagenomics, metatranscriptomics and viromics is one of the main challenges towards a system level understanding of human microbiome.

Latin Americans are highly heterogeneous regarding the type of Native American ancestry. Consideration of specific associations with common diseases may lead to substantial advances in unraveling of disease etiology and disease prevention. Here we investigate possible associations between the type of Native American ancestry and leading causes of death. After an aggregate-data study based on genome-wide genotype data from 1805 admixed Chileans and 639,789 deaths, we validate an identified association with gallbladder cancer relying on individual data from 64 gallbladder cancer patients, with and without a family history, and 170 healthy controls. Native American proportions were markedly underestimated when the two main types of Native American ancestry in Chile, originated from the Mapuche and Aymara indigenous peoples, were combined together. Consideration of the type of Native American ancestry was crucial to identify disease associations. Native American ancestry showed no association with gallbladder cancer mortality (P = 0.26). By contrast, each 1% increase in the Mapuche proportion represented a 3.7% increased mortality risk by gallbladder cancer (95%CI 3.1-4.3%, P = 6×10-27). Individual-data results and extensive sensitivity analyses confirmed the association between Mapuche ancestry and gallbladder cancer. Increasing Mapuche proportions were also associated with an increased mortality due to asthma and, interestingly, with a decreased mortality by diabetes. The mortality due to skin, bladder, larynx, bronchus and lung cancers increased with increasing Aymara proportions. Described methods should be considered in future studies on human population genetics and human health. Complementary individual-based studies are needed to apportion the genetic and non-genetic components of associations identified relying on aggregate-data.

The expression of facial asymmetries has been recurrently related with poverty and/or disadvantaged socioeconomic status. Departing from the developmental instability theory, previous approaches attempted to test the statistical relationship between the stress experienced by individuals grown in poor conditions and an increase in facial and corporal asymmetry. Here we aim to further evaluate such hypothesis on a large sample of admixed Latin Americans individuals by exploring if low socioeconomic status individuals tend to exhibit greater facial fluctuating asymmetry values. To do so, we implement Procrustes analysis of variance and Hierarchical Linear Modelling (HLM) to estimate potential associations between facial fluctuating asymmetry values and socioeconomic status. We report significant relationships between facial fluctuating asymmetry values and age, sex, and genetic ancestry, while socioeconomic status failed to exhibit any strong statistical relationship with facial asymmetry. These results are persistent after the effect of heterozygosity (a proxy for genetic ancestry) is controlled in the model. Our results indicate that, at least on the studied sample, there is no relationship between socioeconomic stress (as intended as low socioeconomic status) and facial asymmetries.

[This corrects the article DOI: 10.1371/journal.pone.0169287.].

Risk of hyperuricemia is modified by genetic and environmental factors. Our aim was to identify factors associated with serum uric acid levels and hyperuricemia in Mexicans. A pilot Genome-wide association study GWAS was performed in a subgroup of participants (n = 411) from the Health Workers Cohort Study (HWCS). Single nucleotide polymorphisms (SNPs) associated with serum uric acid levels were validated in all the HWCS participants (n = 1939) and replicated in independent children (n = 1080) and adult (n = 1073) case-control studies. The meta-analysis of the whole HWCS and replication samples identified three SLC2A9 SNPs: rs1014290 (p = 2.3 × 10−64), rs3775948 (p = 8.2 × 10−64) and rs11722228 (p = 1.1 × 10−17); and an ABCG2 missense SNP, rs2231142 (p = 1.0 × 10−18). Among the non-genetic factors identified, the visceral adiposity index, smoking, the metabolic syndrome and its components (waist circumference, blood pressure, glucose and hyperlipidemia) were associated with increased serum uric acid levels and hyperuricemia (p < 0.05). Among the female HWCS participants, the odds ratio for hyperuricemia was 1.24 (95% CI, 1.01−1.53) per unit increase in soft drink consumption. As reported in other studies, our findings indicate that diet, adiposity and genetic variation contribute to the elevated prevalence of hyperuricemia in Mexico.

Genome-wide association studies (GWAS) have identified copy number variants (CNVs) associated with obesity in chromosomal regions 1p31.1, 10q11.22, 11q11, 16p12.3, and recently 1p21.1, which contains the salivary amylase gene (AMY1). Recent evidence suggests this enzyme may influence gut microbiota composition through carbohydrate (mainly starch) degradation. The role of these CNVs in obesity has been scarcely explored in the Latino population, and thus the aim of our study was to evaluate the association of 1p31.1, 10q11.22, 11q11, 16p12.3 and 1p21.1 CNVs with obesity in 921 Mexican children, to replicate significant associations in 920 Mexican adults, and to analyze the association of AMY1 copy number with gut microbiota in 75 children and 45 adults. Of the five CNVs analyzed, 1q11 CNV was significantly associated with obesity in children, but not in adults. Only AMY1 CNV was significantly associated with obesity in both age groups. Moreover, gut microbiota analyses revealed a positive correlation between AMY1 copy number and Prevotella abundance. This genus has enzymes and gene clusters essential for complex polysaccharide degradation and utilization. To our knowledge, this is the first study to analyze the association of these five CNVs in the Mexican population and to report a correlation between AMY1 CN and gut microbiota in humans.

The current genetic makeup of Latin America has been shaped by a history of extensive admixture between Africans, Europeans and Native Americans, a process taking place within the context of extensive geographic and social stratification. We estimated individual ancestry proportions in a sample of 7,342 subjects ascertained in five countries (Brazil, Chile, Colombia, México and Perú). These individuals were also characterized for a range of physical appearance traits and for self-perception of ancestry. The geographic distribution of admixture proportions in this sample reveals extensive population structure, illustrating the continuing impact of demographic history on the genetic diversity of Latin America. Significant ancestry effects were detected for most phenotypes studied. However, ancestry generally explains only a modest proportion of total phenotypic variation. Genetically estimated and self-perceived ancestry correlate significantly, but certain physical attributes have a strong impact on self-perception and bias self-perception of ancestry relative to genetically estimated ancestry.

The understanding of the complex genotype-phenotype architecture of human pigmentation has clear implications for the evolutionary history of humans, as well as for medical and forensic practices. Although dozens of genes have previously been associated with human skin color, knowledge about this trait remains incomplete. In particular, studies focusing on populations outside the European-North American axis are rare, and, until now, admixed populations have seldom been considered. The present study was designed to help fill this gap. Our objective was to evaluate possible associations of 18 single nucleotide polymorphisms (SNPs), located within nine genes, and one pseudogene with the Melanin Index (MI) in two admixed Brazilian populations (Gaucho, N = 352; Baiano, N = 148) with different histories of geographic and ethnic colonization. Of the total sample, four markers were found to be significantly associated with skin color, but only two (SLC24A5 rs1426654, and SLC45A2 rs16891982) were consistently associated with MI in both samples (Gaucho and Baiano). Therefore, only these 2 SNPs should be preliminarily considered to have forensic significance because they consistently showed the association independently of the admixture level of the populations studied. We do not discard that the other two markers (HERC2 rs1129038 and TYR rs1126809) might be also relevant to admixed samples, but additional studies are necessary to confirm the real importance of these markers for skin pigmentation. Finally, our study shows associations of some SNPs with MI in a modern Brazilian admixed sample, with possible applications in forensic genetics. Some classical genetic markers in Euro-North American populations are not associated with MI in our sample. Our results point out the relevance of considering population differences in selecting an appropriate set of SNPs as phenotype predictors in forensic practice.

BACKGROUND: Several studies have identified multiple obesity-associated loci mainly in European populations. However, their contribution to obesity in other ethnicities such as Mexicans is largely unknown. The aim of this study was to examine 26 obesity-associated single-nucleotide polymorphisms (SNP) in a sample of Mexican mestizos. METHODS: 9 SNPs in biological candidate genes showing replications (PPARG, ADRB3, ADRB2, LEPR, GNB3, UCP3, ADIPOQ, UCP2, and NR3C1), and 17 SNPs in or near genes associated with obesity in first, second and third wave GWAS (INSIG2, FTO, MC4R, TMEM18, FAIM2/BCDIN3, BDNF, SH2B1, GNPDA2, NEGR1, KCTD15, SEC16B/RASAL2, NPC1, SFRF10/ETV5, MAF, PRL, MTCH2, and PTER) were genotyped in 1,156 unrelated Mexican-Mestizos including 683 cases (441 obese class I/II and 242 obese class III) and 473 normal-weight controls. In a second stage we selected 12 of the SNPs showing nominal associations with obesity, to seek associations with quantitative obesity-related traits in 3 cohorts including 1,218 Mexican Mestizo children, 945 Mexican Mestizo adults, and 543 Indigenous Mexican adults. RESULTS: After adjusting for age, sex and admixture, significant associations with obesity were found for 6 genes in the case-control study (ADIPOQ, FTO, TMEM18, INSIG2, FAIM2/BCDIN3 and BDNF). In addition, SH2B1 was associated only with class I/II obesity and MC4R only with class III obesity. SNPs located at or near FAIM2/BCDIN3, TMEM18, INSIG2, GNPDA2 and SEC16B/RASAL2 were significantly associated with BMI and/or WC in the combined analysis of Mexican-mestizo children and adults, and FTO locus was significantly associated with increased BMI in Indigenous Mexican populations. CONCLUSIONS: Our findings replicate the association of 8 obesity-related SNPs with obesity risk in Mexican adults, and confirm the role of some of these SNPs in BMI in Mexican adults and children.

BackgroundVNN1 gene expression levels and the G-137T polymorphism have been associated with high density lipoprotein cholesterol (HDL-C) levels in Mexican American adults. We aim to evaluate the contribution of VNN1 gene expression and the G-137T variant to HDL-C levels and other metabolic traits in Mexican prepubertal children.Methodology/principal findingsVNN1 mRNA expression levels were quantified in peripheral blood leukocytes from 224 unrelated Mexican-Mestizo children aged 6-8 years (107 boys and 117 girls) and were genotyped for the G-137T variant (rs4897612). To account for population stratification, a panel of 10 ancestry informative markers was analyzed. After adjustment for admixture, the TT genotype was significantly associated with lower VNN1 mRNA expression levels (P = 2.9 × 10(-5)), decreased HDL-C levels (β = -6.19, P = 0.028) and with higher body mass index (BMI) z-score (β = 0.48, P = 0.024) in the total sample. In addition, VNN1 expression showed a positive correlation with HDL-C levels (r = 0.220; P = 0.017) and a negative correlation with BMI z-score (r = -0.225; P = 0.015) only in girls.Conclusion/significanceOur data suggest that VNN1 gene expression and the G-137T variant are associated with HDL-C levels in Mexican children, particularly in prepubertal girls.

BackgroundABCA1 genetic variation is known to play a role in HDL-C levels and various studies have also implicated ABCA1 variation in cardiovascular risk. The functional ABCA1/R230C variant is frequent in the Mexican population and has been consistently associated with low HDL-C concentrations. Although it has been associated with other cardiovascular risk factors such as obesity and type 2 diabetes mellitus, it is not known whether it is associated with coronary artery disease (CAD).AimThe purpose of the study was to analyze whether the ABCA1/R230C variant is associated with premature CAD in a case-control association study (GEA or Genetics of Atherosclerotic Disease), and to explore whether BMI modulates the effect of the C230 allele on other metabolic traits using a population-based design.ResultsThe C230 allele was significantly associated with both lower HDL-C levels and a lower risk of premature CAD as compared to controls (OR = 0.566; P(add) = 1.499×10(-5)). In addition, BMI modulated the effect of R230C on body fat distribution, as the correlation between BMI and visceral to subcutaneous adipose tissue (a metric of the propensity to store fat viscerally as compared to subcutaneously) was negative in RR homozygous individuals, but positive in premenopausal women bearing the C230 allele, with a statistically significant interaction (P = 0.005). BMI-R230C interaction was also significant for triglyceride levels in women regardless of their menopausal status (P = 0.036).ConclusionThis is the first study assessing the effect of the R230C/ABCA1 variant in remature CAD. C230 was associated with both decreased HDL-C levels and a lower risk of premature CAD, and gender-specific BMI-R230C interactions were observed for different metabolic traits. These interactions may help explain inconsistencies in associations, and underscore the need to further analyze interactions of this functional and frequent variant with diet, exercise and other environmental factors.

BACKGROUND: Common variants rs6232 and rs6235 in the PCSK1 gene have been associated with obesity in European populations. We aimed to evaluate the contribution of these variants to obesity and related traits in Mexican children and adults. METHODOLOGY/PRINCIPAL FINDINGS: Rs6232 and rs6235 were genotyped in 2382 individuals, 1206 children and 1176 adults. Minor allele frequencies were 0.78% for rs6232 and 19.99% for rs6235. Rs6232 was significantly associated with childhood obesity and adult class III obesity (OR = 3.01 95%CI 1.64-5.53; P = 4 × 10⁻⁴ in the combined analysis). In addition, this SNP was significantly associated with lower fasting glucose levels (P = 0.01) and with increased insulin levels and HOMA-B (P = 0.05 and 0.01, respectively) only in non-obese children. In contrast, rs6235 showed no significant association with obesity or with glucose homeostasis parameters in any group. CONCLUSION/SIGNIFICANCE: Although rs6232 is rare in the Mexican population, it should be considered as an important risk factor for extreme forms of obesity.

Culture and genetics rely on two distinct but not isolated transmission systems. Cultural processes may change the human selective environment and thereby affect which individuals survive and reproduce. Here, we evaluated whether the modes of subsistence in Native American populations and the frequencies of the ABCA1*Arg230Cys polymorphism were correlated. Further, we examined whether the evolutionary consequences of the agriculturally constructed niche in Mesoamerica could be considered as a gene-culture coevolution model. For this purpose, we genotyped 229 individuals affiliated with 19 Native American populations and added data for 41 other Native American groups (n = 1905) to the analysis. In combination with the SNP cluster of a neutral region, this dataset was then used to unravel the scenario involved in 230Cys evolutionary history. The estimated age of 230Cys is compatible with its origin occurring in the American continent. The correlation of its frequencies with the archeological data on Zea pollen in Mesoamerica/Central America, the neutral coalescent simulations, and the F(ST)-based natural selection analysis suggest that maize domestication was the driving force in the increase in the frequencies of 230Cys in this region. These results may represent the first example of a gene-culture coevolution involving an autochthonous American allele.