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

Background

Molecular diagnosis of cystic fibrosis (CF) is challenging in Mexico due to the population's high genetic heterogeneity. To date, 46 pathogenic variants (PVs) have been reported, yielding a detection rate of 77%. We updated the spectrum and frequency of PVs responsible for this disease in mexican patients.

Methods

We extracted genomic DNA from peripheral blood lymphocytes obtained from 297 CF patients and their parents. First, we analyzed the five most frequent PVs in the Mexican population using PCR-mediated site-directed mutagenesis. In patients with at least one identified allele, CFTR sequencing was performed using next-generation sequencing tools and multiplex ligation-dependent probe amplification. For variants not previously classified as pathogenic, we used a combination of in silico prediction, CFTR modeling, and clinical characteristics to determine a genotype-phenotype correlation.

Results

We identified 95 PVs, increasing the detection rate to 87.04%. The most frequent variants were p.(PheF508del) (42.7%), followed by p.(Gly542*) (5.6%), p.(Ser945Leu) (2.9%), p.(Trp1204*) and p.(Ser549Asn) (2.5%), and CFTRdel25-26 and p.(Asn386Ilefs*3) (2.3%). The remaining variants had frequencies of <2.0%, and some were exclusive to one family. We identified 10 novel PVs localized in different exons (frequency range: 0.1-0.8%), all of which produced structural changes, deletions, or duplications in different domains of the protein, resulting in dysfunctional ion flow. The use of different in silico software and American College of Medical Genetics and Genomics (ACMG) and the Association for Molecular Pathology (AMP) criteria allowed us to assume that all of these PVs were pathogenic, causing a severe phenotype.

Conclusions

In a highly heterogeneous population, combinations of different tools are needed to identify the variants responsible for CF and enable the establishment of appropriate strategies for CF diagnosis, prevention, and treatment.

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

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

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.

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.

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.

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.

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.

Mexican Carriers of the HNF1A p.E508K Variant Do Not Experience an Enhanced Response to Sulfonylureas

OBJECTIVE

To assess whether an ethnic-specific variant (p.E508K) in the maturity-onset diabetes of the young (MODY) gene hepatocyte nuclear factor-1α (HNF1A) found in Mexicans is associated with higher sensitivity to sulfonylureas, as documented in patients with MODY3.

RESEARCH DESIGN AND METHODS

We recruited 96 participants (46 variant carriers and 50 age- and sex-matched noncarriers). Response to glipizide (one 2.5-5.0-mg dose), metformin (four 500-mg doses), and an oral glucose challenge was evaluated using a previously validated protocol. Glucose and insulin levels and their areas under the curve (AUCs) were compared between groups.

RESULTS

Carriers of the p.E508K variant had a lower maximum insulin peak during the glipizide challenge as compared with noncarriers with diabetes (P < 0.05). Also, carriers had a lower insulin response after the oral glucose challenge. Following an oral glucose tolerance test in the presence of metformin, carriers of the p.E508K variant with diabetes had a lower maximum insulin peak and total and incremental insulin AUC value as compared with noncarriers with diabetes (P < 0.05). A similar but nonsignificant trend was seen in participants without type 2 diabetes.

CONCLUSIONS

Carriers of variant p.E508K in HNF1A have a reduced insulin response rather than the increased sensitivity to sulfonylureas seen in patients with MODY3.

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.

Type 2 Diabetes Variants Disrupt Function of SLC16A11 through Two Distinct Mechanisms

Type 2 diabetes (T2D) affects Latinos at twice the rate seen in populations of European descent. We recently identified a risk haplotype spanning SLC16A11 that explains ∼20% of the increased T2D prevalence in Mexico. Here, through genetic fine-mapping, we define a set of tightly linked variants likely to contain the causal allele(s). We show that variants on the T2D-associated haplotype have two distinct effects: (1) decreasing SLC16A11 expression in liver and (2) disrupting a key interaction with basigin, thereby reducing cell-surface localization. Both independent mechanisms reduce SLC16A11 function and suggest SLC16A11 is the causal gene at this locus. To gain insight into how SLC16A11 disruption impacts T2D risk, we demonstrate that SLC16A11 is a proton-coupled monocarboxylate transporter and that genetic perturbation of SLC16A11 induces changes in fatty acid and lipid metabolism that are associated with increased T2D risk. Our findings suggest that increasing SLC16A11 function could be therapeutically beneficial for T2D. VIDEO ABSTRACT.

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.

