Lorena Orozco

Analysis of protein-coding genetic variation in 60,706 humans

Large-scale reference data sets of human genetic variation are critical for the medical and functional interpretation of DNA sequence changes. Here we describe the aggregation and analysis of high-quality exome (protein-coding region) DNA sequence data for 60,706 individuals of diverse ancestries generated as part of the Exome Aggregation Consortium (ExAC). This catalogue of human genetic diversity contains an average of one variant every eight bases of the exome, and provides direct evidence for the presence of widespread mutational recurrence. We have used this catalogue to calculate objective metrics of pathogenicity for sequence variants, and to identify genes subject to strong selection against various classes of mutation; identifying 3,230 genes with near-complete depletion of predicted protein-truncating variants, with 72% of these genes having no currently established human disease phenotype. Finally, we demonstrate that these data can be used for the efficient filtering of candidate disease-causing variants, and for the discovery of human ‘knockout’ variants in protein-coding genes.

STAT4 associates with systemic lupus erythematosus through two independent effects that correlate with gene expression and act additively with IRF5 to increase risk

Objectives: To confirm and define the genetic association of STAT4 and systemic lupus erythematosus (SLE), investigate the possibility of correlations with differential splicing and/or expression levels, and genetic interaction with IRF5. Methods: 30 tag SNPs were genotyped in an independent set of Spanish cases and controls. SNPs surviving correction for multiple tests were genotyped in five new sets of cases and controls for replication. STAT4 cDNA was analysed by 5′-RACE PCR and sequencing. Expression levels were measured by quantitative PCR. Results: In the fine mapping, four SNPs were significant after correction for multiple testing, with rs3821236 and rs3024866 as the strongest signals, followed by the previously associated rs7574865, and by rs1467199. Association was replicated in all cohorts. After conditional regression analyses, two major independent signals, represented by SNPs rs3821236 and rs7574865, remained significant across the sets. These SNPs belong to separate haplotype blocks. High levels of STAT4 expression correlated with SNPs rs3821236, rs3024866 (both in the same haplotype block) and rs7574865 but not with other SNPs. Transcription of alternative tissue-specific exons 1, indicating the presence of tissue-specific promoters of potential importance in the expression of STAT4, was also detected. No interaction with associated SNPs of IRF5 was observed using regression analysis. Conclusions: These data confirm STAT4 as a susceptibility gene for SLE and suggest the presence of at least two functional variants affecting levels of STAT4. The results also indicate that the genes STAT4 and IRF5 act additively to increase the risk for SLE.

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.

Tumor necrosis factor-α is a common genetic risk factor for asthma, juvenile rheumatoid arthritis, and systemic lupus erythematosus in a Mexican pediatric population.

There is a great deal of evidence that points to the association of the tumor necrosis factor-alpha (TNF-alpha) gene as a common genetic factor in the pathogenesis of diseases that are caused by inflammatory and/or autoimmune etiologies. Two single nucleotide polymorphisms (SNPs) identified in the TNF-alpha promoter region have been associated with disease susceptibility and severity. We investigated whether -308G/A and -238G/A TNF-alpha polymorphisms were associated with asthma, systemic lupus erythematosus (SLE), and juvenile rheumatoid arthritis (JRA) in a pediatric Mexican population. In a case-control study of 725 patients (asthma: 226, JRA: 171, and SLE: 328) and 400 control subjects, the participants were analyzed using the allelic discrimination technique. The genotype distribution of both TNF-alpha polymorphisms was in Hardy-Weinberg equilibrium in each group. However, there were significant differences in the allele frequency of TNF-alpha-308A between the patients and the healthy controls. This allele was detected in 2.9% of the controls, 6.0% of asthmatic and JRA patients (p = 0.002 and p = 0.0086), and 6.7% of SLE patients (p = 0.00049); statistical significance was maintained after ancestry stratification (asthma: p = 0.0143, JRA: p = 0.0083, and SLE: p = 0.0026). Stratification by gender showed that the risk for the -308A allele in asthma and JRA was greater in females (OR = 4.16, p = 0.0008 and OR = 4.4, p = 0.0002, respectively). The TNF-alpha -238A allele showed an association only with JRA in males (OR = 2.89, p = 0.004). These results support the concept that the TNF-alpha gene is a genetic risk factor for asthma, SLE, and JRA in the pediatric Mexican population.

