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.

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.

Gene-lifestyle interaction and type 2 diabetes: the EPIC interact case-cohort study.

BACKGROUND Understanding of the genetic basis of type 2 diabetes (T2D) has progressed rapidly, but the interactions between common genetic variants and lifestyle risk factors have not been systematically investigated in studies with adequate statistical power. Therefore, we aimed to quantify the combined effects of genetic and lifestyle factors on risk of T2D in order to inform strategies for prevention. METHODS AND FINDINGS The InterAct study includes 12,403 incident T2D cases and a representative sub-cohort of 16,154 individuals from a cohort of 340,234 European participants with 3.99 million person-years of follow-up. We studied the combined effects of an additive genetic T2D risk score and modifiable and non-modifiable risk factors using Prentice-weighted Cox regression and random effects meta-analysis methods. The effect of the genetic score was significantly greater in younger individuals (p for interaction  = 1.20×10-4). Relative genetic risk (per standard deviation [4.4 risk alleles]) was also larger in participants who were leaner, both in terms of body mass index (p for interaction  = 1.50×10-3) and waist circumference (p for interaction  = 7.49×10-9). Examination of absolute risks by strata showed the importance of obesity for T2D risk. The 10-y cumulative incidence of T2D rose from 0.25% to 0.89% across extreme quartiles of the genetic score in normal weight individuals, compared to 4.22% to 7.99% in obese individuals. We detected no significant interactions between the genetic score and sex, diabetes family history, physical activity, or dietary habits assessed by a Mediterranean diet score. CONCLUSIONS The relative effect of a T2D genetic risk score is greater in younger and leaner participants. However, this sub-group is at low absolute risk and would not be a logical target for preventive interventions. The high absolute risk associated with obesity at any level of genetic risk highlights the importance of universal rather than targeted approaches to lifestyle intervention.

Response to Infliximab in Crohn's Disease: Genetic Analysis Supporting Expression Profile.

Substantial proportion of Crohn's disease (CD) patients shows no response or a limited response to treatment with infliximab (IFX) and to identify biomarkers of response would be of great clinical and economic benefit. The expression profile of five genes (S100A8-S100A9, G0S2, TNFAIP6, and IL11) reportedly predicted response to IFX and we aimed at investigating their etiologic role through genetic association analysis. Patients with active CD (350) who received at least three induction doses of IFX were included and classified according to IFX response. A tagging strategy was used to select genetic polymorphisms that cover the variability present in the chromosomal regions encoding the identified genes with altered expression. Following genotyping, differences between responders and nonresponders to IFX were observed in haplotypes of the studied regions: S100A8-S100A9 (rs11205276* G/rs3014866* C/rs724781* C/rs3006488* A; P = 0.05); G0S2 (rs4844486* A/rs1473683* T; P = 0.15); TNFAIP6 (rs11677200* C/rs2342910* A/rs3755480* G/rs10432475* A; P = 0.10); and IL11 (rs1126760* C/rs1042506* G; P = 0.07). These differences were amplified in patients with colonic and ileocolonic location for all but the TNFAIP6 haplotype, which evidenced significant difference in ileal CD patients. Our results support the role of the reported expression signature as predictive of anti-TNF outcome in CD patients and suggest an etiological role of those top-five genes in the IFX response pathway.

Prognostic role of genetic biomarkers in clinical progression of prostate cancer.

The aim of this study was to analyze the use of 12 single-nucleotide polymorphisms in genes ELAC2, RNASEL and MSR1 as biomarkers for prostate cancer (PCa) detection and progression, as well as perform a genetic classification of high-risk patients. A cohort of 451 men (235 patients and 216 controls) was studied. We calculated means of regression analysis using clinical values (stage, prostate-specific antigen, Gleason score and progression) in patients and controls at the basal stage and after a follow-up of 72 months. Significantly different allele frequencies between patients and controls were observed for rs1904577 and rs918 (MSR1 gene) and for rs17552022 and rs5030739 (ELAC2). We found evidence of increased risk for PCa in rs486907 and rs2127565 in variants AA and CC, respectively. In addition, rs627928 (TT-GT), rs486907 (AG) and rs3747531 (CG-CC) were associated with low tumor aggressiveness. Some had a weak linkage, such as rs1904577 and rs2127565, rs4792311 and rs17552022, and rs1904577 and rs918. Our study provides the proof-of-principle that some of the genetic variants (such as rs486907, rs627928 and rs2127565) in genes RNASEL, MSR1 and ELAC2 can be used as predictors of aggressiveness and progression of PCa. In the future, clinical use of these biomarkers, in combination with current ones, could potentially reduce the rate of unnecessary biopsies and specific treatments.

