Multiple sclerosis risk variant HLA-DRB1*1501 associates with high expression of DRB1 gene in different human populations.

The human leukocyte antigen (HLA) DRB1*1501 has been consistently associated with multiple sclerosis (MS) in nearly all populations tested. This points to a specific antigen presentation as the pathogenic mechanism though this does not fully explain the disease association. The identification of expression quantitative trait loci (eQTL) for genes in the HLA locus poses the question of the role of gene expression in MS susceptibility. We analyzed the eQTLs in the HLA region with respect to MS-associated HLA-variants obtained from genome-wide association studies (GWAS). We found that the Tag of DRB1*1501, rs3135388 A allele, correlated with high expression of DRB1, DRB5 and DQB1 genes in a Caucasian population. In quantitative terms, the MS-risk AA genotype carriers of rs3135388 were associated with 15.7-, 5.2- and 8.3-fold higher expression of DQB1, DRB5 and DRB1, respectively, than the non-risk GG carriers. The haplotype analysis of expression-associated variants in a Spanish MS cohort revealed that high expression of DRB1 and DQB1 alone did not contribute to the disease. However, in Caucasian, Asian and African American populations, the DRB1*1501 allele was always highly expressed. In other immune related diseases such as type 1 diabetes, inflammatory bowel disease, ulcerative colitis, asthma and IgA deficiency, the best GWAS-associated HLA SNPs were also eQTLs for different HLA Class II genes. Our data suggest that the DR/DQ expression levels, together with specific structural properties of alleles, seem to be the causal effect in MS and in other immunopathologies rather than specific antigen presentation alone.

Association study of genetic variants of pro-inflammatory chemokine and cytokine genes in systemic lupus erythematosus

BACKGROUND. Several lines of evidence suggest that chemokines and cytokines play an important role in the inflammatory development and progression of systemic lupus erythematosus. The aim of this study was to evaluate the relevance of functional genetic variations of RANTES, IL-8, IL-1alpha, and MCP-1 for systemic lupus erythematosus. METHODS. The study was conducted on 500 SLE patients and 481 ethnically matched healthy controls. Genotyping of polymorphisms in the RANTES, IL-8, IL-1alpha, and MCP-1 genes were performed using a real-time polymerase chain reaction (PCR) system with pre-developed TaqMan allelic discrimination assay. RESULTS. No significant differences between SLE patients and healthy controls were observed when comparing genotype, allele or haplotype frequencies of the RANTES, IL-8, IL-1alpha, and MCP-1 polymorphisms. In addition, no evidence for association with clinical sub-features of SLE was found. CONCLUSION. These results suggest that the tested functional variation of RANTES, IL-8, IL-1alpha, and MCP-1 genes do not confer a relevant role in the susceptibility or severity of SLE in the Spanish population.

Genetic analysis of caveolin-1 and eNOS genes in colorectal cancer.

Caveolae are involved in physical compartmentalization between different groups of signaling events. Its main component, CAV1, modulates different pathways in cellular physiology. The emerging evidence pointing to the role of CAV1 in cancer led us to study whether different alleles of this gene are associated with colorectal cancer (CRC). Since one of the most characterized enzymes regulated by CAV1 is eNOS, we decided to include both genes in this study. We analyzed five SNPs in 360 unrelated CRC patients and 550 controls from the general population. Two of these SNPs were located within eNOS and three within the CAV1 gene. Although haplotype distribution was not associated with CRC, haplotype TiA (CAV1) was associated with familiar forms of CRC (p<0.05). This was especially evident in CRC antecedents and nuclear forms of CRC. If both CG (eNOS) and TiA (CAV1) haplotypes were taken together, this association increased in significance. Thus, we propose that CAV1, either alone or together with eNOS alleles, might modify CRC heritability.

Genome-wide association study of multiple sclerosis confirms a novel locus at 5p13.1.

Multiple Sclerosis (MS) is the most common progressive and disabling neurological condition affecting young adults in the world today. From a genetic point of view, MS is a complex disorder resulting from the combination of genetic and non-genetic factors. We aimed to identify previously unidentified loci conducting a new GWAS of Multiple Sclerosis (MS) in a sample of 296 MS cases and 801 controls from the Spanish population. Meta-analysis of our data in combination with previous GWAS was done. A total of 17 GWAS-significant SNPs, corresponding to three different loci were identified:HLA, IL2RA, and 5p13.1. All three have been previously reported as GWAS-significant. We confirmed our observation in 5p13.1 for rs9292777 using two additional independent Spanish samples to make a total of 4912 MS cases and 7498 controls (ORpooled = 0.84; 95%CI: 0.80-0.89; p = 1.36 × 10-9). This SNP differs from the one reported within this locus in a recent GWAS. Although it is unclear whether both signals are tapping the same genetic association, it seems clear that this locus plays an important role in the pathogenesis of MS.

