Proteins encoded in genomic regions associated with immune-mediated disease physically interact and suggest underlying biology

Genome-wide association studies (GWAS) have defined over 150 genomic regions unequivocally containing variation predisposing to immune-mediated disease. Inferring disease biology from these observations, however, hinges on our ability to discover the molecular processes being perturbed by these risk variants. It has previously been observed that different genes harboring causal mutations for the same Mendelian disease often physically interact. We sought to evaluate the degree to which this is true of genes within strongly associated loci in complex disease. Using sets of loci defined in rheumatoid arthritis (RA) and Crohn's disease (CD) GWAS, we build protein-protein interaction (PPI) networks for genes within associated loci and find abundant physical interactions between protein products of associated genes. We apply multiple permutation approaches to show that these networks are more densely connected than chance expectation. To confirm biological relevance, we show that the components of the networks tend to be expressed in similar tissues relevant to the phenotypes in question, suggesting the network indicates common underlying processes perturbed by risk loci. Furthermore, we show that the RA and CD networks have predictive power by demonstrating that proteins in these networks, not encoded in the confirmed list of disease associated loci, are significantly enriched for association to the phenotypes in question in extended GWAS analysis. Finally, we test our method in 3 non-immune traits to assess its applicability to complex traits in general. We find that genes in loci associated to height and lipid levels assemble into significantly connected networks but did not detect excess connectivity among Type 2 Diabetes (T2D) loci beyond chance. Taken together, our results constitute evidence that, for many of the complex diseases studied here, common genetic associations implicate regions encoding proteins that physically interact in a preferential manner, in line with observations in Mendelian disease.

Estimating the net contribution of interleukin-28B variation to spontaneous hepatitis C virus clearance

The identification of associations between interleukin-28B (IL-28B) variants and the spontaneous clearance of hepatitis C virus (HCV) raises the issues of causality and the net contribution of host genetics to the trait. To estimate more precisely the net effect of IL-28B genetic variation on HCV clearance, we optimized genotyping and compared the host contributions in multiple- and single-source cohorts to control for viral and demographic effects. The analysis included individuals with chronic or spontaneously cleared HCV infections from a multiple-source cohort (n = 389) and a single-source cohort (n = 71). We performed detailed genotyping in the coding region of IL-28B and searched for copy number variations to identify the genetic variant or haplotype carrying the strongest association with viral clearance. This analysis was used to compare the effects of IL-28B variation in the two cohorts. Haplotypes characterized by carriage of the major alleles at IL-28B single-nucleotide polymorphisms (SNPs) were highly overrepresented in individuals with spontaneous clearance versus those with chronic HCV infections (66.1% versus 38.6%, P = 6 × 10(-9) ). The odds ratios for clearance were 2.1 [95% confidence interval (CI) = 1.6-3.0] and 3.9 (95% CI = 1.5-10.2) in the multiple- and single-source cohorts, respectively. Protective haplotypes were in perfect linkage (r(2) = 1.0) with a nonsynonymous coding variant (rs8103142). Copy number variants were not detected. CONCLUSION: We identified IL-28B haplotypes highly predictive of spontaneous HCV clearance. The high linkage disequilibrium between IL-28B SNPs indicates that association studies need to be complemented by functional experiments to identify single causal variants. The point estimate for the genetic effect was higher in the single-source cohort, which was used to effectively control for viral diversity, sex, and coinfections and, therefore, offered a precise estimate of the net host genetic contribution.

LUPA: a European initiative taking advantage of the canine genome architecture for unravelling complex disorders in both human and dogs

The domestic dog offers a unique opportunity to explore the genetic basis of disease, morphology and behaviour. Humans share many diseases with our canine companions, making dogs an ideal model organism for comparative disease genetics. Using newly developed resources, genome-wide association studies in dog breeds are proving to be exceptionally powerful. Towards this aim, veterinarians and geneticists from 12 European countries are collaborating to collect and analyse the DNA from large cohorts of dogs suffering from a range of carefully defined diseases of relevance to human health. This project, named LUPA, has already delivered considerable results. The consortium has collaborated to develop a new high density single nucleotide polymorphism (SNP) array. Mutations for four monogenic diseases have been identified and the information has been utilised to find mutations in human patients. Several complex diseases have been mapped and fine mapping is underway. These findings should ultimately lead to a better understanding of the molecular mechanisms underlying complex diseases in both humans and their best friend.

Genome-wide analysis of rare copy number variations reveals PARK2 as a candidate gene for attention-deficit/hyperactivity disorder

