A comparison of methods for detecting DNA methylation from long-read sequencing of human genomes.

BACKGROUND: Long-read sequencing can enable the detection of base modifications, such as CpG methylation, in single molecules of DNA. The most commonly used methods for long-read sequencing are nanopore developed by Oxford Nanopore Technologies (ONT) and single molecule real-time (SMRT) sequencing developed by Pacific Bioscience (PacBio). In this study, we systematically compare the performance of CpG methylation detection from long-read sequencing. RESULTS: We demonstrate that CpG methylation detection from 7179 nanopore-sequenced DNA samples is highly accurate and consistent with 132 oxidative bisulfite-sequenced (oxBS) samples, isolated from the same blood draws. We introduce quality filters for CpGs that further enhance the accuracy of CpG methylation detection from nanopore-sequenced DNA, while removing at most 30% of CpGs. We evaluate the per-site performance of CpG methylation detection across different genomic features and CpG methylation rates and demonstrate how the latest R10.4 flowcell chemistry and base-calling algorithms improve methylation detection from nanopore sequencing. Additionally, we show how the methylation detection of 50 SMRT-sequenced genomes compares to nanopore sequencing and oxBS. CONCLUSIONS: This study provides the first systematic comparison of CpG methylation detection tools for long-read sequencing methods. We compare two commonly used computational methods for the detection of CpG methylation in a large number of nanopore genomes, including samples sequenced using the latest R10.4 nanopore flowcell chemistry and 50 SMRT sequenced samples. We provide insights into the strengths and limitations of each sequencing method as well as recommendations for standardization and evaluation of tools designed for genome-scale modified base detection using long-read sequencing.

Rare variants with large effects provide functional insights into the pathology of migraine subtypes, with and without aura.

Migraine is a complex neurovascular disease with a range of severity and symptoms, yet mostly studied as one phenotype in genome-wide association studies (GWAS). Here we combine large GWAS datasets from six European populations to study the main migraine subtypes, migraine with aura (MA) and migraine without aura (MO). We identified four new MA-associated variants (in PRRT2, PALMD, ABO and LRRK2) and classified 13 MO-associated variants. Rare variants with large effects highlight three genes. A rare frameshift variant in brain-expressed PRRT2 confers large risk of MA and epilepsy, but not MO. A burden test of rare loss-of-function variants in SCN11A, encoding a neuron-expressed sodium channel with a key role in pain sensation, shows strong protection against migraine. Finally, a rare variant with cis-regulatory effects on KCNK5 confers large protection against migraine and brain aneurysms. Our findings offer new insights with therapeutic potential into the complex biology of migraine and its subtypes.

Large-scale plasma proteomics comparisons through genetics and disease associations.

High-throughput proteomics platforms measuring thousands of proteins in plasma combined with genomic and phenotypic information have the power to bridge the gap between the genome and diseases. Here we performed association studies of Olink Explore 3072 data generated by the UK Biobank Pharma Proteomics Project1 on plasma samples from more than 50,000 UK Biobank participants with phenotypic and genotypic data, stratifying on British or Irish, African and South Asian ancestries. We compared the results with those of a SomaScan v4 study on plasma from 36,000 Icelandic people2, for 1,514 of whom Olink data were also available. We found modest correlation between the two platforms. Although cis protein quantitative trait loci were detected for a similar absolute number of assays on the two platforms (2,101 on Olink versus 2,120 on SomaScan), the proportion of assays with such supporting evidence for assay performance was higher on the Olink platform (72% versus 43%). A considerable number of proteins had genomic associations that differed between the platforms. We provide examples where differences between platforms may influence conclusions drawn from the integration of protein levels with the study of diseases. We demonstrate how leveraging the diverse ancestries of participants in the UK Biobank helps to detect novel associations and refine genomic location. Our results show the value of the information provided by the two most commonly used high-throughput proteomics platforms and demonstrate the differences between them that at times provides useful complementarity.

Evaluation of Large-Scale Proteomics for Prediction of Cardiovascular Events.

