Use of phage display technology to investigate allergen-antibody interactions

Phage display is an advanced technology that can be used to characterize the interactions of antibody with antigen at the molecular level. It provides valuable data when applied to the investigation of IgE interaction with allergens. The aim of this rostrum article is to provide an explanation of the potential of phage display for increasing the understanding of allergen- IgE interaction, the discovery of diagnostic reagents, and the development of novel therapeutics for the treatment of allergic disease. The significance of initial studies that have applied phage display technology in allergy research will be highlighted. Phage display has been used to clone human IgE to timothy grass pollen allergen Phl p 5, to characterize the epitopes for murine and human antibodies to a birch pollen allergen Bet v 1, and to elucidate the epitopes of a murine mAb to the house dust mite allergen Der p 1. The technology has identified peptides that functionally mimic sites of human IgE constant domains and that were used to raise antiserum for blocking binding of IgE to the FcεRI on basophils and subsequent release of histamine. Phage display has also been used to characterize novel peanut and fungal allergens. The method has been used to increase our understanding of the molecular basis of allergen-IgE interactions and to develop clinically relevant reagents with the pharmacologic potential to block the effector phase of allergic reactions. Many advances from these early studies are likely as phage display technology evolves and allergists gain expertise in its research applications.

Genome-wide association analysis identifies 11 risk variants associated with the asthma with hay fever phenotype

Background To date, no genome-wide association study (GWAS) has considered the combined phenotype of asthma with hay fever. Previous analyses of family data from the Tasmanian Longitudinal Health Study provide evidence that this phenotype has a stronger genetic cause than asthma without hay fever. Objective We sought to perform a GWAS of asthma with hay fever to identify variants associated with having both diseases. Methods We performed a meta-analysis of GWASs comparing persons with both physician-diagnosed asthma and hay fever (n = 6,685) with persons with neither disease (n = 14,091). Results At genome-wide significance, we identified 11 independent variants associated with the risk of having asthma with hay fever, including 2 associations reaching this level of significance with allergic disease for the first time: ZBTB10 (rs7009110; odds ratio [OR], 1.14; P = 4 × 10−9) and CLEC16A (rs62026376; OR, 1.17; P = 1 × 10−8). The rs62026376:C allele associated with increased asthma with hay fever risk has been found to be associated also with decreased expression of the nearby DEXI gene in monocytes. The 11 variants were associated with the risk of asthma and hay fever separately, but the estimated associations with the individual phenotypes were weaker than with the combined asthma with hay fever phenotype. A variant near LRRC32 was a stronger risk factor for hay fever than for asthma, whereas the reverse was observed for variants in/near GSDMA and TSLP. Single nucleotide polymorphisms with suggestive evidence for association with asthma with hay fever risk included rs41295115 near IL2RA (OR, 1.28; P = 5 × 10−7) and rs76043829 in TNS1 (OR, 1.23; P = 2 × 10−6). Conclusion By focusing on the combined phenotype of asthma with hay fever, variants associated with the risk of allergic disease can be identified with greater efficiency.

Genome-wide association study to identify genetic determinants of severe asthma

Background: The genetic basis for developing asthma has been extensively studied. However, association studies to date have mostly focused on mild to moderate disease and genetic risk factors for severe asthma remain unclear. Objective: To identify common genetic variants affecting susceptibility to severe asthma. Methods: A genome-wide association study was undertaken in 933 European ancestry individuals with severe asthma based on Global Initiative for Asthma (GINA) criteria 3 or above and 3346 clean controls. After standard quality control measures, the association of 480 889 genotyped single nucleotide polymorphisms (SNPs) was tested. To improve the resolution of the association signals identified, non-genotyped SNPs were imputed in these regions using a dense reference panel of SNP genotypes from the 1000 Genomes Project. Then replication of SNPs of interest was undertaken in a further 231 cases and 1345 controls and a meta-analysis was performed to combine the results across studies. Results: An association was confirmed in subjects with severe asthma of loci previously identified for association with mild to moderate asthma. The strongest evidence was seen for the ORMDL3/GSDMB locus on chromosome 17q12-21 (rs4794820, p=1.03×10 (-8)following meta-analysis) meeting genome-wide significance. Strong evidence was also found for the IL1RL1/IL18R1 locus on 2q12 (rs9807989, p=5.59×10 (-8) following meta-analysis) just below this threshold. No novel loci for susceptibility to severe asthma met strict criteria for genome-wide significance. Conclusions: The largest genome-wide association study of severe asthma to date was carried out and strong evidence found for the association of two previously identified asthma susceptibility loci in patients with severe disease. A number of novel regions with suggestive evidence were also identified warranting further study.

