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.

Differential DNA Methylation from Autistic Children Enriches Evidence for Genes Associated with ASD and New Candidate Genes

The etiology of Autism Spectrum Disorders (ASD) is a result of the interaction between genes and the environment. The study of epigenetic factors that affect gene expression, such as DNA methylation, has become an important area of research in ASD. In recent years, there has been an increasing body of evidence pointing to epigenetic mechanisms that influence brain development, as in the case of ASD, when gene methylation dysregulation is present. Our analysis revealed 853 differentially methylated CpG in ASD patients, affecting 509 genes across the genome. Enrichment analysis showed five related diseases, including autistic disorder and mental disorders, which are particularly significant. In this work, we identified 64 genes that were previously reported in the SFARI gene database, classified according to their impact index. Additionally, we identified new genes that have not been previously reported as candidates with differences in the methylation patterns of Mexican children with ASD.

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.

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.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

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.

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.

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.

MicroRNA delivery systems in glioma therapy and perspectives: A systematic review

Gliomas are the deadliest of all primary brain tumors, and they constitute a serious global health problem. MicroRNAs (miRNAs) are gene expression regulators associated with glioma pathogenesis. Thus, miRNAs represent potential therapeutic agents for treating gliomas. However, miRNAs have not been established as part of the regular clinical armamentarium. This systemic review evaluates current molecular and pre-clinical studies with the aim of defining the most appealing supramolecular platform for administering therapeutic miRNA to patients with gliomas. An integrated analysis suggested that cationic lipid nanoparticles, functionalized with octa-arginine peptides, represent a potentially specific, practical, non-invasive intervention for treating gliomas. This supramolecular platform allows loading both hydrophilic (miRNA) and hydrophobic (anti-tumor drugs, like temozolomide) molecules. This systemic review is the first to describe miRNA delivery systems targeted to gliomas that integrate several types of molecules as active ingredients. Further experimental validation is warranted to confirm the practical value of miRNA delivery systems.

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.

Dysferlinopathy misdiagnosed with juvenile polymyositis in the pre-symptomatic stage of hyperCKemia: a case report and literature review

Background

Dysferlinopathy encompasses a group of rare muscular dystrophies caused by recessive mutations in the DYSF gene. The phenotype ranges from asymptomatic elevated serum creatine kinase (hyperCKemia) to selective and progressive involvement of the proximal and/or distal muscles of the limbs. Bohan and Peter criteria are the most widely used for the diagnosis of polymyositis, but they have limitations and can misclassify muscular dystrophies with inflammation as polymyositis. Most dysferlinopathy patients have muscle biopsies with inflammation and thus are vulnerable to misdiagnosis with polymyositis and inappropriate treatment with steroids and immunosuppressors.

Case presentation

We describe a 14 years-old male patient who was referred for assessment of asymptomatic hyperCKemia (26,372 IU/L). An X-linked dystrophinopathy initially was ruled out by direct genetic testing. Juvenile polymyositis was considered based on muscle biopsy, creatine kinase levels, and electromyography changes. Corticosteroid treatment triggered proximal lower limb muscular weakness, and no full muscular strength recovery was observed after corticosteroid withdrawal. Based on these observations, a limb-girdle muscular dystrophy (LGMD) was suspected, and LGMDR2 was confirmed by whole exome sequencing.

Conclusion

We report a dysferlinopathy patient who was misdiagnosed with juvenile polymyositis and explore in a literature review how common such misdiagnoses are. With diagnosis based only on routine clinicopathological examinations, distinguishing an inflammatory myopathy from dysferlinopathy is quite difficult. We suggest that before establishing a diagnosis of “definite” or “probable” juvenile polymyositis, according to Bohan and Peter or current ACR/EULAR criteria, a muscular dystrophy must first be ruled out.

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.

Rare coding variants in 35 genes associate with circulating lipid levels-A multi-ancestry analysis of 170,000 exomes

