Long-read sequencing technologies substantially overcome the limitations of short-reads but have not been considered as a feasible replacement for population-scale projects, being a combination of too expensive, not scalable enough or too error-prone. Here we develop an efficient and scalable wet lab and computational protocol, Napu, for Oxford Nanopore Technologies long-read sequencing that seeks to address those limitations. We applied our protocol to cell lines and brain tissue samples as part of a pilot project for the National Institutes of Health Center for Alzheimer's and Related Dementias. Using a single PromethION flow cell, we can detect single nucleotide polymorphisms with F1-score comparable to Illumina short-read sequencing. Small indel calling remains difficult within homopolymers and tandem repeats, but achieves good concordance to Illumina indel calls elsewhere. Further, we can discover structural variants with F1-score on par with state-of-the-art de novo assembly methods. Our protocol phases small and structural variants at megabase scales and produces highly accurate, haplotype-specific methylation calls.

Recently, an African ancestry-specific Parkinson disease (PD) risk signal was identified at the gene encoding glucocerebrosidase (GBA1). This variant ( rs3115534 -G) is carried by ~50% of West African PD cases and imparts a dose-dependent increase in risk for disease. The risk variant has varied frequencies across African ancestry groups but is almost absent in European and Asian ancestry populations. GBA1 is a gene of high clinical and therapeutic interest. Damaging biallelic protein-coding variants cause Gaucher disease and monoallelic variants confer risk for PD and dementia with Lewy bodies, likely by reducing the function of glucocerebrosidase. Interestingly, the African ancestry-specific GBA1 risk variant is a noncoding variant, suggesting a different mechanism of action. Using full-length RNA transcript sequencing, we identified partial intron 8 expression in risk variant carriers (G) but not in nonvariant carriers (T). Antibodies targeting the N terminus of glucocerebrosidase showed that this intron-retained isoform is likely not protein coding and subsequent proteomics did not identify a shorter protein isoform, suggesting that the disease mechanism is RNA based. Clustered regularly interspaced short palindromic repeats editing of the reported index variant ( rs3115534 ) revealed that this is the sequence alteration responsible for driving the production of these transcripts containing intron 8. Follow-up analysis of this variant showed that it is in a key intronic branchpoint sequence and, therefore, has important implications in splicing and disease. In addition, when measuring glucocerebrosidase activity, we identified a dose-dependent reduction in risk variant carriers. Overall, we report the functional effect of a GBA1 noncoding risk variant, which acts by interfering with the splicing of functional GBA1 transcripts, resulting in reduced protein levels and reduced glucocerebrosidase activity. This understanding reveals a potential therapeutic target in an underserved and underrepresented population.