Fine mapping and replication of genetic risk Loci in primary sclerosing cholangitis

Background and aims. Primary sclerosing cholangitis (PSC) is a chronic cholestatic liver disease characterized by progressive inflammation and fibrosis of the bile ducts eventually leading to biliary cirrhosis. Recent genetic studies in PSC have identified associations at 2q13, 2q35, 3p21, 4q27, 13q31 and suggestive association at 10p15. The aim of this study was to further characterize and refine the genetic architecture of PSC. Methods. We analyzed previously reported associated SNPs at four of these non-HLA loci and 59 SNPs tagging the IL-2/IL-21 (4q27) and IL2RA (10p15) loci in 992 UK PSC cases and 5162 healthy UK controls. Results. The most associated SNPs identified were rs3197999 (3p21 (MST1), p = 1.9 × 10 -6, OR A vs G = 1.28, 95% CI (1.16-1.42)); rs4147359 (10p15 (IL2RA), p = 2.6 × 10 -4, OR A vs G = 1.20, 95% CI (1.09-1.33)) and rs12511287 (4q27 (IL-2/IL-21), p = 3.0 × 10 -4, OR A vs T = 1.21, 95% CI (1.09-1.35)). In addition, we performed a meta-analysis for selected SNPs using published summary statistics from recent studies. We observed genome-wide significance for rs3197999 (3p21 (MST1), P combined = 3.8 × 10 -12) and rs4147359 (10p15 (IL2RA), P combined = 1.5 × 10 -8). Conclusion. We have for the first time confirmed the association of PSC with genetic variants at 10p15 (IL2RA) locus at genome-wide significance and replicated the associations at MST1 and IL-2/IL-21 loci in a large homogeneous UK population. These results strongly implicate the role of IL-2/IL2RA pathway in PSC and provide further confirmation of MST1 association. © Informa Healthcare.

Using genome-Wwide complex trait analysis to quantify 'missing heritability' in Parkinson's disease

Genome-wide association studies (GWASs) have been successful at identifying single-nucleotide polymorphisms (SNPs) highly associated with common traits; however, a great deal of the heritable variation associated with common traits remains unaccounted for within the genome. Genome-wide complex trait analysis (GCTA) is a statistical method that applies a linear mixed model to estimate phenotypic variance of complex traits explained by genome-wide SNPs, including those not associated with the trait in a GWAS. We applied GCTA to 8 cohorts containing 7096 case and 19 455 control individuals of European ancestry in order to examine the missing heritability present in Parkinson's disease (PD). We meta-analyzed our initial results to produce robust heritability estimates for PD types across cohorts. Our results identify 27% (95% CI 17-38, P = 8.08E - 08) phenotypic variance associated with all types of PD, 15% (95% CI -0.2 to 33, P = 0.09) phenotypic variance associated with early-onset PD and 31% (95% CI 17-44, P = 1.34E - 05) phenotypic variance associated with late-onset PD. This is a substantial increase from the genetic variance identified by top GWAS hits alone (between 3 and 5%) and indicates there are substantially more risk loci to be identified. Our results suggest that although GWASs are a useful tool in identifying the most common variants associated with complex disease, a great deal of common variants of small effect remain to be discovered. © Published by Oxford University Press 2012.

Common variants at 12q14 and 12q24 are associated with hippocampal volume

