Genome-wide association study identifies susceptibility loci for open angle glaucoma at TMCO1 and CDKN2B-AS1

We report a genome-wide association study for open-angle glaucoma (OAG) blindness using a discovery cohort of 590 individuals with severe visual field loss (cases) and 3,956 controls. We identified associated loci at TMCO1 (rs4656461[G] odds ratio (OR) = 1.68, P = 6.1 × 10-10) and CDKN2B-AS1 (rs4977756[A] OR = 1.50, P = 4.7 × 10-9). We replicated these associations in an independent cohort of cases with advanced OAG (rs4656461 P = 0.010; rs4977756 P = 0.042) and two additional cohorts of less severe OAG (rs4656461 combined discovery and replication P = 6.00 × 10-14, OR = 1.51, 95% CI 1.35-1.68; rs4977756 combined P = 1.35 × 10-14, OR = 1.39, 95% CI 1.28-1.51). We show retinal expression of genes at both loci in human ocular tissues. We also show that CDKN2A and CDKN2B are upregulated in the retina of a rat model of glaucoma. © 2011 Nature America, Inc. All rights reserved.

A genome-wide asociation study identifies new psoriasis susceptibility loci and an interaction betwEn HLA-C and ERAP1

To identify new susceptibility loci for psoriasis, we undertOk a genome-wide asociation study of 594,224 SNPs in 2,622 individuals with psoriasis and 5,667 controls. We identified asociations at eight previously unreported genomic loci. Seven loci harbored genes with recognized iMune functions (IL28RA, REL, IFIH1, ERAP1, TRAF3IP2, NFKBIA and TYK2). These asociations were replicated in 9,079 European samples (six loci with a combined P < 5-10 -8 and two loci with a combined P < 5-10-7). We also report compeLing evidence for an interaction betwEn the HLA-C and ERAP1 loci (combined P = 6.95-10-6). ERAP1 plays an important role in MHC claS I peptide proceSing. ERAP1 variants only influenced psoriasis susceptibility in individuals carrying the HLA-C risk aLele. Our findings implicate pathways that integrate epidermal barrier dysfunction with iNate and adaptive iMune dysregulation in psoriasis pathogenesis.

Genome-wide association study of ankylosing spondylitis identifies non-MHC susceptibility loci

To identify susceptibility loci for ankylosing spondylitis, we undertook a genome-wide association study in 2,053 unrelated ankylosing spondylitis cases among people of European descent and 5,140 ethnically matched controls, with replication in an independent cohort of 898 ankylosing spondylitis cases and 1,518 controls. Cases were genotyped with Illumina HumHap370 genotyping chips. In addition to strong association with the major histocompatibility complex (MHC; P 10 800), we found association with SNPs in two gene deserts at 2p15 (rs10865331; combined P = 1.9 × 10 19) and 21q22 (rs2242944; P = 8.3 × 10 20), as well as in the genes ANTXR2 (rs4333130; P = 9.3 × 10 8) and IL1R2 (rs2310173; P = 4.8 × 10 7). We also replicated previously reported associations at IL23R (rs11209026; P = 9.1 × 10 14) and ERAP1 (rs27434; P = 5.3 × 10 12). This study reports four genetic loci associated with ankylosing spondylitis risk and identifies a major role for the interleukin (IL)-23 and IL-1 cytokine pathways in disease susceptibility. © 2010 Nature America, Inc. All rights reserved.

Genome-wide association study of CNVs in 16,000 cases of eight common diseases and 3,000 shared controls

Copy number variants (CNVs) account for a major proportion of human genetic polymorphism and have been predicted to have an important role in genetic susceptibility to common disease. To address this we undertook a large, direct genome-wide study of association between CNVs and eight common human diseases. Using a purpose-designed array we typed 19,000 individuals into distinct copy-number classes at 3,432 polymorphic CNVs, including an estimated 50% of all common CNVs larger than 500 base pairs. We identified several biological artefacts that lead to false-positive associations, including systematic CNV differences between DNAs derived from blood and cell lines. Association testing and follow-up replication analyses confirmed three loci where CNVs were associated with diseaseIRGM for Crohns disease, HLA for Crohns disease, rheumatoid arthritis and type 1 diabetes, and TSPAN8 for type 2 diabetesalthough in each case the locus had previously been identified in single nucleotide polymorphism (SNP)-based studies, reflecting our observation that most common CNVs that are well-typed on our array are well tagged by SNPs and so have been indirectly explored through SNP studies. We conclude that common CNVs that can be typed on existing platforms are unlikely to contribute greatly to the genetic basis of common human diseases. © 2010 Macmillan Publishers Limited. All rights reserved.