Attention-deficit/hyperactivity disorder (ADHD) is a common, highly heritable neurodevelopmental disorder. Genetic loci have not yet been identified by genome-wide association studies. Rare copy number variations (CNVs), such as chromosomal deletions or duplications, have been implicated in ADHD and other neurodevelopmental disorders. To identify rare (frequency 1%) CNVs that increase the risk of ADHD, we performed a whole-genome CNV analysis based on 489 young ADHD patients and 1285 adult population-based controls and identified one significantly associated CNV region. In tests for a global burden of large (>500 kb) rare CNVs, we observed a nonsignificant (P=0.271) 1.126-fold enriched rate of subjects carrying at least one such CNV in the group of ADHD cases. Locus-specific tests of association were used to assess if there were more rare CNVs in cases compared with controls. Detected CNVs, which were significantly enriched in the ADHD group, were validated by quantitative (q)PCR. Findings were replicated in an independent sample of 386 young patients with ADHD and 781 young population-based healthy controls. We identified rare CNVs within the parkinson protein 2 gene (PARK2) with a significantly higher prevalence in ADHD patients than in controls (P=2.8 × 10(-4) after empirical correction for genome-wide testing). In total, the PARK2 locus (chr 6: 162 659 756-162 767 019) harboured three deletions and nine duplications in the ADHD patients and two deletions and two duplications in the controls. By qPCR analysis, we validated 11 of the 12 CNVs in ADHD patients (P=1.2 × 10(-3) after empirical correction for genome-wide testing). In the replication sample, CNVs at the PARK2 locus were found in four additional ADHD patients and one additional control (P=4.3 × 10(-2)). Our results suggest that copy number variants at the PARK2 locus contribute to the genetic susceptibility of ADHD. Mutations and CNVs in PARK2 are known to be associated with Parkinson disease.Molecular Psychiatry advance online publication, 20 November 2012; doi:10.1038/mp.2012.161.

DCLK1 variants are associated across schizophrenia and attention deficit/hyperactivity disorder

Doublecortin and calmodulin like kinase 1 (DCLK1) is implicated in synaptic plasticity and neurodevelopment. Genetic variants in DCLK1 are associated with cognitive traits, specifically verbal memory and general cognition. We investigated the role of DCLK1 variants in three psychiatric disorders that have neuro-cognitive dysfunctions: schizophrenia (SCZ), bipolar affective disorder (BP) and attention deficit/hyperactivity disorder (ADHD). We mined six genome wide association studies (GWASs) that were available publically or through collaboration; three for BP, two for SCZ and one for ADHD. We also genotyped the DCLK1 region in additional samples of cases with SCZ, BP or ADHD and controls that had not been whole-genome typed. In total, 9895 subjects were analysed, including 5308 normal controls and 4,587 patients (1,125 with SCZ, 2,496 with BP and 966 with ADHD). Several DCLK1 variants were associated with disease phenotypes in the different samples. The main effect was observed for rs7989807 in intron 3, which was strongly associated with SCZ alone and even more so when cases with SCZ and ADHD were combined (P-value = 4 × 10(-5) and 4 × 10(-6), respectively). Associations were also observed with additional markers in intron 3 (combination of SCZ, ADHD and BP), intron 19 (SCZ+BP) and the 3'UTR (SCZ+BP). Our results suggest that genetic variants in DCLK1 are associated with SCZ and, to a lesser extent, with ADHD and BP. Interestingly the association is strongest when SCZ and ADHD are considered together, suggesting common genetic susceptibility. Given that DCLK1 variants were previously found to be associated with cognitive traits, these results are consistent with the role of DCLK1 in neurodevelopment and synaptic plasticity.

Alcohol and colorectal cancer: the role of alcohol dehydrogenase 1C polymorphism