Added Value of Deep Sequencing Relative to Population Sequencing in Heavily Pre-Treated HIV-1-Infected Subjects.

OBJECTIVE: To explore the potential of deep HIV-1 sequencing for adding clinically relevant information relative to viral population sequencing in heavily pre-treated HIV-1-infected subjects. METHODS: In a proof-of-concept study, deep sequencing was compared to population sequencing in HIV-1-infected individuals with previous triple-class virological failure who also developed virologic failure to deep salvage therapy including, at least, darunavir, tipranavir, etravirine or raltegravir. Viral susceptibility was inferred before salvage therapy initiation and at virological failure using deep and population sequencing genotypes interpreted with the HIVdb, Rega and ANRS algorithms. The threshold level for mutant detection with deep sequencing was 1%. RESULTS: 7 subjects with previous exposure to a median of 15 antiretrovirals during a median of 13 years were included. Deep salvage therapy included darunavir, tipranavir, etravirine or raltegravir in 4, 2, 2 and 5 subjects, respectively. Self-reported treatment adherence was adequate in 4 and partial in 2; one individual underwent treatment interruption during follow-up. Deep sequencing detected all mutations found by population sequencing and identified additional resistance mutations in all but one individual, predominantly after virological failure to deep salvage therapy. Additional genotypic information led to consistent decreases in predicted susceptibility to etravirine, efavirenz, nucleoside reverse transcriptase inhibitors and indinavir in 2, 1, 2 and 1 subject, respectively. Deep sequencing data did not consistently modify the susceptibility predictions achieved with population sequencing for darunavir, tipranavir or raltegravir. CONCLUSIONS: In this subset of heavily pre-treated individuals, deep sequencing improved the assessment of genotypic resistance to etravirine, but did not consistently provide additional information on darunavir, tipranavir or raltegravir susceptibility. These data may inform the design of future studies addressing the clinical value of minority drug-resistant variants in treatment-experienced subjects.

A Genome-Wide Association Study of Upper Aerodigestive Tract Cancers Conducted within the INHANCE Consortium

Genome-wide association studies (GWAS) have been successful in identifying common genetic variation involved in susceptibility to etiologically complex disease. We conducted a GWAS to identify common genetic variation involved in susceptibility to upper aero-digestive tract (UADT) cancers. Genome-wide genotyping was carried out using the Illumina HumanHap300 beadchips in 2,091 UADT cancer cases and 3,513 controls from two large European multi-centre UADT cancer studies, as well as 4,821 generic controls. The 19 top-ranked variants were investigated further in an additional 6,514 UADT cancer cases and 7,892 controls of European descent from an additional 13 UADT cancer studies participating in the INHANCE consortium. Five common variants presented evidence for significant association in the combined analysis (p≤5×10−7). Two novel variants were identified, a 4q21 variant (rs1494961, p = 1×10−8) located near DNA repair related genes HEL308 and FAM175A (or Abraxas) and a 12q24 variant (rs4767364, p = 2×10−8) located in an extended linkage disequilibrium region that contains multiple genes including the aldehyde dehydrogenase 2 (ALDH2) gene. Three remaining variants are located in the ADH gene cluster and were identified previously in a candidate gene study involving some of these samples. The association between these three variants and UADT cancers was independently replicated in 5,092 UADT cancer cases and 6,794 controls non-overlapping samples presented here (rs1573496-ADH7, p = 5×10−8; rs1229984-ADH1B, p = 7×10−9; and rs698-ADH1C, p = 0.02). These results implicate two variants at 4q21 and 12q24 and further highlight three ADH variants in UADT cancer susceptibility.