Blood pressure loci identified with a gene-centric array.

Raised blood pressure (BP) is a major risk factor for cardiovascular disease. Previous studies have identified 47 distinct genetic variants robustly associated with BP, but collectively these explain only a few percent of the heritability for BP phenotypes. To find additional BP loci, we used a bespoke gene-centric array to genotype an independent discovery sample of 25,118 individuals that combined hypertensive case-control and general population samples. We followed up four SNPs associated with BP at our p < 8.56 × 10(-7) study-specific significance threshold and six suggestively associated SNPs in a further 59,349 individuals. We identified and replicated a SNP at LSP1/TNNT3, a SNP at MTHFR-NPPB independent (r(2) = 0.33) of previous reports, and replicated SNPs at AGT and ATP2B1 reported previously. An analysis of combined discovery and follow-up data identified SNPs significantly associated with BP at p < 8.56 × 10(-7) at four further loci (NPR3, HFE, NOS3, and SOX6). The high number of discoveries made with modest genotyping effort can be attributed to using a large-scale yet targeted genotyping array and to the development of a weighting scheme that maximized power when meta-analyzing results from samples ascertained with extreme phenotypes, in combination with results from nonascertained or population samples. Chromatin immunoprecipitation and transcript expression data highlight potential gene regulatory mechanisms at the MTHFR and NOS3 loci. These results provide candidates for further study to help dissect mechanisms affecting BP and highlight the utility of studying SNPs and samples that are independent of those studied previously even when the sample size is smaller than that in previous studies.

Common variants of the liver fatty acid binding protein gene influence the risk of type 2 diabetes and insulin resistance in Spanish population

SUMMARY The main objective was to evaluate the association between SNPs and haplotypes of the FABP1-4 genes and type 2 diabetes, as well as its interaction with fat intake, in one general Spanish population. The association was replicated in a second population in which HOMA index was also evaluated. METHODS 1217 unrelated individuals were selected from a population-based study [Hortega study: 605 women; mean age 54 y; 7.8% with type 2 diabetes]. The replication population included 805 subjects from Segovia, a neighboring region of Spain (446 females; mean age 52 y; 10.3% with type 2 diabetes). DM2 mellitus was defined in a similar way in both studies. Fifteen SNPs previously associated with metabolic traits or with potential influence in the gene expression within the FABP1-4 genes were genotyped with SNPlex and tested. Age, sex and BMI were used as covariates in the logistic regression model. RESULTS One polymorphism (rs2197076) and two haplotypes of the FABP-1 showed a strong association with the risk of DM2 in the original population. This association was further confirmed in the second population as well as in the pooled sample. None of the other analyzed variants in FABP2, FABP3 and FABP4 genes were associated. There was not a formal interaction between rs2197076 and fat intake. A significant association between the rs2197076 and the haplotypes of the FABP1 and HOMA-IR was also present in the replication population. CONCLUSIONS The study supports the role of common variants of the FABP-1 gene in the development of type 2 diabetes in Caucasians.

Comprehensive analysis of NRG1 common and rare variants in Hirschsprung patients.

Hirschsprung disease (HSCR, OMIM 142623) is a developmental disorder characterized by the absence of ganglion cells along variable lengths of the distal gastrointestinal tract, which results in tonic contraction of the aganglionic gut segment and functional intestinal obstruction. The RET proto-oncogene is the major gene for HSCR with differential contributions of its rare and common, coding and noncoding mutations to the multifactorial nature of this pathology. Many other genes have been described to be associated with the pathology, as NRG1 gene (8p12), encoding neuregulin 1, which is implicated in the development of the enteric nervous system (ENS), and seems to contribute by both common and rare variants. Here we present the results of a comprehensive analysis of the NRG1 gene in the context of the disease in a series of 207 Spanish HSCR patients, by both mutational screening of its coding sequence and evaluation of 3 common tag SNPs as low penetrance susceptibility factors, finding some potentially damaging variants which we have functionally characterized. All of them were found to be associated with a significant reduction of the normal NRG1 protein levels. The fact that those mutations analyzed alter NRG1 protein would suggest that they would be related with HSCR disease not only in Chinese but also in a Caucasian population, which reinforces the implication of NRG1 gene in this pathology.