Importance: Whether protein risk scores derived from a single plasma sample could be useful for risk assessment for atherosclerotic cardiovascular disease (ASCVD), in conjunction with clinical risk factors and polygenic risk scores, is uncertain. Objective: To develop protein risk scores for ASCVD risk prediction and compare them to clinical risk factors and polygenic risk scores in primary and secondary event populations. Design, Setting, and Participants: The primary analysis was a retrospective study of primary events among 13 540 individuals in Iceland (aged 40-75 years) with proteomics data and no history of major ASCVD events at recruitment (study duration, August 23, 2000 until October 26, 2006; follow-up through 2018). We also analyzed a secondary event population from a randomized, double-blind lipid-lowering clinical trial (2013-2016), consisting of individuals with stable ASCVD receiving statin therapy and for whom proteomic data were available for 6791 individuals. Exposures: Protein risk scores (based on 4963 plasma protein levels and developed in a training set in the primary event population); polygenic risk scores for coronary artery disease and stroke; and clinical risk factors that included age, sex, statin use, hypertension treatment, type 2 diabetes, body mass index, and smoking status at the time of plasma sampling. Main Outcomes and Measures: Outcomes were composites of myocardial infarction, stroke, and coronary heart disease death or cardiovascular death. Performance was evaluated using Cox survival models and measures of discrimination and reclassification that accounted for the competing risk of non-ASCVD death. Results: In the primary event population test set (4018 individuals [59.0% women]; 465 events; median follow-up, 15.8 years), the protein risk score had a hazard ratio (HR) of 1.93 per SD (95% CI, 1.75 to 2.13). Addition of protein risk score and polygenic risk scores significantly increased the C index when added to a clinical risk factor model (C index change, 0.022 [95% CI, 0.007 to 0.038]). Addition of the protein risk score alone to a clinical risk factor model also led to a significantly increased C index (difference, 0.014 [95% CI, 0.002 to 0.028]). Among White individuals in the secondary event population (6307 participants; 432 events; median follow-up, 2.2 years), the protein risk score had an HR of 1.62 per SD (95% CI, 1.48 to 1.79) and significantly increased C index when added to a clinical risk factor model (C index change, 0.026 [95% CI, 0.011 to 0.042]). The protein risk score was significantly associated with major adverse cardiovascular events among individuals of African and Asian ancestries in the secondary event population. Conclusions and Relevance: A protein risk score was significantly associated with ASCVD events in primary and secondary event populations. When added to clinical risk factors, the protein risk score and polygenic risk score both provided statistically significant but modest improvement in discrimination.

DanMAC5: a browser of aggregated sequence variants from 8,671 whole genome sequenced Danish individuals.

OBJECTIVES: Allele counts of sequence variants obtained by whole genome sequencing (WGS) often play a central role in interpreting the results of genetic and genomic research. However, such variant counts are not readily available for individuals in the Danish population. Here, we present a dataset with allele counts for sequence variants (single nucleotide variants (SNVs) and indels) identified from WGS of 8,671 (5,418 females) individuals from the Danish population. The data resource is based on WGS data from three independent research projects aimed at assessing genetic risk factors for cardiovascular, psychiatric, and headache disorders. To enable the sharing of information on sequence variation in Danish individuals, we created summarized statistics on allele counts from anonymized data and made them available through the European Genome-phenome Archive (EGA, https://identifiers.org/ega. DATASET: EGAD00001009756 ) and in a dedicated browser, DanMAC5 (available at www.danmac5.dk ). The summary level data and the DanMAC5 browser provide insight into the allelic spectrum of sequence variants segregating in the Danish population, which is important in variant interpretation. DATA DESCRIPTION: Three WGS datasets with an average coverage of 30x were processed independently using the same quality control pipeline. Subsequently, we summarized, filtered, and merged allele counts to create a high-quality summary level dataset of sequence variants.

The population genomic legacy of the second plague pandemic.

Human populations have been shaped by catastrophes that may have left long-lasting signatures in their genomes. One notable example is the second plague pandemic that entered Europe in ca. 1,347 CE and repeatedly returned for over 300 years, with typical village and town mortality estimated at 10%-40%.1 It is assumed that this high mortality affected the gene pools of these populations. First, local population crashes reduced genetic diversity. Second, a change in frequency is expected for sequence variants that may have affected survival or susceptibility to the etiologic agent (Yersinia pestis).2 Third, mass mortality might alter the local gene pools through its impact on subsequent migration patterns. We explored these factors using the Norwegian city of Trondheim as a model, by sequencing 54 genomes spanning three time periods: (1) prior to the plague striking Trondheim in 1,349 CE, (2) the 17th-19th century, and (3) the present. We find that the pandemic period shaped the gene pool by reducing long distance immigration, in particular from the British Isles, and inducing a bottleneck that reduced genetic diversity. Although we also observe an excess of large FST values at multiple loci in the genome, these are shaped by reference biases introduced by mapping our relatively low genome coverage degraded DNA to the reference genome. This implies that attempts to detect selection using ancient DNA (aDNA) datasets that vary by read length and depth of sequencing coverage may be particularly challenging until methods have been developed to account for the impact of differential reference bias on test statistics.

Population-level deficit of homozygosity unveils CPSF3 as an intellectual disability syndrome gene.

Predicting the pathogenicity of biallelic missense variants can be challenging. Here, we use a deficit of observed homozygous carriers of missense variants, versus an expected number in a set of 153,054 chip-genotyped Icelanders, to identify potentially pathogenic genotypes. We follow three missense variants with a complete deficit of homozygosity and find that their pathogenic effect in homozygous state ranges from severe childhood disease to early embryonic lethality. One of these variants is in CPSF3, a gene not previously linked to disease. From a set of clinically sequenced Icelanders, and by sequencing archival samples targeted through the Icelandic genealogy, we find four homozygous carriers. Additionally, we find two homozygous carriers of Mexican descent of another missense variant in CPSF3. All six homozygous carriers of missense variants in CPSF3 show severe intellectual disability, seizures, microcephaly, and abnormal muscle tone. Here, we show how the absence of certain homozygous genotypes from a large population set can elucidate causes of previously unexplained recessive diseases and early miscarriage.

Large-scale integration of the plasma proteome with genetics and disease.

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) < 1%). We tested plasma protein levels for association with 373 diseases and other traits and identified 257,490 associations. We integrated pQTL and genetic associations with diseases and other traits and found that 12% of 45,334 lead associations in the GWAS Catalog are with variants in high linkage disequilibrium with pQTL. We identified 938 genes encoding potential drug targets with variants that influence levels of possible biomarkers. Combining proteomics, genomics and transcriptomics, we provide a valuable resource that can be used to improve understanding of disease pathogenesis and to assist with drug discovery and development.

