Meta-Analysis of Genome-Wide Association Studies Identifies Six New Loci for Serum Calcium Concentrations

Calcium is vital to the normal functioning of multiple organ systems and its serum concentration is tightly regulated. Apart from CASR, the genes associated with serum calcium are largely unknown. We conducted a genome-wide association meta-analysis of 39,400 individuals from 17 population-based cohorts and investigated the 14 most strongly associated loci in ≤21,679 additional individuals. Seven loci (six new regions) in association with serum calcium were identified and replicated. Rs1570669 near CYP24A1 (P = 9.1E-12), rs10491003 upstream of GATA3 (P = 4.8E-09) and rs7481584 in CARS (P = 1.2E-10) implicate regions involved in Mendelian calcemic disorders: Rs1550532 in DGKD (P = 8.2E-11), also associated with bone density, and rs7336933 near DGKH/KIAA0564 (P = 9.1E-10) are near genes that encode distinct isoforms of diacylglycerol kinase. Rs780094 is in GCKR. We characterized the expression of these genes in gut, kidney, and bone, and demonstrate modulation of gene expression in bone in response to dietary calcium in mice. Our results shed new light on the genetics of calcium homeostasis.     

Genome-Wide Association Scan Meta-Analysis Identifies Three Loci Influencing Adiposity and Fat Distribution

To identify genetic loci influencing central obesity and fat distribution, we performed a meta-analysis of 16 genome-wide association studies (GWAS, N=38,580) informative for adult waist circumference (WC) and waist–hip ratio (WHR). We selected 26 SNPs for follow-up, for which the evidence of association with measures of central adiposity (WC and/or WHR) was strong and disproportionate to that for overall adiposity or height. Follow-up studies in a maximum of 70,689 individuals identified two loci strongly associated with measures of central adiposity; these map near TFAP2B (WC, P=1.9×10⁻¹¹) and MSRA (WC, P=8.9×10⁻⁹). A third locus, near LYPLAL1, was associated with WHR in women only (P=2.6×10⁻⁸). The variants near TFAP2B appear to influence central adiposity through an effect on overall obesity/fat-mass, whereas LYPLAL1 displays a strong female-only association with fat distribution. By focusing on anthropometric measures of central obesity and fat distribution, we have identified three loci implicated in the regulation of human adiposity.     

Meta-Analysis of Genome-Wide Scans for Human Adult Stature Identifies Novel Loci and Associations with Measures of Skeletal Frame Size

Recent genome-wide (GW) scans have identified several independent loci affecting human stature, but their contribution through the different skeletal components of height is still poorly understood. We carried out a genome-wide scan in 12,611 participants, followed by replication in an additional 7,187 individuals, and identified 17 genomic regions with GW-significant association with height. Of these, two are entirely novel (rs11809207 in CATSPER4, combined P-value=6.1×10⁻⁸ and rs910316 in TMED10, P-value=1.4×10⁻⁷) and two had previously been described with weak statistical support (rs10472828 in NPR3, P-value=3×10⁻⁷ and rs849141 in JAZF1, P-value=3.2×10⁻¹¹). One locus (rs1182188 at GNA12) identifies the first height eQTL. We also assessed the contribution of height loci to the upper- (trunk) and lower-body (hip axis and femur) skeletal components of height. We find evidence for several loci associated with trunk length (including rs6570507 in GPR126, P-value=4×10⁻⁵ and rs6817306 in LCORL, P-value=4×10⁻⁴), hip axis length (including rs6830062 at LCORL, P-value=4.8×10⁻⁴ and rs4911494 at UQCC, P-value=1.9×10⁻⁴), and femur length (including rs710841 at PRKG2, P-value=2.4×10⁻⁵ and rs10946808 at HIST1H1D, P-value=6.4×10⁻⁶). Finally, we used conditional analyses to explore a possible differential contribution of the height loci to these different skeletal size measurements. In addition to validating four novel loci controlling adult stature, our study represents the first effort to assess the contribution of genetic loci to three skeletal components of height. Further statistical tests in larger numbers of individuals will be required to verify if the height loci affect height preferentially through these subcomponents of height.     