Fine mapping and functional analysis of the multiple sclerosis risk gene CD6.

CD6 has recently been identified and validated as risk gene for multiple sclerosis (MS), based on the association of a single nucleotide polymorphism (SNP), rs17824933, located in intron 1. CD6 is a cell surface scavenger receptor involved in T-cell activation and proliferation, as well as in thymocyte differentiation. In this study, we performed a haptag SNP screen of the CD6 gene locus using a total of thirteen tagging SNPs, of which three were non-synonymous SNPs, and replicated the recently reported GWAS SNP rs650258 in a Spanish-Basque collection of 814 controls and 823 cases. Validation of the six most strongly associated SNPs was performed in an independent collection of 2265 MS patients and 2600 healthy controls. We identified association of haplotypes composed of two non-synonymous SNPs [rs11230563 (R225W) and rs2074225 (A257V)] in the 2(nd) SRCR domain with susceptibility to MS (P max(T) permutation = 1×10(-4)). The effect of these haplotypes on CD6 surface expression and cytokine secretion was also tested. The analysis showed significantly different CD6 expression patterns in the distinct cell subsets, i.e. - CD4(+) naïve cells, P = 0.0001; CD8(+) naïve cells, P<0.0001; CD4(+) and CD8(+) central memory cells, P = 0.01 and 0.05, respectively; and natural killer T (NKT) cells, P = 0.02; with the protective haplotype (RA) showing higher expression of CD6. However, no significant changes were observed in natural killer (NK) cells, effector memory and terminally differentiated effector memory T cells. Our findings reveal that this new MS-associated CD6 risk haplotype significantly modifies expression of CD6 on CD4(+) and CD8(+) T cells.

The TT genotype of the STAT4 rs7574865 polymorphism is associated with high disease activity and disability in patients with early arthritis

BACKGROUND The number of copies of the HLA-DRB1 shared epitope, and the minor alleles of the STAT4 rs7574865 and the PTPN22 rs2476601 polymorphisms have all been linked with an increased risk of developing rheumatoid arthritis. In the present study, we investigated the effects of these genetic variants on disease activity and disability in patients with early arthritis. METHODOLOGY AND RESULTS We studied 640 patients with early arthritis (76% women; median age, 52 years), recording disease-related variables every 6 months during a 2-year follow-up. HLA-DRB1 alleles were determined by PCR-SSO, while rs7574865 and rs2476601 were genotyped with the Taqman 5' allelic discrimination assay. Multivariate analysis was performed using generalized estimating equations for repeated measures. After adjusting for confounding variables such as gender, age and ACPA, the TT genotype of rs7574865 in STAT4 was associated with increased disease activity (DAS28) as compared with the GG genotype (β coefficient [95% confidence interval] = 0.42 [0.01-0.83], p = 0.044). Conversely, the presence of the T allele of rs2476601 in PTPN22 was associated with diminished disease activity during follow-up in a dose-dependent manner (CT genotype = -0.27 [-0.56- -0.01], p = 0.042; TT genotype = -0.68 [-1.64- -0.27], p = 0.162). After adjustment for gender, age and disease activity, homozygosity for the T allele of rs7574865 in STAT4 was associated with greater disability as compared with the GG genotype. CONCLUSIONS Our data suggest that patients with early arthritis who are homozygous for the T allele of rs7574865 in STAT4 may develop a more severe form of the disease with increased disease activity and disability.

Genome-wide mRNA expression correlates of viral control in CD4+ T-cells from HIV-1-infected individuals.