Association of TLR7 copy number variation with susceptibility to childhood-onset systemic lupus erythematosus in Mexican population

Objective Variations in gene copy number (CNV) have been recognised as a hereditable source of susceptibility in human complex diseases. Recent studies have shown that Tlr7 gene dosage has a significant contribution in the autoimmune-enhancing effect in mouse models of systemic lupus erythematosus (SLE). A study was therefore performed to investigate whether CNVs in TLR7 contribute to the genetic component of childhood-onset SLE. Methods A case–control association study was performed in 328 Mexican children with SLE and 403 healthy controls. Determination of CNVs of TLR7 was achieved by real-time PCR using the ΔΔCt method. Expression levels of TLR7 and interferon α (IFNα) were determined in 23 patients. In addition, a stratification analysis was performed to investigate the association of TLR7 gene copy number (CN) with lupus nephritis. Results A significant increase was found in the relative TLR7 gene CN in females patients with SLE compared with female controls (p<0.0001). However, logistic regression analysis by gender showed a higher OR (OR 6.61, p=0.005) in male patients with >1 copy of TLR7 than in female patients with >2 copies (OR 3.07, p<0.0001). This association was not observed with lupus nephritis. TLR7 mRNA levels correlated significantly with TLR7 CN and with IFNα mRNA levels. Conclusion These results show that an increase in TLR7 CN is a risk factor for childhood-onset SLE and provide new evidence for a role for X-linked gene dosage in SLE susceptibility. There is also evidence to suggest that TLR7 may be involved in the pathogenesis of SLE through the induction of IFNα.

Impact of genetic ancestry and sociodemographic status on the clinical expression of systemic lupus erythematosus in American Indian-European populations.

OBJECTIVE American Indian-Europeans, Asians, and African Americans have an excess morbidity from systemic lupus erythematosus (SLE) and a higher prevalence of lupus nephritis than do Caucasians. The aim of this study was to analyze the relationship between genetic ancestry and sociodemographic characteristics and clinical features in a large cohort of American Indian-European SLE patients. METHODS A total of 2,116 SLE patients of American Indian-European origin and 4,001 SLE patients of European descent for whom we had clinical data were included in the study. Genotyping of 253 continental ancestry-informative markers was performed on the Illumina platform. Structure and Admixture software were used to determine genetic ancestry proportions of each individual. Logistic regression was used to test the association between genetic ancestry and sociodemographic and clinical characteristics. Odds ratios (ORs) were calculated with 95% confidence intervals (95% CIs). RESULTS The average American Indian genetic ancestry of 2,116 SLE patients was 40.7%. American Indian genetic ancestry conferred increased risks of renal involvement (P < 0.0001, OR 3.50 [95% CI 2.63- 4.63]) and early age at onset (P < 0.0001). American Indian ancestry protected against photosensitivity (P < 0.0001, OR 0.58 [95% CI 0.44-0.76]), oral ulcers (P < 0.0001, OR 0.55 [95% CI 0.42-0.72]), and serositis (P < 0.0001, OR 0.56 [95% CI 0.41-0.75]) after adjustment for age, sex, and age at onset. However, age and sex had stronger effects than genetic ancestry on malar rash, discoid rash, arthritis, and neurologic involvement. CONCLUSION In general, American Indian genetic ancestry correlates with lower sociodemographic status and increases the risk of developing renal involvement and SLE at an earlier age.

Association of PDCD1 polymorphisms with childhood-onset systemic lupus erythematosus