Genome-Wide Association Study Identifies Novel Pharmacogenomic Loci For Therapeutic Response to Montelukast in Asthma

Background

Genome-wide association study (GWAS) is a powerful tool to identify novel pharmacogenetic single nucleotide polymorphisms (SNPs). Leukotriene receptor antagonists (LTRAs) are a major class of asthma medications, and genetic factors contribute to variable responses to these drugs. We used GWAS to identify novel SNPs associated with the response to the LTRA, montelukast, in asthmatics.

Methods

Using genome-wide genotype and phenotypic data available from American Lung Association - Asthma Clinical Research Center (ALA-ACRC) cohorts, we evaluated 8-week change in FEV₁ related to montelukast administration in a discovery population of 133 asthmatics. The top 200 SNPs from the discovery GWAS were then tested in 184 additional samples from two independent cohorts.

Results

Twenty-eight SNP associations from the discovery GWAS were replicated. Of these, rs6475448 achieved genome-wide significance (combined P = 1.97 x 10^-09), and subjects from all four studies who were homozygous for rs6475448 showed increased ΔFEV₁ from baseline in response to montelukast.

Conclusions

Through GWAS, we identified a novel pharmacogenomic locus related to improved montelukast response in asthmatics.     

Genome-Wide Linkage Scan Identifies Two Novel Genetic Loci for Coronary Artery Disease: In GeneQuest Families

Coronary artery disease (CAD) is the leading cause of death worldwide. Recent genome-wide association studies (GWAS) identified >50 common variants associated with CAD or its complication myocardial infarction (MI), but collectively they account for <20% of heritability, generating a phenomena of “missing heritability”. Rare variants with large effects may account for a large portion of missing heritability. Genome-wide linkage studies of large families and follow-up fine mapping and deep sequencing are particularly effective in identifying rare variants with large effects. Here we show results from a genome-wide linkage scan for CAD in multiplex GeneQuest families with early onset CAD and MI. Whole genome genotyping was carried out with 408 markers that span the human genome by every 10 cM and linkage analyses were performed using the affected relative pair analysis implemented in GENEHUNTER. Affected only nonparametric linkage (NPL) analysis identified two novel CAD loci with highly significant evidence of linkage on chromosome 3p25.1 (peak NPL  = 5.49) and 3q29 (NPL  = 6.84). We also identified four loci with suggestive linkage on 9q22.33, 9q34.11, 17p12, and 21q22.3 (NPL  = 3.18–4.07). These results identify novel loci for CAD and provide a framework for fine mapping and deep sequencing to identify new susceptibility genes and novel variants associated with risk of CAD.     

Genome-Wide Linkage Analysis Identifies Loci for Testicle and Ovary Traits in Chickens

Development of testes or ovaries is critical to chicken breeders. Understanding the genetic mechanisms influencing the development of the testes and ovaries could enhance selection efforts which target reproductive traits. The linkage analysis was conducted within an F2 population derived from Beijing-You chickens and a commercial broiler line. The results have identified one quantitative trait loci (QTL, designated T1) for bilateral testicular weight (TW) and the percentage of TW to carcass weight, and five QTLs (designated O1–O5) for ovary weight (follicle-free, OW) and the percentage of OW to carcass weight. For the testes traits, QTL T1 is located between 6.55 and 8.56 Mb on GGA13. Especially, the gene gamma-amino butyric acid A receptor, alpha 1 (GABRA1) located near the T1 peak. For ovarian traits, QTL O2 was located at 29.31 Mb on GGA7. G protein-coupled receptor 39 (GPR39) present at the O2 peak was expressed at higher levels within the reproductive tract. It is also involved in the regulation of several reproductive functions. Other QTL peaks and the genes’ function in the ovary and testes need to be evaluated. The QTLs and the genes identified in this study could provide valuable information for establishing reproductive traits in chickens, and need further investigation.     

