Scywalker: scalable end-to-end data analysis workflow for long-read single-cell transcriptome sequencing.

MOTIVATION: Existing nanopore single-cell data analysis tools showed severe limitations in handling current data sizes. RESULTS: We introduce scywalker, an innovative and scalable package developed to comprehensively analyze long-read sequencing data of full-length single-cell or single-nuclei cDNA. We developed novel scalable methods for cell barcode demultiplexing and single-cell isoform calling and quantification and incorporated these in an easily deployable package. Scywalker streamlines the entire analysis process, from sequenced fragments in FASTQ format to demultiplexed pseudobulk isoform counts, into a single command suitable for execution on either server or cluster. Scywalker includes data quality control, cell type identification, and an interactive report. Assessment of datasets from the human brain, Arabidopsis leaves, and previously benchmarked data from mixed cell lines demonstrate excellent correlation with short-read analyses at both the cell-barcoding and gene quantification levels. At the isoform level, we show that scywalker facilitates the direct identification of cell-type-specific expression of novel isoforms. AVAILABILITY AND IMPLEMENTATION: Scywalker is available on github.com/derijkp/scywalker under the GNU General Public License (GPL) and at https://zenodo.org/records/13359438/files/scywalker-0.108.0-Linux-x86_64.tar.gz.

Suppression of the Arabidopsis cinnamoyl-CoA reductase 1-6 intronic T-DNA mutation by epigenetic modification.

Arabidopsis (Arabidopsis thaliana) transfer DNA (T-DNA) insertion collections are popular resources for fundamental plant research. Cinnamoyl-CoA reductase 1 (CCR1) catalyzes an essential step in the biosynthesis of the cell wall polymer lignin. Accordingly, the intronic T-DNA insertion mutant ccr1-6 has reduced lignin levels and shows a stunted growth phenotype. Here, we report restoration of the ccr1-6 mutant phenotype and CCR1 expression levels after a genetic cross with a UDP-glucosyltransferase 72e1 (ugt72e1),-e2,-e3 T-DNA mutant. We discovered that the phenotypic recovery was not dependent on the UGT72E family loss of function but due to an epigenetic phenomenon called trans T-DNA suppression. Via trans T-DNA suppression, the gene function of an intronic T-DNA mutant was restored after the introduction of an additional T-DNA sharing identical sequences, leading to heterochromatinization and splicing out of the T-DNA-containing intron. Consequently, the suppressed ccr1-6 allele was named epiccr1-6. Long-read sequencing revealed that epiccr1-6, not ccr1-6, carries dense cytosine methylation over the full length of the T-DNA. We showed that the SAIL T-DNA in the UGT72E3 locus could trigger the trans T-DNA suppression of the GABI-Kat T-DNA in the CCR1 locus. Furthermore, we scanned the literature for other potential cases of trans T-DNA suppression in Arabidopsis and found that 22% of the publications matching our query report on double or higher-order T-DNA mutants that meet the minimal requirements for trans T-DNA suppression. These combined observations indicate that intronic T-DNA mutants need to be used with caution since methylation of intronic T-DNA might derepress gene expression and can thereby confound results.

Increased prime edit rates in KCNQ2 and SCN1A via single nicking all-in-one plasmids.

BACKGROUND: Prime editing (PE) is the most recent gene editing technology able to introduce targeted alterations to the genome, including single base pair changes, small insertions, and deletions. Several improvements to the PE machinery have been made in the past few years, and these have been tested in a range of model systems including immortalized cell lines, stem cells, and animal models. While double nicking RNA (dncRNA) PE systems PE3 and PE5 currently show the highest editing rates, they come with reduced accuracy as undesired indels or SNVs arise at edited loci. Here, we aimed to improve single ncRNA (sncRNA) systems PE2 and PE4max by generating novel all-in-one (pAIO) plasmids driven by an EF-1α promoter, which is especially suitable for human-induced pluripotent stem cell (hiPSC) models. RESULTS: pAIO-EF1α-PE2 and pAIO-EF1α-PE4max were used to edit the voltage gated potassium channel gene KCNQ2 and voltage gated sodium channel gene SCN1A. Two clinically relevant mutations were corrected using pAIO-EF1α-PE2 including the homozygous truncating SCN1A R612* variant in HEK293T cells and the heterozygous gain-of-function KCNQ2 R201C variant in patient-derived hiPSC. We show that sncRNA PE yielded detectable editing rates in hiPSC ranging between 6.4% and 9.8%, which was further increased to 41% after a GFP-based fluorescence-activated cell sorting (FACS) cell sorting step. Furthermore, we show that selecting the high GFP expressing population improved editing efficiencies up to 3.2-fold compared to the low GFP expressing population, demonstrating that not only delivery but also the number of copies of the PE enzyme and/or pegRNA per cell are important for efficient editing. Edit rates were not improved when an additional silent protospacer-adjacent motif (PAM)-removing alteration was introduced in hiPSC at the target locus. Finally, there were no genome-wide off-target effects using pAIO-EF1α-PE2 and no off-target editing activity near the edit locus highlighting the accuracy of snc prime editors. CONCLUSION: Taken together, our study shows an improved efficacy of EF-1α driven sncRNA pAIO-PE plasmids in hiPSC reaching high editing rates, especially after FACS sorting. Optimizing these sncRNA PE systems is of high value when considering future therapeutic in vivo use, where accuracy will be extremely important.