Differences between germline genomes of monozygotic twins.

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. CpG>TpG 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.

Lifelong Reduction in LDL (Low-Density Lipoprotein) Cholesterol due to a Gain-of-Function Mutation in LDLR.

BACKGROUND: Loss-of-function mutations in the LDL (low-density lipoprotein) receptor gene (LDLR) cause elevated levels of LDL cholesterol and premature cardiovascular disease. To date, a gain-of-function mutation in LDLR with a large effect on LDL cholesterol levels has not been described. Here, we searched for sequence variants in LDLR that have a large effect on LDL cholesterol levels. METHODS: We analyzed whole-genome sequencing data from 43 202 Icelanders. Single-nucleotide polymorphisms and structural variants including deletions, insertions, and duplications were genotyped using whole-genome sequencing-based data. LDL cholesterol associations were carried out in a sample of >100 000 Icelanders with genetic information (imputed or whole-genome sequencing). Molecular analyses were performed using RNA sequencing and protein expression assays in Epstein-Barr virus-transformed lymphocytes. RESULTS: We discovered a 2.5-kb deletion (del2.5) overlapping the 3' untranslated region of LDLR in 7 heterozygous carriers from a single family. Mean level of LDL cholesterol was 74% lower in del2.5 carriers than in 101 851 noncarriers, a difference of 2.48 mmol/L (96 mg/dL; P=8.4×10-8). Del2.5 results in production of an alternative mRNA isoform with a truncated 3' untranslated region. The truncation leads to a loss of target sites for microRNAs known to repress translation of LDLR. In Epstein-Barr virus-transformed lymphocytes derived from del2.5 carriers, expression of alternative mRNA isoform was 1.84-fold higher than the wild-type isoform (P=0.0013), and there was 1.79-fold higher surface expression of the LDL receptor than in noncarriers (P=0.0086). We did not find a highly penetrant detrimental impact of lifelong very low levels of LDL cholesterol due to del2.5 on health of the carriers. CONCLUSIONS: Del2.5 is the first reported gain-of-function mutation in LDLR causing a large reduction in LDL cholesterol. These data point to a role for alternative polyadenylation of LDLR mRNA as a potent regulator of LDL receptor expression in humans.

A PRPH splice-donor variant associates with reduced sural nerve amplitude and risk of peripheral neuropathy.

Nerve conduction (NC) studies generate measures of peripheral nerve function that can reveal underlying pathology due to axonal loss, demyelination or both. We perform a genome-wide association study of sural NC amplitude and velocity in 7045 Icelanders and find a low-frequency splice-donor variant in PRPH (c.996+1G>A; MAF = 1.32%) associating with decreased NC amplitude but not velocity. PRPH encodes peripherin, an intermediate filament (IF) protein involved in cytoskeletal development and maintenance of neurons. Through RNA and protein studies, we show that the variant leads to loss-of-function (LoF), as when over-expressed in a cell line devoid of other IFs, it does not allow formation of the normal filamentous structure of peripherin, yielding instead punctate protein inclusions. Recall of carriers for neurological assessment confirms that from an early age, homozygotes have significantly lower sural NC amplitude than non-carriers and are at risk of a mild, early-onset, sensory-negative, axonal polyneuropathy.