A Genome-Wide Association Study Identifies Potential Susceptibility Loci for Hirschsprung Disease

Hirschsprung disease (HSCR) is a congenital and heterogeneous disorder characterized by the absence of intramural nervous plexuses along variable lengths of the hindgut. Although RET is a well-established risk factor, a recent genome-wide association study (GWAS) of HSCR has identified NRG1 as an additional susceptibility locus. To discover additional risk loci, we performed a GWAS of 123 sporadic HSCR patients and 432 unaffected controls using a large-scale platform with coverage of over 1 million polymorphic markers. The result was that our study replicated the findings of RET-CSGALNACT2-RASGEF1A genomic region (_rawP = 5.69×10⁻¹⁹ before a Bonferroni correction; _corrP = 4.31×10⁻¹³ after a Bonferroni correction) and NRG1 as susceptibility loci. In addition, this study identified SLC6A20 (_adjP = 2.71×10⁻⁶), RORA (_adjP = 1.26×10⁻⁵), and ABCC9 (_adjP = 1.86×10⁻⁵) as new potential susceptibility loci under adjusting the already known loci on the RET-CSGALNACT2-RASGEF1A and NRG1 regions, although none of the SNPs in these genes passed the Bonferroni correction. In further subgroup analysis, the RET-CSGALNACT2-RASGEF1A genomic region was observed to have different significance levels among subgroups: short-segment (S-HSCR, _corrP = 1.71×10⁻⁵), long-segment (L-HSCR, _corrP = 6.66×10⁻⁴), and total colonic aganglionosis (TCA, _corrP>0.05). This differential pattern in the significance level suggests that other genomic loci or mechanisms may affect the length of aganglionosis in HSCR subgroups during enteric nervous system (ENS) development. Although functional evaluations are needed, our findings might facilitate improved understanding of the mechanisms of HSCR pathogenesis.     

Genome-Wide Scan Informed by Age-Related Disease Identifies Loci for Exceptional Human Longevity

We developed a new statistical framework to find genetic variants associated with extreme longevity. The method, informed GWAS (iGWAS), takes advantage of knowledge from large studies of age-related disease in order to narrow the search for SNPs associated with longevity. To gain support for our approach, we first show there is an overlap between loci involved in disease and loci associated with extreme longevity. These results indicate that several disease variants may be depleted in centenarians versus the general population. Next, we used iGWAS to harness information from 14 meta-analyses of disease and trait GWAS to identify longevity loci in two studies of long-lived humans. In a standard GWAS analysis, only one locus in these studies is significant (APOE/TOMM40) when controlling the false discovery rate (FDR) at 10%. With iGWAS, we identify eight genetic loci to associate significantly with exceptional human longevity at FDR < 10%. We followed up the eight lead SNPs in independent cohorts, and found replication evidence of four loci and suggestive evidence for one more with exceptional longevity. The loci that replicated (FDR < 5%) included APOE/TOMM40 (associated with Alzheimer’s disease), CDKN2B/ANRIL (implicated in the regulation of cellular senescence), ABO (tags the O blood group), and SH2B3/ATXN2 (a signaling gene that extends lifespan in Drosophila and a gene involved in neurological disease). Our results implicate new loci in longevity and reveal a genetic overlap between longevity and age-related diseases and traits, including coronary artery disease and Alzheimer’s disease. iGWAS provides a new analytical strategy for uncovering SNPs that influence extreme longevity, and can be applied more broadly to boost power in other studies of complex phenotypes.     

Genome-Wide Meta-Analysis of Myopia and Hyperopia Provides Evidence for Replication of 11 Loci