Aging is associated with reductions in hippocampal volume that are accelerated by Alzheimer's disease and vascular risk factors. Our genome-wide association study (GWAS) of dementia-free persons (n = 9,232) identified 46 SNPs at four loci with P values of <4.0 × 10 -7. In two additional samples (n = 2,318), associations were replicated at 12q14 within MSRB3-WIF1 (discovery and replication; rs17178006; P = 5.3 × 10 -11) and at 12q24 near HRK-FBXW8 (rs7294919; P = 2.9 × 10 -11). Remaining associations included one SNP at 2q24 within DPP4 (rs6741949; P = 2.9 × 10 -7) and nine SNPs at 9p33 within ASTN2 (rs7852872; P = 1.0 × 10 -7); along with the chromosome 12 associations, these loci were also associated with hippocampal volume (P < 0.05) in a third younger, more heterogeneous sample (n = 7,794). The SNP in ASTN2 also showed suggestive association with decline in cognition in a largely independent sample (n = 1,563). These associations implicate genes related to apoptosis (HRK), development (WIF1), oxidative stress (MSR3B), ubiquitination (FBXW8) and neuronal migration (ASTN2), as well as enzymes targeted by new diabetes medications (DPP4), indicating new genetic influences on hippocampal size and possibly the risk of cognitive decline and dementia.

Genome-wide association study identifies five new schizophrenia loci

We examined the role of common genetic variation in schizophrenia in a genome-wide association study of substantial size: a stage 1 discovery sample of 21,856 individuals of European ancestry and a stage 2 replication sample of 29,839 independent subjects. The combined stage 1 and 2 analysis yielded genome-wide significant associations with schizophrenia for seven loci, five of which are new (1p21.3, 2q32.3, 8p23.2, 8q21.3 and 10q24.32-q24.33) and two of which have been previously implicated (6p21.32-p22.1 and 18q21.2). The strongest new finding (P = 1.6 × 10 -11) was with rs1625579 within an intron of a putative primary transcript for MIR137 (microRNA 137), a known regulator of neuronal development. Four other schizophrenia loci achieving genome-wide significance contain predicted targets of MIR137, suggesting MIR137-mediated dysregulation as a previously unknown etiologic mechanism in schizophrenia. In a joint analysis with a bipolar disorder sample (16,374 affected individuals and 14,044 controls), three loci reached genome-wide significance: CACNA1C (rs4765905, P = 7.0 × 10 -9), ANK3 (rs10994359, P = 2.5 × 10 -8) and the ITIH3-ITIH4 region (rs2239547, P = 7.8 × 10 -9).

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 ulcerative colitis identifies three new susceptibility loci, including the HNF4A region

Ulcerative colitis is a common form of inflammatory bowel disease with a complex etiology. As part of the Wellcome Trust Case Control Consortium 2, we performed a genome-wide association scan for ulcerative colitis in 2,361 cases and 5,417 controls. Loci showing evidence of association at P 1 × 10 5 were followed up by genotyping in an independent set of 2,321 cases and 4,818 controls. We find genome-wide significant evidence of association at three new loci, each containing at least one biologically relevant candidate gene, on chromosomes 20q13 (HNF4A; P = 3.2 × 10 17), 16q22 (CDH1 and CDH3; P = 2.8 × 10 8) and 7q31 (LAMB1; P = 3.0 × 10 8). Of note, CDH1 has recently been associated with susceptibility to colorectal cancer, an established complication of longstanding ulcerative colitis. The new associations suggest that changes in the integrity of the intestinal epithelial barrier may contribute to the pathogenesis of ulcerative colitis. © 2009 Nature America, Inc. All rights reserved.

Synthesis, X-ray structures, and spectroscopic and electrochemical properties of (mu-oxo)bis(mu-carboxylato)diruthenium complexes having six imidazole bases as terminal ligands

