Genetic variation near IRS1 associates with reduced adiposity and an impaired metabolic profile.

Genome-wide association studies have identified 32 loci influencing body mass index, but this measure does not distinguish lean from fat mass. To identify adiposity loci, we meta-analyzed associations between ∼2.5 million SNPs and body fat percentage from 36,626 individuals and followed up the 14 most significant (P < 10(-6)) independent loci in 39,576 individuals. We confirmed a previously established adiposity locus in FTO (P = 3 × 10(-26)) and identified two new loci associated with body fat percentage, one near IRS1 (P = 4 × 10(-11)) and one near SPRY2 (P = 3 × 10(-8)). Both loci contain genes with potential links to adipocyte physiology. Notably, the body-fat-decreasing allele near IRS1 is associated with decreased IRS1 expression and with an impaired metabolic profile, including an increased visceral to subcutaneous fat ratio, insulin resistance, dyslipidemia, risk of diabetes and coronary artery disease and decreased adiponectin levels. Our findings provide new insights into adiposity and insulin resistance.

Defining the role of common variation in the genomic and biological architecture of adult human height.

Using genome-wide data from 253,288 individuals, we identified 697 variants at genome-wide significance that together explained one-fifth of the heritability for adult height. By testing different numbers of variants in independent studies, we show that the most strongly associated ∼2,000, ∼3,700 and ∼9,500 SNPs explained ∼21%, ∼24% and ∼29% of phenotypic variance. Furthermore, all common variants together captured 60% of heritability. The 697 variants clustered in 423 loci were enriched for genes, pathways and tissue types known to be involved in growth and together implicated genes and pathways not highlighted in earlier efforts, such as signaling by fibroblast growth factors, WNT/β-catenin and chondroitin sulfate-related genes. We identified several genes and pathways not previously connected with human skeletal growth, including mTOR, osteoglycin and binding of hyaluronic acid. Our results indicate a genetic architecture for human height that is characterized by a very large but finite number (thousands) of causal variants.

Large-scale gene-centric analysis identifies novel variants for coronary artery disease.

Coronary artery disease (CAD) has a significant genetic contribution that is incompletely characterized. To complement genome-wide association (GWA) studies, we conducted a large and systematic candidate gene study of CAD susceptibility, including analysis of many uncommon and functional variants. We examined 49,094 genetic variants in ∼2,100 genes of cardiovascular relevance, using a customised gene array in 15,596 CAD cases and 34,992 controls (11,202 cases and 30,733 controls of European descent; 4,394 cases and 4,259 controls of South Asian origin). We attempted to replicate putative novel associations in an additional 17,121 CAD cases and 40,473 controls. Potential mechanisms through which the novel variants could affect CAD risk were explored through association tests with vascular risk factors and gene expression. We confirmed associations of several previously known CAD susceptibility loci (eg, 9p21.3:p<10(-33); LPA:p<10(-19); 1p13.3:p<10(-17)) as well as three recently discovered loci (COL4A1/COL4A2, ZC3HC1, CYP17A1:p<5×10(-7)). However, we found essentially null results for most previously suggested CAD candidate genes. In our replication study of 24 promising common variants, we identified novel associations of variants in or near LIPA, IL5, TRIB1, and ABCG5/ABCG8, with per-allele odds ratios for CAD risk with each of the novel variants ranging from 1.06-1.09. Associations with variants at LIPA, TRIB1, and ABCG5/ABCG8 were supported by gene expression data or effects on lipid levels. Apart from the previously reported variants in LPA, none of the other ∼4,500 low frequency and functional variants showed a strong effect. Associations in South Asians did not differ appreciably from those in Europeans, except for 9p21.3 (per-allele odds ratio: 1.14 versus 1.27 respectively; P for heterogeneity = 0.003). This large-scale gene-centric analysis has identified several novel genes for CAD that relate to diverse biochemical and cellular functions and clarified the literature with regard to many previously suggested genes.

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.

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.

ELOVL6 Genetic Variation Is Related to Insulin Sensitivity: A New Candidate Gene in Energy Metabolism

