Vincenzo Forgetta

An atlas of genetic influences on human blood metabolites.

Genome-wide association scans with high-throughput metabolic profiling provide unprecedented insights into how genetic variation influences metabolism and complex disease. Here we report the most comprehensive exploration of genetic loci influencing human metabolism thus far, comprising 7,824 adult individuals from 2 European population studies. We report genome-wide significant associations at 145 metabolic loci and their biochemical connectivity with more than 400 metabolites in human blood. We extensively characterize the resulting in vivo blueprint of metabolism in human blood by integrating it with information on gene expression, heritability and overlap with known loci for complex disorders, inborn errors of metabolism and pharmacological targets. We further developed a database and web-based resources for data mining and results visualization. Our findings provide new insights into the role of inherited variation in blood metabolic diversity and identify potential new opportunities for drug development and for understanding disease.

Allelic Variation in TLR4 Is Linked to Susceptibility to Salmonella enterica Serovar Typhimurium Infection in Chickens

ABSTRACT Toll-like receptor 4 (TLR4) is part of a group of evolutionarily conserved pattern recognition receptors involved in the activation of the immune system in response to various pathogens and in the innate defense against infection. We describe here the cloning and characterization of the avian orthologue of mammalian TLR4. Chicken TLR4 encodes a 843-amino-acid protein that contains a leucine-rich repeat extracellular domain, a short transmembrane domain typical of type I transmembrane proteins, and a Toll-interleukin-1R signaling domain characteristic of all TLR proteins. The chicken TLR4 protein shows 46% identity (64% similarity) to human TLR4 and 41% similarity to other TLR family members. Northern blot analysis reveals that TLR4 is expressed at approximately the same level in all tissues tested, including brain, thymus, kidney, intestine, muscle, liver, lung, bursa of Fabricius, heart, and spleen. The probe detected only one transcript of ca. 4.4 kb in length for all tissues except muscle where the size of TLR4 mRNA was ca. 9.6 kb. We have mapped TLR4 to microchromosome E41W17 in a region harboring the gene for tenascin C and known to be well conserved between the chicken and mammalian genomes. This region of the chicken genome was shown previously to harbor a Salmonella susceptibility locus. By using linkage analysis, TLR4 was shown to be linked to resistance to infection with Salmonella enterica serovar Typhimurium in chickens (likelihood ratio test of 10.2, P = 0.00138), suggesting a role of TLR4 in the host response of chickens to Salmonella infection.

Vitamin D and Risk of Multiple Sclerosis: A Mendelian Randomization Study

Background Observational studies have demonstrated an association between decreased vitamin D level and risk of multiple sclerosis (MS); however, it remains unclear whether this relationship is causal. We undertook a Mendelian randomization (MR) study to evaluate whether genetically lowered vitamin D level influences the risk of MS. Methods and Findings We identified single nucleotide polymorphisms (SNPs) associated with 25-hydroxyvitamin D (25OHD) level from SUNLIGHT, the largest (n = 33,996) genome-wide association study to date for vitamin D. Four SNPs were genome-wide significant for 25OHD level (p-values ranging from 6 × 10−10 to 2 × 10−109), and all four SNPs lay in, or near, genes strongly implicated in separate mechanisms influencing 25OHD. We then ascertained their effect on 25OHD level in 2,347 participants from a population-based cohort, the Canadian Multicentre Osteoporosis Study, and tested the extent to which the 25OHD-decreasing alleles explained variation in 25OHD level. We found that the count of 25OHD-decreasing alleles across these four SNPs was strongly associated with lower 25OHD level (n = 2,347, F-test statistic = 49.7, p = 2.4 × 10−12). Next, we conducted an MR study to describe the effect of genetically lowered 25OHD on the odds of MS in the International Multiple Sclerosis Genetics Consortium study, the largest genetic association study to date for MS (including up to 14,498 cases and 24,091 healthy controls). Alleles were weighted by their relative effect on 25OHD level, and sensitivity analyses were performed to test MR assumptions. MR analyses found that each genetically determined one-standard-deviation decrease in log-transformed 25OHD level conferred a 2.0-fold increase in the odds of MS (95% CI: 1.7–2.5; p = 7.7 × 10−12; I 2 = 63%, 95% CI: 0%–88%). This result persisted in sensitivity analyses excluding SNPs possibly influenced by population stratification or pleiotropy (odds ratio [OR] = 1.7, 95% CI: 1.3–2.2; p = 2.3 × 10−5; I 2 = 47%, 95% CI: 0%–85%) and including only SNPs involved in 25OHD synthesis or metabolism (ORsynthesis = 2.1, 95% CI: 1.6–2.6, p = 1 × 10−9; ORmetabolism = 1.9, 95% CI: 1.3–2.7, p = 0.002). While these sensitivity analyses decreased the possibility that pleiotropy may have biased the results, residual pleiotropy is difficult to exclude entirely. Conclusions A genetically lowered 25OHD level is strongly associated with increased susceptibility to MS. Whether vitamin D sufficiency can delay, or prevent, MS onset merits further investigation in long-term randomized controlled trials.