BACKGROUND: Chronic alcohol consumption is a risk factor for colorectal cancer. Animal experiments as well as genetic linkage studies in Japanese individuals with inactive acetaldehyde dehydrogenase leading to elevated acetaldehyde concentrations following ethanol ingestion support the hypothesis that acetaldehyde may be responsible for this carcinogenic effect of alcohol. In Caucasians, a polymorphism of alcohol dehydrogenase 1C (ADH1C) exists resulting in different acetaldehyde concentrations following ethanol oxidation. METHODS: To evaluate whether the association between alcohol consumption and colorectal tumor development is modified by ADH1C polymorphism, we recruited 173 individuals with colorectal tumors diagnosed by colonoscopy and 788 control individuals without colorectal tumors. Genotyping was performed using genomic DNA extracted from whole blood followed by polymerase chain reaction. RESULTS: Genotype ADH1C*1/1 was more frequent in patients with alcohol-associated colorectal neoplasia compared to patients without cancers in the multivariate model controlling for age, gender, and alcohol intake (odds ratio = 1.674, 95% confidence interval = 1.110-2.524, 2-sided p from Wald test = 0.0139). In addition, the joint test of the genetic effect and interaction between ADH1C genotype and alcohol intake (2-sided p = 0.0007) indicated that the difference in ADH1C*1 polymorphisms between controls and colorectal neoplasia is strongly influenced by the alcohol consumption and that only individuals drinking more than 30 g ethanol per day with the genotype ADH1C*1/1 had an increased risk for colorectal tumors. CONCLUSIONS: These data identify ADH1C homozygosity as a genetic risk marker for colorectal tumors in individuals consuming more than 30 g alcohol per day and emphasize the role of acetaldehyde as a carcinogenic agent in alcohol-related colorectal carcinogenesis.

Contribution of genome-wide significant single-nucleotide polymorphisms and antiretroviral therapy to dyslipidemia in HIV-infected individuals: a longitudinal study

BACKGROUND: HIV-infected individuals have an increased risk of myocardial infarction. Antiretroviral therapy (ART) is regarded as a major determinant of dyslipidemia in HIV-infected individuals. Previous genetic studies have been limited by the validity of the single-nucleotide polymorphisms (SNPs) interrogated and by cross-sectional design. Recent genome-wide association studies have reliably associated common SNPs to dyslipidemia in the general population. METHODS AND RESULTS: We validated the contribution of 42 SNPs (33 identified in genome-wide association studies and 9 previously reported SNPs not included in genome-wide association study chips) and of longitudinally measured key nongenetic variables (ART, underlying conditions, sex, age, ethnicity, and HIV disease parameters) to dyslipidemia in 745 HIV-infected study participants (n=34 565 lipid measurements; median follow-up, 7.6 years). The relative impact of SNPs and ART to lipid variation in the study population and their cumulative influence on sustained dyslipidemia at the level of the individual were calculated. SNPs were associated with lipid changes consistent with genome-wide association study estimates. SNPs explained up to 7.6% (non-high-density lipoprotein cholesterol), 6.2% (high-density lipoprotein cholesterol), and 6.8% (triglycerides) of lipid variation; ART explained 3.9% (non-high-density lipoprotein cholesterol), 1.5% (high-density lipoprotein cholesterol), and 6.2% (triglycerides). An individual with the most dyslipidemic antiretroviral and genetic background had an approximately 3- to 5-fold increased risk of sustained dyslipidemia compared with an individual with the least dyslipidemic therapy and genetic background. CONCLUSIONS: In the HIV-infected population treated with ART, the weight of the contribution of common SNPs and ART to dyslipidemia was similar. When selecting an ART regimen, genetic information should be considered in addition to the dyslipidemic effects of ART agents.

Common noncoding UMOD gene variants induce salt-sensitive hypertension and kidney damage by increasing uromodulin expression

Hypertension and chronic kidney disease (CKD) are complex traits representing major global health problems1,2. Multiple genome-wide association studies have identified common variants in the promoter of the UMOD gene3–9, which encodes uromodulin, the major protein secreted in normal urine, that cause independent susceptibility to CKD and hypertension. Despite compelling genetic evidence for the association between UMOD risk variants and disease susceptibility in the general population, the underlying biological mechanism is not understood. Here, we demonstrate that UMOD risk variants increased UMOD expression in vitro and in vivo. Uromodulin overexpression in transgenic mice led to salt-sensitive hypertension and to the presence of age-dependent renal lesions similar to those observed in elderly individuals homozygous for UMOD promoter risk variants. The link between uromodulin and hypertension is due to activation of the renal sodium cotransporter NKCC2. We demonstrated the relevance of this mechanism in humans by showing that pharmacological inhibition of NKCC2 was more effective in lowering blood pressure in hypertensive patients who are homozygous for UMOD promoter risk variants than in other hypertensive patients. Our findings link genetic susceptibility to hypertension and CKD to the level of uromodulin expression and uromodulin’s effect on salt reabsorption in the kidney. These findings point to uromodulin as a therapeutic target for lowering blood pressure and preserving renal function.

Genetic and microbial factors modulating the ubiquitin proteasome system in inflammatory bowel disease