Large-scale gene sequencing studies for complex traits have the potential to identify causal genes with therapeutic implications. We performed gene-based association testing of blood lipid levels with rare (minor allele frequency < 1%) predicted damaging coding variation by using sequence data from >170,000 individuals from multiple ancestries: 97,493 European, 30,025 South Asian, 16,507 African, 16,440 Hispanic/Latino, 10,420 East Asian, and 1,182 Samoan. We identified 35 genes associated with circulating lipid levels; some of these genes have not been previously associated with lipid levels when using rare coding variation from population-based samples. We prioritize 32 genes in array-based genome-wide association study (GWAS) loci based on aggregations of rare coding variants; three (EVI5, SH2B3, and PLIN1) had no prior association of rare coding variants with lipid levels. Most of our associated genes showed evidence of association among multiple ancestries. Finally, we observed an enrichment of gene-based associations for low-density lipoprotein cholesterol drug target genes and for genes closest to GWAS index single-nucleotide polymorphisms (SNPs). Our results demonstrate that gene-based associations can be beneficial for drug target development and provide evidence that the gene closest to the array-based GWAS index SNP is often the functional gene for blood lipid levels.

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13023-021-01924-z.

Determinants of penetrance and variable expressivity in monogenic metabolic conditions across 77,184 exomes

Hundreds of thousands of genetic variants have been reported to cause severe monogenic diseases, but the probability that a variant carrier develops the disease (termed penetrance) is unknown for virtually all of them. Additionally, the clinical utility of common polygenetic variation remains uncertain. Using exome sequencing from 77,184 adult individuals (38,618 multi-ancestral individuals from a type 2 diabetes case-control study and 38,566 participants from the UK Biobank, for whom genotype array data were also available), we apply clinical standard-of-care gene variant curation for eight monogenic metabolic conditions. Rare variants causing monogenic diabetes and dyslipidemias display effect sizes significantly larger than the top 1% of the corresponding polygenic scores. Nevertheless, penetrance estimates for monogenic variant carriers average 60% or lower for most conditions. We assess epidemiologic and genetic factors contributing to risk prediction in monogenic variant carriers, demonstrating that inclusion of polygenic variation significantly improves biomarker estimation for two monogenic dyslipidemias.

Reconstruction of ancient microbial genomes from the human gut

Loss of gut microbial diversity^1–6 in industrial populations is associated with chronic diseases⁷, 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.

Ancient microbiomes from palaeofaeces are more similar to non-industrialized than industrialized human gut microbiomes regardless of geography, but 39% of their de novo reconstructed genomes represent previously undescribed microbial species.

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

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

MON-708 Characterizing DNA Methylation Signatures in Adipose Tissue from Metabolic Impaired Asymptomatic Individuals

Obesity remains as a global epidemic characterized by progressive metabolic dysregulation in glucose homeostasis. Along with a genetic association in the development of T2D, epigenetic regulation has been suggested as a significant contributor in altered gene expression. Recent studies have described DNA methylation changes in insulin-sensitive tissues involved in T2D pathogenesis, however epigenetic dynamics on early stages to metabolic alterations is still unclear.

We investigated potential DNA methylation signatures in 34 asymptomatic individuals from the GEMM family study. We compared differentially methylated CpG sites (DMC: B value>0 and delta Beta >|10%|; Infinium EPIC array) from subcutaneous adipose tissue (SCAT) in different groups of individuals according to BMI (kg/m²) and HbA1c (%) levels as follow: Group A Control (C): n=9, 22.0±1.9 kg/m², 4.8±0.3%; Group B Overweight (OW) with normal HbA1c: n=6, 27.8±1.6 kg/m², 5±0.2%; Group C Obese (OB) with normal HbA1c: n=6, 34.6±4.2 kg/m², 5.2±0.2%; Group D Prediabetes (PD): n=7, 31.1±5.7 kg/m², 5.9±0.2% and Group E T2D: n=6, 30.6±7.3 kg/m², 7.2±0.9%.