Rare variants in IFFO1, DTNB, NLRC3 and SLC22A10 associate with Alzheimer's disease CSF profile of neuronal injury and inflammation.

Alzheimer's disease (AD) biomarkers represent several neurodegenerative processes, such as synaptic dysfunction, neuronal inflammation and injury, as well as amyloid pathology. We performed an exome-wide rare variant analysis of six AD biomarkers (ß-amyloid, total/phosphorylated tau, NfL, YKL-40, and Neurogranin) to discover genes associated with these markers. Genetic and biomarker information was available for 480 participants from two studies: EMIF-AD and ADNI. We applied a principal component (PC) analysis to derive biomarkers combinations, which represent statistically independent biological processes. We then tested whether rare variants in 9576 protein-coding genes associate with these PCs using a Meta-SKAT test. We also tested whether the PCs are intermediary to gene effects on AD symptoms with a SMUT test. One PC loaded on NfL and YKL-40, indicators of neuronal injury and inflammation. Four genes were associated with this PC: IFFO1, DTNB, NLRC3, and SLC22A10. Mediation tests suggest, that these genes also affect dementia symptoms via inflammation/injury. We also observed an association between a PC loading on Neurogranin, a marker for synaptic functioning, with GABBR2 and CASZ1, but no mediation effects. The results suggest that rare variants in IFFO1, DTNB, NLRC3, and SLC22A10 heighten susceptibility to neuronal injury and inflammation, potentially by altering cytoskeleton structure and immune activity disinhibition, resulting in an elevated dementia risk. GABBR2 and CASZ1 were associated with synaptic functioning, but mediation analyses suggest that the effect of these two genes on synaptic functioning is not consequential for AD development.

Whole-exome rare-variant analysis of Alzheimer's disease and related biomarker traits.

INTRODUCTION: Despite increasing evidence of a role of rare genetic variation in the risk of Alzheimer's disease (AD), limited attention has been paid to its contribution to AD-related biomarker traits indicative of AD-relevant pathophysiological processes. METHODS: We performed whole-exome gene-based rare-variant association studies (RVASs) of 17 AD-related traits on whole-exome sequencing (WES) data generated in the European Medical Information Framework for Alzheimer's Disease Multimodal Biomarker Discovery (EMIF-AD MBD) study (n = 450) and whole-genome sequencing (WGS) data from ADNI (n = 808). RESULTS: Mutation screening revealed a novel probably pathogenic mutation (PSEN1 p.Leu232Phe). Gene-based RVAS revealed the exome-wide significant contribution of rare coding variation in RBKS and OR7A10 to cognitive performance and protection against left hippocampal atrophy, respectively. DISCUSSION: The identification of these novel gene-trait associations offers new perspectives into the role of rare coding variation in the distinct pathophysiological processes culminating in AD, which may lead to identification of novel therapeutic and diagnostic targets.

Mitochondrial GpC and CpG DNA Hypermethylation Cause Metabolic Stress-Induced Mitophagy and Cholestophagy.