Complexes [Ru2O(O2CR)(2)(1-MeIm)(6)](ClO4)(2) (la-c), [Ru2O(O2CR)(2)(ImH)(6)](ClO4)(2) (2a,b), and [Ru2O(O2CR)(2)(4-MeImH)(6)](ClO4)(2) (3a,b) with a (mu-oxo)bis(mu-carboxylato)diruthenium(III) core have been prepared by reacting Ru2Cl(O2CR)(4) with the corresponding imidazole base, viz. 1-methylimidazole (1-MeIm), imidazole (ImH), and 4-methylimidazole (4-MeImH) in methanol, followed by treatment with NaClO4 in water (R: Me, a; C6H4-p-OMe, b; C6H4-p-Me, c). Diruthenium(III,IV) complexes [Ru2O(O2CR)(2)(1-MeIm)(6)](ClO4)(3) (R: Me, 4a; C6H4-p-OMe, 4b; C6H4-p-Me, 4c) have been prepared by one-electron oxidation of 1 in MeCN with K2S2O8 in water. Complexes la, 2a . 3H(2)O, and 4a . 1.5H(2)O have been structurally characterized. Crystal data for the complexes are as follows: la, orthorhombic, P2(1)2(1)2(1), a = 7.659(3) Angstrom, b = 22.366(3) Angstrom, c = 23.688(2) Angstrom, V = 4058(2) Angstrom(3), Z = 4, R = 0.0475, and R-w = 0.0467 for 2669 reflections with F-o > 2 sigma(F-o); 2a . 3H(2)O, triclinic, , a = 13.735(3) Angstrom, b = 14.428(4) Angstrom, c = 20.515(8) Angstrom, alpha = 87.13(3)degrees, beta = 87.61(3)degrees, gamma = 63.92(2)degrees, V = 3646(2) Angstrom(3), Z = 4, R = 0.0485 and R-w = 0.0583 for 10 594 reflections with F-o > 6 sigma(F-o); 4a . 1.5H(2)O triclinic, , a = 11.969(3) Angstrom, b = 12.090(6) Angstrom, c = 17.421(3) Angstrom, alpha = 108.93(2)degrees, beta = 84.42(2)degrees, gamma = 105.97(2)degrees, V = 2292(1) Angstrom(3), Z = 2, R = 0.0567, and R-w = 0.0705 for 6775 reflections with F-o > 6 sigma(F-o). The complexes have a diruthenium unit held by an oxo and two carboxylate ligands, and the imidazole ligands occupy the terminal sites of the core. The Ru-Ru distance and the Ru-O-oxo-Ru angle in la and 2a . 3H(2)O are 3.266(1), 3.272(1) Angstrom and 122.4(4), 120.5(2)degrees, while in 4a . 1.5H(2)O these values are 3.327(1) Angstrom and 133.6(2)degrees. The diruthenium(III) complexes 1-3 are blue in color and they exhibit an intense visible band in the range 560-575 nm. The absorption is charge transfer in nature involving the Ru(III)-d pi and O-oxo-p pi orbitals. The diruthenium(III,IV) complexes are red in color and show an intense band near 500 nm. The diruthenium(III) core readily gets oxidized with K2S2O8 forming quantitatively the diruthenium(III,IV) complex. The visible spectral record of the conversion shows an isosbestic point at 545 nm for 1 and at 535 nm for 2 and 3. Protonation of the oxide bridge by HClO4 in methanol yields the [Ru-2(mu-OH)(mu-O2CR)(2)](3+) core. The hydroxo species shows a visible band al 550 nm. The pK(a) value for la is 2.45. The protonated species are unstable. The 1-MeIm species converts to the diruthenium(III,IV) core, while the imidazole complex converts to [Ru(ImH)(6)](3+) and some uncharacterized products. Complex [Ru(ImH)(6)](ClO4)(3) has been structurally characterized. The diruthenium(III) complexes are essentially diamagnetic and show characteristic H-1 NMR spectra indicating the presence of the dimeric structure in solution. The diruthenium(III,IV) complexes are paramagnetic and display rhombic EPR spectral features. Complexes 1-3 are redox active. Complex 1 shows the one-electron reversible Ru-2(III)/(RuRuIV)-Ru-III, one-electron quasireversible (RuRuIV)-Ru-III/Ru-2(IV), and two-electron quasireversible Ru-2(III)/Ru-2(II) couples near 0.4, 1.5, and -1.0 V vs SCE In MeCN-0.1 M TBAP, respectively, in the cyclic and differential pulse voltammetric studies. Complexes 2 and 3 exhibit only reversible Ru-2(III)/(RuRuIV)-Ru-III and the quasireversible (RuRuIV)-Ru-III/Ru-2(IV) couples near 0.4 and 1.6 V vs SCE, respectively, The observation of a quasireversible one-step two-electron transfer reduction process in 1 is significant considering its relevance to the rapid and reversible Fe-2(III)/Fe-2(II) redox process known for the tribridged diiron core in the oxy and deoxy forms of hemerythrin.