OBJECTIVE: Altered microbiota composition, changes in immune responses and impaired intestinal barrier functions are observed in IBD. Most of these features are controlled by proteases and their inhibitors to maintain gut homeostasis. Unrestrained or excessive proteolysis can lead to pathological gastrointestinal conditions. The aim was to validate the identified protease IBD candidates from a previously performed systematic review through a genetic association study and functional follow-up. DESIGN: We performed a genetic association study in a large multicentre cohort of patients with Crohn's disease (CD) and UC from five European IBD referral centres in a total of 2320 CD patients, 2112 UC patients and 1796 healthy controls. Subsequently, we did an extensive functional assessment of the candidate genes to explore their causality in IBD pathogenesis. RESULTS: Ten single nucleotide polymorphisms (SNPs) in four genes were significantly associated with CD: CYLD, USP40, APEH and USP3. CYLD was the most significant gene with the intronically located rs12324931 the strongest associated SNP (pFDR=1.74e-17, OR=2.24 (1.83 to 2.74)). Five SNPs in four genes were significantly associated with UC: USP40, APEH, DAG1 and USP3. CYLD, as well as some of the other associated genes, is part of the ubiquitin proteasome system (UPS). We therefore determined if the IBD-associated adherent-invasive Escherichia coli (AIEC) can modulate the UPS functioning. Infection of intestinal epithelial cells with the AIEC LF82 reference strain modulated the UPS turnover by reducing poly-ubiquitin conjugate accumulation, increasing 26S proteasome activities and decreasing protein levels of the NF-κB regulator CYLD. This resulted in IκB-α degradation and NF-κB activation. This activity was very important for the pathogenicity of AIEC since decreased CYLD resulted in increased ability of AIEC LF82 to replicate intracellularly. CONCLUSIONS: Our results reveal the UPS, and CYLD specifically, as an important contributor to IBD pathogenesis, which is favoured by both genetic and microbial factors.

IL28B expression depends on a novel TT/-G polymorphism which improves HCV clearance prediction

Approximately 3% of the world population is chronically infected with the hepatitis C virus (HCV), with potential development of cirrhosis and hepatocellular carcinoma. Despite the availability of new antiviral agents, treatment remains suboptimal. Genome-wide association studies (GWAS) identified rs12979860, a polymorphism nearby IL28B, as an important predictor of HCV clearance. We report the identification of a novel TT/-G polymorphism in the CpG region upstream of IL28B, which is a better predictor of HCV clearance than rs12979860. By using peripheral blood mononuclear cells (PBMCs) from individuals carrying different allelic combinations of the TT/-G and rs12979860 polymorphisms, we show that induction of IL28B and IFN-γ–inducible protein 10 (IP-10) mRNA relies on TT/-G, but not rs12979860, making TT/-G the only functional variant identified so far. This novel step in understanding the genetic regulation of IL28B may have important implications for clinical practice, as the use of TT/G genotyping instead of rs12979860 would improve patient management.

New generation genomic platforms in investigation of complex diseases and BEN.