Genome-wide association study identifies eight loci associated with blood pressure

Elevated blood pressure is a common, heritable cause of cardiovascular disease worldwide. To date, identification of common genetic variants influencing blood pressure has proven challenging. We tested 2.5 million genotyped and imputed SNPs for association with systolic and diastolic blood pressure in 34,433 subjects of European ancestry from the Global BPgen consortium and followed up findings with direct genotyping (N 71,225 European ancestry, N 12,889 Indian Asian ancestry) and in silico comparison (CHARGE consortium, N = 29,136). We identified association between systolic or diastolic blood pressure and common variants in eight regions near the CYP17A1 (P = 7 × 10 24), CYP1A2 (P = 1 × 10 23), FGF5 (P = 1 × 10 21), SH2B3 (P = 3 × 10 18), MTHFR (P = 2 × 10 13), c10orf107 (P = 1 × 10 9), ZNF652 (P = 5 × 10 9) and PLCD3 (P = 1 × 10 8) genes. All variants associated with continuous blood pressure were associated with dichotomous hypertension. These associations between common variants and blood pressure and hypertension offer mechanistic insights into the regulation of blood pressure and may point to novel targets for interventions to prevent cardiovascular disease.

Genome-wide association studies and musculoskeletal diseases

Bone and joint diseases are major causes of morbidity and mortality worldwide, and their prevalence is increasing as the average population age increases. Most common musculoskeletal diseases show significant heritability, and few have treatments that prevent disease or can induce true treatment-free, disease-free remission. Furthermore, despite valiant efforts of hypothesis-driven research, our understanding of the etiopathogenesis of these conditions is, with few exceptions, at best moderate. Therefore, there has been a long-standing interest in genetics research in musculoskeletal disease as a hypothesis-free method for investigating disease etiopathogenesis. Important contributions have been made through the identification of monogenic causes of disease, but the holy grail of human genetics research has been the identification of the genes responsible for common diseases. The development of genome-wide association (GWA) studies has revolutionized this field, and led to an explosion in the number of genes identified that are definitely involved in musculoskeletal disease pathogenesis. However, this approach will not identify all common disease genes, and although the current progress is exciting and proves the potential of this research discipline, other approaches will be required to identify many of the types of genetic variation likely to be involved.

Completion of the genome sequence of Lettuce necrotic yellows virus, type species of the genus Cytorhabdovirus

We completed the genome sequence of Lettuce necrotic yellows virus (LNYV) by determining the nucleotide sequences of the 4a (putative phosphoprotein), 4b, M (matrix protein), G (glycoprotein) and L (polymerase) genes. The genome consists of 12,807 nucleotides and encodes six genes in the order 3′ leader-N-4a(P)-4b-M-G-L-5′ trailer. Sequences were derived from clones of a cDNA library from LNYV genomic RNA and from fragments amplified using reverse transcription-polymerase chain reaction. The 4a protein has a low isoelectric point characteristic for rhabdovirus phosphoproteins. The 4b protein has significant sequence similarities with the movement proteins of capillo- and trichoviruses and may be involved in cell-to-cell movement. The putative G protein sequence contains a predicted 25 amino acids signal peptide and endopeptidase cleavage site, three predicted glycosylation sites and a putative transmembrane domain. The deduced L protein sequence shows similarities with the L proteins of other plant rhabdoviruses and contains polymerase module motifs characteristic for RNA-dependent RNA polymerases of negative-strand RNA viruses. Phylogenetic analysis of this motif among rhabdoviruses placed LNYV in a group with other sequenced cytorhabdoviruses, most closely related to Strawberry crinkle virus.

Decoration of Pasteurella multocida lipopolysaccharide with phosphocholine is important for virulence

Phosphocholine (PCho) is an important substituent of surface structures expressed by a number of bacterial pathogens. Its role in virulence has been investigated in several species, in which it has been shown to play a role in bacterial adhesion to mucosal surfaces, in resistance to antimicrobial peptides, or in sensitivity to complement-mediated killing. The lipopolysaccharide (LPS) structure of Pasteurella multocida strain Pm70, whose genome sequence is known, has recently been determined and does not contain PCho. However, LPS structures from the closely related, virulent P. multocida strains VP161 and X-73 were shown to contain PCho on their terminal galactose sugar residues. To determine if PCho was involved in the virulence of P. multocida, we used subtractive hybridization of the VP161 genome against the Pm70 genome to identify a four-gene locus (designated pcgDABC) which we show is required for the addition of the PCho residues to LPS. The proteins predicted to be encoded by pcgABC showed identity to proteins involved in choline uptake, phosphorylation, and nucleotide sugar activation of PCho. We constructed a P. multocida VP161 pcgC mutant and demonstrated that this strain produces LPS that lacks PCho on the terminal galactose residues. This pcgC mutant displayed reduced in vivo growth in a chicken infection model and was more sensitive to the chicken antimicrobial peptide fowlicidin-1 than the wild-type P. multocida strain