BACKGROUND: The elongase of long chain fatty acids family 6 (ELOVL6) is an enzyme that specifically catalyzes the elongation of saturated and monounsaturated fatty acids with 12, 14 and 16 carbons. ELOVL6 is expressed in lipogenic tissues and it is regulated by sterol regulatory element binding protein 1 (SREBP-1). OBJECTIVE: We investigated whether ELOVL6 genetic variation is associated with insulin sensitivity in a population from southern Spain. DESIGN: We undertook a prospective, population-based study collecting phenotypic, metabolic, nutritional and genetic information. Measurements were made of weight and height and the body mass index (BMI) was calculated. Insulin resistance was measured by homeostasis model assessment. The type of dietary fat was assessed from samples of cooking oil taken from the participants' kitchens and analyzed by gas chromatography. Five SNPs of the ELOVL6 gene were analyzed by SNPlex. RESULTS: Carriers of the minor alleles of the SNPs rs9997926 and rs6824447 had a lower risk of having high HOMA_IR, whereas carriers of the minor allele rs17041272 had a higher risk of being insulin resistant. An interaction was detected between the rs6824447 polymorphism and the intake of oil in relation with insulin resistance, such that carriers of this minor allele who consumed sunflower oil had lower HOMA_IR than those who did not have this allele (P = 0.001). CONCLUSIONS: Genetic variations in the ELOVL6 gene were associated with insulin sensitivity in this population-based study.

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.

Analysis of TNFAIP3, a feedback inhibitor of nuclear factor-κB and the neighbor intergenic 6q23 region in rheumatoid arthritis susceptibility

INTRODUCTION Genome-wide association studies of rheumatoid arthritis (RA) have identified an association of the disease with a 6q23 region devoid of genes. TNFAIP3, an RA candidate gene, flanks this region, and polymorphisms in both the TNFAIP3 gene and the intergenic region are associated with systemic lupus erythematosus. We hypothesized that there is a similar association with RA, including polymorphisms in TNFAIP3 and the intergenic region. METHODS To test this hypothesis, we selected tag-single nucleotide polymorphisms (SNPs) in both loci. They were analyzed in 1,651 patients with RA and 1,619 control individuals of Spanish ancestry. RESULTS Weak evidence of association was found both in the 6q23 intergenic region and in the TNFAIP3 locus. The rs582757 SNP and a common haplotype in the TNFAIP3 locus exhibited association with RA. In the intergenic region, two SNPs were associated, namely rs609438 and rs13207033. The latter was only associated in patients with anti-citrullinated peptide antibodies. Overall, statistical association was best explained by the interdependent contribution of SNPs from the two loci TNFAIP3 and the 6q23 intergenic region. CONCLUSIONS Our data are consistent with the hypothesis that several RA genetic factors exist in the 6q23 region, including polymorphisms in the TNFAIP3 gene, like that previously described for systemic lupus erythematosus.

Genetic polymorphisms of RANTES, IL1-A, MCP-1 and TNF-A genes in patients with prostate cancer

BACKGROUND Inflammation has been implicated as an etiological factor in several human cancers, including prostate cancer. Allelic variants of the genes involved in inflammatory pathways are logical candidates as genetic determinants of prostate cancer risk. The purpose of this study was to investigate whether single nucleotide polymorphisms of genes that lead to increased levels of pro-inflammatory cytokines and chemokines are associated with an increased prostate cancer risk. METHODS A case-control study design was used to test the association between prostate cancer risk and the polymorphisms TNF-A-308 A/G (rs 1800629), RANTES-403 G/A (rs 2107538), IL1-A-889 C/T (rs 1800587) and MCP-1 2518 G/A (rs 1024611) in 296 patients diagnosed with prostate cancer and in 311 healthy controls from the same area. RESULTS Diagnosis of prostate cancer was significantly associated with TNF-A GA + AA genotype (OR, 1.61; 95% CI, 1.09-2.64) and RANTES GA + AA genotype (OR, 1.44; 95% CI, 1.09-2.38). A alleles in TNF-A and RANTES influenced prostate cancer susceptibility and acted independently of each other in these subjects. No epistatic effect was found for the combination of different polymorphisms studied. Finally, no overall association was found between prostate cancer risk and IL1-A or MCP-1 polymorphisms. CONCLUSION Our results and previously published findings on genes associated with innate immunity support the hypothesis that polymorphisms in proinflammatory genes may be important in prostate cancer development.

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.

Mutation analysis of genes that control the G1/S cell cycle in melanoma: TP53, CDKN1A, CDKN2A, and CDKN2B.