Reversible Insertion of Methyl Isothiocyanate into Copper(I)Aryloxides

MeNCS undergoes insertion into the copper(I)-aryloxide bond to form [N-methylimino(aryloxy)methanethiolato]-copper(I) complexes. This insertion occurs in the absence of ancillary ligands unlike the analogous insertion of PhNCS. The reaction with 4-methylphenoxide results in the formation of hexakis[[N-methylimino(4-methylphenoxy) methanethiolato]copper(I)] (1), which has been characterized by X-ray crystallography. Crystal data for 1: hexagonal , a = 12.365(3) Angstrom, c = 36.734(16) Angstrom, gamma = 120 degrees, Z = 3, V = 4863(3) Angstrom(3), R = 0.0306. Reactions of 2,6-dimethyl- and 4-chlorophenoxides also result in analogous copper(I) complexes 2 and 3. Addition of stochiometric amounts of PPh(3) to the oligomeric complexes typically results in the extrusion of MeNCS. The ease of extrusion is dependent on the substituents on the aryloxide, and this deinsertion is accelerated by water. However, the extrusion reaction is slow enough in the case of the N-methylimino(2,6-dimethylphenoxy)-methanethiolate complex and the isolation of an intermediate monomeric product bis(triphenylphosphine)[N-methylimino(2, 6-dimethylphenoxy)methanethiolato] copper(I) (4) is possible. Crystal data for 4: triclinic , a = 10.088(2) Angstrom, b = 11.302(1) Angstrom, c = 17.990(2) Angstrom, alpha = 94.06(1)degrees, beta = 95.22(2)degrees, gamma = 103.94(1)degrees, Z = 2, V = 1974.4(7) Angstrom(3), R = 0.0361. In the presence of of PPh(3), the insertion reaction becomes reversible. This allows the exchange of the heterocumulene MeNCS or the aryloxy group in these molecules with another heterocumulene or a phenol, respectively, when catalytic amounts of PPh(3) are added. Oligomers with exchanged heterocumulmes and phenols could be characterized by independent synthesis.

Generalizability of established prostate cancer risk variants in men of African ancestry

Genome-wide association studies have identified more than 80 risk variants for prostate cancer, mainly in European or Asian populations. The generalizability of these variants in other racial/ethnic populations needs to be understood before the loci can be used widely in risk modeling. In our study, we examined 82 previously reported risk variants in 4,853 prostate cancer cases and 4,678 controls of African ancestry. We performed association testing for each variant using logistic regression adjusted for age, study and global ancestry. Of the 82 known risk variants, 68 (83%) had effects that were directionally consistent in their association with prostate cancer risk and 30 (37%) were significantly associated with risk at p < 0.05, with the most statistically significant variants being rs116041037 (p = 3.7 × 10(-26) ) and rs6983561 (p = 1.1 × 10(-16) ) at 8q24, as well as rs7210100 (p = 5.4 × 10(-8) ) at 17q21. By exploring each locus in search of better markers, the number of variants that captured risk in men of African ancestry (p < 0.05) increased from 30 (37%) to 44 (54%). An aggregate score comprised of these 44 markers was strongly associated with prostate cancer risk [per-allele odds ratio (OR) = 1.12, p = 7.3 × 10(-98) ]. In summary, the consistent directions of effects for the vast majority of variants in men of African ancestry indicate common functional alleles that are shared across populations. Further exploration of these susceptibility loci is needed to identify the underlying biologically relevant variants to improve prostate cancer risk modeling in populations of African ancestry.