There is great interindividual variability in HIV-1 viral setpoint after seroconversion, some of which is known to be due to genetic differences among infected individuals. Here, our focus is on determining, genome-wide, the contribution of variable gene expression to viral control, and to relate it to genomic DNA polymorphism. RNA was extracted from purified CD4+ T-cells from 137 HIV-1 seroconverters, 16 elite controllers, and 3 healthy blood donors. Expression levels of more than 48,000 mRNA transcripts were assessed by the Human-6 v3 Expression BeadChips (Illumina). Genome-wide SNP data was generated from genomic DNA using the HumanHap550 Genotyping BeadChip (Illumina). We observed two distinct profiles with 260 genes differentially expressed depending on HIV-1 viral load. There was significant upregulation of expression of interferon stimulated genes with increasing viral load, including genes of the intrinsic antiretroviral defense. Upon successful antiretroviral treatment, the transcriptome profile of previously viremic individuals reverted to a pattern comparable to that of elite controllers and of uninfected individuals. Genome-wide evaluation of cis-acting SNPs identified genetic variants modulating expression of 190 genes. Those were compared to the genes whose expression was found associated with viral load: expression of one interferon stimulated gene, OAS1, was found to be regulated by a SNP (rs3177979, p = 4.9E-12); however, we could not detect an independent association of the SNP with viral setpoint. Thus, this study represents an attempt to integrate genome-wide SNP signals with genome-wide expression profiles in the search for biological correlates of HIV-1 control. It underscores the paradox of the association between increasing levels of viral load and greater expression of antiviral defense pathways. It also shows that elite controllers do not have a fully distinctive mRNA expression pattern in CD4+ T cells. Overall, changes in global RNA expression reflect responses to viral replication rather than a mechanism that might explain viral control.

Sensitivity of genome-wide-association signals to phenotyping strategy: the PROP-TAS2R38 taste association as a benchmark.

Natural genetic variation can have a pronounced influence on human taste perception, which in turn may influence food preference and dietary choice. Genome-wide association studies represent a powerful tool to understand this influence. To help optimize the design of future genome-wide-association studies on human taste perception we have used the well-known TAS2R38-PROP association as a tool to determine the relative power and efficiency of different phenotyping and data-analysis strategies. The results show that the choice of both data collection and data processing schemes can have a very substantial impact on the power to detect genotypic variation that affects chemosensory perception. Based on these results we provide practical guidelines for the design of future GWAS studies on chemosensory phenotypes. Moreover, in addition to the TAS2R38 gene past studies have implicated a number of other genetic loci to affect taste sensitivity to PROP and the related bitter compound PTC. None of these other locations showed genome-wide significant associations in our study. To facilitate further, target-gene driven, studies on PROP taste perception we provide the genome-wide list of p-values for all SNPs genotyped in the current study.

A genome-wide association search for type 2 diabetes genes in African Americans.

African Americans are disproportionately affected by type 2 diabetes (T2DM) yet few studies have examined T2DM using genome-wide association approaches in this ethnicity. The aim of this study was to identify genes associated with T2DM in the African American population. We performed a Genome Wide Association Study (GWAS) using the Affymetrix 6.0 array in 965 African-American cases with T2DM and end-stage renal disease (T2DM-ESRD) and 1029 population-based controls. The most significant SNPs (n = 550 independent loci) were genotyped in a replication cohort and 122 SNPs (n = 98 independent loci) were further tested through genotyping three additional validation cohorts followed by meta-analysis in all five cohorts totaling 3,132 cases and 3,317 controls. Twelve SNPs had evidence of association in the GWAS (P<0.0071), were directionally consistent in the Replication cohort and were associated with T2DM in subjects without nephropathy (P<0.05). Meta-analysis in all cases and controls revealed a single SNP reaching genome-wide significance (P<2.5×10(-8)). SNP rs7560163 (P = 7.0×10(-9), OR (95% CI) = 0.75 (0.67-0.84)) is located intergenically between RND3 and RBM43. Four additional loci (rs7542900, rs4659485, rs2722769 and rs7107217) were associated with T2DM (P<0.05) and reached more nominal levels of significance (P<2.5×10(-5)) in the overall analysis and may represent novel loci that contribute to T2DM. We have identified novel T2DM-susceptibility variants in the African-American population. Notably, T2DM risk was associated with the major allele and implies an interesting genetic architecture in this population. These results suggest that multiple loci underlie T2DM susceptibility in the African-American population and that these loci are distinct from those identified in other ethnic populations.

The systemic lupus erythematosus IRF5 risk haplotype is associated with systemic sclerosis