The INSIG2 rs7566605 polymorphism is not associated with body mass index and breast cancer risk

Background The single nucleotide polymorphism rs7566605, located in the promoter of the INSIG2 gene, has been the subject of a strong scientific effort aimed to elucidate its possible association with body mass index (BMI). The first report showing that rs7566605 could be associated with body fatness was a genome-wide association study (GWAS) which used BMI as the primary phenotype. Many follow-up studies sought to validate the association of rs7566605 with various markers of obesity, with several publications reporting inconsistent findings. BMI is considered to be one of the measures of choice to evaluate body fatness and there is evidence that body fatness is related with an increased risk of breast cancer (BC). Methods we tested in a large-scale association study (3,973 women, including 1,269 invasive BC cases and 2,194 controls), nested within the EPIC cohort, the involvement of rs7566605 as predictor of BMI and BC risk. Results and Conclusions In this study we were not able to find any statistically significant association between this SNP and BMI, nor did we find any significant association between the SNP and an increased risk of breast cancer overall and by subgroups of age, or menopausal status.

Comprehensive analysis of RET common and rare variant patientsin a series of Spanish Hirschsprung patients confirms a synergistic effect of both kinds of events

Background. RET is the major gene associated to Hirschsprung disease (HSCR) with differential contributions of its rare and common, coding and noncoding mutations to the multifactorial nature of this pathology. In the present study, we have performed a comprehensive study of our HSCR series evaluating the involvement of both RET rare variants (RVs) and common variants (CVs) in the context of the disease. Methods. RET mutational screening was performed by dHPLC and direct sequencing for the identification of RVs. In addition Taqman technology was applied for the genotyping of 3 RET CVs previously associated to HSCR, including a variant lying in an enhancer domain within RET intron 1 (rs2435357). Statistical analyses were performed using the SPSS v.17.0 to analyze the distribution of the variants. Results. Our results confirm the strongest association to HSCR for the "enhancer" variant, and demonstrate a significantly higher impact of it in male versus female patients. Integration of the RET RVs and CVs analysis showed that in 91.66% of cases with both kinds of mutational events, the enhancer allele is in trans with the allele bearing the RET RV. Conclusions. A gender effect exists on both the transmission and distribution of rare coding and common HSCR causing mutations. In addition, these RET CVs and RVs seem to act in a synergistic way leading to HSCR phenotype.

Association analysis of urotensin II gene (UTS2) and flanking regions with biochemical parameters related to insulin resistance.

Background. Urotensin II (UII) is a potent vasoconstrictor peptide, which signals through a G-protein coupled receptor (GPCR) known as GPR14 or urotensin receptor (UTR). UII exerts a broad spectrum of actions in several systems such as vascular cell, heart muscle or pancreas, where it inhibits insulin release. Objective. Given the reported role of UII in insulin secretion, we have performed a genetic association analysis of the UTS2 gene and flanking regions with biochemical parameters related to insulin resistance (fasting glucose, glucose 2 hours after a glucose overload, fasting insulin and insulin resistance estimated as HOMA). Results and Conclusions. We have identified several polymorphisms associated with the analysed clinical traits, not only at the UTS2 gene, but also in thePER3 gene, located upstream from UTS2. Our results are compatible with a role for UII in glucose homeostasis and diabetes although we cannot rule out the possibility that PER3 gene may underlie the reported associations.

Prevalence and resistance mutations of non-B HIV-1 subtypes among immigrants in Southern Spain along the decade 2000-2010.