Metabolic dysfunction-associated steatotic liver disease (MASLD) is characterized by a constant accumulation of lipids in the liver. This hepatic lipotoxicity is associated with a dysregulation of the first step in lipid catabolism, known as beta oxidation, which occurs in the mitochondrial matrix. Eventually, this dysregulation will lead to mitochondrial dysfunction. To evaluate the possible involvement of mitochondrial DNA methylation in this lipid metabolic dysfunction, we investigated the functional metabolic effects of mitochondrial overexpression of CpG (MSssI) and GpC (MCviPI) DNA methyltransferases in relation to gene expression and (mito)epigenetic signatures. Overall, the results show that mitochondrial GpC and, to a lesser extent, CpG methylation increase bile acid metabolic gene expression, inducing the onset of cholestasis through mito-nuclear epigenetic reprogramming. Moreover, both increase the expression of metabolic nuclear receptors and thereby induce basal overactivation of mitochondrial respiration. The latter promotes mitochondrial swelling, favoring lipid accumulation and metabolic-stress-induced mitophagy and autophagy stress responses. In conclusion, both mitochondrial GpC and CpG methylation create a metabolically challenging environment that induces mitochondrial dysfunction, which may contribute to the progression of MASLD.

Genomic Analysis of a Strain Collection Containing Multidrug-, Extensively Drug-, Pandrug-, and Carbapenem-Resistant Modern Clinical Isolates of Acinetobacter baumannii.

In this study, we characterize a new collection that comprises multidrug-resistant (MDR), extensively drug-resistant (XDR), pandrug-resistant (PDR), and carbapenem-resistant modern clinical isolates of Acinetobacter baumannii collected from hospitals through national microbiological surveillance in Belgium. Bacterial isolates (n = 43) were subjected to whole-genome sequencing (WGS), combining Illumina (MiSeq) and Nanopore (MinION) technologies, from which high-quality genomes (chromosome and plasmids) were de novo assembled. Antimicrobial susceptibility testing was performed along with genome analyses, which identified intrinsic and acquired resistance determinants along with their genetic environments and vehicles. Furthermore, the bacterial isolates were compared to the most prevalent A. baumannii sequence type 2 (ST2) (Pasteur scheme) genomes available from the BIGSdb database. Of the 43 strains, 40 carried determinants of resistance to carbapenems; blaOXA-23 (n = 29) was the most abundant acquired antimicrobial resistance gene, with 39 isolates encoding at least two different types of OXA enzymes. According to the Pasteur scheme, the majority of the isolates were globally disseminated clones of ST2 (n = 25), while less frequent sequence types included ST636 (n = 6), ST1 (n = 4), ST85 and ST78 (n = 2 each), and ST604, ST215, ST158, and ST10 (n = 1 each). Using the Oxford typing scheme, we identified 22 STs, including two novel types (ST2454 and ST2455). While the majority (26/29) of blaOXA-23 genes were chromosomally carried, all blaOXA-72 genes were plasmid borne. Our results show the presence of high-risk clones of A. baumannii within Belgian health care facilities with frequent occurrences of genes encoding carbapenemases, highlighting the crucial need for constant surveillance.

Structural variants identified by Oxford Nanopore PromethION sequencing of the human genome.

We sequenced the genome of the Yoruban reference individual NA19240 on the long-read sequencing platform Oxford Nanopore PromethION for evaluation and benchmarking of recently published aligners and germline structural variant calling tools, as well as a comparison with the performance of structural variant calling from short-read sequencing data. The structural variant caller Sniffles after NGMLR or minimap2 alignment provides the most accurate results, but additional confidence or sensitivity can be obtained by a combination of multiple variant callers. Sensitive and fast results can be obtained by minimap2 for alignment and a combination of Sniffles and SVIM for variant identification. We describe a scalable workflow for identification, annotation, and characterization of tens of thousands of structural variants from long-read genome sequencing of an individual or population. By discussing the results of this well-characterized reference individual, we provide an approximation of what can be expected in future long-read sequencing studies aiming for structural variant identification.

Loss of DPP6 in neurodegenerative dementia: a genetic player in the dysfunction of neuronal excitability.