We found 43 overlapping genes with shared pathways in all groups, mainly those related to metabolism and adipogenesis. We also documented particular altered methylated genes, in each group (OW: 386, OB:1005, PD:76 and T2D:189). Pathway enrichment analysis in OB and T2D was mainly related to glucose metabolism, while in OW and PD was NOTCH signaling. All groups displayed a consistent hypermethylation in RARA, ESR1 and NCOR2, well known genes involved in lipid metabolism. Additionally, we describe for the first time, a progression toward hypomethylation in ARHGAP15 and MTAP, related with an impaired metabolic status. Otherwise, analysis of overlapping CpG sites revealed a consistently hypermethylated state in OW (86.42%), OB (86.48%) and PD (51.72%), in contrast with the hypomethylation state (56.3%) observed in the T2D group, previously observed elsewhere (1).

In conclusion, comparison of methylation in SCAT obtained from OW, OB, PD and T2D individuals, display potential pathways and DMC signatures specific in each group. Common novel overlapping genes in global DNA methylation profiles of SCAT, were also observed.

Reference: (1) Barajas-Olmos et al., BMC Med Genet. 2018 Feb 21;19(1):1-8.

Nothing to Disclose: FE, FB, AM, EH, GEMM, ER, RB, LO.

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.

Variation in Actionable Pharmacogenetic Markers in Natives and Mestizos From Mexico

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

Exome sequencing of 20,791 cases of type 2 diabetes and 24,440 controls

Protein-coding genetic variants that strongly affect disease risk can yield relevant clues to disease pathogenesis. Here we report exome-sequencing analyses of 20,791 individuals with type 2 diabetes (T2D) and 24,440 non-diabetic control participants from 5 ancestries. We identify gene-level associations of rare variants (with minor allele frequencies of less than 0.5%) in 4 genes at exome-wide significance, including a series of more than 30 SLC30A8 alleles that conveys protection against T2D, and in 12 gene sets, including those corresponding to T2D drug targets (P = 6.1 × 10-3) and candidate genes from knockout mice (P = 5.2 × 10-3). Within our study, the strongest T2D gene-level signals for rare variants explain at most 25% of the heritability of the strongest common single-variant signals, and the gene-level effect sizes of the rare variants that we observed in established T2D drug targets will require 75,000-185,000 sequenced cases to achieve exome-wide significance. We propose a method to interpret these modest rare-variant associations and to incorporate these associations into future target or gene prioritization efforts.

Next-generation sequencing for identifying a novel/de novo pathogenic variant in a Mexican patient with cystic fibrosis: a case report

Background

Mexico is among the countries showing the highest heterogeneity of CFTR variants. However, no de novo variants have previously been reported in Mexican patients with cystic fibrosis (CF).

Case presentation

Here, we report the first case of a novel/de novo variant in a Mexican patient with CF. Our patient was an 8-year-old male who had exhibited the clinical onset of CF at one month of age, with steatorrhea, malabsorption, poor weight gain, anemia, and recurrent respiratory tract infections. Complete sequencing of the CFTR gene by next generation sequencing (NGS) revealed two different variants in trans, including the previously reported CF-causing variant c.3266G > A (p.Trp1089*, W1089*), that was inherited from the mother, and the novel/de novo CFTR variant c.1762G > T (p.Glu588*).

Conclusion

Our results demonstrate the efficiency of targeted NGS for making a rapid and precise diagnosis in patients with clinically suspected CF. This method can enable the provision of accurate genetic counselling, and improve our understanding of the molecular basis of genetic diseases.

Comparing signals of natural selection between three Indigenous North American populations

While many studies have highlighted human adaptations to diverse environments worldwide, genomic studies of natural selection in Indigenous populations in the Americas have been absent from this literature until very recently. Since humans first entered the Americas some 20,000 years ago, they have settled in many new environments across the continent. This diversity of environments has placed variable selective pressures on the populations living in each region, but the effects of these pressures have not been extensively studied to date. To help fill this gap, we collected genome-wide data from three Indigenous North American populations from different geographic regions of the continent (Alaska, southeastern United States, and central Mexico). We identified signals of natural selection in each population and compared signals across populations to explore the differences in selective pressures among the three regions sampled. We find evidence of adaptation to cold and high-latitude environments in Alaska, while in the southeastern United States and central Mexico, pathogenic environments seem to have created important selective pressures. This study lays the foundation for additional functional and phenotypic work on possible adaptations to varied environments during the history of population diversification in the Americas.