A genome-wide association study of tick burden and milk composition in cattle.

To study the genetic basis of tick burden and milk production and their interrelationship, we collected a sample of 1961 cattle with multiple tick counts from northern Australia of which 973 had dairy production data in the Australian Dairy Herd Information Service database. We calculated heritabilities, genetic and phenotypic correlations for these traits and showed a negative relationship between tick counts and milk and milk component yield. Tests of polymorphisms of four genes associated with milk yield, ABCG2, DGAT1, GHR and PRLR, showed no statistically significant effect on tick burden but highly significant associations to milk component yield in these data and we confirmed separate effects for GHR and PRLR on bovine chromosome 20. To begin to identify some of the molecular genetic bases for these traits, we genotyped a sample of 189 of these cattle for 7397 single nucleotide polymorphisms in a genome-wide association study. Although the allele effects for adjusted milk fat and protein yield were highly correlated (r = 0.66), the correlations of allele effects of these milk component yields and tick burden were small (|r| <= 0.10). These results agree in general with the phenotypic correlations between tick counts and milk component yield and suggest that selection on markers for tick burden or milk component yield may have no undesirable effect on the other trait.

Genome-wide mapping of cystitis due to Streptococcus agalactiae and Escherichia coli in mice identifies a unique bladder transcriptome that signifies pathogen-specific antimicrobial defense against urinary tract infection

The most common causes of urinary tract infections (UTIs) are Gram-negative pathogens such as Escherichia coli; however, Gram-positive organisms including Streptococcus agalactiae, or group B streptococcus (GBS), also cause UTI. In GBS infection, UTI progresses to cystitis once the bacteria colonize bladder, but the host responses triggered in the bladder immediately following infection are largely unknown. Here, we used genome-wide expression profiling to map the bladder transcriptome of GBS UTI in mice infected transurethrally with uropathogenic GBS that was cultured from a 35 year-old women with cystitis. RNA from bladders was applied to Affymetrix Gene-1.0ST microarrays; qRT-PCR was used to analyze selected gene responses identified in array datasets. A surprisingly small significant gene list of 172 genes was identified at 24h; this compared to 2507 genes identified in a side-by-side comparison with uropathogenic E. coli (UPEC). No genes exhibited significantly altered expression at 2h in GBS-infected mice according to arrays despite high bladder bacterial loads at this early time point. The absence of a marked early host response to GBS juxtaposed with broad-based bladder responses activated by UPEC at 2h. Bioinformatics analyses including integrative systems-level network mapping revealed multiple activated biological pathways in the GBS cystitis transcriptome that regulate leukocyte activation, inflammation, apoptosis, and cytokine-chemokine biosynthesis. These findings define a novel, minimalistic type of bladder host response triggered by GBS UTI, which comprises collective antimicrobial pathways that differ dramatically from those activated by UPEC. Overall, this study emphasizes the unique nature of bladder immune activation mechanisms triggered by distinct uropathogens.

Genome-wide association study identifies novel genetic variants contributing to variation in blood metabolite levels

Metabolites are small molecules involved in cellular metabolism, which can be detected in biological samples using metabolomic techniques. Here we present the results of genome-wide association and meta-analyses for variation in the blood serum levels of 129 metabolites as measured by the Biocrates metabolomic platform. In a discovery sample of 7,478 individuals of European descent, we find 4,068 genome- and metabolome-wide significant (Z-test, P<1.09 × 10−9) associations between single-nucleotide polymorphisms (SNPs) and metabolites, involving 59 independent SNPs and 85 metabolites. Five of the fifty-nine independent SNPs are new for serum metabolite levels, and were followed-up for replication in an independent sample (N=1,182). The novel SNPs are located in or near genes encoding metabolite transporter proteins or enzymes (SLC22A16, ARG1, AGPS and ACSL1) that have demonstrated biomedical or pharmaceutical importance. The further characterization of genetic influences on metabolic phenotypes is important for progress in biological and medical research.

Genome-wide meta-analysis identifies five new susceptibility loci for cutaneous malignant melanoma