Emerging evidence suggested a converging mechanism in neurodegenerative brain diseases (NBD) involving early neuronal network dysfunctions and alterations in the homeostasis of neuronal firing as culprits of neurodegeneration. In this study, we used paired-end short-read and direct long-read whole genome sequencing to investigate an unresolved autosomal dominant dementia family significantly linked to 7q36. We identified and validated a chromosomal inversion of ca. 4 Mb, segregating on the disease haplotype and disrupting the coding sequence of dipeptidyl-peptidase 6 gene (DPP6). DPP6 resequencing identified significantly more rare variants-nonsense, frameshift, and missense-in early-onset Alzheimer's disease (EOAD, p value = 0.03, OR = 2.21 95% CI 1.05-4.82) and frontotemporal dementia (FTD, p = 0.006, OR = 2.59, 95% CI 1.28-5.49) patient cohorts. DPP6 is a type II transmembrane protein with a highly structured extracellular domain and is mainly expressed in brain, where it binds to the potassium channel Kv4.2 enhancing its expression, regulating its gating properties and controlling the dendritic excitability of hippocampal neurons. Using in vitro modeling, we showed that the missense variants found in patients destabilize DPP6 and reduce its membrane expression (p < 0.001 and p < 0.0001) leading to a loss of protein. Reduced DPP6 and/or Kv4.2 expression was also detected in brain tissue of missense variant carriers. Loss of DPP6 is known to cause neuronal hyperexcitability and behavioral alterations in Dpp6-KO mice. Taken together, the results of our genomic, genetic, expression and modeling analyses, provided direct evidence supporting the involvement of DPP6 loss in dementia. We propose that loss of function variants have a higher penetrance and disease impact, whereas the missense variants have a variable risk contribution to disease that can vary from high to low penetrance. Our findings of DPP6, as novel gene in dementia, strengthen the involvement of neuronal hyperexcitability and alteration in the homeostasis of neuronal firing as a disease mechanism to further investigate.

The European Reference Genome Atlas: piloting a decentralised approach to equitable biodiversity genomics.

A genomic database of all Earth's eukaryotic species could contribute to many scientific discoveries; however, only a tiny fraction of species have genomic information available. In 2018, scientists across the world united under the Earth BioGenome Project (EBP), aiming to produce a database of high-quality reference genomes containing all ~1.5 million recognized eukaryotic species. As the European node of the EBP, the European Reference Genome Atlas (ERGA) sought to implement a new decentralised, equitable and inclusive model for producing reference genomes. For this, ERGA launched a Pilot Project establishing the first distributed reference genome production infrastructure and testing it on 98 eukaryotic species from 33 European countries. Here we outline the infrastructure and explore its effectiveness for scaling high-quality reference genome production, whilst considering equity and inclusion. The outcomes and lessons learned provide a solid foundation for ERGA while offering key learnings to other transnational, national genomic resource projects and the EBP.

Phenotypic Characterization and Heterogeneity among Modern Clinical Isolates of Acinetobacter baumannii.

Acinetobacter baumannii is an opportunistic pathogenic bacterium prioritized by WHO and CDC because of its increasing antibiotic resistance. Heterogeneity among strains represents the hallmark of A. baumannii bacteria. We wondered to what extent extensively used strains, so-called reference strains, reflect the dynamic nature and intrinsic heterogeneity of these bacteria. We analyzed multiple phenotypic traits of 43 nonredundant, modern, and multidrug-resistant, extensively drug-resistant, and pandrug-resistant clinical isolates and broadly used strains of A. baumannii. Comparison of these isolates at the genetic and phenotypic levels confirmed a high degree of heterogeneity. Importantly, we observed that a significant portion of modern clinical isolates strongly differs from several historically established strains in the light of colony morphology, cellular density, capsule production, natural transformability, and in vivo virulence. The significant differences between modern clinical isolates of A. baumannii and established strains could hamper the study of A. baumannii, especially concerning its virulence and resistance mechanisms. Hence, we propose a variable collection of modern clinical isolates that are characterized at the genetic and phenotypic levels, covering a wide range of the phenotypic spectrum, with six different macrocolony type groups, from avirulent to hypervirulent phenotypes, and with naturally noncapsulated to hypermucoid strains, with intermediate phenotypes as well. Strain-specific mechanistic observations remain interesting per se, and established "reference" strains have undoubtedly been shown to be very useful to study basic mechanisms of A. baumannii biology. However, any study based on a specific strain of A. baumannii should be compared to modern and clinically relevant isolates. IMPORTANCE Acinetobacter baumannii is a bacterium prioritized by the CDC and WHO because of its increasing antibiotic resistance, leading to treatment failures. The hallmark of this pathogen is the high heterogeneity observed among isolates, due to a very dynamic genome. In this context, we tested if a subset of broadly used isolates, considered "reference" strains, was reflecting the genetic and phenotypic diversity found among currently circulating clinical isolates. We observed that the so-called reference strains do not cover the whole diversity of the modern clinical isolates. While formerly established strains successfully generated a strong base of knowledge in the A. baumannii field and beyond, our study shows that a rational choice of strain, related to a specific biological question, should be taken into consideration. Any data obtained with historically established strains should also be compared to modern and clinically relevant isolates, especially concerning drug screening, resistance, and virulence contexts.

