Search

The authors have withdrawn this manuscript because this paper was posted prematurely in advance of a UK Biobank Pharma Proteomics Project consortium effort. Therefore, the authors do not wish this work to be cited as reference for the project. If you have any questions, please contact the corresponding author

We describe the analysis of whole genome sequences (WGS) of 150,119 individuals from the UK biobank (UKB). This constitutes a set of high quality variants, including 585,040,410 SNPs, representing 7.0% of all possible human SNPs, and 58,707,036 indels. The large set of variants allows us to characterize selection based on sequence variation within a population through a Depletion Rank (DR) score for windows along the genome. DR analysis shows that coding exons represent a small fraction of regions in the genome subject to strong sequence conservation. We define three cohorts within the UKB, a large British Irish cohort (XBI) and smaller African (XAF) and South Asian (XSA) cohorts. A haplotype reference panel is provided that allows reliable imputation of most variants carried by three or more sequenced individuals. We identified 895,055 structural variants and 2,536,688 microsatellites, groups of variants typically excluded from large scale WGS studies. Using this formidable new resource, we provide several examples of trait associations for rare variants with large effects not found previously through studies based on exome sequencing and/or imputation.

The spread of SARS-CoV-2 is dependent on several factors, both biological and behavioral. The effectiveness of various non-pharmaceutical interventions can largely be attributed to changes in human behavior, but quantifying this effect remains challenging. Reconstructing the transmission tree of the third wave of SARS-CoV-2 infections in Iceland using contact tracing and viral sequence data from 2522 cases enables us to compare the infectiousness of distinct groups of persons directly. We find that people diagnosed outside of quarantine are 89% more infectious than those diagnosed while in quarantine, and infectiousness decreases as a function of the time spent in quarantine. Furthermore, we find that people of working age, 16-66 years old, are 47% more infectious than those outside that age range. Lastly, the transmission tree enables us to model the effect that given population prevalence of vaccination would have had on the third wave had they been administered before that time using several different strategies. We find that vaccinating in order of ascending age or uniformly at random would have prevented more infections per vaccination than vaccinating in order of descending age.

A pressing concern in the SARS-CoV-2 epidemic and other viral outbreaks, is the extent to which the containment measures are halting the viral spread. A straightforward way to assess this is to tally the active cases and the recovered ones throughout the epidemic. Here, we show how epidemic control can be assessed with molecular information during a well characterized epidemic in Iceland. We demonstrate how the viral concentration decreased in those newly diagnosed as the epidemic transitioned from exponential growth phase to containment phase. The viral concentration in the cases identified in population screening decreased faster than in those symptomatic and considered at high risk and that were targeted by the healthcare system. The viral concentration persists in recovering individuals as we found that half of the cases are still positive after two weeks. We demonstrate that accumulation of mutations in SARS-CoV-2 genome can be exploited to track the rate of new viral generations throughout the different phases of the epidemic, where the accumulation of mutations decreases as the transmission rate decreases in the containment phase. Overall, the molecular signatures of SARS-CoV-2 infections contain valuable epidemiological information that can be used to assess the effectiveness of containment measures., The concentration of SARS-CoV-2 changes during an individual’s infection, and mutations accumulate as viruses are transmitted between people. Here, the authors use data from Iceland to demonstrate how this information can be exploited at the population-level to determine the phase of the epidemic.