A Mendelian randomization study of the effect of type-2 diabetes on coronary heart disease

In observational studies, type-2 diabetes (T2D) is associated with an increased risk of coronary heart disease (CHD), yet interventional trials have shown no clear effect of glucose-lowering on CHD. Confounding may have therefore influenced these observational estimates. Here we use Mendelian randomization to obtain unconfounded estimates of the influence of T2D and fasting glucose (FG) on CHD risk. Using multiple genetic variants associated with T2D and FG, we find that risk of T2D increases CHD risk (odds ratio (OR)=1.11 (1.05–1.17), per unit increase in odds of T2D, P=8.8 × 10−5; using data from 34,840/114,981 T2D cases/controls and 63,746/130,681 CHD cases/controls). FG in non-diabetic individuals tends to increase CHD risk (OR=1.15 (1.00–1.32), per mmol·per l, P=0.05; 133,010 non-diabetic individuals and 63,746/130,681 CHD cases/controls). These findings provide evidence supporting a causal relationship between T2D and CHD and suggest that long-term trials may be required to discern the effects of T2D therapies on CHD risk.

Whole-genome sequence-based analysis of thyroid function

Normal thyroid function is essential for health, but its genetic architecture remains poorly understood. Here, for the heritable thyroid traits thyrotropin (TSH) and free thyroxine (FT4), we analyse whole-genome sequence data from the UK10K project (N=2,287). Using additional whole-genome sequence and deeply imputed data sets, we report meta-analysis results for common variants (MAF≥1%) associated with TSH and FT4 (N=16,335). For TSH, we identify a novel variant in SYN2 (MAF=23.5%, P=6.15 × 10−9) and a new independent variant in PDE8B (MAF=10.4%, P=5.94 × 10−14). For FT4, we report a low-frequency variant near B4GALT6/SLC25A52 (MAF=3.2%, P=1.27 × 10−9) tagging a rare TTR variant (MAF=0.4%, P=2.14 × 10−11). All common variants explain ≥20% of the variance in TSH and FT4. Analysis of rare variants (MAF<1%) using sequence kernel association testing reveals a novel association with FT4 in NRG1. Our results demonstrate that increased coverage in whole-genome sequence association studies identifies novel variants associated with thyroid function.

Fourteen-Genome Comparison Identifies DNA Markers for Severe-Disease-Associated Strains of Clostridium difficile

ABSTRACT Clostridium difficile is a common cause of infectious diarrhea in hospitalized patients. A severe and increased incidence of C. difficile infection (CDI) is associated predominantly with the NAP1 strain; however, the existence of other severe-disease-associated (SDA) strains and the extensive genetic diversity across C. difficile complicate reliable detection and diagnosis. Comparative genome analysis of 14 sequenced genomes, including those of a subset of NAP1 isolates, allowed the assessment of genetic diversity within and between strain types to identify DNA markers that are associated with severe disease. Comparative genome analysis of 14 isolates, including five publicly available strains, revealed that C. difficile has a core genome of 3.4 Mb, comprising ∼3,000 genes. Analysis of the core genome identified candidate DNA markers that were subsequently evaluated using a multistrain panel of 177 isolates, representing more than 50 pulsovars and 8 toxinotypes. A subset of 117 isolates from the panel had associated patient data that allowed assessment of an association between the DNA markers and severe CDI. We identified 20 candidate DNA markers for species-wide detection and 10,683 single nucleotide polymorphisms (SNPs) associated with the predominant SDA strain (NAP1). A species-wide detection candidate marker, the sspA gene, was found to be the same across 177 sequenced isolates and lacked significant similarity to those of other species. Candidate SNPs in genes CD1269 and CD1265 were found to associate more closely with disease severity than currently used diagnostic markers, as they were also present in the toxin A-negative and B-positive (A-B+) strain types. The genetic markers identified illustrate the potential of comparative genomics for the discovery of diagnostic DNA-based targets that are species specific or associated with multiple SDA strains.

Identification of 153 new loci associated with heel bone mineral density and functional involvement of GPC6 in osteoporosis

Osteoporosis is a common disease diagnosed primarily by measurement of bone mineral density (BMD). We undertook a genome-wide association study (GWAS) in 142,487 individuals from the UK Biobank to identify loci associated with BMD as estimated by quantitative ultrasound of the heel. We identified 307 conditionally independent single-nucleotide polymorphisms (SNPs) that attained genome-wide significance at 203 loci, explaining approximately 12% of the phenotypic variance. These included 153 previously unreported loci, and several rare variants with large effect sizes. To investigate the underlying mechanisms, we undertook (1) bioinformatic, functional genomic annotation and human osteoblast expression studies; (2) gene-function prediction; (3) skeletal phenotyping of 120 knockout mice with deletions of genes adjacent to lead independent SNPs; and (4) analysis of gene expression in mouse osteoblasts, osteocytes and osteoclasts. The results implicate GPC6 as a novel determinant of BMD, and also identify abnormal skeletal phenotypes in knockout mice associated with a further 100 prioritized genes.