Mitochondrial DNA methylation in metabolic associated fatty liver disease.

Introduction: Hepatic lipid accumulation and mitochondrial dysfunction are hallmarks of metabolic associated fatty liver disease (MAFLD), yet molecular parameters underlying MAFLD progression are not well understood. Differential methylation within the mitochondrial DNA (mtDNA) has been suggested to be associated with dysfunctional mitochondria, also during progression to Metabolic Steatohepatitis (MeSH). This study further investigates whether mtDNA methylation is associated with hepatic lipid accumulation and MAFLD. Methods: HepG2 cells were constructed to stably express mitochondria-targeted viral and prokaryotic cytosine DNA methyltransferases (mtM.CviPI or mtM.SssI for GpC or CpG methylation, respectively). A catalytically inactive variant (mtM.CviPI-Mut) was constructed as a control. Mouse and human patients' samples were also investigated. mtDNA methylation was assessed by pyro- or nanopore sequencing. Results and discussion: Differentially induced mtDNA hypermethylation impaired mitochondrial gene expression and metabolic activity in HepG2-mtM.CviPI and HepG2-mtM.SssI cells and was associated with increased lipid accumulation, when compared to the controls. To test whether lipid accumulation causes mtDNA methylation, HepG2 cells were subjected to 1 or 2 weeks of fatty acid treatment, but no clear differences in mtDNA methylation were detected. In contrast, hepatic Nd6 mitochondrial gene body cytosine methylation and Nd6 gene expression were increased in mice fed a high-fat high cholesterol diet (HFC for 6 or 20 weeks), when compared to controls, while mtDNA content was unchanged. For patients with simple steatosis, a higher ND6 methylation was confirmed using Methylation Specific PCR, but no additional distinctive cytosines could be identified using pyrosequencing. This study warrants further investigation into a role for mtDNA methylation in promoting mitochondrial dysfunction and impaired lipid metabolism in MAFLD.

Null mutations in progranulin cause ubiquitin-positive frontotemporal dementia linked to chromosome 17q21.

Frontotemporal dementia (FTD) with ubiquitin-immunoreactive neuronal inclusions (both cytoplasmic and nuclear) of unknown nature has been linked to a chromosome 17q21 region (FTDU-17) containing MAPT (microtubule-associated protein tau). FTDU-17 patients have consistently been shown to lack a tau-immunoreactive pathology, a feature characteristic of FTD with parkinsonism linked to mutations in MAPT (FTDP-17). Furthermore, in FTDU-17 patients, mutations in MAPT and genomic rearrangements in the MAPT region have been excluded by both genomic sequencing and fluorescence in situ hybridization on mechanically stretched chromosomes. Here we demonstrate that FTDU-17 is caused by mutations in the gene coding for progranulin (PGRN), a growth factor involved in multiple physiological and pathological processes including tumorigenesis. Besides the production of truncated PGRN proteins due to premature stop codons, we identified a mutation within the splice donor site of intron 0 (IVS0 + 5G > C), indicating loss of the mutant transcript by nuclear degradation. The finding was made within an extensively documented Belgian FTDU-17 founder family. Transcript and protein analyses confirmed the absence of the mutant allele and a reduction in the expression of PGRN. We also identified a mutation (c.3G > A) in the Met1 translation initiation codon, indicating loss of PGRN due to lack of translation of the mutant allele. Our data provide evidence that PGRN haploinsufficiency leads to neurodegeneration because of reduced PGRN-mediated neuronal survival. Furthermore, in a Belgian series of familial FTD patients, PGRN mutations were 3.5 times more frequent than mutations in MAPT, underscoring a principal involvement of PGRN in FTD pathogenesis.

Scarless excision of an insertion sequence restores capsule production and virulence in Acinetobacter baumannii.