Most of the non-B HIV-1 subtypes are predominant in Sub-Saharan Africa and India although they have been found worldwide. In the last decade, immigration from these areas has increased considerably in Spain. The objective of this study was to evaluate the prevalence of non-B subtypes circulating in a cohort of HIV-1-infected immigrants in Seville, Southern Spain and to identify drug resistance-associated mutations. METHODS: Complete protease and first 220 codons of the reverse transcriptase coding regions were amplified and sequenced by population sequencing. HIV-1 subtypes were determined using Stanford University Drug Resistance Database, and phylogenetic analysis was performed comparing multiple reported sequences. Drug resistance mutations were defined according to the International AIDS Society-USA. RESULTS: From 2000 to 2010 a total of 1,089 newly diagnosed HIV-1-infected patients were enrolled in our cohort. Of these, 121 were immigrants, of which 98 had ethical approval and informed consent to include in our study. Twenty-nine immigrants (29/98, 29.6%) were infected with non-B subtypes, of which 15/29 (51.7%) were CRF02-AG, mostly from Sub-Saharan Africa, and 2/29 (6.9%) were CRF01-AE from Eastern Europe. A, C, F, J and G subtypes from Eastern Europe, Central-South America and Sub-Saharan Africa were also present. Some others harboured recombinant forms CRF02-AG/CRF01-AE, CRF2-AG/G and F/B, B/C, and K/G, in PR and RT-coding regions. Patients infected with non-B subtypes showed a high frequency of minor protease inhibitor resistance mutations, M36I, L63P, and K20R/I. Only one patient, CRF02_AG, showed major resistance mutation L90M. Major RT inhibitor resistance mutations K70R and A98G were present in one patient with subtype G, L100I in one patient with CRF01_AE, and K103N in another patient with CRF01_AE. Three patients had other mutations such as V118I, E138A and V90I. CONCLUSIONS: The circulation of non-B subtypes has significantly increased in Southern Spain during the last decade, with 29.6% prevalence, in association with demographic changes among immigrants. This could be an issue in the treatment and management of these patients. Resistance mutations have been detected in these patients with a prevalence of 7% among treatment-naïve patients compared with the 21% detected among patients under HAART or during treatment interruption.

TRAIL/TRAIL Receptor System and Susceptibility to Multiple Sclerosis

The TNF-related apoptosis inducing ligand (TRAIL)/TRAIL receptor system participates in crucial steps in immune cell activation or differentiation. It is able to inhibit proliferation and activation of T cells and to induce apoptosis of neurons and oligodendrocytes, and seems to be implicated in autoimmune diseases. Thus, TRAIL and TRAIL receptor genes are potential candidates for involvement in susceptibility to multiple sclerosis (MS). To test whether single-nucleotide polymorphisms (SNPs) in the human genes encoding TRAIL, TRAILR-1, TRAILR-2, TRAILR-3 and TRAILR-4 are associated with MS susceptibility, we performed a candidate gene case-control study in the Spanish population. 59 SNPs in the TRAIL and TRAIL receptor genes were analysed in 628 MS patients and 660 controls, and validated in an additional cohort of 295 MS patients and 233 controls. Despite none of the SNPs withstood the highly conservative Bonferroni correction, three SNPs showing uncorrected p values<0.05 were successfully replicated: rs4894559 in TRAIL gene, p = 9.8×10(-4), OR = 1.34; rs4872077, in TRAILR-1 gene, p = 0.005, OR = 1.72; and rs1001793 in TRAILR-2 gene, p = 0.012, OR = 0.84. The combination of the alleles G/T/A in these SNPs appears to be associated with a reduced risk of developing MS (p = 2.12×10(-5), OR = 0.59). These results suggest that genes of the TRAIL/TRAIL receptor system exerts a genetic influence on MS.

Functional and transcriptional induction of aquaporin-1 gene by hypoxia; analysis of promoter and role of Hif-1α.