A regulatory single nucleotide polymorphism (SNP) PD1.3G/A located on programmed cell death 1 (PDCD1) gene, was shown to be involved in susceptibility to systemic lupus erythematosus (SLE) in Swedish, European American, and Mexican cases. However, association to childhood-onset SLE has not been analyzed. The aim of this study was to investigate the association of PDCD1 polymorphisms and haplotypes with susceptibility to childhood-onset SLE in Mexican population. Three PDCD1 SNPs, PD1.3G/A, PD1.5C/T, PD1.6G/A, were analyzed in 250 childhood-onset SLE Mexican patients and 355 healthy controls in a case–control association study. Polymorphisms were genotyped by TaqMan technology. Stratification analysis was performed on the SLE cohort to investigate the SNP association with renal disorder. In addition, haplotypes were constructed with these three SNPs. The PD1.3A allele was significantly associated to childhood-onset SLE (P=0.0019, odds ratio (OR) 2.73, 95% confidence interval (95% CI) 1.35–5.56). The other PDCD1 SNPs did not show association. A total of 155 patients (62%) had nephritis, and no association was observed with PDCD1 SNPs. The ACG haplotype (PD1.3A, PD1.5C, PD1.6G) included almost all PD1.3A alleles, and it was more frequent in SLE patients (5.5%) than in controls (2.1%) (P=0.003; OR 2.73, 95% CI 1.37–5.46). The haplotype structure in Mexican controls was significantly different from those reported in Spanish and Swedish. Our results support association of the PD1.3A SNP to susceptibility of childhood-onset SLE in Mexican population and does not show association to lupus nephritis in this age group.

Genome-Wide Association Study in an Amerindian Ancestry Population Reveals Novel Systemic Lupus Erythematosus Risk Loci and the Role of European Admixture.

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease with a strong genetic component. We undertook the present work to perform the first genome‐wide association study on individuals from the Americas who are enriched for Native American heritage.

Spectrum of CFTR mutations in Mexican cystic fibrosis patients: identification of five novel mutations (W1098C, 846delT, P750L, 4160insGGGG and 297-1G-->A).

Abstract. We have analyzed 97 CF unrelated Mexican families for mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) gene. Our initial screening for 12 selected CFTR mutations led to mutation detection in 56.66% of the tested chromosomes. In patients with at least one unknown mutation after preliminary screening, an extensive analysis of the CFTR gene by single stranded conformation polymorphism (SSCP) or by multiplex heteroduplex (mHET) analysis was performed. A total of 34 different mutations representing 74.58% of the CF chromosomes were identified, including five novel CFTR mutations: W1098C, P750L, 846delT, 4160insGGGG and 297–1G→A. The level of detection of the CF mutations in Mexico is still lower than that observed in other populations with a relatively low frequency of the ΔF508 mutation, mainly from southern Europe. The CFTR gene analysis described here clearly demonstrated the high heterogeneity of our CF population, which could be explained by the complex ethnic composition of the Mexican population, in particular by the strong impact of the genetic pool from southern European countries.

Genetic Variants Associated With Quantitative Glucose Homeostasis Traits Translate to Type 2 Diabetes in Mexican Americans: The GUARDIAN (Genetics Underlying Diabetes in Hispanics) Consortium

Insulin sensitivity, insulin secretion, insulin clearance, and glucose effectiveness exhibit strong genetic components, although few studies have examined their genetic architecture or influence on type 2 diabetes (T2D) risk. We hypothesized that loci affecting variation in these quantitative traits influence T2D. We completed a multicohort genome-wide association study to search for loci influencing T2D-related quantitative traits in 4,176 Mexican Americans. Quantitative traits were measured by the frequently sampled intravenous glucose tolerance test (four cohorts) or euglycemic clamp (three cohorts), and random-effects models were used to test the association between loci and quantitative traits, adjusting for age, sex, and admixture proportions (Discovery). Analysis revealed a significant (P < 5.00 × 10−8) association at 11q14.3 (MTNR1B) with acute insulin response. Loci with P < 0.0001 among the quantitative traits were examined for translation to T2D risk in 6,463 T2D case and 9,232 control subjects of Mexican ancestry (Translation). Nonparametric meta-analysis of the Discovery and Translation cohorts identified significant associations at 6p24 (SLC35B3/TFAP2A) with glucose effectiveness/T2D, 11p15 (KCNQ1) with disposition index/T2D, and 6p22 (CDKAL1) and 11q14 (MTNR1B) with acute insulin response/T2D. These results suggest that T2D and insulin secretion and sensitivity have both shared and distinct genetic factors, potentially delineating genomic components of these quantitative traits that drive the risk for T2D.