A Loss-of-Function Splice Acceptor Variant in IGF2 Is Protective for Type 2 Diabetes

Type 2 diabetes (T2D) affects more than 415 million people worldwide, and its costs to the health care system continue to rise. To identify common or rare genetic variation with potential therapeutic implications for T2D, we analyzed and replicated genome-wide protein coding variation in a total of 8,227 individuals with T2D and 12,966 individuals without T2D of Latino descent. We identified a novel genetic variant in the IGF2 gene associated with ∼20% reduced risk for T2D. This variant, which has an allele frequency of 17% in the Mexican population but is rare in Europe, prevents splicing between IGF2 exons 1 and 2. We show in vitro and in human liver and adipose tissue that the variant is associated with a specific, allele-dosage-dependent reduction in the expression of IGF2 isoform 2. In individuals who do not carry the protective allele, expression of IGF2 isoform 2 in adipose is positively correlated with both incidence of T2D and increased plasma glycated hemoglobin in individuals without T2D, providing support that the protective effects are mediated by reductions in IGF2 isoform 2. Broad phenotypic examination of carriers of the protective variant revealed no association with other disease states or impaired reproductive health. These findings suggest that reducing IGF2 isoform 2 expression in relevant tissues has potential as a new therapeutic strategy for T2D, even beyond the Latin American population, with no major adverse effects on health or reproduction.

Heterogenous Distribution of MTHFR Gene Variants among Mestizos and Diverse Amerindian Groups from Mexico

Methylenetetrahydrofolate reductase (MTHFR) is a key enzyme in folate metabolism. Folate deficiency has been related to several conditions, including neural tube defects (NTDs) and cardiovascular diseases. Hence, MTHFR genetic variants have been studied worldwide, particularly the C677T and A1298C. We genotyped the C677T and A1298C MTHFR polymorphisms in Mexican Amerindians (MAs), from the largest sample included in a genetic study (n = 2026, from 62 ethnic groups), and in a geographically-matched Mexican Mestizo population (MEZ, n = 638). The 677T allele was most frequent in Mexican individuals, particularly in MAs. The frequency of this allele in both MAs and MEZs was clearly enriched in the South region of the country, followed by the Central East and South East regions. In contrast, the frequency of the 1298C risk allele in Mexicans was one of the lowest in the world. Both in MAs and MEZs the variants 677T and 1298C displayed opposite allele frequency gradients from southern to northern Mexico. Our findings suggest that in Mestizos the 677T allele was derived from Amerindians while the 1298C allele was a European contribution. Some subgroups showed an allele frequency distribution that highlighted their genetic diversity. Notably, the distribution of the frequency of the 677T allele was consistent with that of the high incidence of NTDs reported in MEZ.

Transcriptome landmarks of the functional maturity of rat beta-cells, from lactation to adulthood

Research on the postnatal development of pancreatic beta-cells has become an important subject in recent years. Understanding the mechanisms that govern beta-cell postnatal maturation could bring new opportunities to therapeutic approaches for diabetes. The weaning period consists of a critical postnatal window for structural and physiologic maturation of rat beta-cells. To investigate transcriptome changes involved in the maturation of beta-cells neighboring this period, we performed microarray analysis in fluorescence-activated cell-sorted (FACS) beta-cell-enriched populations. Our results showed a variety of gene sets including those involved in the integration of metabolism, modulation of electrical activity, and regulation of the cell cycle that play important roles in the maturation process. These observations were validated using reverse hemolytic plaque assay, electrophysiological recordings, and flow cytometry analysis. Moreover, we suggest some unexplored pathways such as sphingolipid metabolism, insulin-vesicle trafficking, regulation of transcription/transduction by miRNA-30, trafficking proteins, and cell cycle proteins that could play important roles in the process mentioned above for further investigation.