(Full text is available at http://www.manu.edu.mk/prilozi). New generation genomic platforms enable us to decipher the complex genetic basis of complex diseases and Balkan Endemic Nephropathy (BEN) at a high-throughput basis. They give valuable information about predisposing Single Nucleotide Polymorphisms (SNPs), Copy Number Variations (CNVs) or Loss of Heterozygosity (LOH) (using SNP-array) and about disease-causing mutations along the whole sequence of candidate-genes (using Next Generation Sequencing). This information could be used for screening of individuals in risk families and moving the main medicine stream to the prevention. They also might have an impact on more effective treatment. Here we discuss these genomic platforms and report some applications of SNP-array technology in a case with familial nephrotic syndrome. Key words: complex diseases, genome wide association studies, SNP, genomic arrays, next generation sequ-encing.

Two Loci on Chromosome 5 Are Associated with Serum IgE Levels in Labrador Retrievers

Crosslinking of immunoglobulin E antibodies (IgE) bound at the surface of mast cells and subsequent mediator release is considered the most important trigger for allergic reactions. Therefore, the genetic control of IgE levels is studied in the context of allergic diseases, such as asthma, atopic rhinitis, or atopic dermatitis (AD). We performed genome-wide association studies in 161 Labrador Retrievers with regard to total and allergen-specific immunoglobulin E (IgE) levels. We identified a genome-wide significant association on CFA 5 with the antigen-specific IgE responsiveness to Acarus siro. We detected a second genome-wide significant association with respect to the antigen-specific IgE responsiveness to Tyrophagus putrescentiae at a different locus on chromosome 5. A. siro and T. putrescentiae both belong to the family Acaridae and represent so-called storage or forage mites. These forage mites are discussed as major allergen sources in canine AD. No obvious candidate gene for the regulation of IgE levels is located under the two association signals. Therefore our studies offer a chance of identifying a novel mechanism controlling the host's IgE response.

The brown coat colour of Coppernecked goats is associated with a non-synonymous variant at the TYRP1 locus on chromosome 8

The recent development of a goat SNP genotyping microarray enables genome-wide association studies in this important livestock species. We investigated the genetic basis of the black and brown coat colour in Valais Blacknecked and Coppernecked goats. A genome-wide association analysis using goat SNP50 BeadChip genotypes of 22 cases and 23 controls allowed us to map the locus for the brown coat colour to goat chromosome 8. The TYRP1 gene is located within the associated chromosomal region, and TYRP1 variants cause similar coat colour phenotypes in different species. We thus considered TYRP1 as a strong positional and functional candidate. We resequenced the caprine TYRP1 gene by Sanger and Illumina sequencing and identified two non-synonymous variants, p.Ile478Thr and p.Gly496Asp, that might have a functional impact on the TYRP1 protein. However, based on the obtained pedigree and genotype data, the brown coat colour in these goats is not due to a single recessive loss-of-function allele. Surprisingly, the genotype distribution and the pedigree data suggest that the (496) Asp allele might possibly act in a dominant manner. The (496) Asp allele was present in 77 of 81 investigated Coppernecked goats and did not occur in black goats. This strongly suggests heterogeneity underlying the brown coat colour in Coppernecked goats. Functional experiments or targeted matings will be required to verify the unexpected preliminary findings.

Evaluation of genome-wide loci of iron metabolism in hereditary hemochromatosis identifies PCSK7 as a host risk factor of liver cirrhosis

Genome-wide association studies (GWAS) have revealed genetic determinants of iron metabolism, but correlation of these with clinical phenotypes is pending. Homozygosity for HFE C282Y is the predominant genetic risk factor for hereditary hemochromatosis (HH) and may cause liver cirrhosis. However, this genotype has a low penetrance. Thus, detection of yet unknown genetic markers that identify patients at risk of developing severe liver disease is necessary for better prevention. Genetic loci associated with iron metabolism (TF, TMPRSS6, PCSK7, TFR2 and Chr2p14) in recent GWAS and liver fibrosis (PNPLA3) in recent meta-analysis were analyzed for association with either liver cirrhosis or advanced fibrosis in 148 German HFE C282Y homozygotes. Replication of associations was sought in additional 499 Austrian/Swiss and 112 HFE C282Y homozygotes from Sweden. Only variant rs236918 in the PCSK7 gene (proprotein convertase subtilisin/kexin type 7) was associated with cirrhosis or advanced fibrosis (P = 1.02 × 10(-5)) in the German cohort with genotypic odds ratios of 3.56 (95% CI 1.29-9.77) for CG heterozygotes and 5.38 (95% CI 2.39-12.10) for C allele carriers. Association between rs236918 and cirrhosis was confirmed in Austrian/Swiss HFE C282Y homozygotes (P = 0.014; ORallelic = 1.82 (95% CI 1.12-2.95) but not in Swedish patients. Post hoc combined analyses of German/Swiss/Austrian patients with available liver histology (N = 244, P = 0.00014, ORallelic = 2.84) and of males only (N = 431, P = 2.17 × 10(-5), ORallelic = 2.54) were consistent with the premier finding. Association between rs236918 and cirrhosis was not confirmed in alcoholic cirrhotics, suggesting specificity of this genetic risk factor for HH. PCSK7 variant rs236918 is a risk factor for cirrhosis in HH patients homozygous for the HFE C282Y mutation.

Disease-modifying genes and monogenic disorders: experience in cystic fibrosis

The mechanisms responsible for the determination of phenotypes are still not well understood; however, it has become apparent that modifier genes must play a considerable role in the phenotypic heterogeneity of Mendelian disorders. Significant advances in genetic technologies and molecular medicine allow huge amounts of information to be generated from individual samples within a reasonable time frame. This review focuses on the role of modifier genes using the example of cystic fibrosis, the most common lethal autosomal recessive disorder in the white population, and discusses the advantages and limitations of candidate gene approaches versus genome-wide association studies. Moreover, the implications of modifier gene research for other monogenic disorders, as well as its significance for diagnostic, prognostic, and therapeutic approaches are summarized. Increasing insight into modifying mechanisms opens up new perspectives, dispelling the idea of genetic disorders being caused by one single gene.