BACKGROUND The role of genes involved in the control of progression from the G1 to the S phase of the cell cycle in melanoma tumors in not fully known. The aim of our study was to analyse mutations in TP53, CDKN1A, CDKN2A, and CDKN2B genes in melanoma tumors and melanoma cell lines METHODS We analysed 39 primary and metastatic melanomas and 9 melanoma cell lines by single-stranded conformational polymorphism (SSCP). RESULTS The single-stranded technique showed heterozygous defects in the TP53 gene in 8 of 39 (20.5%) melanoma tumors: three new single point mutations in intronic sequences (introns 1 and 2) and exon 10, and three new single nucleotide polymorphisms located in introns 1 and 2 (C to T transition at position 11701 in intron 1; C insertion at position 11818 in intron 2; and C insertion at position 11875 in intron 2). One melanoma tumor exhibited two heterozygous alterations in the CDKN2A exon 1 one of which was novel (stop codon, and missense mutation). No defects were found in the remaining genes. CONCLUSION These results suggest that these genes are involved in melanoma tumorigenesis, although they may be not the major targets. Other suppressor genes that may be informative of the mechanism of tumorigenesis in skin melanomas should be studied.

Methods for testing association between uncertain genotypes and quantitative traits.

Interpretability and power of genome-wide association studies can be increased by imputing unobserved genotypes, using a reference panel of individuals genotyped at higher marker density. For many markers, genotypes cannot be imputed with complete certainty, and the uncertainty needs to be taken into account when testing for association with a given phenotype. In this paper, we compare currently available methods for testing association between uncertain genotypes and quantitative traits. We show that some previously described methods offer poor control of the false-positive rate (FPR), and that satisfactory performance of these methods is obtained only by using ad hoc filtering rules or by using a harsh transformation of the trait under study. We propose new methods that are based on exact maximum likelihood estimation and use a mixture model to accommodate nonnormal trait distributions when necessary. The new methods adequately control the FPR and also have equal or better power compared to all previously described methods. We provide a fast software implementation of all the methods studied here; our new method requires computation time of less than one computer-day for a typical genome-wide scan, with 2.5 M single nucleotide polymorphisms and 5000 individuals.

The T-381C SNP in BNP gene may be modestly associated with type 2 diabetes: an updated meta-analysis in 49 279 subjects.

A recent study reported an association between the brain natriuretic peptide (BNP) promoter T-381C polymorphism (rs198389) and protection against type 2 diabetes (T2D). As replication in several studies is mandatory to confirm genetic results, we analyzed the T-381C polymorphism in seven independent case-control cohorts and in 291 T2D-enriched pedigrees totalling 39 557 subjects of European origin. A meta-analysis of the seven case-control studies (n = 39 040) showed a nominal protective effect [odds ratio (OR) = 0.86 (0.79-0.94), P = 0.0006] of the CC genotype on T2D risk, consistent with the previous study. By combining all available data (n = 49 279), we further confirmed a modest contribution of the BNP T-381C polymorphism for protection against T2D [OR = 0.86 (0.80-0.92), P = 1.4 x 10(-5)]. Potential confounders such as gender, age, obesity status or family history were tested in 4335 T2D and 4179 normoglycemic subjects and they had no influence on T2D risk. This study provides further evidence of a modest contribution of the BNP T-381C polymorphism in protection against T2D and illustrates the difficulty of unambiguously proving modest-sized associations even with large sample sizes.

Genome-wide linkage and association analyses to identify genes influencing adiponectin levels: the GEMS Study.

Adiponectin has a variety of metabolic effects on obesity, insulin sensitivity, and atherosclerosis. To identify genes influencing variation in plasma adiponectin levels, we performed genome-wide linkage and association scans of adiponectin in two cohorts of subjects recruited in the Genetic Epidemiology of Metabolic Syndrome Study. The genome-wide linkage scan was conducted in families of Turkish and southern European (TSE, n = 789) and Northern and Western European (NWE, N = 2,280) origin. A whole genome association (WGA) analysis (500K Affymetrix platform) was carried out in a set of unrelated NWE subjects consisting of approximately 1,000 subjects with dyslipidemia and 1,000 overweight subjects with normal lipids. Peak evidence for linkage occurred at chromosome 8p23 in NWE subjects (lod = 3.10) and at chromosome 3q28 near ADIPOQ, the adiponectin structural gene, in TSE subjects (lod = 1.70). In the WGA analysis, the single-nucleotide polymorphisms (SNPs) most strongly associated with adiponectin were rs3774261 and rs6773957 (P < 10(-7)). These two SNPs were in high linkage disequilibrium (r(2) = 0.98) and located within ADIPOQ. Interestingly, our fourth strongest region of association (P < 2 x 10(-5)) was to an SNP within CDH13, whose protein product is a newly identified receptor for high-molecular-weight species of adiponectin. Through WGA analysis, we confirmed previous studies showing SNPs within ADIPOQ to be strongly associated with variation in adiponectin levels and further observed these to have the strongest effects on adiponectin levels throughout the genome. We additionally identified a second gene (CDH13) possibly influencing variation in adiponectin levels. The impact of these SNPs on health and disease has yet to be determined.