Systemic sclerosis (SSc) is a fibrotic autoimmune disease in which the genetic component plays an important role. One of the strongest SSc association signals outside the human leukocyte antigen (HLA) region corresponds to interferon (IFN) regulatory factor 5 (IRF5), a major regulator of the type I IFN pathway. In this study we aimed to evaluate whether three different haplotypic blocks within this locus, which have been shown to alter the protein function influencing systemic lupus erythematosus (SLE) susceptibility, are involved in SSc susceptibility and clinical phenotypes. For that purpose, we genotyped one representative single-nucleotide polymorphism (SNP) of each block (rs10488631, rs2004640, and rs4728142) in a total of 3,361 SSc patients and 4,012 unaffected controls of Caucasian origin from Spain, Germany, The Netherlands, Italy and United Kingdom. A meta-analysis of the allele frequencies was performed to analyse the overall effect of these IRF5 genetic variants on SSc. Allelic combination and dependency tests were also carried out. The three SNPs showed strong associations with the global disease (rs4728142: P  = 1.34×10(-8), OR  = 1.22, CI 95%  = 1.14-1.30; rs2004640: P  = 4.60×10(-7), OR  = 0.84, CI 95%  = 0.78-0.90; rs10488631: P  = 7.53×10(-20), OR  = 1.63, CI 95%  = 1.47-1.81). However, the association of rs2004640 with SSc was not independent of rs4728142 (conditioned P  = 0.598). The haplotype containing the risk alleles (rs4728142*A-rs2004640*T-rs10488631*C: P  = 9.04×10(-22), OR  = 1.75, CI 95%  = 1.56-1.97) better explained the observed association (likelihood P-value  = 1.48×10(-4)), suggesting an additive effect of the three haplotypic blocks. No statistical significance was observed in the comparisons amongst SSc patients with and without the main clinical characteristics. Our data clearly indicate that the SLE risk haplotype also influences SSc predisposition, and that this association is not sub-phenotype-specific.

Two functional variants of IRF5 influence the development of macular edema in patients with non-anterior uveitis.

OBJECTIVE Interferon (IFN) signaling plays a crucial role in autoimmunity. Genetic variation in interferon regulatory factor 5 (IRF5), a major regulator of the type I interferon induction, has been associated with risk of developing several autoimmune diseases. In the current study we aimed to evaluate whether three sets of correlated IRF5 genetic variants, independently associated with SLE and with different functional roles, are involved in uveitis susceptibility and its clinical subphenotypes. METHODS Three IRF5 polymorphisms, rs2004640, rs2070197 and rs10954213, representative of each group, were genotyped using TaqMan® allelic discrimination assays in a total of 263 non-anterior uveitis patients and 724 healthy controls of Spanish origin. RESULTS A clear association between two of the three analyzed genetic variants, rs2004640 and rs10954213, and the absence of macular edema was observed in the case/control analysis (P FDR =5.07E-03, OR=1.48, CI 95%=1.14-1.92 and P FDR =3.37E-03, OR=1.54, CI 95%=1.19-2.01, respectively). Consistently, the subphenotype analysis accordingly with the presence/absence of this clinical condition also reached statistical significance (rs2004640: P=0.037, OR=0.69, CI 95%=0.48-0.98; rs10954213: P=0.030, OR=0.67, CI 95%=0.47-0.96), thus suggesting that both IRF5 genetic variants are specifically associated with the lack of macular edema in uveitis patients. CONCLUSION Our results clearly showed for the first time that two functional genetic variants of IRF5 may play a role in the development of macular edema in non-anterior uveitis patients. Identifying genetic markers for macular edema could lead to the possibility of developing novel treatments or preventive therapies.

Evaluation of the IL2/IL21, IL2RA and IL2RB genetic variants influence on the endogenous non-anterior uveitis genetic predisposition.

BACKGROUND Recently, different genetic variants located within the IL2/IL21 genetic region as well as within both IL2RA and IL2RB loci have been associated to multiple autoimmune disorders. We aimed to investigate for the first time the potential influence of the IL2/IL21, IL2RA and IL2RB most associated polymorphisms with autoimmunity on the endogenous non-anterior uveitis genetic predisposition. METHODS A total of 196 patients with endogenous non-anterior uveitis and 760 healthy controls, all of them from Caucasian population, were included in the current study. The IL2/IL21 (rs2069762, rs6822844 and rs907715), IL2RA (2104286, rs11594656 and rs12722495) and IL2RB (rs743777) genetic variants were genotyped using TaqMan® allelic discrimination assays. RESULTS A statistically significant difference was found for the rs6822844 (IL2/IL21 region) minor allele frequency in the group of uveitis patients compared with controls (P(-value)=0.02, OR=0.64 CI 95%=0.43-0.94) although the significance was lost after multiple testing correction. Furthermore, no evidence of association with uveitis was detected for the analyzed genetic variants of the IL2RA or IL2RB loci. CONCLUSION Our results indicate that analyzed IL2/IL21, IL2RA and IL2RB polymorphisms do not seem to play a significant role on the non-anterior uveitis genetic predisposition although further studies are needed in order to clear up the influence of these loci on the non-anterior uveitis susceptibility.