Thirteen common susceptibility loci have been reproducibly associated with cutaneous malignant melanoma (CMM). We report the results of an international 2-stage meta-analysis of CMM genome-wide association studies (GWAS). This meta-analysis combines 11 GWAS (5 previously unpublished) and a further three stage 2 data sets, totaling 15,990 CMM cases and 26,409 controls. Five loci not previously associated with CMM risk reached genome-wide significance (P < 5 × 10−8), as did 2 previously reported but unreplicated loci and all 13 established loci. Newly associated SNPs fall within putative melanocyte regulatory elements, and bioinformatic and expression quantitative trait locus (eQTL) data highlight candidate genes in the associated regions, including one involved in telomere biology.

A mega-analysis of genome-wide association studies for major depressive disorder

Prior genome-wide association studies (GWAS) of major depressive disorder (MDD) have met with limited success. We sought to increase statistical power to detect disease loci by conducting a GWAS mega-analysis for MDD. In the MDD discovery phase, we analyzed more than 1.2 million autosomal and X chromosome single-nucleotide polymorphisms (SNPs) in 18 759 independent and unrelated subjects of recent European ancestry (9240 MDD cases and 9519 controls). In the MDD replication phase, we evaluated 554 SNPs in independent samples (6783 MDD cases and 50 695 controls). We also conducted a cross-disorder meta-analysis using 819 autosomal SNPs with P<0.0001 for either MDD or the Psychiatric GWAS Consortium bipolar disorder (BIP) mega-analysis (9238 MDD cases/8039 controls and 6998 BIP cases/7775 controls). No SNPs achieved genome-wide significance in the MDD discovery phase, the MDD replication phase or in pre-planned secondary analyses (by sex, recurrent MDD, recurrent early-onset MDD, age of onset, pre-pubertal onset MDD or typical-like MDD from a latent class analyses of the MDD criteria). In the MDD-bipolar cross-disorder analysis, 15 SNPs exceeded genome-wide significance (P<5 x 10(-8)), and all were in a 248 kb interval of high LD on 3p21.1 (chr3:52 425 083-53 822 102, minimum P=5.9 x 10(-9) at rs2535629). Although this is the largest genome-wide analysis of MDD yet conducted, its high prevalence means that the sample is still underpowered to detect genetic effects typical for complex traits. Therefore, we were unable to identify robust and replicable findings. We discuss what this means for genetic research for MDD. The 3p21.1 MDD-BIP finding should be interpreted with caution as the most significant SNP did not replicate in MDD samples, and genotyping in independent samples will be needed to resolve its status.

Genetic relationship between five psychiatric disorders estimated from genome-wide SNPs

Most psychiatric disorders are moderately to highly heritable. The degree to which genetic variation is unique to individual disorders or shared across disorders is unclear. To examine shared genetic etiology, we use genome-wide genotype data from the Psychiatric Genomics Consortium (PGC) for cases and controls in schizophrenia, bipolar disorder, major depressive disorder, autism spectrum disorders (ASD) and attention-deficit/hyperactivity disorder (ADHD). We apply univariate and bivariate methods for the estimation of genetic variation within and covariation between disorders. SNPs explained 17-29% of the variance in liability. The genetic correlation calculated using common SNPs was high between schizophrenia and bipolar disorder (0.68 +/- 0.04 s.e.), moderate between schizophrenia and major depressive disorder (0.43 +/- 0.06 s.e.), bipolar disorder and major depressive disorder (0.47 +/- 0.06 s.e.), and ADHD and major depressive disorder (0.32 +/- 0.07 s.e.), low between schizophrenia and ASD (0.16 +/- 0.06 s.e.) and non-significant for other pairs of disorders as well as between psychiatric disorders and the negative control of Crohn's disease. This empirical evidence of shared genetic etiology for psychiatric disorders can inform nosology and encourages the investigation of common pathophysiologies for related disorders.

A genome-wide association study of sleep habits and insomnia

Several aspects of sleep behavior such as timing, duration and quality have been demonstrated to be heritable. To identify common variants that influence sleep traits in the population, we conducted a genome-wide association study of six sleep phenotypes assessed by questionnaire in a sample of 2,323 individuals from the Australian Twin Registry. Genotyping was performed on the Illumina 317, 370, and 610K arrays and the SNPs in common between platforms were used to impute non-genotyped SNPs. We tested for association with more than 2,000,000 common polymorphisms across the genome. While no SNPs reached the genome-wide significance threshold, we identified a number of associations in plausible candidate genes. Most notably, a group of SNPs in the third intron of the CACNA1C gene ranked as most significant in the analysis of sleep latency (P = 1.3 x 10(-)(6)). We attempted to replicate this association in an independent sample from the Chronogen Consortium (n = 2,034), but found no evidence of association (P = 0.73). We have identified several other suggestive associations that await replication in an independent sample. We did not replicate the results from previous genome-wide analyses of self-reported sleep phenotypes after correction for multiple testing.