Refractive error (RE) is a complex, multifactorial disorder characterized by a mismatch between the optical power of the eye and its axial length that causes object images to be focused off the retina. The two major subtypes of RE are myopia (nearsightedness) and hyperopia (farsightedness), which represent opposite ends of the distribution of the quantitative measure of spherical refraction. We performed a fixed effects meta-analysis of genome-wide association results of myopia and hyperopia from 9 studies of European-derived populations: AREDS, KORA, FES, OGP-Talana, MESA, RSI, RSII, RSIII and ERF. One genome-wide significant region was observed for myopia, corresponding to a previously identified myopia locus on 8q12 (p = 1.25×10⁻⁸), which has been reported by Kiefer et al. as significantly associated with myopia age at onset and Verhoeven et al. as significantly associated to mean spherical-equivalent (MSE) refractive error. We observed two genome-wide significant associations with hyperopia. These regions overlapped with loci on 15q14 (minimum p value = 9.11×10⁻¹¹) and 8q12 (minimum p value 1.82×10⁻¹¹) previously reported for MSE and myopia age at onset. We also used an intermarker linkage- disequilibrium-based method for calculating the effective number of tests in targeted regional replication analyses. We analyzed myopia (which represents the closest phenotype in our data to the one used by Kiefer et al.) and showed replication of 10 additional loci associated with myopia previously reported by Kiefer et al. This is the first replication of these loci using myopia as the trait under analysis. “Replication-level” association was also seen between hyperopia and 12 of Kiefer et al.''s published loci. For the loci that show evidence of association to both myopia and hyperopia, the estimated effect of the risk alleles were in opposite directions for the two traits. This suggests that these loci are important contributors to variation of refractive error across the distribution.     

Genome-Wide Association Study in BRCA1 Mutation Carriers Identifies Novel Loci Associated with Breast and Ovarian Cancer Risk

BRCA1-associated breast and ovarian cancer risks can be modified by common genetic variants. To identify further cancer risk-modifying loci, we performed a multi-stage GWAS of 11,705 BRCA1 carriers (of whom 5,920 were diagnosed with breast and 1,839 were diagnosed with ovarian cancer), with a further replication in an additional sample of 2,646 BRCA1 carriers. We identified a novel breast cancer risk modifier locus at 1q32 for BRCA1 carriers (rs2290854, P = 2.7×10⁻⁸, HR = 1.14, 95% CI: 1.09–1.20). In addition, we identified two novel ovarian cancer risk modifier loci: 17q21.31 (rs17631303, P = 1.4×10⁻⁸, HR = 1.27, 95% CI: 1.17–1.38) and 4q32.3 (rs4691139, P = 3.4×10⁻⁸, HR = 1.20, 95% CI: 1.17–1.38). The 4q32.3 locus was not associated with ovarian cancer risk in the general population or BRCA2 carriers, suggesting a BRCA1-specific association. The 17q21.31 locus was also associated with ovarian cancer risk in 8,211 BRCA2 carriers (P = 2×10⁻⁴). These loci may lead to an improved understanding of the etiology of breast and ovarian tumors in BRCA1 carriers. Based on the joint distribution of the known BRCA1 breast cancer risk-modifying loci, we estimated that the breast cancer lifetime risks for the 5% of BRCA1 carriers at lowest risk are 28%–50% compared to 81%–100% for the 5% at highest risk. Similarly, based on the known ovarian cancer risk-modifying loci, the 5% of BRCA1 carriers at lowest risk have an estimated lifetime risk of developing ovarian cancer of 28% or lower, whereas the 5% at highest risk will have a risk of 63% or higher. Such differences in risk may have important implications for risk prediction and clinical management for BRCA1 carriers.     

A Genome-Wide Association Study of Psoriasis and Psoriatic Arthritis Identifies New Disease Loci