The plasma proteome can help bridge the gap between the genome and diseases. Here we describe genome-wide association studies (GWASs) of plasma protein levels measured with 4,907 aptamers in 35,559 Icelanders. We found 18,084 associations between sequence variants and levels of proteins in plasma (protein quantitative trait loci; pQTL), of which 19% were with rare variants (minor allele frequency (MAF)

Publisher's version (útgefin grein), Aims: To explore whether variability in dietary cholesterol and phytosterol absorption impacts the risk of coronary artery disease (CAD) using as instruments sequence variants in the ABCG5/8 genes, key regulators of intestinal absorption of dietary sterols. Methods and results: We examined the effects of ABCG5/8 variants on non-high-density lipoprotein (non-HDL) cholesterol (N up to 610 532) and phytosterol levels (N = 3039) and the risk of CAD in Iceland, Denmark, and the UK Biobank (105 490 cases and 844 025 controls). We used genetic scores for non-HDL cholesterol to determine whether ABCG5/8 variants confer greater risk of CAD than predicted by their effect on non-HDL cholesterol. We identified nine rare ABCG5/8 coding variants with substantial impact on non-HDL cholesterol. Carriers have elevated phytosterol levels and are at increased risk of CAD. Consistent with impact on ABCG5/8 transporter function in hepatocytes, eight rare ABCG5/8 variants associate with gallstones. A genetic score of ABCG5/8 variants predicting 1 mmol/L increase in non-HDL cholesterol associates with two-fold increase in CAD risk [odds ratio (OR) = 2.01, 95% confidence interval (CI) 1.75-2.31, P = 9.8 × 10-23] compared with a 54% increase in CAD risk (OR = 1.54, 95% CI 1.49-1.59, P = 1.1 × 10-154) associated with a score of other non-HDL cholesterol variants predicting the same increase in non-HDL cholesterol (P for difference in effects = 2.4 × 10-4). Conclusions: Genetic variation in cholesterol absorption affects levels of circulating non-HDL cholesterol and risk of CAD. Our results indicate that both dietary cholesterol and phytosterols contribute directly to atherogenesis., This work was supported by deCODE genetics/Amgen. S.B. is supported by the Novo Nordisk Foundation (NNF14CC0001 and NNF17OC0027594). F.W.A. is supported by University College London Hospital National Institute for Health Research Biomedical Research Centre (BRC288A).

BACKGROUNDLimited data exist on how SARS-CoV-2 enters and spreads in the general population.METHODSWe used two strategies for SARS-CoV-2 testing: targeted testing of high-risk individuals (n=4,551) and a population screening (n=5,502). We sequenced SARS-CoV-2 from 340 individuals.RESULTSOn March 22 2020, 528 had tested positive for SARS-CoV-2 in the targeted testing (11.6%) and 50 in the population screening (0.9%); approximately 0.2% of the Icelandic population. Large fractions of positives had travelled outside Iceland (38.4% and 34.0%). Fewer under 10 years old were positive than those older: 2.8% vs. 12.3% for targeted testing (P=1.6e-9) and 0.0% vs. 1.0% for population screening (P=0.031). Fewer females were positive in the targeted testing than males (9.5% vs. 14.6%, P=6.8e-9). SARS-CoV-2 came from eight clades, seven A clades and one B clade. The clade composition differed between the testing groups and changed with time. In the early targeted testing, 65.0% of clades were A2a1 and A2a2 derived from Italian and Austrian skiing areas, but in the later targeted testing went down to 30.6% and were overtaken by A1a and A2a, the most common clades in the population screening.CONCLUSIONSARS-CoV-2 has spread widely in Iceland outside of the high-risk groups. Several strains cause these infections and their relative contribution changed rapidly. Children and females are less vulnerable than adults and males. To contain the pandemic we must increase the scope of the testing.

De novo mutations (DNMs) cause a large proportion of severe rare diseases of childhood. DNMs that occur early may result in mosaicism of both somatic and germ cells. Such early mutations can cause recurrence of disease. We scanned 1,007 sibling pairs from 251 families and identified 878 DNMs shared by siblings (ssDNMs) at 448 genomic sites. We estimated DNM recurrence probability based on parental mosaicism, sharing of DNMs among siblings, parent-of-origin, mutation type and genomic position. We detected 57.2% of ssDNMs in the parental blood. The recurrence probability of a DNM decreases by 2.27% per year for paternal DNMs and 1.78% per year for maternal DNMs. Maternal ssDNMs are more likely to be T>C mutations than paternal ssDNMs, and less likely to be C>T mutations. Depending on the properties of the DNM, the recurrence probability ranges from 0.011% to 28.5%. We have launched an online calculator to allow estimation of DNM recurrence probability for research purposes.

Publisher's version (útgefin grein) Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations., Mutations in genes encoding subunits of the phagocyte NADPH oxidase complex are recognized to cause chronic granulomatous disease (CGD), a severe primary immunodeficiency. Here we describe how deficiency of CYBC1, a previously uncharacterized protein in humans (C17orf62), leads to reduced expression of NADPH oxidase’s main subunit (gp91phox) and results in CGD. Analyzing two brothers diagnosed with CGD we identify a homozygous loss-of-function mutation, p.Tyr2Ter, in CYBC1. Imputation of p.Tyr2Ter into 155K chipgenotyped Icelanders reveals six additional homozygotes, all with signs of CGD, manifesting as colitis, rare infections, or a severely impaired PMA-induced neutrophil oxidative burst. Homozygosity for p.Tyr2Ter consequently associates with inflammatory bowel disease (IBD) in Iceland (P = 8.3 × 10−8; OR = 67.6), as well as reduced height (P = 3.3 × 10−4; −8.5 cm). Overall, we find that CYBC1 deficiency results in CGD characterized by colitis and a distinct profile of infections indicative of macrophage dysfunction., We wish to thank the family of the two probands, as well as all the other individuals who participated in the study and whose contribution made this work possible.

The characterization of mutational processes that generate sequence diversity in the human genome is of paramount importance both to medical genetics and to evolutionary studies. To understand how the age and sex of transmitting parents affect de novo mutations, here we sequence 1,548 Icelanders, their parents, and, for a subset of 225, at least one child, to 35× genome-wide coverage. We find 108,778 de novo mutations, both single nucleotide polymorphisms and indels, and determine the parent of origin of 42,961. The number of de novo mutations from mothers increases by 0.37 per year of age (95% CI 0.32-0.43), a quarter of the 1.51 per year from fathers (95% CI 1.45-1.57). The number of clustered mutations increases faster with the mother's age than with the father's, and the genomic span of maternal de novo mutation clusters is greater than that of paternal ones. The types of de novo mutation from mothers change substantially with age, with a 0.26% (95% CI 0.19-0.33%) decrease in cytosine-phosphate-guanine to thymine-phosphate-guanine (CpG>TpG) de novo mutations and a 0.33% (95% CI 0.28-0.38%) increase in C>G de novo mutations per year, respectively. Remarkably, these age-related changes are not distributed uniformly across the genome. A striking example is a 20 megabase region on chromosome 8p, with a maternal C>G mutation rate that is up to 50-fold greater than the rest of the genome. The age-related accumulation of maternal non-crossover gene conversions also mostly occurs within these regions. Increased sequence diversity and linkage disequilibrium of C>G variants within regions affected by excess maternal mutations indicate that the underlying mutational process has persisted in humans for thousands of years. Moreover, the regional excess of C>G variation in humans is largely shared by chimpanzees, less by gorillas, and is almost absent from orangutans. This demonstrates that sequence diversity in humans results from evolving interactions between age, sex, mutation type, and genomic location.

Contains fulltext : 174179.pdf (Publisher’s version ) (Closed access) We performed a genome-wide association study of total hip replacements, based on variants identified through whole-genome sequencing, which included 4,657 Icelandic patients and 207,514 population controls. We discovered two rare signals that strongly associate with osteoarthritis total hip replacement: a missense variant, c.1141G>C (p.Asp369His), in the COMP gene (allelic frequency = 0.026%, P = 4.0 x 10-12, odds ratio (OR) = 16.7) and a frameshift mutation, rs532464664 (p.Val330Glyfs*106), in the CHADL gene that associates through a recessive mode of inheritance (homozygote frequency = 0.15%, P = 4.5 x 10-18, OR = 7.71). On average, c.1141G>C heterozygotes and individuals homozygous for rs532464664 had their hip replacement operation 13.5 years and 4.9 years earlier than others (P = 0.0020 and P = 0.0026), respectively. We show that the full-length CHADL transcript is expressed in cartilage. Furthermore, the premature stop codon introduced by the CHADL frameshift mutation results in nonsense-mediated decay of the mutant transcripts.

Common sequence variants at the haptoglobin gene (HP) have been associated with blood lipid levels. Through whole-genome sequencing of 8,453 Icelanders, we discovered a splice donor founder mutation in HP (NM_001126102.1:c.190 + 1G > C, minor allele frequency = 0.56%). This mutation occurs on the HP1 allele of the common copy number variant in HP and leads to a loss of function of HP1. It associates with lower levels of haptoglobin (P = 2.1 × 10-54), higher levels of non-high density lipoprotein cholesterol (β = 0.26 mmol/l, P = 2.6 × 10-9) and greater risk of coronary artery disease (odds ratio = 1.30, 95% confidence interval: 1.10-1.54, P = 0.0024). Through haplotype analysis and with RNA sequencing, we provide evidence of a causal relationship between one of the two haptoglobin isoforms, namely Hp1, and lower levels of non-HDL cholesterol. Furthermore, we show that the HP1 allele associates with various other quantitative biological traits.

Despite the important role that monozygotic twins have played in genetics research, little is known about their genomic differences. Here we show that monozygotic twins differ on average by 5.2 early developmental mutations and that approximately 15% of monozygotic twins have a substantial number of these early developmental mutations specific to one of them. Using the parents and offspring of twins, we identified pre-twinning mutations. We observed instances where a twin was formed from a single cell lineage in the pre-twinning cell mass and instances where a twin was formed from several cell lineages. CpGTpG mutations increased in frequency with embryonic development, coinciding with an increase in DNA methylation. Our results indicate that allocations of cells during development shapes genomic differences between monozygotic twins.

Background: The ATP-binding cassette (ABC) G5 and G8 transporters, encoded by the ABCG5 and ABCG8 (ABCG5/8) genes, have a key role in controlling intestinal absorption of cholesterol and other sterols, and in promoting their hepatobiliary secretion. Here, we use novel rare ABCG5/8 variants, together with common variants, to explore whether variability in sterol absorption impacts the risk of coronary artery disease (CAD). Methods: We examined 41 ABCG5/8 coding variants for association with non-high density lipoprotein (non-HDL) cholesterol in data from Iceland and the Global Lipid Genetics Consortium (N up to 472,392 combined). We tested variants for association with CAD in up to 67,130 cases and 700,075 controls from Iceland and the UK Biobank. Findings: We identified several rare ABCG5/8 coding variants with substantial impact on circulating levels of non-HDL cholesterol, and demonstrate that heterozygous carriers are at increased risk of CAD and other related diseases. Consistent with causal effects of other sterols such as phytosterols, on atherosclerotic disease we demonstrate that the degree of CAD risk conferred by both common and rare ABCG5/8 variants is greater than predicted based on their effect on non-HDL cholesterol levels. Interpretation: Our data highlight intestinal absorption as an important mechanism in the regulation of non-HDL /LDL cholesterol levels and CAD risk. Further, the results support the notion that other sterols such as phytosterols contribute to the pathogenesis of CAD. Funding Statement: The study was funded by deCODE genetics/Amgen Inc. Declaration of Interests: The authors who are affiliated with deCODE genetics/AMGEN, declare competing financial interests as employees. The remaining authors declare no competing financial interests. Ethics Approval Statement: The study was approved by The Icelandic Data Protection Authority and the National Bioethics Committee of Iceland (approvals no. VSNb201508003/03.01, VSNb2015010033/03.12 and VSNb2014100020/03.12 with amendments). All participating subjects donating biological samples signed informed consents. Personal identities of the participants and biological samples were encrypted with a third-party system approved and monitored by the Icelandic Data Protection Authority.

The ages of puberty, first sexual intercourse and first birth signify the onset of reproductive ability, behaviour and success, respectively. These sequenced events have behavioural, physiological and health significance, and may also influence overall reproductive fitness. In a genome-wide association study of 125,667 white men and women aged 40-69 in the UK Biobank Study, we identify 38 sequence variants with association P-values

Spinal muscular atrophies (SMAs) are a heterogeneous group of disorders characterized by muscular atrophy, weakness, and hypotonia due to suspected lower motor neuron degeneration (LMND). In a large cohort of 3,465 individuals suspected with SMA submitted for SMN1 testing to our routine diagnostic laboratory, 48.8% carried a homozygous SMN1 deletion, 2.8% a subtle mutation, and an SMN1 deletion, whereas 48.4% remained undiagnosed. Recently, several other genes implicated in SMA/LMND have been reported. Despite several efforts to establish a diagnostic algorithm for non-5q-SMA (SMA without deletion or point mutations in SMN1 [5q13.2]), data from large-scale studies are not available. We tested the clinical utility of targeted sequencing in non-5q-SMA by developing two different gene panels. We first analyzed 30 individuals with a small panel including 62 genes associated with LMND using IonTorrent-AmpliSeq target enrichment. Then, additional 65 individuals were tested with a broader panel encompassing up to 479 genes implicated in neuromuscular diseases (NMDs) with Agilent-SureSelect target enrichment. The NMD panel provided a higher diagnostic yield (33%) than the restricted LMND panel (13%). Nondiagnosed cases were further subjected to exome or genome sequencing. Our experience supports the use of gene panels covering a broad disease spectrum for diseases that are highly heterogeneous and clinically difficult to differentiate.

De novo mutations (DNMs) cause a large proportion of severe rare diseases of childhood. DNMs that occur early may result in mosaicism of both somatic and germ cells. Such early mutations can cause recurrence of disease. We scanned 1,007 sibling pairs from 251 families and identified 878 DNMs shared by siblings (ssDNMs) at 448 genomic sites. We estimated DNM recurrence probability based on parental mosaicism, sharing of DNMs among siblings, parent-of-origin, mutation type and genomic position. We detected 57.2% of ssDNMs in the parental blood. The recurrence probability of a DNM decreases by 2.27% per year for paternal DNMs and 1.78% per year for maternal DNMs. Maternal ssDNMs are more likely to be TC mutations than paternal ssDNMs, and less likely to be CT mutations. Depending on the properties of the DNM, the recurrence probability ranges from 0.011% to 28.5%. We have launched an online calculator to allow estimation of DNM recurrence probability for research purposes.

Background: Rare missense mutations in the gene encoding coatomer subunit alpha (COPA) have recently been shown to cause autoimmune interstitial lung, joint and kidney disease, also known as COPA syndrome, under a dominant mode of inheritance. Case presentation: Here we describe an Icelandic family with three affected individuals over two generations with a rare clinical presentation of lung and joint disease and a histological diagnosis of follicular bronchiolitis. We performed whole-genome sequencing (WGS) of the three affected as well as three unaffected members of the family, and searched for rare genotypes associated with disease using 30,067 sequenced Icelanders as a reference population. We assessed all coding and splicing variants, prioritizing variants in genes known to cause interstitial lung disease. We detected a heterozygous missense mutation, p.Glu241Lys, in the COPA gene, private to the affected family members. The mutation occurred de novo in the paternal germline of the index case and was absent from 30,067 Icelandic genomes and 141,353 individuals from the genome Aggregation Database (gnomAD). The mutation occurs within the conserved and functionally important WD40 domain of the COPA protein. Conclusions: This is the second report of the p.Glu241Lys mutation in COPA, indicating the recurrent nature of the mutation. The mutation was reported to co-segregate with COPA syndrome in a large family from the USA with five affected members, and classified as pathogenic. The two separate occurrences of the p.Glu241Lys mutation in cases and its absence from a large number of sequenced genomes confirms its role in the pathogenesis of the COPA syndrome. Keywords: COPA syndrome, Lung disease, Arthritis, Immune dysregulation, Case report

Background Epileptic encephalopathies are a group of childhood epilepsies that display high phenotypic and genetic heterogeneity. The recent, extensive use of next-generation sequencing has identified a large number of genes in epileptic encephalopathies, including UBA5 in which biallelic mutations were first described as pathogenic in 2016 (Colin E et al., Am J Hum Genet 99(3):695-703, 2016. Muona M et al., Am J Hum Genet 99(3):683-694, 2016). UBA5 encodes an activating enzyme for a post-translational modification mechanism known as ufmylation, and is the first gene from the ufmylation pathway that is linked to disease. Case presentation We sequenced the genomes of two sisters with early-onset epileptic encephalopathy along with their unaffected parents in an attempt to find a genetic cause for their condition. The sisters, born in 2004 and 2006, presented with infantile spasms at six months of age, which later progressed to recurrent, treatment-resistant seizures. We detected a compound heterozygous genotype in UBA5 in the sisters, a genotype not seen elsewhere in an Icelandic reference set of 30,067 individuals nor in public databases. One of the mutations, c.684G > A, is a paternally inherited exonic splicing mutation, occuring at the last nucleotide of exon 7 of UBA5. The mutation is predicted to disrupt the splice site, resulting in loss-of-function of one allele of UBA5. The second mutation is a maternally inherited missense mutation, p.Ala371Thr, previously reported as pathogenic when in compound heterozygosity with a loss-of-function mutation in UBA5 and is believed to produce a hypomorphic allele. Supportive of this, we have identified three adult Icelanders homozygous for the p.Ala371Thr mutation who show no signs of neurological disease. Conclusions We describe compound heterozygous mutations in the UBA5 gene in two sisters with early-onset epileptic encephalopathy. To our knowledge, this is the first description of mutations in UBA5 since the initial discovery that pathogenic biallelic variants in the gene cause early-onset epileptic encephalopathy. We further provide confirmatory evidence that p.Ala371Thr is a hypomorphic mutation, by presenting three adult homozygotes who show no signs of neurological disease. Electronic supplementary material The online version of this article (10.1186/s12881-017-0466-8) contains supplementary material, which is available to authorized users.

Understanding of sequence diversity is the cornerstone of analysis of genetic disorders, population genetics, and evolutionary biology. Here, we present an update of our sequencing set to 15,220 Icelanders who we sequenced to an average genome-wide coverage of 34X. We identified 39,020,168 autosomal variants passing GATK filters: 31,079,378 SNPs and 7,940,790 indels. Calling de novo mutations (DNMs) is a formidable challenge given the high false positive rate in sequencing datasets relative to the mutation rate. Here we addressed this issue by using segregation of alleles in three-generation families. Using this transmission assay, we controlled the false positive rate and identified 108,778 high quality DNMs. Furthermore, we used our extended family structure and read pair tracing of DNMs to a panel of phased SNPs, to determine the parent of origin of 42,961 DNMs.

IL-33 is a tissue-derived cytokine that induces and amplifies eosinophilic inflammation and has emerged as a promising new drug target for asthma and allergic disease. Common variants at IL33 and IL1RL1, encoding the IL-33 receptor ST2, associate with eosinophil counts and asthma. Through whole-genome sequencing and imputation into the Icelandic population, we found a rare variant in IL33 (NM_001199640:exon7:c.487-1G>C (rs146597587-C), allele frequency = 0.65%) that disrupts a canonical splice acceptor site before the last coding exon. It is also found at low frequency in European populations. rs146597587-C associates with lower eosinophil counts (β = -0.21 SD, P = 2.5×10–16, N = 103,104), and reduced risk of asthma in Europeans (OR = 0.47; 95%CI: 0.32, 0.70, P = 1.8×10–4, N cases = 6,465, N controls = 302,977). Heterozygotes have about 40% lower total IL33 mRNA expression than non-carriers and allele-specific analysis based on RNA sequencing and phased genotypes shows that only 20% of the total expression is from the mutated chromosome. In half of those transcripts the mutation causes retention of the last intron, predicted to result in a premature stop codon that leads to truncation of 66 amino acids. The truncated IL-33 has normal intracellular localization but neither binds IL-33R/ST2 nor activates ST2-expressing cells. Together these data demonstrate that rs146597587-C is a loss of function mutation and support the hypothesis that IL-33 haploinsufficiency protects against asthma., The Dutch asthma study was supported by the Netherlands Asthma Foundation grant AF (AF 95.09, AF 98.48, AF 3.2.02.51 and AF 3.2.07.015) and a grant from the University Medical Center Groningen. The Vlagtwedde-Vlaardingen cohort study was supported by the Ministry of Health and Environmental Hygiene of the Netherlands and the Netherlands Asthma Fund (grant 187) and the Netherlands Asthma Fund grant no. 3.2.02.51, the Stichting Astma Bestrijding, BBMRI-NL (Complementiation project), and the European Respiratory Society COPD research award 2011 to H.M. Boezen. COPSAC is funded by private and public research funds all listed on www.copsac.com. The Lundbeck Foundation; Danish State Budget; Danish Council for Strategic Research; Danish Council for Independent Research and The Capital Region Research Foundation have provided core support for COPSAC. The Denmark-1 study has been supported with an unrestricted grant from AstraZeneca and ALK-Abello, Denmark. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Data supporting the findings of this study are available within the article and its Supplementary Information files. Summary level data of markers tested for association is described in Scientific data as of 2015 with the identifier http://dx.doi.org/10.1038/sdata.2015.11 (ref. 64). Whole-genome sequencing summary data are available at the European Variant Archive (EVA, https://www.ebi.ac.uk/eva/) under the accession code PRJEB8636. GTEx data (http://www.gtexportal.org) accessed 8 August 2016; LD Score Database (ftp://atguftp.mgh.harvard.edu/brendan/1k_eur_r2_hm3snps_se_weights.RDS) accessed 23 June 2015; GWAS catalogue (https://www.ebi.ac.u/gwas/home) accessed 8 September 2016., Lumbar disc herniation (LDH) is common and often debilitating. Microdiscectomy of herniated lumbar discs (LDHsurg) is performed on the most severe cases to resolve the resulting sciatica. Here we perform a genome-wide association study on 4,748 LDHsurg cases and 282,590 population controls and discover 37 highly correlated markers associating with LDHsurg at 8q24.21 (between CCDC26 and GSDMC), represented by rs6651255[C] (OR=0.81; P=5.6 × 10−12) with a stronger effect among younger patients than older. As rs6651255[C] also associates with height, we performed a Mendelian randomization analysis using height polygenic risk scores as instruments to estimate the effect of height on LDHsurg risk, and found that the marker's association with LDHsurg is much greater than predicted by its effect on height. In light of presented findings, we speculate that the effect of rs6651255 on LDHsurg is driven by susceptibility to developing severe and persistent sciatica upon LDH., The financial support from the European Commission to the NeuroPain project (FP7#HEALTH-2013-602891-2) and the National Institutes of Health (R01DE022905) is acknowledged.

Through whole-genome sequencing of 2,630 Icelanders and imputation into 11,114 Icelandic cases and 267,140 controls followed by testing in Danish and Iranian samples, we discovered 4 previously unreported variants affecting risk of type 2 diabetes (T2D). A low-frequency (1.47%) variant in intron 1 of CCND2, rs76895963[G], reduces risk of T2D by half (odds ratio (OR) = 0.53, P = 5.0 x 10(-21)) and is correlated with increased CCND2 expression. Notably, this variant is also associated with both greater height and higher body mass index (1.17 cm per allele, P = 5.5 x 10(-12) and 0.56 kg/m(2) per allele, P = 6.5 x 10(-7), respectively). In addition, two missense variants in PAM, encoding p. Asp563Gly (frequency of 4.98%) and p. Ser539Trp (frequency of 0.65%), confer moderately higher risk of T2D (OR = 1.23, P = 3.9 x 10(-10) and OR = 1.47, P = 1.7 x 10(-5), respectively), and a rare (0.20%) frameshift variant in PDX1, encoding p.Gly218Alafs*12, associates with high risk of T2D (OR = 2.27, P = 7.3 x 10(-7)).

Mutation of the DNA molecule is one of the most fundamental processes in biology. In this study, we use 283 parent-offspring trios to estimate the rate of mutation for both single nucleotide variants (SNVs) and short length variants (indels) in humans and examine the mutation process. We found 17812 SNVs, corresponding to a mutation rate of 1.29 × 10−8 per position per generation (PPPG) and 1282 indels corresponding to a rate of 9.29 × 10−10 PPPG. We estimate that around 3% of human de novo SNVs are part of a multi-nucleotide mutation (MNM), with 558 (3.1%) of mutations positioned less than 20kb from another mutation in the same individual (median distance of 525bp). The rate of de novo mutations is greater in late replicating regions (p = 8.29 × 10−19) and nearer recombination events (p = 0.0038) than elsewhere in the genome., Author Summary In each generation new genetic variants are introduced by mutations. In this study we use whole genome sequence data from Icelandic families to directly observe such new mutations. Our estimate of the mutation rate implies that a newborn with 30-year-old parents will on average carry 75 new SNV mutations and 6 new short indel mutations. We observe that some of the found mutations occur much closer together than would be expected by chance. Our analysis shows that mutational hotspots cannot explain this clustering, instead the clustering mutations are likely created by a single mutational event. We observe a different composition for mutations that cluster very close together compared to more distant clustering mutations. This suggests that there is likely more than one type of underlying mutational mechanism creating the multi nucleotide mutation events. We furthermore observe a higher mutation rate near recombination events but find that this effect cannot explain the large number of clustering mutations that we observe.

Creatine kinase (CK) and lactate dehydrogenase (LDH) are widely used markers of tissue damage. To search for sequence variants influencing serum levels of CK and LDH, 28.3 million sequence variants identified through whole-genome sequencing of 2,636 Icelanders were imputed into 63,159 and 98,585 people with CK and LDH measurements, respectively. Here we describe 13 variants associating with serum CK and 16 with LDH levels, including four that associate with both. Among those, 15 are non-synonymous variants and 12 have a minor allele frequency below 5%. We report sequence variants in genes encoding the enzymes being measured (CKM and LDHA), as well as in genes linked to muscular (ANO5) and immune/inflammatory function (CD163/CD163L1, CSF1, CFH, HLA-DQB1, LILRB5, NINJ1 and STAB1). A number of the genes are linked to the mononuclear/phagocyte system and clearance of enzymes from the serum. This highlights the variety in the sources of normal diversity in serum levels of enzymes., Creatine kinase (CK) and lactate dehydrogenase (LDH) are biomarkers of tissue damages including myopathy and myocardial infarction. Here, Patrick Sulem and colleagues perform a genome-wide association study to identify common and rare genetic variants that associates with serum CK or LDH levels.

To access publisher's full text version of this article click on the hyperlink at the bottom of the page Mycobacterium tuberculosis infections cause 9 million new tuberculosis cases and 1.5 million deaths annually. To identify variants conferring risk of tuberculosis, we tested 28.3 million variants identified through whole-genome sequencing of 2,636 Icelanders for association with tuberculosis (8,162 cases and 277,643 controls), pulmonary tuberculosis (PTB) and M. tuberculosis infection. We found association of three variants in the region harboring genes encoding the class II human leukocyte antigens (HLAs): rs557011[T] (minor allele frequency (MAF) = 40.2%), associated with M. tuberculosis infection (odds ratio (OR) = 1.14, P = 3.1 × 10(-13)) and PTB (OR = 1.25, P = 5.8 × 10(-12)), and rs9271378[G] (MAF = 32.5%), associated with PTB (OR = 0.78, P = 2.5 × 10(-12))-both located between HLA-DQA1 and HLA-DRB1-and a missense variant encoding p.Ala210Thr in HLA-DQA1 (MAF = 19.1%, rs9272785), associated with M. tuberculosis infection (P = 9.3 × 10(-9), OR = 1.14). We replicated association of these variants with PTB in samples of European ancestry from Russia and Croatia (P < 5.9 × 10(-4)). These findings show that the HLA class II region contributes to genetic risk of tuberculosis, possibly through reduced presentation of protective M. tuberculosis antigens to T cells. US National Institute of Allergy and Infectious Diseases HHSN266200400064C UK Wellcome Trust 088838/Z/09/Z 095198/Z/10/Z info:eu-repo/grantAgreement/EC/FP7/201483 European Research Council Starting 260477 Royal Society UF0763346 RG090638 Wellcome Trust Senior Research fellowship National Institute for Health Research (NIHR) Cambridge Biomedical Research Centre

Pyrroloquinoline quinone [4,5-dihydro-4,5-dioxo-1H-pyrrolo[2,3-f]quinoline-2,7,9-tricarboxylic acid (PQQ)] is a bacterial vitamin that serves as a cofactor in numerous alcohol dehydrogenases. Its biosynthesis in Klebsiella pneumoniae is facilitated by six genes, pqqABCDEF, and proceeds by an unknown pathway. The protein encoded by pqqC catalyzes the final step of PQQ formation, which involves a ring closure and an overall eight-electron oxidation of 3a-(2-amino-2-carboxyethyl)-4,5-dioxo-4,5,6,7,8,9-hexahydroquinoline-7,9-dicarboxylic acid (AHQQ) in the absence of a redox-active metal or cofactor. A recent crystal structure has implicated numerous PQQ-PqqC interactions [Magnusson et al. (2004) Proc. Natl. Acad. Sci. U.S.A. 101, 7913-7918]. To investigate the mechanism of the PqqC reaction, the active site residue His84 has been mutated to H84A and H84N, and the kinetic and spectroscopic properties have been compared to each other and the wild-type enzyme using aerobic and anaerobic conditions. Both mutants form PQQ under aerobic conditions with rate constants of 0.09 min-1 and 0.056 min-1 relative to 0.34 min-1 for the wild-type enzyme. In addition to the initial E-AHQQ complex (532-536 nm) and the product E-PQQ complex (346-366 nm), a number of spectral intermediates are observed between 316 and 344 nm. The anaerobic reaction is particularly informative, showing that while mixing of H84N with AHQQ leads to a 344 nm intermediate, this is unable to proceed to a final 318 nm species; by contrast H84A forms the 344 nm species as a precursor to the 318 nm species. In the context of the proposed chemical mechanism for PqqC [Magnusson et al. (2004) Proc. Natl. Acad. Sci. U.S.A. 101, 7913-7918], we assign the 344 nm intermediate to a quinoid species and the 318 nm intermediate to an initial quinol species. The proposed role of H84 is as a proton donor to the oxyanion of the quinoid species such that subsequent C-H bond cleavage can occur to form a monoanionic quinol. In the absence of a proton donor (as occurs in H84N), the normal reaction path is precluded as this would require formation of an unstable, dianionic species. Unlike H84N, H84A appears to be small enough to allow entry of active site water, which is postulated to adopt the role of active site proton donor.

We have accumulated considerable data on the genetic makeup of the Icelandic population by sequencing the whole genomes of 2,636 Icelanders to depth of at least 10X and by chip genotyping 101,584 more. The sequencing was done with Illumina technology. The median sequencing depth was 20X and 909 individuals were sequenced to a depth of at least 30X. We found 20 million single nucleotide polymorphisms (SNPs) and 1.5 million insertions/deletions (indels) that passed stringent quality control. Almost all the common SNPs (derived allele frequency (DAF) over 2%) that we identified in Iceland have been observed by either dbSNP (build 137) or the Exome Sequencing Project (ESP) while only 60 and 20% of rare (DAF

A common characteristic of AdoCbl1-dependent enzymes is the generation and utilization of enzyme-bound organic radicals (1–3). The Co–C5′ bond of AdoCbl is relatively weak having a bond dissociation energy of ~ 30 kcal mol−1 (4). A prevalent mechanism for AdoCbl-dependent enzymes begins with the homolytic cleavage of the Co–C5′ bond to give cob(II)alamin and the 5′-deoxyadenosyl radical (1). The 5′-deoxyadenosyl radical being an unstabilized, primary alkyl radical is expected to have a short lifetime, low concentration, and a high reactivity towards H atom abstraction. Thus far, these properties of the 5′-deoxyadenosyl radical have precluded direct observation of the radical by EPR spectroscopy. Cleavage of the Co–C5′ bond is kinetically coupled to H atom abstraction from a substrate molecule in the reactions of methylmalonyl-CoA mutase (5), glutamate mutase (6), ethanolamine ammonia-lyase (7), and diol dehydrase (8). These kinetic experiments have highlighted the transient nature of the 5′-deoxyadenosyl radical species. In the case of the AdoCbl-dependent Class II ribonucleoside triphosphate reductase the 5′-deoxyadenosyl radical abstracts an H atom from the thiol of Cys408 of the protein (9). An allosteric effector-dependent epimerization of chiral [5′-2H1]AdoCbl by the C408A and C408S mutant forms of the enzyme requires the transient existence of the 5′-deoxyadenosyl radical (10). DDH (diol dehydrase, EC:4.2.1.28) uses 2- and 3-carbon vicinal diols as substrates and catalyzes their transformation into aldehydes (8, 11). The minimal mechanism for DDH with 1,2-propanediol as the substrate is depicted in Scheme 1. Reversible cleavage of the Co–C5′ bond of AdoCbl generates the 5′-deoxyadenosyl radical a, which abstracts a hydrogen atom from C1 of 1,2-propanediol to produce 5′-deoxyadenosine b and the substrate radical c. The substrate radical undergoes isomerization to the product-related radical d by an unknown mechanism. Hypothetical mechanisms for the radical rearrangement include the formation of a radical cation (12) or a radical anion intermediate (13). Concerted mechanisms having transition state stabilization by partial proton transfer or deprotonation of the spectator −OH group have also been discussed (14, 15). Crystal structures of DDH reveal that the vicinyl −OH groups of 1,2-propanediol are coordinated to K+, and K+ may be involved in catalysis (16, 17). Magnetic interactions between active site radicals and the nuclear spins of 39K and 203,205Tl have, however, escaped detection by EPR and ENDOR (18). Moreover, site directed mutagenesis has revealed essential roles of Glu-α170 and Asp-α335 (19) in putative H-bonding contacts with the vicinal −OH groups of the substrate. Hydrogen transfer from 5′-deoxyadenosine b to the radical d regenerates the 5′-deoxyadenosyl radical, which recombines with cob(II)alamin. The enzyme catalyzes the dehydration of the gem-diol to complete the reaction (20). Reactions of the DDH-AdoCbl complex with substrates or inhibitors elicit EPR signals corresponding to organic radicals spin coupled to the low spin Co2+ of cob(II)alamin (18, 21–26). Scheme 1 AnAdoCbl, shown in Scheme 2, is an analog of AdoCbl. AnAdoCbl has a Co–C5′ bond dissociation energy which is ~ 6 kcal mol−1 weaker than the corresponding bond in AdoCbl (27). The weakened Co–C5′ bond in anAdoCbl is attributed to allylic stabilization in the anhydroadenosyl radical relative to the 5′-deoxyadenosyl radical. Structures of the two radicals are compared in Scheme 3. AnAdoCbl is a functional coenzyme for DDH, displaying a turnover rate 0.02% of that with AdoCbl as coenzyme (28). Binding of anAdoCbl to DDH induces Co–C5′ bond cleavage. The Co–C5′ bond cleavage is observed spectrophotometrically as the formation of cob(II)alamin either in the presence or absence of the substrate (28). However, the enzyme undergoes slow inactivation as observed by the formation of cob(III)alamin, irrespective of the presence of the substrate. Inactivation is a consequence of electron transfer from cob(II)alamin to the anhydroadenosyl radical and subsequent quenching of the allylic anion by a solvent derived proton (28). Scheme 2 Scheme 3 Cleavage of the Co–C5′ bond of anAdoCbl generates the anhydroadenosyl radical and cob(II)alamin in the active site of DDH. The present paper presents EPR evidence for the presence of strongly-coupled spin-triplets involving the anhydroadenosyl radical and the low spin Co2+ of cob(II)alamin in complexes of DDH with the cofactor analog in the presence and absence of substrate. Analysis of the spin-triplet EPR spectra by simulation and interpretation of the ZFS in terms of molecular geometries is presented.

Loss-of-function mutations cause many mendelian diseases. Here we aimed to create a catalog of autosomal genes that are completely knocked out in humans by rare loss-of-function mutations. We sequenced the whole genomes of 2,636 Icelanders and imputed the sequence variants identified in this set into 101,584 additional chip-genotyped and phased Icelanders. We found a total of 6,795 autosomal loss-of-function SNPs and indels in 4,924 genes. Of the genotyped Icelanders, 7.7% are homozygotes or compound heterozygotes for loss-of-function mutations with a minor allele frequency (MAF) below 2% in 1,171 genes (complete knockouts). Genes that are highly expressed in the brain are less often completely knocked out than other genes. Homozygous loss-of-function offspring of two heterozygous parents occurred less frequently than expected (deficit of 136 per 10,000 transmissions for variants with MAF

Here we describe the insights gained from sequencing the whole genomes of 2,636 Icelanders to a median depth of 20×. We found 20 million SNPs and 1.5 million insertions-deletions (indels). We describe the density and frequency spectra of sequence variants in relation to their functional annotation, gene position, pathway and conservation score. We demonstrate an excess of homozygosity and rare protein-coding variants in Iceland. We imputed these variants into 104,220 individuals down to a minor allele frequency of 0.1% and found a recessive frameshift mutation in MYL4 that causes early-onset atrial fibrillation, several mutations in ABCB4 that increase risk of liver diseases and an intronic variant in GNAS associating with increased thyroid-stimulating hormone levels when maternally inherited. These data provide a study design that can be used to determine how variation in the sequence of the human genome gives rise to human diversity.

The mechanism of propagation of the radical center between the cofactor, substrate, and product in the adenosylcobalamin- (AdoCbl) dependent reaction of ethanolamine ammonia-lyase has been probed by pulsed electron nuclear double resonance (ENDOR) spectroscopy. The radical of S-2-aminopropanol, which appears in the steady state of the reaction, was used in ENDOR experiments to determine the nuclear spin transition frequencies of (2)H introduced from either deuterated substrate or deuterated coenzyme and of (13)C introduced into the ribosyl moiety of AdoCbl. A (2)H doublet (1.4 MHz splitting) was observed centered about the Larmor frequency of (2)H. Identical ENDOR frequencies were observed for (2)H irrespective of its mode of introduction into the complex. A (13)C doublet ENDOR signal was observed from samples prepared with [U-(13)C-ribosyl]-AdoCbl. The (13)C coupling tensor obtained from the ENDOR powder pattern shows that the (13)C has scalar as well as dipole-dipole coupling to the unpaired electron located at C1 of S-2-aminopropanol. The dipole-dipole coupling is consistent with a distance of 3.4+/-0.2 A between C1 of the radical and C5' of the labeled cofactor component. These results establish that the C5' carbon of the 5'-deoxyadenosyl radical moves approximately 7 A from its position as part of AdoCbl to a position where it is in contact with C1 of the substrate which lies approximately 12 A from the Co(2+) of cob(II)alamin. These findings are also consistent with the contention that 5'-deoxyadenosine is the sole mediator of hydrogen transfers in ethanolamine ammonia-lyase.

The question of how coenzyme B12-dependent enzymes facilitate the cleavage of the Co−C bond of the cofactor is of interest. We have synthesized an analogue of 5‘-deoxyadenosylcobalamin (AdoCbl1) designed to stabilize the 5‘-deoxyadenosyl radical (5‘-deoxyadenosine-5‘-yl) that is produced upon homolysis of the Co−C bond. By replacement of the upper axial ligand of AdoCbl by a 3‘,4‘-anhydro-5‘-deoxyadenosyl moiety, the radical formed on the nucleoside analogue is stabilized by allylic delocalization. The compound, 5‘-deoxy- 3‘,4‘-anhydroadenosylcobalamin (3‘,4‘-anAdoCbl) was synthesized by chemical and enzymatic methods. The final step was coupling of cob(I)alamin and 3‘,4‘-anhydroATP catalyzed by CobA, an ATP:corrinoid adenosyltransferase. 3‘,4‘-anAdoCbl displays interesting properties. The compound has not been purified to homogeneity due to its thermal and oxygen sensitivity. It was characterized by UV−vis spectroscopy, ESI-MS, and NMR spectroscopy. The bond dissociation energy of the Co−C bond of the an...

Biogenesis of pyrroloquinoline quinone (PQQ) in Klebsiella pneumoniae requires the expression of six genes (pqqA-F). One of these genes (pqqE) encodes a 43 kDa protein (PqqE) that plays a role in the initial steps in PQQ formation [Veletrop, J. S., et al. (1995) J. Bacteriol. 177, 5088-5098]. PqqE contains two highly conserved cysteine motifs at the N- and C-termini, with the N-terminal motif comprised of a CX(3)CX(2)C consensus sequence that is unique to a family of proteins known as radical S-adenosyl-l-methionine (SAM) enzymes [Sofia, H. J., et al. (2001) Nucleic Acids Res. 29, 1097-1106]. PqqE from K. pneumoniae was cloned into Escherichia coli and expressed as the native protein and with an N-terminal His(6) tag. Anaerobic expression and purification of the His(6)-tagged PqqE results in an enzyme with a brownish-red hue indicative of Fe-S cluster formation. Spectroscopic and physical analyses indicate that PqqE contains a mixture of Fe-S clusters, with the predominant form of the enzyme containing two [4Fe-4S] clusters. PqqE isolated anaerobically yields an active enzyme capable of cleaving SAM to methionine and 5'-deoxyadenosine in an uncoupled reaction (k(obs) = 0.011 +/- 0.001 min(-1)). In this reaction, the 5'-deoxyadenosyl radical either abstracts a hydrogen atom from a solvent accessible position in the enzyme or obtains a proton and electron from buffer. The putative PQQ substrate PqqA has not yet been shown to be modified by PqqE, implying that PqqA must be modified before becoming the substrate for PqqE and/or that another protein in the biosynthetic pathway is critical for the initial steps in PQQ biogenesis.

Pyrroloquinoline quinone (PQQ) is a tricyclic o-quinone, which serves as a cofactor in several enzyme-catalyzed redox reactions in certain bacteria. PQQ is also important for human health, and its role as a vitamin in mammals has recently been suggested. Although much is known about the function of enzymes that use PQQ as cofactor, relatively little is known about the biosynthesis of this coenzyme. Six gene products in Klebsiella pneumoniae (PqqA-F) are involved in PQQ biosynthesis, and PqqC has been shown to catalyze the last step in the pathway. The chemical structure of the substrate for PqqC has remained elusive and has hampered our understanding of the nature of this reaction. In this report we describe the purification and structure of the substrate as deduced by a number of spectroscopic and chemical methods. The substrate is 3a-(2-amino-2-carboxyethyl)-4,5-dioxo-4,5,6,7,8,9-hexahydroquinoline-7,9-dicarboxylic acid-a fully reduced derivative of PQQ, which has not undergone ring cyclization. These results show that PqqC catalyzes a novel reaction, which involves ring closure and an amazing eight-electron oxidation of the substrate.

The biosynthesis of pyrroloquinoline quinone (PQQ), a vitamin and redox cofactor of quinoprotein dehydrogenases, is facilitated by an unknown pathway that requires the expression of six genes, pqqA to -F. PqqC, the protein encoded by pqqC, catalyzes the final step in the pathway in a reaction that involves ring cyclization and eight-electron oxidation of 3a-(2-amino-2-carboxyethyl)-4,5-dioxo-4,5,6,7,8,9-hexahydroquinoline-7,9-dicarboxylic-acid to PQQ. Herein, we describe the crystal structures of PqqC and its complex with PQQ and determine the stoichiometry of H2O2 formation and O2 uptake during the reaction. The PqqC structure(s) reveals a compact seven-helix bundle that provides the scaffold for a positively charged active site cavity. Product binding induces a large conformational change, which results in the active site recruitment of amino acid side chains proposed to play key roles in the catalytic mechanism. PqqC is unusual in that it transfers redox equivalents to molecular oxygen without the assistance of a redox active metal or cofactor. The structure of the enzyme-product complex shows additional electron density next to R179 and C5 of PQQ, which can be modeled as O2 or H2O2, indicating a site for oxygen binding. We propose a reaction sequence that involves base-catalyzed cyclization and a series of quinone-quinol tautomerizations that are followed by cycles of O2/H2O2-mediated oxidations. Pyrroloquinoline quinone [4,5-dihydro-4,5-dioxo-1H-pyrrolo-[2,3-f]quinoline-2,7,9-tricarboxylic acid; PQQ (Fig. 1)] is an aromatic, tricyclic ortho-quinone that serves as the redox cofactor for several bacterial dehydrogenases. Among the best known examples are methanol dehydrogenase and glucose dehydrogenase (1, 2). PQQ belongs to the family of quinone cofactors that has been recognized as the third class of redox cofactors following pyridine nucleotide- and flavin-dependent cofactors (3). Although plants and animals do not produce PQQ themselves, PQQ has invoked considerable interest because of its presence in human milk and its remarkable antioxidant properties (4–6). Recently, the first potential eukaryotic PQQ-dependent enzyme [aminoadipic 6-semialdehyde-dehydrogenase (AASDH; U26)] has been identified, indicating that PQQ may function as a vitamin in mammals as well (7)., R.S. was supported by Austrian Science Fund Fellowship J2209-B04, J.P.K. by National Institutes of Health Grant GM39296, and O.T.M. by a postdoctoral fellowship from the Miller Institute for Basic Research in Science, University of California, Berkeley.

S-3',4'-anhydroadenosyl-l-methionine is an analogue of the S-adenosyl-l-methionine coenzyme. Here we report on a rapid solvent exchange of the methylene protons at the 5'-position of this analogue. The rate of H/D exchange was measured by nuclear magnetic resonance spectroscopy under buffered conditions in deuterium oxide. The reaction is specific base catalyzed and displays a second-order rate constant of 2 x 10(4) M(-1) s(-1), which corresponds to a rate enhancement of 10(12) compared to solvent exchange of alpha-methylene protons in acyclic, aliphatic sulfonium ions. No other carbon bonded hydrogens in the molecule exchange with solvent under the experimental conditions. Allylic stabilization of a carbanionic-like transition state for the solvent exchange process can account for these results. Solvent exchange under these mild conditions provides a simple way to prepare a 5'-2H-labeled form of the coenzyme analogue.

3',4'-Anhydroadenosylcobalamin (anAdoCbl) is an analogue of the adenosylcobalamin (AdoCbl) coenzyme (Magnusson, O.Th., and Frey, P. A. (2000) J. Am. Chem. Soc. 122, 8807-8813). This compound supports activity for diol dehydrase at 0.02% of that observed with AdoCbl. In a side reaction, however, anAdoCbl induces suicide inactivation by an electron-transfer mechanism. Homolytic cleavage of the Co-C bond of anAdoCbl at the active site of diol dehydrase was observed by spectrophotometric detection of cob(II)alamin. Anaerobic conversion of enzyme bound cob(II)alamin to cob(III)alamin, both in the absence and presence of substrate, indicates that the coenzyme derived 5'-deoxy-3',4'-anhydroadenosine-5'-yl serves as the oxidizing agent. This hypothesis is supported by the stoichiometric formation of 3',5'-dideoxyadenosine-4',5'-ene as the nucleoside cleavage product, as determined by high-performance liquid chromatography, mass spectrometry, and nuclear magnetic resonance spectroscopy. Experiments performed in deuterium oxide show that a single solvent exchangeable proton is incorporated into the product. These data are consistent with the intermediate formation of a transient allylic anion formed after one electron transfer from cob(II)alamin to the allylic 5'-deoxy-3',4'-anhydroadenosyl radical. Selective protonation at C3' was demonstrated by spectroscopic characterization of the purified product. This study provides an example of suicide inactivation of a radical enzyme brought about by a side reaction of an analogue of the radical intermediate.

An allylic analogue of the 5'-deoxyadenosyl radical has been characterized at the active site of lysine 2,3-aminomutase (LAM) by electron paramagnetic resonance (EPR) spectroscopy. The anhydroadenosyl radical, 5'-deoxy-3',4'-anhydroadenosine-5'-yl, is a surrogate of the less stable 5'-deoxyadenosyl radical, which has never been observed but has been postulated to be a radical intermediate in the catalytic cycles of a number of enzymes. An earlier communication [Magnusson, O.Th., Reed, G. H., and Frey, P. A. (1999) J. Am. Chem. Soc. 121, 9764-9765] included the initial spectroscopic identification at 77 K of the radical, which is formed upon replacement of S-adenosylmethionine by S-3',4'-anhydroadenosylmethionine as a coenzyme for LAM. The electron paramagnetic resonance spectrum of the radical changes dramatically between 77 and 4.5 K. This unusual temperature dependence is attributed to a spin-spin interaction between the radical and thermally populated, higher spin states of the [4Fe-4S]+2 center, which is diamagnetic at 4.5 K. The EPR spectra of the radical at 4.5 K have been analyzed using isotopic substitutions and simulations. Analysis of the nuclear hyperfine splitting shows that the unpaired spin is distributed equally between C5'- and C3'- as expected for an allylic radical. Hyperfine splitting from the beta-proton at C-2'(H) shows that the dihedral angle to the p(z)-orbital at C-3' is approximately 37 degrees. This conformation is in good agreement with a structural model of the radical. The rate of formation of the allylic radical shows that it is kinetically competent as an intermediate. Measurements of 2H kinetic isotope effects indicate that with lysine as the substrate, the rate-limiting steps follow initial reductive cleavage of the coenzyme analogue.