We identify a new mechanism mediating capsule production and virulence in the WHO and CDC priority ESKAPE pathogen Acinetobacter baumannii. Non-capsulated and avirulent bacteria can revert into a capsulated and virulent state upon scarless excision of an ISAba13 insertion sequence under stress conditions. Reversion events fully restore capsule production and in vivo virulence. This increases our knowledge about A. baumannii genome dynamics, and the regulation of capsule production, virulence and resistance.

Long-Read Sequencing to Unravel Complex Structural Variants of CEP78 Leading to Cone-Rod Dystrophy and Hearing Loss.

Inactivating variants as well as a missense variant in the centrosomal CEP78 gene have been identified in autosomal recessive cone-rod dystrophy with hearing loss (CRDHL), a rare syndromic inherited retinal disease distinct from Usher syndrome. Apart from this, a complex structural variant (SV) implicating CEP78 has been reported in CRDHL. Here we aimed to expand the genetic architecture of typical CRDHL by the identification of complex SVs of the CEP78 region and characterization of their underlying mechanisms. Approaches used for the identification of the SVs are shallow whole-genome sequencing (sWGS) combined with quantitative polymerase chain reaction (PCR) and long-range PCR, or ExomeDepth analysis on whole-exome sequencing (WES) data. Targeted or whole-genome nanopore long-read sequencing (LRS) was used to delineate breakpoint junctions at the nucleotide level. For all SVs cases, the effect of the SVs on CEP78 expression was assessed using quantitative PCR on patient-derived RNA. Apart from two novel canonical CEP78 splice variants and a frameshifting single-nucleotide variant (SNV), two SVs affecting CEP78 were identified in three unrelated individuals with CRDHL: a heterozygous total gene deletion of 235 kb and a partial gene deletion of 15 kb in a heterozygous and homozygous state, respectively. Assessment of the molecular consequences of the SVs on patient's materials displayed a loss-of-function effect. Delineation and characterization of the 15-kb deletion using targeted LRS revealed the previously described complex CEP78 SV, suggestive of a recurrent genomic rearrangement. A founder haplotype was demonstrated for the latter SV in cases of Belgian and British origin, respectively. The novel 235-kb deletion was delineated using whole-genome LRS. Breakpoint analysis showed microhomology and pointed to a replication-based underlying mechanism. Moreover, data mining of bulk and single-cell human and mouse transcriptional datasets, together with CEP78 immunostaining on human retina, linked the CEP78 expression domain with its phenotypic manifestations. Overall, this study supports that the CEP78 locus is prone to distinct SVs and that SV analysis should be considered in a genetic workup of CRDHL. Finally, it demonstrated the power of sWGS and both targeted and whole-genome LRS in identifying and characterizing complex SVs in patients with ocular diseases.

CHMP2B C-truncating mutations in frontotemporal lobar degeneration are associated with an aberrant endosomal phenotype in vitro.

The charged multivesicular body protein 2B gene (CHMP2B) was recently associated with frontotemporal lobar degeneration (FTLD) linked to chromosome 3 in a Danish FTLD family (FTD-3). In this family, a mutation in the acceptor splice site of exon 6 produced two aberrant transcripts predicting two C-truncated CHMP2B proteins due to a read through of intron 5 (p.Met178ValfsX2) and a cryptic splicing event within exon 6 (p.Met178LeufsX30). Extensive mutation analysis of CHMP2B in Belgian patients (N = 146) identified one nonsense mutation in exon 5 (c.493C>T) in a familial FTLD patient, predicting a C-truncated protein p.Gln165X analogous to the Danish mutant proteins. Overexpression of Belgian p.Gln165X in human neuroblastoma SK-N-SH cells showed the formation of large, aberrant endosomal structures that were highly similar to those observed for Danish p.Met178ValfsX2. Together, these data suggest that C-truncating mutations in CHMP2B might underlie the pathogenic mechanism in FTLD by disturbing endosome function. We also describe a missense mutation in exon 5 of CHMP2B (p.Asn143Ser) in a familial patient with cortical basal degeneration. However, the pathogenic character of this mutation remains elusive.

NanoSatellite: accurate characterization of expanded tandem repeat length and sequence through whole genome long-read sequencing on PromethION.