A genome-wide association study was performed to identify genetic factors involved in susceptibility to psoriasis (PS) and psoriatic arthritis (PSA), inflammatory diseases of the skin and joints in humans. 223 PS cases (including 91 with PSA) were genotyped with 311,398 single nucleotide polymorphisms (SNPs), and results were compared with those from 519 Northern European controls. Replications were performed with an independent cohort of 577 PS cases and 737 controls from the U.S., and 576 PSA patients and 480 controls from the U.K.. Strongest associations were with the class I region of the major histocompatibility complex (MHC). The most highly associated SNP was rs10484554, which lies 34.7 kb upstream from HLA-C (P = 7.8×10⁻¹¹, GWA scan; P = 1.8×10⁻³⁰, replication; P = 1.8×10⁻³⁹, combined; U.K. PSA: P = 6.9×10⁻¹¹). However, rs2395029 encoding the G2V polymorphism within the class I gene HCP5 (combined P = 2.13×10⁻²⁶ in U.S. cases) yielded the highest ORs with both PS and PSA (4.1 and 3.2 respectively). This variant is associated with low viral set point following HIV infection and its effect is independent of rs10484554. We replicated the previously reported association with interleukin 23 receptor and interleukin 12B (IL12B) polymorphisms in PS and PSA cohorts (IL23R: rs11209026, U.S. PS, P = 1.4×10⁻⁴; U.K. PSA: P = 8.0×10⁻⁴; IL12B:rs6887695, U.S. PS, P = 5×10⁻⁵ and U.K. PSA, P = 1.3×10⁻³) and detected an independent association in the IL23R region with a SNP 4 kb upstream from IL12RB2 (P = 0.001). Novel associations replicated in the U.S. PS cohort included the region harboring lipoma HMGIC fusion partner (LHFP) and conserved oligomeric golgi complex component 6 (COG6) genes on chromosome 13q13 (combined P = 2×10⁻⁶ for rs7993214; OR = 0.71), the late cornified envelope gene cluster (LCE) from the Epidermal Differentiation Complex (PSORS4) (combined P = 6.2×10⁻⁵ for rs6701216; OR 1.45) and a region of LD at 15q21 (combined P = 2.9×10⁻⁵ for rs3803369; OR = 1.43). This region is of interest because it harbors ubiquitin-specific protease-8 whose processed pseudogene lies upstream from HLA-C. This region of 15q21 also harbors the gene for SPPL2A (signal peptide peptidase like 2a) which activates tumor necrosis factor alpha by cleavage, triggering the expression of IL12 in human dendritic cells. We also identified a novel PSA (and potentially PS) locus on chromosome 4q27. This region harbors the interleukin 2 (IL2) and interleukin 21 (IL21) genes and was recently shown to be associated with four autoimmune diseases (Celiac disease, Type 1 diabetes, Grave''s disease and Rheumatoid Arthritis).     

A Genome-Wide Association Study Identifies Five Loci Influencing Facial Morphology in Europeans

Inter-individual variation in facial shape is one of the most noticeable phenotypes in humans, and it is clearly under genetic regulation; however, almost nothing is known about the genetic basis of normal human facial morphology. We therefore conducted a genome-wide association study for facial shape phenotypes in multiple discovery and replication cohorts, considering almost ten thousand individuals of European descent from several countries. Phenotyping of facial shape features was based on landmark data obtained from three-dimensional head magnetic resonance images (MRIs) and two-dimensional portrait images. We identified five independent genetic loci associated with different facial phenotypes, suggesting the involvement of five candidate genes—PRDM16, PAX3, TP63, C5orf50, and COL17A1—in the determination of the human face. Three of them have been implicated previously in vertebrate craniofacial development and disease, and the remaining two genes potentially represent novel players in the molecular networks governing facial development. Our finding at PAX3 influencing the position of the nasion replicates a recent GWAS of facial features. In addition to the reported GWA findings, we established links between common DNA variants previously associated with NSCL/P at 2p21, 8q24, 13q31, and 17q22 and normal facial-shape variations based on a candidate gene approach. Overall our study implies that DNA variants in genes essential for craniofacial development contribute with relatively small effect size to the spectrum of normal variation in human facial morphology. This observation has important consequences for future studies aiming to identify more genes involved in the human facial morphology, as well as for potential applications of DNA prediction of facial shape such as in future forensic applications.     