Genetically decreased vitamin D and risk of Alzheimer disease

Objective: To test whether genetically decreased vitamin D levels are associated with Alzheimer disease (AD) using mendelian randomization (MR), a method that minimizes bias due to confounding or reverse causation. Methods: We selected single nucleotide polymorphisms (SNPs) that are strongly associated with 25-hydroxyvitamin D (25OHD) levels (p < 5 × 10−8) from the Study of Underlying Genetic Determinants of Vitamin D and Highly Related Traits (SUNLIGHT) Consortium (N = 33,996) to act as instrumental variables for the MR study. We measured the effect of each of these SNPs on 25OHD levels in the Canadian Multicentre Osteoporosis Study (CaMos; N = 2,347) and obtained the corresponding effect estimates for each SNP on AD risk from the International Genomics of Alzheimers Project (N = 17,008 AD cases and 37,154 controls). To produce MR estimates, we weighted the effect of each SNP on AD by its effect on 25OHD and meta-analyzed these estimates using a fixed-effects model to provide a summary effect estimate. Results: The SUNLIGHT Consortium identified 4 SNPs to be genome-wide significant for 25OHD, which described 2.44% of the variance in 25OHD in CaMos. All 4 SNPs map to genes within the vitamin D metabolic pathway. MR analyses demonstrated that a 1-SD decrease in natural log–transformed 25OHD increased AD risk by 25% (odds ratio 1.25, 95% confidence interval 1.03–1.51, p = 0.021). After sensitivity analysis in which we removed SNPs possibly influenced by pleiotropy and population stratification, the results were largely unchanged. Conclusions: Our results provide evidence supporting 25OHD as a causal risk factor for AD. These findings provide further rationale to understand the effect of vitamin D supplementation on cognition and AD risk in randomized controlled trials.

Mendelian randomisation applied to drug development in cardiovascular disease: a review

Despite increased expenditure, productivity of the pharmaceutical industry has decreased and currently 90% of developed molecules entering phase II and phase III clinical trials fail to gain regulatory approval. Most of these failures are due to lack of therapeutic efficacy rather than lack of safety, suggesting that drug development failures may often be due to poor drug target validation. Currently, drug targets are largely validated using in vitro assays and animal models which may not translate well to human disease. Emerging methods from human genetics, such as Mendelian randomisation (MR), can enable the validation of putative biomarker drug targets in humans prior to the initiation of clinical trials. MR studies can provide evidence as to whether genetically determined levels of a biomarker influence disease aetiology, enabling investigators to infer whether the biomarker is causal. We review the extent to which MR techniques may be helpful in biomarker validation by assessing the concordance between the results from MR studies and phase III clinical trials for lipid therapy in cardiovascular disease. Our findings show that concordance is highest when MR provides evidence suggesting that a biomarker is not causal. In contrast, there are many examples of clinical trials that still failed despite targeting confirmed causal biomarkers. We discuss why such trials may not succeed, despite evidence for causality in MR studies, and outline important limitations when using MR for biomarker validation in drug development. Nonetheless, given the current inefficiencies in drug development, MR methods have potential to improve the success rate of drug development and ultimately the delivery of new molecules to clinical care.

Sequencing of the Dutch elm disease fungus genome using the Roche/454 GS-FLX Titanium System in a comparison of multiple genomics core facilities.

As part of the DNA Sequencing Research Group of the Association of Biomolecular Resource Facilities, we have tested the reproducibility of the Roche/454 GS-FLX Titanium System at five core facilities. Experience with the Roche/454 system ranged from <10 to >340 sequencing runs performed. All participating sites were supplied with an aliquot of a common DNA preparation and were requested to conduct sequencing at a common loading condition. The evaluation of sequencing yield and accuracy metrics was assessed at a single site. The study was conducted using a laboratory strain of the Dutch elm disease fungus Ophiostoma novo-ulmi strain H327, an ascomycete, vegetatively haploid fungus with an estimated genome size of 30-50 Mb. We show that the Titanium System is reproducible, with some variation detected in loading conditions, sequencing yield, and homopolymer length accuracy. We demonstrate that reads shorter than the theoretical minimum length are of lower overall quality and not simply truncated reads. The O. novo-ulmi H327 genome assembly is 31.8 Mb and is comprised of eight chromosome-length linear scaffolds, a circular mitochondrial conti of 66.4 kb, and a putative 4.2-kb linear plasmid. We estimate that the nuclear genome encodes 8613 protein coding genes, and the mitochondrion encodes 15 genes and 26 tRNAs.