Technological limitations have hindered the large-scale genetic investigation of tandem repeats in disease. We show that long-read sequencing with a single Oxford Nanopore Technologies PromethION flow cell per individual achieves 30× human genome coverage and enables accurate assessment of tandem repeats including the 10,000-bp Alzheimer's disease-associated ABCA7 VNTR. The Guppy "flip-flop" base caller and tandem-genotypes tandem repeat caller are efficient for large-scale tandem repeat assessment, but base calling and alignment challenges persist. We present NanoSatellite, which analyzes tandem repeats directly on electric current data and improves calling of GC-rich tandem repeats, expanded alleles, and motif interruptions.

Founder mutation p.R1441C in the leucine-rich repeat kinase 2 gene in Belgian Parkinson's disease patients.

We determined the prevalence of mutations in two major functional domains of the leucine-rich repeat kinase 2 gene (LRRK2) in Belgian Parkinson's disease (PD) patients (N=304) of which 18.1% were familial PD patients. Ten patients were heterozygous for five different missense mutations (3.29%) of whom six carried the same mutation p.R1441C (1.97%). All six p.R1441C carriers were familial PD patients explaining 10.7% of familial PD in the Belgian patient group. Moreover, they shared a common disease haplotype of 21 consecutive markers in a region of 438 kb, suggesting that they are distant descendants of a single common ancestor. Clinically, p.R1441C carriers had typical levodopa-responsive parkinsonism with tremor as the most common presenting feature. Their age at onset was highly variable and ranged from 39 to 73 years, suggesting the influence of modifying factors. The remaining four patients were heterozygous each for a novel missense mutation located in the Roc or kinase domain. The pathogenic nature of these mutations remains to be determined, though we have genetic evidence that at least some represent rare but benign variants rather than causal mutations. The latter observation indicates that prudence is needed in diagnostic testing of LRRK2 in PD patients. Functional data should underlie a conclusion on the pathogenic nature of some mutations that have not been conclusively linked to disease.

Mutations other than null mutations producing a pathogenic loss of progranulin in frontotemporal dementia.

Null mutations in the progranulin gene (GRN, PGRN) were recently identified as the causal mechanism underlying frontotemporal dementia (FTD) with ubiquitin-positive brain pathology linked to chromosome 17 (FTDU-17). In a Belgian and French FTD series comprising 332 patients, we reported 13 PGRN null mutations which were mainly nonsense and frameshift mutations resulting in premature stop codons. Here we report in the same patient series three missense mutations of which two (c.743C>T, p.Pro248Leu and c.1294C>T, p.Arg432Cys) were predicted in silico to severely affect protein folding and/or processing leading to PGRN protein haploinsufficiency. In addition, we observed three sequence variations in the 5' regulatory region that might potentially affect PGRN transcription activity. Our findings extend the mutation spectrum in PGRN leading to loss of functional PGRN as the basis for FTD.

Progranulin null mutations in both sporadic and familial frontotemporal dementia.

Frontotemporal dementia (FTD) is the second most frequent type of neurodegenerative dementias. Mutations in the progranulin gene (GRN, PGRN) were recently identified in FTDU-17, an FTD subtype characterized by ubiquitin-immunoreactive inclusions and linkage to chromosome 17q21. We looked for PGRN mutations in a large series of 210 FTD patients (52 familial, 158 sporadic) to accurately evaluate the frequency of PGRN mutations in both sporadic and familial FTD, and FTD with associated motoneuron disease (FTD-MND), as well as to study the clinical phenotype of patients with a PGRN mutation. We identified nine novel PGRN null mutations in 10 index patients. The relative frequency of PGRN null mutations in FTD was 4.8% (10/210) and 12.8% (5/39) in pure familial forms. Interestingly, 5/158 (3.2%) apparently sporadic FTD patients carried a PGRN mutation, suggesting the possibility of de novo mutations or incomplete penetrance. In contrast, none of the 43 patients with FTD-MND had PGRN mutations, supporting that FTDU-17 and FTD-MND are genetically distinct. The clinical phenotype of PGRN mutation carriers was particular because of the wide range in onset age and the frequent occurrence of early apraxia (50%), visual hallucinations (30%), and parkinsonism (30%) during the course of the disease. This study supports that PGRN null mutations represent a more frequent cause of FTD than MAPT mutations (4.8% vs. 2.9%) but are not responsible for FTD-MND. It also demonstrates that half of the patients with a PGRN mutation in our series had no apparent family history of dementia. Taking this into account, genetic testing should be now considered more systematically, even in patients without obvious familial history of FTD.