Genome-Wide Association Study Link Novel Loci to Endometriosis

Endometriosis is a common gynecological condition with complex etiology defined by the presence of endometrial glands and stroma outside the womb. Endometriosis is a common cause of both cyclic and chronic pelvic pain, reduced fertility, and reduced quality-of-life. Diagnosis and treatment of endometriosis is, on average, delayed by 7–10 years from the onset of symptoms. Absence of a timely and non-invasive diagnostic tool is presently the greatest barrier to the identification and treatment of endometriosis. Twin and family studies have documented an increased relative risk in families. To identify genetic factors that contribute to endometriosis we conducted a two-stage genome-wide association study (GWAS) of a European cohort including 2,019 surgically confirmed endometriosis cases and 14,471 controls. Three of the SNPs we identify associated at P<5×10⁻⁸ in our combined analysis belong to two loci: LINC00339-WNT4 on 1p36.12 (rs2235529; P = 8.65×10⁻⁹, OR = 1.29, CI = 1.18–1.40) and RND3-RBM43 on 2q23.3 (rs1519761; P = 4.70×10⁻⁸, OR = 1.20, Cl = 1.13–1.29, and rs6757804; P = 4.05×10⁻⁸, OR = 1.20, Cl = 1.13–1.29). Using an adjusted Bonferoni significance threshold of 4.51×10⁻⁷ we identify two additional loci in our meta-analysis that associate with endometriosis:, RNF144B-ID4 on 6p22.3 (rs6907340; P = 2.19×10⁻⁷, OR = 1.20, Cl = 1.12–1.28), and HNRNPA3P1-LOC100130539 on 10q11.21 (rs10508881; P = 4.08×10⁻⁷, OR = 1.19, Cl = 1.11–1.27). Consistent with previously suggested associations to WNT4 our study implicate a 150 kb region around WNT4 that also include LINC00339 and CDC42. A univariate analysis of documented infertility, age at menarche, and family history did not show allelic association with these SNP markers. Clinical data from patients in our study reveal an average delay in diagnosis of 8.4 years and confirm a strong correlation between endometriosis severity and infertility (n = 1182, P<0.001, OR = 2.18). This GWAS of endometriosis was conducted with high diagnostic certainty in cases, and with stringent handling of population substructure. Our findings broaden the understanding of the genetic factors that play a role in endometriosis.     

Genome-Wide RNAi Screen Identifies Novel Host Proteins Required for Alphavirus Entry

The enveloped alphaviruses include important and emerging human pathogens such as Chikungunya virus and Eastern equine encephalitis virus. Alphaviruses enter cells by clathrin-mediated endocytosis, and exit by budding from the plasma membrane. While there has been considerable progress in defining the structure and function of the viral proteins, relatively little is known about the host factors involved in alphavirus infection. We used a genome-wide siRNA screen to identify host factors that promote or inhibit alphavirus infection in human cells. Fuzzy homologue (FUZ), a protein with reported roles in planar cell polarity and cilia biogenesis, was required for the clathrin-dependent internalization of both alphaviruses and the classical endocytic ligand transferrin. The tetraspanin membrane protein TSPAN9 was critical for the efficient fusion of low pH-triggered virus with the endosome membrane. FUZ and TSPAN9 were broadly required for infection by the alphaviruses Sindbis virus, Semliki Forest virus, and Chikungunya virus, but were not required by the structurally-related flavivirus Dengue virus. Our results highlight the unanticipated functions of FUZ and TSPAN9 in distinct steps of alphavirus entry and suggest novel host proteins that may serve as targets for antiviral therapy.     