A loss-of-function variant in ALOX15 protects against nasal polyps and chronic rhinosinusitis.

Nasal polyps (NP) are lesions on the nasal and paranasal sinus mucosa and are a risk factor for chronic rhinosinusitis (CRS). We performed genome-wide association studies on NP and CRS in Iceland and the UK (using UK Biobank data) with 4,366 NP cases, 5,608 CRS cases, and >700,000 controls. We found 10 markers associated with NP and 2 with CRS. We also tested 210 markers reported to associate with eosinophil count, yielding 17 additional NP associations. Of the 27 NP signals, 7 associate with CRS and 13 with asthma. Most notably, a missense variant in ALOX15 that causes a p.Thr560Met alteration in arachidonate 15-lipoxygenase (15-LO) confers large genome-wide significant protection against NP (P = 8.0 × 10-27, odds ratio = 0.32; 95% confidence interval = 0.26, 0.39) and CRS (P = 1.1 × 10-8, odds ratio = 0.64; 95% confidence interval = 0.55, 0.75). p.Thr560Met, carried by around 1 in 20 Europeans, was previously shown to cause near total loss of 15-LO enzymatic activity. Our findings identify 15-LO as a potential target for therapeutic intervention in NP and CRS.

Multiple transmissions of de novo mutations in families.

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.

A homozygous loss-of-function mutation leading to CYBC1 deficiency causes chronic granulomatous disease.

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 chip-genotyped 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.

Transcriptional profiling and biomarker identification reveal tissue specific effects of expanded ataxin-3 in a spinocerebellar ataxia type 3 mouse model.

BACKGROUND: Spinocerebellar ataxia type 3 (SCA3) is a progressive neurodegenerative disorder caused by expansion of the polyglutamine repeat in the ataxin-3 protein. Expression of mutant ataxin-3 is known to result in transcriptional dysregulation, which can contribute to the cellular toxicity and neurodegeneration. Since the exact causative mechanisms underlying this process have not been fully elucidated, gene expression analyses in brains of transgenic SCA3 mouse models may provide useful insights. METHODS: Here we characterised the MJD84.2 SCA3 mouse model expressing the mutant human ataxin-3 gene using a multi-omics approach on brain and blood. Gene expression changes in brainstem, cerebellum, striatum and cortex were used to study pathological changes in brain, while blood gene expression and metabolites/lipids levels were examined as potential biomarkers for disease. RESULTS: Despite normal motor performance at 17.5 months of age, transcriptional changes in brain tissue of the SCA3 mice were observed. Most transcriptional changes occurred in brainstem and striatum, whilst cerebellum and cortex were only modestly affected. The most significantly altered genes in SCA3 mouse brain were Tmc3, Zfp488, Car2, and Chdh. Based on the transcriptional changes, α-adrenergic and CREB pathways were most consistently altered for combined analysis of the four brain regions. When examining individual brain regions, axon guidance and synaptic transmission pathways were most strongly altered in striatum, whilst brainstem presented with strongest alterations in the pi-3 k cascade and cholesterol biosynthesis pathways. Similar to other neurodegenerative diseases, reduced levels of tryptophan and increased levels of ceramides, di- and triglycerides were observed in SCA3 mouse blood. CONCLUSIONS: The observed transcriptional changes in SCA3 mouse brain reveal parallels with previous reported neuropathology in patients, but also shows brain region specific effects as well as involvement of adrenergic signalling and CREB pathway changes in SCA3. Importantly, the transcriptional changes occur prior to onset of motor- and coordination deficits.

Ancient genomes from Iceland reveal the making of a human population.

Opportunities to directly study the founding of a human population and its subsequent evolutionary history are rare. Using genome sequence data from 27 ancient Icelanders, we demonstrate that they are a combination of Norse, Gaelic, and admixed individuals. We further show that these ancient Icelanders are markedly more similar to their source populations in Scandinavia and the British-Irish Isles than to contemporary Icelanders, who have been shaped by 1100 years of extensive genetic drift. Finally, we report evidence of unequal contributions from the ancient founders to the contemporary Icelandic gene pool. These results provide detailed insights into the making of a human population that has proven extraordinarily useful for the discovery of genotype-phenotype associations.

Compound heterozygous mutations in UBA5 causing early-onset epileptic encephalopathy in two sisters.

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.

Whole genome characterization of sequence diversity of 15,220 Icelanders.

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.