Association of a low-frequency variant in HNF1A with type 2 diabetes in a Latino population

IMPORTANCE

Latino populations have one of the highest prevalences of type 2 diabetes worldwide.

OBJECTIVES

To investigate the association between rare protein-coding genetic variants and prevalence of type 2 diabetes in a large Latino population and to explore potential molecular and physiological mechanisms for the observed relationships.

DESIGN, SETTING, AND PARTICIPANTS

Whole-exome sequencing was performed on DNA samples from 3756 Mexican and US Latino individuals (1794 with type 2 diabetes and 1962 without diabetes) recruited from 1993 to 2013. One variant was further tested for allele frequency and association with type 2 diabetes in large multiethnic data sets of 14,276 participants and characterized in experimental assays.

MAIN OUTCOME AND MEASURES

Prevalence of type 2 diabetes. Secondary outcomes included age of onset, body mass index, and effect on protein function.

RESULTS

A single rare missense variant (c.1522G>A [p.E508K]) was associated with type 2 diabetes prevalence (odds ratio [OR], 5.48; 95% CI, 2.83-10.61; P = 4.4 × 10(-7)) in hepatocyte nuclear factor 1-α (HNF1A), the gene responsible for maturity onset diabetes of the young type 3 (MODY3). This variant was observed in 0.36% of participants without type 2 diabetes and 2.1% of participants with it. In multiethnic replication data sets, the p.E508K variant was seen only in Latino patients (n = 1443 with type 2 diabetes and 1673 without it) and was associated with type 2 diabetes (OR, 4.16; 95% CI, 1.75-9.92; P = .0013). In experimental assays, HNF-1A protein encoding the p.E508K mutant demonstrated reduced transactivation activity of its target promoter compared with a wild-type protein. In our data, carriers and noncarriers of the p.E508K mutation with type 2 diabetes had no significant differences in compared clinical characteristics, including age at onset. The mean (SD) age for carriers was 45.3 years (11.2) vs 47.5 years (11.5) for noncarriers (P = .49) and the mean (SD) BMI for carriers was 28.2 (5.5) vs 29.3 (5.3) for noncarriers (P = .19).

CONCLUSIONS AND RELEVANCE

Using whole-exome sequencing, we identified a single low-frequency variant in the MODY3-causing gene HNF1A that is associated with type 2 diabetes in Latino populations and may affect protein function. This finding may have implications for screening and therapeutic modification in this population, but additional studies are required.

Knockdown of dystrophin Dp71 impairs PC12 cells cycle: localization in the spindle and cytokinesis structures implies a role for Dp71 in cell division

The function of dystrophin Dp71 in neuronal cells remains to be established. Previously, we revealed the involvement of this protein in both nerve growth factor (NGF)-induced neuronal differentiation and cell adhesion by isolation and characterization of PC12 neuronal cells with depleted levels of Dp71. In this work, a novel phenotype of Dp71-knockdown cells was characterized, which is their delayed growth rate. Cell cycle analyses revealed an altered behavior of Dp71-depleted cells, which consists of a delay in G0/G1 transition and an increase in apoptosis during nocodazole-induced mitotic arrest. Dp71 associates with lamin B1 and β-dystroglycan, proteins involved in aspects of the cell division cycle; therefore, we compared the distribution of Dp71 with that of lamin B1 and β-dystroglycan in PC12 cells at mitosis and cytokinesis by means of immunofluorescence and confocal microscopy analysis. All of these three proteins exhibited a similar immunostaining pattern, localized at mitotic spindle, cleavage furrow, and midbody. It is noteworthy that a drastic decreased staining in mitotic spindle, cleavage furrow, and midbody was observed for both lamin B1 and β-dystroglycan in Dp71-depleted cells. Furthermore, we demonstrated the interaction of Dp71 with lamin B1 in PC12 cells by immunoprecipitation and pull-down assays, and importantly, we revealed that knockdown of Dp71 expression caused a marked reduction in lamin B1 levels and altered localization of the nuclear envelope protein emerin. Our data indicate that Dp71 is a component of the mitotic spindle and cytokinesis multi-protein apparatuses that might modulate the cell division cycle by affecting lamin B1 and β-dystroglycan levels.