Genome-Wide Association Study Identifies a Novel Canine Glaucoma Locus

Glaucoma is an optic neuropathy and one of the leading causes of blindness. Its hereditary forms are classified into primary closed-angle (PCAG), primary open-angle (POAG) and primary congenital glaucoma (PCG). Although many loci have been mapped in human, only a few genes have been identified that are associated with the development of glaucoma and the genetic basis of the disease remains poorly understood. Glaucoma has also been described in many dog breeds, including Dandie Dinmont Terriers (DDT) in which it is a late-onset (>7 years) disease. We designed clinical and genetic studies to better define the clinical features of glaucoma in the DDT and to identify the genetic cause. Clinical diagnosis was based on ophthalmic examinations of the affected dogs and 18 additionally investigated unaffected DDTs. We collected DNA from over 400 DTTs and a genome wide association study was performed in a cohort of 23 affected and 23 controls, followed by a fine mapping, a replication study and candidate gene sequencing. The clinical study suggested that ocular abnormalities including abnormal iridocorneal angles and pectinate ligament dysplasia are common (50% and 72%, respectively) in the breed and the disease resembles human PCAG. The genetic study identified a novel 9.5 Mb locus on canine chromosome 8 including the 1.6 Mb best associated region (p = 1.63×10⁻¹⁰, OR = 32 for homozygosity). Mutation screening in five candidate genes did not reveal any causative variants. This study indicates that although ocular abnormalities are common in DDTs, the genetic risk for glaucoma is conferred by a novel locus on CFA8. The canine locus shares synteny to a region in human chromosome 14q, which harbors several loci associated with POAG and PCG. Our study reveals a new locus for canine glaucoma and ongoing molecular studies will likely help to understand the genetic etiology of the disease.     

Integration of Genome-Wide SNP Data and Gene-Expression Profiles Reveals Six Novel Loci and Regulatory Mechanisms for Amino Acids and Acylcarnitines in Whole Blood

Profiling amino acids and acylcarnitines in whole blood spots is a powerful tool in the laboratory diagnosis of several inborn errors of metabolism. Emerging data suggests that altered blood levels of amino acids and acylcarnitines are also associated with common metabolic diseases in adults. Thus, the identification of common genetic determinants for blood metabolites might shed light on pathways contributing to human physiology and common diseases. We applied a targeted mass-spectrometry-based method to analyze whole blood concentrations of 96 amino acids, acylcarnitines and pathway associated metabolite ratios in a Central European cohort of 2,107 adults and performed genome-wide association (GWA) to identify genetic modifiers of metabolite concentrations. We discovered and replicated six novel loci associated with blood levels of total acylcarnitine, arginine (both on chromosome 6; rs12210538, rs17657775), propionylcarnitine (chromosome 10; rs12779637), 2-hydroxyisovalerylcarnitine (chromosome 21; rs1571700), stearoylcarnitine (chromosome 1; rs3811444), and aspartic acid traits (chromosome 8; rs750472). Based on an integrative analysis of expression quantitative trait loci in blood mononuclear cells and correlations between gene expressions and metabolite levels, we provide evidence for putative causative genes: SLC22A16 for total acylcarnitines, ARG1 for arginine, HLCS for 2-hydroxyisovalerylcarnitine, JAM3 for stearoylcarnitine via a trans-effect at chromosome 1, and PPP1R16A for aspartic acid traits. Further, we report replication and provide additional functional evidence for ten loci that have previously been published for metabolites measured in plasma, serum or urine.In conclusion, our integrative analysis of SNP, gene-expression and metabolite data points to novel genetic factors that may be involved in the regulation of human metabolism. At several loci, we provide evidence for metabolite regulation via gene-expression and observed overlaps with GWAS loci for common diseases. These results form a strong rationale for subsequent functional and disease-related studies.