Aquaporin-1 (AQP1) is a water channel that is highly expressed in tissues with rapid O(2) transport. It has been reported that this protein contributes to gas permeation (CO(2), NO and O(2)) through the plasma membrane. We show that hypoxia increases Aqp1 mRNA and protein levels in tissues, namely mouse brain and lung, and in cultured cells, the 9L glioma cell line. Stopped-flow light-scattering experiments confirmed an increase in the water permeability of 9L cells exposed to hypoxia, supporting the view that hypoxic Aqp1 up-regulation has a functional role. To investigate the molecular mechanisms underlying this regulatory process, transcriptional regulation was studied by transient transfections of mouse endothelial cells with a 1297 bp 5' proximal Aqp1 promoter-luciferase construct. Incubation in hypoxia produced a dose- and time-dependent induction of luciferase activity that was also obtained after treatments with hypoxia mimetics (DMOG and CoCl(2)) and by overexpressing stabilized mutated forms of HIF-1α. Single mutations or full deletions of the three putative HIF binding domains present in the Aqp1 promoter partially reduced its responsiveness to hypoxia, and transfection with Hif-1α siRNA decreased the in vitro hypoxia induction of Aqp1 mRNA and protein levels. Our results indicate that HIF-1α participates in the hypoxic induction of AQP1. However, we also demonstrate that the activation of Aqp1 promoter by hypoxia is complex and multifactorial and suggest that besides HIF-1α other transcription factors might contribute to this regulatory process. These data provide a conceptual framework to support future research on the involvement of AQP1 in a range of pathophysiological conditions, including edema, tumor growth, and respiratory diseases.

Frequency of Fabry disease in male and female haemodialysis patients in Spain

BACKGROUND: Fabry disease (FD), an X-linked lysosomal storage disorder, is caused by a reduced activity of the lysosomal enzyme alpha-galactosidase A. The disorder ultimately leads to organ damage (including renal failure) in males and females. However, heterozygous females usually present a milder phenotype with a later onset and a slower progression. METHODS: A combined enzymatic and genetic strategy was used, measuring the activity of alpha-galactosidase A and genotyping the alpha-galactosidase A gene (GLA) in dried blood samples (DBS) of 911 patients undergoing haemodialysis in centers across Spain. RESULTS: GLA alterations were found in seven unrelated patients (4 males and 3 females). Two novel mutations (p.Gly346AlafsX347 and p.Val199GlyfsX203) were identified as well as a previously described mutation, R118C. The R118C mutation was present in 60% of unrelated patients with GLA causal mutations. The D313Y alteration, considered by some authors as a pseudo-deficiency allele, was also found in two out of seven patients. CONCLUSIONS: Excluding the controversial D313Y alteration, FD presents a frequency of one in 182 individuals (0.55%) within this population of males and females undergoing haemodialysis. Moreover, our findings suggest that a number of patients with unexplained and atypical symptoms of renal disease may have FD. Screening programmes for FD in populations of individuals presenting severe kidney dysfunction, cardiac alterations or cerebrovascular disease may lead to the diagnosis of FD in those patients, the study of their families and eventually the implementation of a specific therapy.

Novel MLPA procedure using self-designed probes allows comprehensive analysis for CNVs of the genes involved in Hirschsprung disease.

Background: Hirschsprung disease is characterized by the absence of intramural ganglion cells in the enteric plexuses, due to a fail during enteric nervous system formation. Hirschsprung has a complex genetic aetiology and mutations in several genes have been related to the disease. There is a clear predominance of missense/nonsense mutations in these genes whereas copy number variations (CNVs) have been seldom described, probably due to the limitations of conventional techniques usually employed for mutational analysis. In this study, we have looked for CNVs in some of the genes related to Hirschsprung (EDNRB, GFRA1, NRTN and PHOX2B) using the Multiple Ligation-dependent Probe Amplification (MLPA) approach. Methods: CNVs screening was performed in 208 HSCR patients using a self-designed set of MLPA probes, covering the coding region of those genes. Results: A deletion comprising the first 4 exons in GFRA1 gene was detected in 2 sporadic HSCR patients and in silico approaches have shown that the critical translation initiation signal in the mutant gene was abolished. In this study, we have been able to validate the reliability of this technique for CNVs screening in HSCR. Conclusions: The implemented MLPA based technique presented here allows CNV analysis of genes involved in HSCR that have not been not previously evaluated. Our results indicate that CNVs could be implicated in the pathogenesis of HSCR, although they seem to be an uncommon molecular cause of HSCR.