A note on acyclic edge coloring of complete bipartite graphs

An acyclic edge coloring of a graph is a proper edge coloring such that there are no bichromatic (2-colored) cycles. The acyclic chromatic index of a graph is the minimum number k such that there is an acyclic edge coloring using k colors and is denoted by a'(G). Let Delta = Delta(G) denote the maximum degree of a vertex in a graph G. A complete bipartite graph with n vertices on each side is denoted by K-n,K-n. Alon, McDiarmid and Reed observed that a'(K-p-1,K-p-1) = p for every prime p. In this paper we prove that a'(K-p,K-p) <= p + 2 = Delta + 2 when p is prime. Basavaraju, Chandran and Kummini proved that a'(K-n,K-n) >= n + 2 = Delta + 2 when n is odd, which combined with our result implies that a'(K-p,K-p) = p + 2 = Delta + 2 when p is an odd prime. Moreover we show that if we remove any edge from K-p,K-p, the resulting graph is acyclically Delta + 1 = p + 1-edge-colorable. (C) 2009 Elsevier B.V. All rights reserved.

Nutrient uptake in tropical turfgrasses growing in winter in southern Queensland

The effects of fertilisers on 8 tropical turfgrasses growing in 100-L bags of sand were studied over winter in Murrumba Downs, just north of Brisbane in southern Queensland (latitude 27.4°S, longitude 153.1°E). The species used were: Axonopus compressus (broad-leaf carpetgrass), Cynodon dactylon (bermudagrass 'Winter Green') and C. dactylon x C. transvaalensis hybrid ('Tifgreen'), Digitaria didactyla (Queensland blue couch), Paspalum notatum (bahiagrass '38824'), Stenotaphrum secundatum (buffalograss 'Palmetto'), Eremochloa ophiuroides (centipedegrass 'Centec') and Zoysia japonica (zoysiagrass 'ZT-11'). Control plots were fertilised with complete fertilisers every month from May to September (72 kg N/ha, 31 kg P/ha, 84 kg K/ha, 48 kg S/ha, 30 kg Ca/ha and 7.2 kg Mg/ha), and unfertilised plots received no fertiliser. Carpetgrass and standard bermudagrass were the most sensitive species to nutrient supply, with lower shoot dry weights in the unfertilised plots (shoots mowed to thatch level) compared with the fertilised plots in June. There were lower shoot dry weights in the unfertilised plots in July for all species, except for buffalograss, centipedegrass and zoysiagrass, and lower shoot dry weights in the unfertilised plots in August for all species, except for centipedegrass. At the end of the experiment in September, unfertilised plots were 11% of the shoot dry weights of fertilised plots, with all species affected. Mean shoot nitrogen concentrations fell from 3.2 to 1.7% in the unfertilised plots from May to August, below the sufficiency range for turfgrasses (2.8-3.5%). There were also declines in P (0.45-0.36%), K (2.4-1.5%), S (0.35-0.25%), Mg (0.24-0.18%) and B (9-6 mg/kg), which were all in the sufficiency range. The shoots in the control plots took up the following levels (kg/ha.month) of nutrients: N, 10.0-27.0; P, 1.6-4.0; K, 8.2-19.8; S, 1.0-4.2; Ca, 1.1-3.3; and Mg, 0.8-2.2, compared with applications (kg/ha.month) of: N, 72; P, 31; K, 84; S, 48; Ca, 30; and Mg, 7.2, indicating a recovery of 14-38% for N, 5-13% for P, 10-24% for K, 2-9% for S, 4-11% for Ca and 11-30% for Mg. These results suggest that buffalograss, centipedegrass and zoysiagrass are less sensitive to low nutrient supply than carpetgrass, bermudagrass, blue couch and bahiagrass. Data on nutrient uptake showed that the less sensitive species required only half or less of the nitrogen required to maintain the growth of the other grasses, indicating potential savings for turf managers in fertiliser costs and the environment in terms of nutrients entering waterways.

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.