Sequencing of human genomes with nanopore technology

Advances in long-read single-cell transcriptomics

Long-read single-cell transcriptomics (scRNA-Seq) is revolutionizing the way we profile heterogeneity in disease. Traditional short-read scRNA-Seq methods are limited in their ability to provide complete transcript coverage, resolve isoforms, and identify novel transcripts. The scRNA-Seq protocols developed for long-read sequencing platforms overcome these limitations by enabling the characterization of full-length transcripts. Long-read scRNA-Seq techniques initially suffered from comparatively poor accuracy compared to short read scRNA-Seq. However, with improvements in accuracy, accessibility, and cost efficiency, long-reads are gaining popularity in the field of scRNA-Seq. This review details the advances in long-read scRNA-Seq, with an emphasis on library preparation protocols and downstream bioinformatics analysis tools.

Molecular diversity in isocitrate dehydrogenase-wild-type glioblastoma

In the dynamic landscape of glioblastoma, the 2021 World Health Organization Classification of Central Nervous System tumours endeavoured to establish biological homogeneity, yet isocitrate dehydrogenase-wild-type (IDH-wt) glioblastoma persists as a tapestry of clinical and molecular diversity. Intertumoural heterogeneity in IDH-wt glioblastoma presents a formidable challenge in treatment strategies. Recent strides in genetics and molecular biology have enhanced diagnostic precision, revealing distinct subtypes and invasive patterns that influence survival in patients with IDH-wt glioblastoma. Genetic and molecular biomarkers, such as the overexpression of neurofibromin 1, phosphatase and tensin homolog and/or cyclin-dependent kinase inhibitor 2A, along with specific immune cell abundance and neurotransmitters, correlate with favourable outcomes. Conversely, increased expression of epidermal growth factor receptor tyrosine kinase, platelet-derived growth factor receptor alpha and/or vascular endothelial growth factor receptor, coupled with the prevalence of glioma stem cells, tumour-associated myeloid cells, regulatory T cells and exhausted effector cells, signifies an unfavourable prognosis. The methylation status of O6-methylguanine-DNA methyltransferase and the influence of microenvironmental factors and neurotransmitters further shape treatment responses. Understanding intertumoural heterogeneity is complemented by insights into intratumoural dynamics and cellular interactions within the tumour microenvironment. Glioma stem cells and immune cell composition significantly impact progression and outcomes, emphasizing the need for personalized therapies targeting pro-tumoural signalling pathways and resistance mechanisms. A successful glioblastoma management demands biomarker identification, combination therapies and a nuanced approach considering intratumoural variability. These advancements herald a transformative era in glioblastoma comprehension and treatment.

A universal molecular control for DNA, mRNA and protein expression

The expression of genes encompasses their transcription into mRNA followed by translation into protein. In recent years, next-generation sequencing and mass spectrometry methods have profiled DNA, RNA and protein abundance in cells. However, there are currently no reference standards that are compatible across these genomic, transcriptomic and proteomic methods, and provide an integrated measure of gene expression. Here, we use synthetic biology principles to engineer a multi-omics control, termed pREF, that can act as a universal molecular standard for next-generation sequencing and mass spectrometry methods. The pREF sequence encodes 21 synthetic genes that can be in vitro transcribed into spike-in mRNA controls, and in vitro translated to generate matched protein controls. The synthetic genes provide qualitative controls that can measure sensitivity and quantitative accuracy of DNA, RNA and peptide detection. We demonstrate the use of pREF in metagenome DNA sequencing and RNA sequencing experiments and evaluate the quantification of proteins using mass spectrometry. Unlike previous spike-in controls, pREF can be independently propagated and the synthetic mRNA and protein controls can be sustainably prepared by recipient laboratories using common molecular biology techniques. Together, this provides a universal synthetic standard able to integrate genomic, transcriptomic and proteomic methods.

HaplotagLR: An efficient and configurable utility for haplotagging long reads

Understanding the functional effects of sequence variation is crucial in genomics. Individual human genomes contain millions of variants that contribute to phenotypic variability and disease risks at the population level. Because variants rarely act in isolation, we must consider potential interactions of neighboring variants to accurately predict functional effects. We can accomplish this using haplotagging, which matches sequencing reads to their parental haplotypes using alleles observed at known heterozygous variants. However, few published tools for haplotagging exist and these share several technical and usability-related shortcomings that limit applicability, in particular a lack of insight or control over error rates, and lack of key metrics on the underlying sources of haplotagging error. Here we present HaplotagLR: a user-friendly tool that haplotags long sequencing reads based on a multinomial model and existing phased variant lists. HaplotagLR is user-configurable and includes a basic error model to control the empirical FDR in its output. We show that HaplotagLR outperforms the leading haplotagging method in simulated datasets, especially at high levels of specificity, and displays 7% greater sensitivity in haplotagging real data. HaplotagLR advances both the immediate utility of haplotagging and paves the way for further improvements to this important method.

Exact global alignment using A* with chaining seed heuristic and match pruning

MOTIVATION

Sequence alignment has been at the core of computational biology for half a century. Still, it is an open problem to design a practical algorithm for exact alignment of a pair of related sequences in linear-like time.

RESULTS

We solve exact global pairwise alignment with respect to edit distance by using the A* shortest path algorithm. In order to efficiently align long sequences with high divergence, we extend the recently proposed seed heuristic with match chaining, gap costs, and inexact matches. We additionally integrate the novel match pruning technique and diagonal transition to improve the A* search. We prove the correctness of our algorithm, implement it in the A*PA aligner, and justify our extensions intuitively and empirically.

On random sequences of divergence d=4% and length n, the empirical runtime of A*PA scales near-linearly with length (best fit n1.06, n≤107 bp). A similar scaling remains up to d=12% (best fit n1.24, n≤107 bp). For n=107 bp and d=4%, A*PA reaches >500× speedup compared to the leading exact aligners Edlib and BiWFA. The performance of A*PA is highly influenced by long gaps. On long (n>500kb) ONT reads of a human sample it efficiently aligns sequences with d<10%, leading to 3× median speedup compared to Edlib and BiWFA. When the sequences come from different human samples, A*PA performs 1.7× faster than Edlib and BiWFA.

AVAILABILITY AND IMPLEMENTATION

github.com/RagnarGrootKoerkamp/astar-pairwise-aligner.

Blastocystis genetic diversity in animal and human samples from different departments of Colombia using complete sequencing of the 18S rRNA gene (SSU rRNA) by Oxford Nanopore Technologies (ONT)

Blastocystis is an intestinal microeukaryote that has raised attention due to its wide distribution in animals and humans. The risk of zoonotic circulation primarily arises from close contact with infected animals. Therefore, the following study aimed to evaluate the diversity and frequency of Blastocystis subtypes in Colombian human and animal samples using complete sequencing of the 18S rRNA gene. For this purpose, 341 human stool samples and 277 animal fecal samples (from cattle, sheep, goat, pigs, cats, and dogs), were collected from different Colombian regions and analyzed using PCR-based detection and full-length 18S SSU rRNA gene Next-Generation Sequencing (NGS). Among the 618 samples from both hosts, humans and animals, the results revealed widespread Blastocystis frequency, with 48.09% (n = 164) in humans and 31.4% (n = 87) detection in animals. Dogs, cats, sheep, pigs, and wild animals tested positive, aligning with global prevalence patterns. Also, 29 human samples and 23 animal samples were sequenced using ONT technology from which 11 long-read unique sequences were generated and cluster with their compared reference sequences. The subtype distribution varied within hosts, detecting ST1 and ST3 in both human and animal samples. Subtypes ST5, ST10, ST14, ST15, ST21, ST24, ST25 and ST26 were limited to animals hosts, some of which are considered to have zoonotic potential. On the other hand, ST2 was found exclusively in human samples from Bolivar region. Mixed infections occurred in both animal and humans, 60.86% and 27.58% respectively. Moreover, to our knowledge, this is the first study in Colombia identifying ST15 in pigs and ST25 in sheep. The subtypes (STs) identified in this study indicate that certain animals may serve as reservoirs with the potential for zoonotic transmission. The identification of zoonotic subtypes highlights the use of Next Generation Sequencing as the depth and resolution of the sequences increases providing insights into STs of medical and veterinarian significance. It also reveals the coexistence of diverse subtypes among hosts. Further research is essential for understanding transmission dynamics, health implications, and detection strategies for Blastocystis occurrence in animals and humans, mainly associated to the role of animals as reservoirs and their close interaction with humans.

Molecular characterization of G6PD mutations identifies new mutations and a high frequency of intronic variants in Thai females

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is an X-linked enzymopathy caused by mutations in the G6PD gene. A medical concern associated with G6PD deficiency is acute hemolytic anemia induced by certain foods, drugs, and infections. Although phenotypic tests can correctly identify hemizygous males, as well as homozygous and compound heterozygous females, heterozygous females with a wide range of G6PD activity may be misclassified as normal. This study aimed to develop multiplex high-resolution melting (HRM) analyses to enable the accurate detection of G6PD mutations, especially among females with heterozygous deficiency. Multiplex HRM assays were developed to detect six G6PD variants, i.e., G6PD Gaohe (c.95A>G), G6PD Chinese-4 (c.392G>T), G6PD Mahidol (c.487G>A), G6PD Viangchan (c.871G>A), G6PD Chinese-5 (c.1024C>T), and G6PD Union (c.1360C>T) in two reactions. The assays were validated and then applied to genotype G6PD mutations in 248 Thai females. The sensitivity of the HRM assays developed was 100% [95% confidence interval (CI): 94.40%-100%] with a specificity of 100% (95% CI: 88.78%-100%) for detecting these six mutations. The prevalence of G6PD deficiency was estimated as 3.63% (9/248) for G6PD deficiency and 31.05% (77/248) for intermediate deficiency by phenotypic assay. The developed HRM assays identified three participants with normal enzyme activity as heterozygous for G6PD Viangchan. Interestingly, a deletion in intron 5 nucleotide position 637/638 (c.486-34delT) was also detected by the developed HRM assays. G6PD genotyping revealed a total of 12 G6PD genotypes, with a high prevalence of intronic variants. Our results suggested that HRM analysis-based genotyping is a simple and reliable approach for detecting G6PD mutations, and could be used to prevent the misdiagnosis of heterozygous females by phenotypic assay. This study also sheds light on the possibility of overlooking intronic variants, which could affect G6PD expression and contribute to enzyme deficiency.

Structural and non-coding variants increase the diagnostic yield of clinical whole genome sequencing for rare diseases

BACKGROUND

Whole genome sequencing is increasingly being used for the diagnosis of patients with rare diseases. However, the diagnostic yields of many studies, particularly those conducted in a healthcare setting, are often disappointingly low, at 25-30%. This is in part because although entire genomes are sequenced, analysis is often confined to in silico gene panels or coding regions of the genome.

METHODS

We undertook WGS on a cohort of 122 unrelated rare disease patients and their relatives (300 genomes) who had been pre-screened by gene panels or arrays. Patients were recruited from a broad spectrum of clinical specialties. We applied a bioinformatics pipeline that would allow comprehensive analysis of all variant types. We combined established bioinformatics tools for phenotypic and genomic analysis with our novel algorithms (SVRare, ALTSPLICE and GREEN-DB) to detect and annotate structural, splice site and non-coding variants.

RESULTS

Our diagnostic yield was 43/122 cases (35%), although 47/122 cases (39%) were considered solved when considering novel candidate genes with supporting functional data into account. Structural, splice site and deep intronic variants contributed to 20/47 (43%) of our solved cases. Five genes that are novel, or were novel at the time of discovery, were identified, whilst a further three genes are putative novel disease genes with evidence of causality. We identified variants of uncertain significance in a further fourteen candidate genes. The phenotypic spectrum associated with RMND1 was expanded to include polymicrogyria. Two patients with secondary findings in FBN1 and KCNQ1 were confirmed to have previously unidentified Marfan and long QT syndromes, respectively, and were referred for further clinical interventions. Clinical diagnoses were changed in six patients and treatment adjustments made for eight individuals, which for five patients was considered life-saving.

CONCLUSIONS

Genome sequencing is increasingly being considered as a first-line genetic test in routine clinical settings and can make a substantial contribution to rapidly identifying a causal aetiology for many patients, shortening their diagnostic odyssey. We have demonstrated that structural, splice site and intronic variants make a significant contribution to diagnostic yield and that comprehensive analysis of the entire genome is essential to maximise the value of clinical genome sequencing.

Toward Cytogenomics: Technical Assessment of Long-Read Nanopore Whole-Genome Sequencing for Detecting Large Chromosomal Alterations in Mantle Cell Lymphoma

The current advances and success of next-generation sequencing hold the potential for the transition of cancer cytogenetics toward comprehensive cytogenomics. However, the conventional use of short reads impedes the resolution of chromosomal aberrations. Thus, this study evaluated the detection and reproducibility of extensive copy number alterations and chromosomal translocations using long-read Oxford Nanopore Technologies whole-genome sequencing compared with short-read Illumina sequencing. Using the mantle cell lymphoma cell line Granta-519, almost 99% copy-number reproducibility at the 100-kilobase resolution between replicates was demonstrated, with 98% concordance to Illumina. Collectively, the performance of copy number calling from 1.5 million to 7.5 million long reads was comparable to 1 billion Illumina-based reads (50× coverage). Expectedly, the long-read resolution of canonical translocation t(11;14)(q13;q32) was superior, with a sequence similarity of 89% to the already published CCND1/IGH junction (9× coverage), spanning up to 69 kilobases. The cytogenetic profile of Granta-519 was in general agreement with the literature and karyotype, although several differences remained unresolved. In conclusion, contemporary long-read sequencing is primed for future cytogenomics or sequencing-guided cytogenetics. The combined strength of long- and short-read sequencing is apparent, where the high-precision junctional mapping complements and splits paired-end reads. The potential is emphasized by the flexible single-sample genomic data acquisition of Oxford Nanopore Technologies with the high resolution of allelic imbalances using Illumina short-read sequencing.

Block Aligner: an adaptive SIMD-accelerated aligner for sequences and position-specific scoring matrices

MOTIVATION

Efficiently aligning sequences is a fundamental problem in bioinformatics. Many recent algorithms for computing alignments through Smith-Waterman-Gotoh dynamic programming (DP) exploit Single Instruction Multiple Data (SIMD) operations on modern CPUs for speed. However, these advances have largely ignored difficulties associated with efficiently handling complex scoring matrices or large gaps (insertions or deletions).

RESULTS

We propose a new SIMD-accelerated algorithm called Block Aligner for aligning nucleotide and protein sequences against other sequences or position-specific scoring matrices. We introduce a new paradigm that uses blocks in the DP matrix that greedily shift, grow, and shrink. This approach allows regions of the DP matrix to be adaptively computed. Our algorithm reaches over 5-10 times faster than some previous methods while incurring an error rate of less than 3% on protein and long read datasets, despite large gaps and low sequence identities.

AVAILABILITY AND IMPLEMENTATION

Our algorithm is implemented for global, local, and X-drop alignments. It is available as a Rust library (with C bindings) at https://github.com/Daniel-Liu-c0deb0t/block-aligner.

Impact of genetic structural variants in factor XI deficiency: identification, accurate characterization, and inferred mechanism by long-read sequencing

BACKGROUND

Congenital factor XI (FXI) deficiency is a probably underestimated coagulopathy that confers antithrombotic protection. Characterization of genetic defects in F11 is mainly focused on the identification of single-nucleotide variants and small insertion/deletions because they represent up to 99% of the alterations accounting for factor deficiency, with only 3 gross gene defects of structural variants (SVs) having been described.

OBJECTIVES

To identify and characterize the SVs affecting F11.

METHODS

The study was performed in 93 unrelated subjects with FXI deficiency recruited in Spanish hospitals over a period of 25 years (1997-2022). F11 was analyzed by next-generation sequencing, multiplex ligand probe amplification, and long-read sequencing.

RESULTS

Our study identified 30 different genetic variants. Interestingly, we found 3 SVs, all heterozygous: a complex duplication affecting exons 8 and 9, a tandem duplication of exon 14, and a large deletion affecting the whole gene. Nucleotide resolution obtained by long-read sequencing revealed Alu repetitive elements involved in all breakpoints. The large deletion was probably generated de novo in the paternal allele during gametogenesis, and despite affecting 30 additional genes, no syndromic features were described.

CONCLUSION

SVs may account for a high proportion of F11 genetic defects implicated in the molecular pathology of congenital FXI deficiency. These SVs, likely caused by a nonallelic homologous recombination involving repetitive elements, are heterogeneous in both type and length and may be de novo. These data support the inclusion of methods to detect SVs in this disorder, with long-read-based methods being the most appropriate because they detect all SVs and achieve adequate nucleotide resolution.

Diagnostics and analysis of SARS-CoV-2: current status, recent advances, challenges and perspectives

The disastrous spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has induced severe public healthcare issues and weakened the global economy significantly. Although SARS-CoV-2 infection is not as fatal as the initial outbreak, many infected victims suffer from long COVID. Therefore, rapid and large-scale testing is critical in managing patients and alleviating its transmission. Herein, we review the recent advances in techniques to detect SARS-CoV-2. The sensing principles are detailed together with their application domains and analytical performances. In addition, the advantages and limits of each method are discussed and analyzed. Besides molecular diagnostics and antigen and antibody tests, we also review neutralizing antibodies and emerging SARS-CoV-2 variants. Further, the characteristics of the mutational locations in the different variants with epidemiological features are summarized. Finally, the challenges and possible strategies are prospected to develop new assays to meet different diagnostic needs. Thus, this comprehensive and systematic review of SARS-CoV-2 detection technologies may provide insightful guidance and direction for developing tools for the diagnosis and analysis of SARS-CoV-2 to support public healthcare and effective long-term pandemic management and control.

Dwell Time Prolongation and Identification of Single Nucleotides Passing through a Solid-State Nanopore by Using Ammonium Sulfate Aqueous Solution

The ionic current blockades when poly(dT)₆₀ or dNTPs passed through SiN nanopores in an aqueous solution containing (NH₄)₂SO₄ were investigated. The dwell time of poly(dT)₆₀ in the nanopores in an aqueous solution containing (NH₄)₂SO₄ was significantly longer compared to that in an aqueous solution that did not contain (NH₄)₂SO₄. This dwell time prolongation effect due to the aqueous solution containing (NH₄)₂SO₄ was also confirmed when dCTP passed through the nanopores. In addition, when the nanopores were fabricated via dielectric breakdown in the aqueous solution containing (NH₄)₂SO₄, the dwell time prolongation effect for dCTP still occurred even after the aqueous solution was displaced with the aqueous solution without (NH₄)₂SO₄. Furthermore, we measured the ionic current blockades when the four types of dNTPs passed through the same nanopore, and the four types of dNTPs could be statistically identified according to their current blockade values.

Accurate sequencing of DNA motifs able to form alternative (non-B) structures

Approximately 13% of the human genome at certain motifs have the potential to form noncanonical (non-B) DNA structures (e.g., G-quadruplexes, cruciforms, and Z-DNA), which regulate many cellular processes but also affect the activity of polymerases and helicases. Because sequencing technologies use these enzymes, they might possess increased errors at non-B structures. To evaluate this, we analyzed error rates, read depth, and base quality of Illumina, Pacific Biosciences (PacBio) HiFi, and Oxford Nanopore Technologies (ONT) sequencing at non-B motifs. All technologies showed altered sequencing success for most non-B motif types, although this could be owing to several factors, including structure formation, biased GC content, and the presence of homopolymers. Single-nucleotide mismatch errors had low biases in HiFi and ONT for all non-B motif types but were increased for G-quadruplexes and Z-DNA in all three technologies. Deletion errors were increased for all non-B types but Z-DNA in Illumina and HiFi, as well as only for G-quadruplexes in ONT. Insertion errors for non-B motifs were highly, moderately, and slightly elevated in Illumina, HiFi, and ONT, respectively. Additionally, we developed a probabilistic approach to determine the number of false positives at non-B motifs depending on sample size and variant frequency, and applied it to publicly available data sets (1000 Genomes, Simons Genome Diversity Project, and gnomAD). We conclude that elevated sequencing errors at non-B DNA motifs should be considered in low-read-depth studies (single-cell, ancient DNA, and pooled-sample population sequencing) and in scoring rare variants. Combining technologies should maximize sequencing accuracy in future studies of non-B DNA.

Critically short telomeres derepress retrotransposons to promote genome instability in embryonic stem cells

Telomeres, at the ends of chromosomes, protect chromosomes from fusion and preserve genomic stability. However, the molecular mechanisms underlying telomere attrition-induced genome instability remain to be understood. We systematically analyzed the expression of retrotransposons and performed genomic sequencing of different cell and tissue types with telomeres of varying lengths due to telomerase deficiency. We found that critically short telomeres altered retrotransposon activity to promote genomic instability in mouse embryonic stem cells, as evidenced by elevated numbers of single nucleotide variants, indels and copy number variations (CNVs). Transpositions of retrotransposons such as LINE1 resulting from the short telomeres can also be found in these genomes with elevated number of mutations and CNVs. Retrotransposon activation is linked to increased chromatin accessibility, and reduced heterochromatin abundance correlates with short telomeres. Re-elongation of telomeres upon recovery of telomerase partly represses retrotransposons and heterochromatin accumulation. Together, our findings suggest a potential mechanism by which telomeres maintain genomic stability by suppressing chromatin accessibility and retrotransposon activity.

Unraveling the genomic complexity of sylvatic mosquitoes in changing Neotropical environments

Sylvatic New World mosquitoes (e.g. Old-growth Forest species) can transmit viruses among non-human primates. This could be a continuous source of viral cycling and spillover events from animals to humans, particularly in changing environments. However, most species of Neotropical sylvatic mosquitoes (genera Aedes, Haemagogus, and Sabethes), which include vector and non-vector species, currently lack genomic resources because there is no reliable and accurate approach for creating de novo reference genomes for these insects. This is a major knowledge gap in the biology of these mosquitoes, restricting our ability to predict and mitigate the emergence and spread of novel arboviruses in Neotropical regions. We discuss recent advances and potential solutions for generating hybrid de novo assemblies from vector and non-vector species using pools of consanguineous offspring. We also discussed research opportunities likely to emerge from these genomic resources.

Revolving ATPase motors as asymmetrical hexamers in translocating lengthy dsDNA via conformational changes and electrostatic interactions in phi29, T7, herpesvirus, mimivirus, E. coli, and Streptomyces

Investigations of the parallel architectures of biomotors in both prokaryotic and eukaryotic systems suggest a similar revolving mechanism in the use of ATP to drive translocation of the lengthy double-stranded (ds)DNA genomes. This mechanism is exemplified by the dsDNA packaging motor of bacteriophage phi29 that operates through revolving but not rotating dsDNA to "Push through a one-way valve". This unique and novel revolving mechanism discovered in phi29 DNA packaging motor was recently reported in other systems including the dsDNA packaging motor of herpesvirus, the dsDNA ejecting motor of bacteriophage T7, the plasmid conjugation machine TraB in Streptomyces, the dsDNA translocase FtsK of gram-negative bacteria, and the genome-packaging motor in mimivirus. These motors exhibit an asymmetrical hexameric structure for transporting the genome via an inch-worm sequential action. This review intends to delineate the revolving mechanism from a perspective of conformational changes and electrostatic interactions. In phi29, the positively charged residues Arg-Lys-Arg in the N-terminus of the connector bind the negatively charged interlocking domain of pRNA. ATP binding to an ATPase subunit induces the closed conformation of the ATPase. The ATPase associates with an adjacent subunit to form a dimer facilitated by the positively charged arginine finger. The ATP-binding induces a positive charging on its DNA binding surface via an allostery mechanism and thus the higher affinity for the negatively charged dsDNA. ATP hydrolysis induces an expanded conformation of the ATPase with a lower affinity for dsDNA due to the change of the surface charge, but the (ADP+Pi)-bound subunit in the dimer undergoes a conformational change that repels dsDNA. The positively charged lysine rings of the connector attract dsDNA stepwise and periodically to keep its revolving motion along the channel wall, thus maintaining the one-way translocation of dsDNA without reversal and sliding out. The finding of the presence of the asymmetrical hexameric architectures of many ATPases that use the revolving mechanism may provide insights into the understanding of translocation of the gigantic genomes including chromosomes in complicated systems without coiling and tangling to speed up dsDNA translocation and save energy.

Nanopores: synergy from DNA sequencing to industrial filtration - small holes with big impact

Nanopores in thin membranes play important roles in science and industry. Single nanopores have provided a step-change in portable DNA sequencing and understanding nanoscale transport while multipore membranes facilitate food processing and purification of water and medicine. Despite the unifying use of nanopores, the fields of single nanopores and multipore membranes differ - to varying degrees - in terms of materials, fabrication, analysis, and applications. Such a partial disconnect hinders scientific progress as important challenges are best resolved together. This Viewpoint suggests how synergistic crosstalk between the two fields can provide considerable mutual benefits in fundamental understanding and the development of advanced membranes. We first describe the main differences including the atomistic definition of single pores compared to the less defined conduits in multipore membranes. We then outline steps to improve communication between the two fields such as harmonizing measurements and modelling of transport and selectivity. The resulting insight is expected to improve the rational design of porous membranes. The Viewpoint concludes with an outlook of other developments that can be best achieved by collaboration across the two fields to advance the understanding of transport in nanopores and create next-generation porous membranes tailored for sensing, filtration, and other applications.

Current Advances in Nanotechnology for the Next Generation of Sequencing (NGS)

This communication aims at discussing strategies based on developments from nanotechnology focused on the next generation of sequencing (NGS). In this regard, it should be noted that even in the advanced current situation of many techniques and methods accompanied with developments of technology, there are still existing challenges and needs focused on real samples and low concentrations of genomic materials. The approaches discussed/described adopt spectroscopical techniques and new optical setups. PCR bases are introduced to understand the role of non-covalent interactions by discussing about Nobel prizes related to genomic material detection. The review also discusses colorimetric methods, polymeric transducers, fluorescence detection methods, enhanced plasmonic techniques such as metal-enhanced fluorescence (MEF), semiconductors, and developments in metamaterials. In addition, nano-optics, challenges linked to signal transductions, and how the limitations reported in each technique could be overcome are considered in real samples. Accordingly, this study shows developments where optical active nanoplatforms generate signal detection and transduction with enhanced performances and, in many cases, enhanced signaling from single double-stranded deoxyribonucleic acid (DNA) interactions. Future perspectives on miniaturized instrumentation, chips, and devices aimed at detecting genomic material are analyzed. However, the main concept in this report derives from gained insights into nanochemistry and nano-optics. Such concepts could be incorporated into other higher-sized substrates and experimental and optical setups.

Optimal gap-affine alignment in O(s) space

MOTIVATION

Pairwise sequence alignment remains a fundamental problem in computational biology and bioinformatics. Recent advances in genomics and sequencing technologies demand faster and scalable algorithms that can cope with the ever-increasing sequence lengths. Classical pairwise alignment algorithms based on dynamic programming are strongly limited by quadratic requirements in time and memory. The recently proposed wavefront alignment algorithm (WFA) introduced an efficient algorithm to perform exact gap-affine alignment in O(ns) time, where s is the optimal score and n is the sequence length. Notwithstanding these bounds, WFA's O(s2) memory requirements become computationally impractical for genome-scale alignments, leading to a need for further improvement.

RESULTS

In this article, we present the bidirectional WFA algorithm, the first gap-affine algorithm capable of computing optimal alignments in O(s) memory while retaining WFA's time complexity of O(ns). As a result, this work improves the lowest known memory bound O(n) to compute gap-affine alignments. In practice, our implementation never requires more than a few hundred MBs aligning noisy Oxford Nanopore Technologies reads up to 1 Mbp long while maintaining competitive execution times.

AVAILABILITY AND IMPLEMENTATION

All code is publicly available at https://github.com/smarco/BiWFA-paper.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

LRphase: an efficient method for assigning haplotype identity to long reads

Understanding the functional effects of sequence variation is among the primary goals of contemporary genomics. Individual human genomes contain millions of variants which are thought to contribute to phenotypic variability and differential disease risks at the population level. However, because variants rarely act in isolation, we cannot accurately predict functional effects without first considering the potential effects of other interacting variants on the same chromosome. This information can be obtained by phasing the read data from sequencing experiments. However, no standalone tools are available to simply phase reads based on known haplotypes. Here we present LRphase: a user-friendly utility for simple phasing of long sequencing reads.

Advances in sequencing technologies for amyotrophic lateral sclerosis research

Amyotrophic lateral sclerosis (ALS) is caused by upper and lower motor neuron loss and has a fairly rapid disease progression, leading to fatality in an average of 2-5 years after symptom onset. Numerous genes have been implicated in this disease; however, many cases remain unexplained. Several technologies are being used to identify regions of interest and investigate candidate genes. Initial approaches to detect ALS genes include, among others, linkage analysis, Sanger sequencing, and genome-wide association studies. More recently, next-generation sequencing methods, such as whole-exome and whole-genome sequencing, have been introduced. While those methods have been particularly useful in discovering new ALS-linked genes, methodological advances are becoming increasingly important, especially given the complex genetics of ALS. Novel sequencing technologies, like long-read sequencing, are beginning to be used to uncover the contribution of repeat expansions and other types of structural variation, which may help explain missing heritability in ALS. In this review, we discuss how popular and/or upcoming methods are being used to discover ALS genes, highlighting emerging long-read sequencing platforms and their role in aiding our understanding of this challenging disease.

CNN6mA: Interpretable neural network model based on position-specific CNN and cross-interactive network for 6mA site prediction

N6-methyladenine (6mA) plays a critical role in various epigenetic processing including DNA replication, DNA repair, silencing, transcription, and diseases such as cancer. To understand such epigenetic mechanisms, 6 mA has been detected by high-throughput technologies on a genome-wide scale at single-base resolution, together with conventional methods such as immunoprecipitation, mass spectrometry and capillary electrophoresis, but these experimental approaches are time-consuming and laborious. To complement these problems, we have developed a CNN-based 6 mA site predictor, named CNN6mA, which proposed two new architectures: a position-specific 1-D convolutional layer and a cross-interactive network. In the position-specific 1-D convolutional layer, position-specific filters with different window sizes were applied to an inquiry sequence instead of sharing the same filters over all positions in order to extract the position-specific features at different levels. The cross-interactive network explored the relationships between all the nucleotide patterns within the inquiry sequence. Consequently, CNN6mA outperformed the existing state-of-the-art models in many species and created the contribution score vector that intelligibly interpret the prediction mechanism. The source codes and web application in CNN6mA are freely accessible at https://github.com/kuratahiroyuki/CNN6mA.git and http://kurata35.bio.kyutech.ac.jp/CNN6mA/, respectively.

PBSIM3: a simulator for all types of PacBio and ONT long reads

Long-read sequencers, such as Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT) sequencers, have improved their read length and accuracy, thereby opening up unprecedented research. Many tools and algorithms have been developed to analyze long reads, and rapid progress in PacBio and ONT has further accelerated their development. Together with the development of high-throughput sequencing technologies and their analysis tools, many read simulators have been developed and effectively utilized. PBSIM is one of the popular long-read simulators. In this study, we developed PBSIM3 with three new functions: error models for long reads, multi-pass sequencing for high-fidelity read simulation and transcriptome sequencing simulation. Therefore, PBSIM3 is now able to meet a wide range of long-read simulation requirements.

Czy plastik może rozpocząć nową erę w archiwizacji danych?

With the rapid development of information technology, many aspects of our lives are undergoing a digital transformation. An increasing number of users are going online every year, and constantly improving artificial intelligence is gaining popularity, which leads to the growing production of information. Nowadays, information is usually stored in data centres, which will be forced to increase their space with the constant flow of new bits of information. Together with the increase in their space, energy consumption and associated maintenance costs are escalating. In 2021, global data centre power consumption was 220–320 TWh, which is about 0.9–1.3% of global power consumption. Continuous power supply for database operations is responsible for about 1% of total carbon dioxide emissions. Furthermore, it has already been reported that with the exponentially growing amount of data, in about 20 years, the amount of silicon for microprocessors will no longer be sufficient to store all the information. Therefore, scientists are looking for alternatives to the currently used data storage solutions and are developing new technologies using chemical molecules. Recently, even plastic has been explored as a data carrier. In this work, we present examples of new technologies for data storage in polymers. We have discussed polymers as data carriers in comparison with currently used solutions and deliberated whether plastic can become a future material for information archiving.

Genotype Complements the Phenotype: Identification of the Pathogenicity of an LMNA Splice Variant by Nanopore Long-Read Sequencing in a Large DCM Family

Dilated cardiomyopathy (DCM) is a common cause of heart failure (HF) and is of familial origin in 20−40% of cases. Genetic testing by next-generation sequencing (NGS) has yielded a definite diagnosis in many cases; however, some remain elusive. In this study, we used a combination of NGS, human-induced pluripotent-stem-cell-derived cardiomyocytes (iPSC-CMs) and nanopore long-read sequencing to identify the causal variant in a multi-generational pedigree of DCM. A four-generation family with familial DCM was investigated. Next-generation sequencing (NGS) was performed on 22 family members. Skin biopsies from two affected family members were used to generate iPSCs, which were then differentiated into iPSC-CMs. Short-read RNA sequencing was used for the evaluation of the target gene expression, and long-read RNA nanopore sequencing was used to evaluate the relevance of the splice variants. The pedigree suggested a highly penetrant, autosomal dominant mode of inheritance. The phenotype of the family was suggestive of laminopathy, but previous genetic testing using both Sanger and panel sequencing only yielded conflicting evidence for LMNA p.R644C (rs142000963), which was not fully segregated. By re-sequencing four additional affected family members, further non-coding LMNA variants could be detected: rs149339264, rs199686967, rs201379016, and rs794728589. To explore the roles of these variants, iPSC-CMs were generated. RNA sequencing showed the LMNA expression levels to be significantly lower in the iPSC-CMs of the LMNA variant carriers. We demonstrated a dysregulated sarcomeric structure and altered calcium homeostasis in the iPSC-CMs of the LMNA variant carriers. Using targeted nanopore long-read sequencing, we revealed the biological significance of the variant c.356+1G>A, which generates a novel 5′ splice site in exon 1 of the cardiac isomer of LMNA, causing a nonsense mRNA product with almost complete RNA decay and haploinsufficiency. Using novel molecular analysis and nanopore technology, we demonstrated the pathogenesis of the rs794728589 (c.356+1G>A) splice variant in LMNA. This study highlights the importance of precise diagnostics in the clinical management and workup of cardiomyopathies.

Improvement of large copy number variant detection by whole genome nanopore sequencing

INTRODUCTION

Whole-genome sequencing using nanopore technologies can uncover structural variants, which are DNA rearrangements larger than 50 base pairs. Nanopore technologies can also characterize their boundaries with single-base accuracy, owing to the kilobase-long reads that encompass either full variants or their junctions. Other methods, such as next-generation short read sequencing or PCR assays, are limited in their capabilities to detect or characterize structural variants. However, the existing software for nanopore sequencing data analysis still reports incomplete variant sets, which also contain erroneous calls, a considerable obstacle for the molecular diagnosis or accurate genotyping of populations.

METHODS

We compared multiple factors affecting variant calling, such as reference genome version, aligner (minimap2, NGMLR, and lra) choice, and variant caller combinations (Sniffles, CuteSV, SVIM, and NanoVar), to find the optimal group of tools for calling large (>50 kb) deletions and duplications, using data from seven patients exhibiting gross gene defects on SERPINC1 and from a reference variant set as the control. The goal was to obtain the most complete, yet reasonably specific group of large variants using a single cell of PromethION sequencing, which yielded lower depth coverage than short-read sequencing. We also used a custom method for the statistical analysis of the coverage value to refine the resulting datasets.

RESULTS

We found that for large deletions and duplications (>50 kb), the existing software performed worse than for smaller ones, in terms of both sensitivity and specificity, and newer tools had not improved this. Our novel software, disCoverage, could polish variant callers' results, improving specificity by up to 62% and sensitivity by 15%, the latter requiring other data or samples.

CONCLUSION

We analyzed the current situation of >50-kb copy number variants with nanopore sequencing, which could be improved. The methods presented in this work could help to identify the known deletions and duplications in a set of patients, while also helping to filter out erroneous calls for these variants, which might aid the efforts to characterize a not-yet well-known fraction of genetic variability in the human genome.

Epigenetic tumor heterogeneity in the era of single-cell profiling with nanopore sequencing

Nanopore sequencing has brought the technology to the next generation in the science of sequencing. This is achieved through research advancing on: pore efficiency, creating mechanisms to control DNA translocation, enhancing signal-to-noise ratio, and expanding to long-read ranges. Heterogeneity regarding epigenetics would be broad as mutations in the epigenome are sensitive to cause new challenges in cancer research. Epigenetic enzymes which catalyze DNA methylation and histone modification are dysregulated in cancer cells and cause numerous heterogeneous clones to evolve. Detection of this heterogeneity in these clones plays an indispensable role in the treatment of various cancer types. With single-cell profiling, the nanopore sequencing technology could provide a simple sequence at long reads and is expected to be used soon at the bedside or doctor's office. Here, we review the advancements of nanopore sequencing and its use in the detection of epigenetic heterogeneity in cancer.

Satellite DNAs and human sex chromosome variation

Satellite DNAs are present on every chromosome in the cell and are typically enriched in repetitive, heterochromatic parts of the human genome. Sex chromosomes represent a unique genomic and epigenetic context. In this review, we first report what is known about satellite DNA biology on human X and Y chromosomes, including repeat content and organization, as well as satellite variation in typical euploid individuals. Then, we review sex chromosome aneuploidies that are among the most common types of aneuploidies in the general population, and are better tolerated than autosomal aneuploidies. This is demonstrated also by the fact that aging is associated with the loss of the X, and especially the Y chromosome. In addition, supernumerary sex chromosomes enable us to study general processes in a cell, such as analyzing heterochromatin dosage (i.e. additional Barr bodies and long heterochromatin arrays on Yq) and their downstream consequences. Finally, genomic and epigenetic organization and regulation of satellite DNA could influence chromosome stability and lead to aneuploidy. In this review, we argue that the complete annotation of satellite DNA on sex chromosomes in human, and especially in centromeric regions, will aid in explaining the prevalence and the consequences of sex chromosome aneuploidies.

Genotyping of familial Mediterranean fever gene (MEFV)-Single nucleotide polymorphism-Comparison of Nanopore with conventional Sanger sequencing

Background

Through continuous innovation and improvement, Nanopore sequencing has become a powerful technology. Because of its fast processing time, low cost, and ability to generate long reads, this sequencing technique would be particularly suitable for clinical diagnostics. However, its raw data accuracy is inferior in contrast to other sequencing technologies. This constraint still results in limited use of Nanopore sequencing in the field of clinical diagnostics and requires further validation and IVD certification.

Methods

We evaluated the performance of latest Nanopore sequencing in combination with a dedicated data-analysis pipeline for single nucleotide polymorphism (SNP) genotyping of the familial Mediterranean fever gene (MEFV) by amplicon sequencing of 47 clinical samples. Mutations in MEFV are associated with Mediterranean fever, a hereditary periodic fever syndrome. Conventional Sanger sequencing, which is commonly applied in clinical genetic diagnostics, was used as a reference method.

Results

Nanopore sequencing enabled the sequencing of 10 target regions within MEFV with high read depth (median read depth 7565x) in all samples and identified a total of 435 SNPs in the whole sample collective, of which 29 were unique. Comparison of both sequencing workflows showed a near perfect agreement with no false negative calls. Precision, Recall, and F1-Score of the Nanopore sequencing workflow were > 0.99, respectively.

Conclusions

These results demonstrated the great potential of current Nanopore sequencing for application in clinical diagnostics, at least for SNP genotyping by amplicon sequencing. Other more complex applications, especially structural variant identification, require further in-depth clinical validation.

Multiplex structural variant detection by whole-genome mapping and nanopore sequencing

Identification of structural variants (SVs) breakpoints is important in studying mutations, mutagenic causes, and functional impacts. Next-generation sequencing and whole-genome optical mapping are extensively used in SV discovery and characterization. However, multiple platforms and computational approaches are needed for comprehensive analysis, making it resource-intensive and expensive. Here, we propose a strategy combining optical mapping and cas9-assisted targeted nanopore sequencing to analyze SVs. Optical mapping can economically and quickly detect SVs across a whole genome but does not provide sequence-level information or precisely resolve breakpoints. Furthermore, since only a subset of all SVs is known to affect biology, we attempted to type a subset of all SVs using targeted nanopore sequencing. Using our approach, we resolved the breakpoints of five deletions, five insertions, and an inversion, in a single experiment.

Global Increase in Breast Cancer Incidence: Risk Factors and Preventive Measures

Breast cancer is a global cause for concern owing to its high incidence around the world. The alarming increase in breast cancer cases emphasizes the management of disease at multiple levels. The management should start from the beginning that includes stringent cancer screening or cancer registry to effective diagnostic and treatment strategies. Breast cancer is highly heterogeneous at morphology as well as molecular levels and needs different therapeutic regimens based on the molecular subtype. Breast cancer patients with respective subtype have different clinical outcome prognoses. Breast cancer heterogeneity emphasizes the advanced molecular testing that will help on-time diagnosis and improved survival. Emerging fields such as liquid biopsy and artificial intelligence would help to under the complexity of breast cancer disease and decide the therapeutic regimen that helps in breast cancer management. In this review, we have discussed various risk factors and advanced technology available for breast cancer diagnosis to combat the worst breast cancer status and areas that need to be focused for the better management of breast cancer.

Genetic Surveillance of Five SARS-CoV-2 Clinical Samples in Henan Province Using Nanopore Sequencing

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has rapidly spread and poses a major threat to public health worldwide. The whole genome sequencing plays a crucial role in virus surveillance and evolutionary analysis. In this study, five genome sequences of SARS-CoV-2 were obtained from nasopharyngeal swab samples from Zhengzhou, China. Following RNA extraction and cDNA synthesis, multiplex PCR was performed with two primer pools to produce the overlapped amplicons of ~1,200 bp. The viral genomes were obtained with 96% coverage using nanopore sequencing. Forty-five missense nucleotide mutations were identified; out of these, 5 mutations located at Nsp2, Nsp3, Nsp14, and ORF10 genes occurred with a <0.1% frequency in the global dataset. On the basis of mutation profiles, five genomes were clustered into two sublineages (B.1.617.2 and AY.31) or subclades (21A and 21I). The phylogenetic analysis of viral genomes from several regions of China and Myanmar revealed that five patients had different viral transmission chains. Taken together, we established a nanopore sequencing platform for genetic surveillance of SARS-CoV-2 and identified the variants circulating in Zhengzhou during August 2021. Our study provided crucial support for government policymaking and prevention and control of COVID-19.

Comparison of Illumina and Oxford Nanopore Sequencing Technologies for Pathogen Detection from Clinical Matrices Using Molecular Inversion Probes

Next-generation sequencing is rapidly finding footholds in numerous microbiological fields, including infectious disease diagnostics. Here, we describe a molecular inversion probe panel for the identification of bacterial, viral, and parasitic pathogens. We describe the ability of Illumina and Oxford Nanopore Technologies (ONT) to sequence small amplicons originating from this panel for the identification of pathogens in complex matrices. The panel correctly classified 31 bacterial pathogens directly from positive blood culture bottles with a genus-level concordance of 96.7% and 90.3% on the Illumina and ONT platforms, respectively. Both sequencing platforms detected 18 viral and parasitic organisms directly from mock clinical samples of plasma and whole blood at concentrations of 10⁴ PFU/mL with few exceptions. In general, Illumina sequencing exhibited greater read counts with lower percent mapped reads; however, this resulted in no effect on limits of detection compared with ONT sequencing. Mock clinical evaluation of the probe panel on the Illumina and ONT platforms resulted in positive predictive values of 0.91 and 0.88 and negative predictive values of 1 and 1 from de-identified human chikungunya virus samples compared with gold standard quantitative RT-PCR. Overall, these data show that molecular inversion probes are an adaptable technology capable of pathogen detection from complex sample matrices on current next-generation sequencing platforms.

Uncovering Essential Tremor Genetics: The Promise of Long-Read Sequencing

Long-read sequencing (LRS) technologies have been recently introduced to overcome intrinsic limitations of widely-used next-generation sequencing (NGS) technologies, namely the sequencing limited to short-read fragments (150–300 base pairs). Since its introduction, LRS has permitted many successes in unraveling hidden mutational mechanisms. One area in clinical neurology in need of rethinking as it applies to genetic mechanisms is essential tremor (ET). This disorder, among the most common in neurology, is a syndrome often exhibiting an autosomal dominant pattern of inheritance whose large phenotypic spectrum suggest a multitude of genetic etiologies. Exome sequencing has revealed the genetic etiology only in rare ET families (FUS, SORT1, SCN4A, NOS3, KCNS2, HAPLN4/BRAL2, and USP46). We hypothesize that a reason for this shortcoming may be non-classical genetic mechanism(s) underpinning ET, among them trinucleotide, tetranucleotide, or pentanucleotide repeat disorders. In support of this hypothesis, trinucleotide (e.g., GGC repeats in NOTCH2NLC) and pentanucleotide repeat disorders (e.g., ATTTC repeats in STARD7) have been revealed as pathogenic in patients with a past history of what has come to be referred to as “ET plus,” bilateral hand tremor associated with epilepsy and/or leukoencephalopathy. A systematic review of LRS in neurodegenerative disorders showed that 10 of the 22 (45%) genetic etiologies ascertained by LRS include tremor in their phenotypic spectrum, suggesting that future clinical applications of LRS for tremor disorders may uncover genetic subtypes of familial ET that have eluded NGS, particularly those with associated leukoencephalopathy or family history of epilepsy. LRS provides a pathway for potentially uncovering novel genes and genetic mechanisms, helping narrow the large proportion of “idiopathic” ET.

High-fidelity biosensing of dNTPs and nucleic acids by controllable subnanometer channel PaMscS

Current tools for dNTP analysis mainly rely on expensive fluorescent labeling, mass spectrometry or electrochemistry. Single-molecule assay by protein nanopores with an internal diameter of ca. 1-3.6 nm provides a useful tool for dNTP sensing. However, the most commonly used protein nanopores require additional modifications to enable dNTP detection. In this study, the PaMscS channel (mechanosensitive channel of small conductance from Pseudomonas aeruginosa) embedded in the bilayer lipid membrane (BLM) of E. coli polar lipid extract was applied as a nanopore for single molecular sensing. Two mutants of PaMscS nanopores on the side portal region (PaMscS W130A and PaMscS K180R) were selected for direct dNTP or pyrophosphoric acid (PPi) detection without aptamer or protein modification. Notably, the PaMscS mutant pore can be adjusted by regulation of osmolarity differences, which is crucial for the optimal detection of specific molecules. In addition, we established a PaMscS-based diagnosis method for the rapid sensing of disease-associated nucleic acids by monitoring the consumption of dNTPs, with 86% specificity and 100% sensitivity among 22 clinical samples. This protein nanopore, without aptamer or modification, paves a new way for dNTPs, PPi direct sensing and nucleic acid detection with low cost but high versatility.

ECNano: A cost-effective workflow for target enrichment sequencing and accurate variant calling on 4800 clinically significant genes using a single MinION flowcell

BACKGROUND

The application of long-read sequencing using the Oxford Nanopore Technologies (ONT) MinION sequencer is getting more diverse in the medical field. Having a high sequencing error of ONT and limited throughput from a single MinION flowcell, however, limits its applicability for accurate variant detection. Medical exome sequencing (MES) targets clinically significant exon regions, allowing rapid and comprehensive screening of pathogenic variants. By applying MES with MinION sequencing, the technology can achieve a more uniform capture of the target regions, shorter turnaround time, and lower sequencing cost per sample.

METHOD

We introduced a cost-effective optimized workflow, ECNano, comprising a wet-lab protocol and bioinformatics analysis, for accurate variant detection at 4800 clinically important genes and regions using a single MinION flowcell. The ECNano wet-lab protocol was optimized to perform long-read target enrichment and ONT library preparation to stably generate high-quality MES data with adequate coverage. The subsequent variant-calling workflow, Clair-ensemble, adopted a fast RNN-based variant caller, Clair, and was optimized for target enrichment data. To evaluate its performance and practicality, ECNano was tested on both reference DNA samples and patient samples.

RESULTS

ECNano achieved deep on-target depth of coverage (DoC) at average > 100× and > 98% uniformity using one MinION flowcell. For accurate ONT variant calling, the generated reads sufficiently covered 98.9% of pathogenic positions listed in ClinVar, with 98.96% having at least 30× DoC. ECNano obtained an average read length of 1000 bp. The long reads of ECNano also covered the adjacent splice sites well, with 98.5% of positions having ≥ 30× DoC. Clair-ensemble achieved > 99% recall and accuracy for SNV calling. The whole workflow from wet-lab protocol to variant detection was completed within three days.

CONCLUSION

We presented ECNano, an out-of-the-box workflow comprising (1) a wet-lab protocol for ONT target enrichment sequencing and (2) a downstream variant detection workflow, Clair-ensemble. The workflow is cost-effective, with a short turnaround time for high accuracy variant calling in 4800 clinically significant genes and regions using a single MinION flowcell. The long-read exon captured data has potential for further development, promoting the application of long-read sequencing in personalized disease treatment and risk prediction.

New Insights for Biosensing: Lessons from Microbial Defense Systems

Microorganisms have gained defense systems during the lengthy process of evolution over millions of years. Such defense systems can protect them from being attacked by invading species (e.g., CRISPR-Cas for establishing adaptive immune systems and nanopore-forming toxins as virulence factors) or enable them to adapt to different conditions (e.g., gas vesicles for achieving buoyancy control). These microorganism defense systems (MDS) have inspired the development of biosensors that have received much attention in a wide range of fields including life science research, food safety, and medical diagnosis. This Review comprehensively analyzes biosensing platforms originating from MDS for sensing and imaging biological analytes. We first describe a basic overview of MDS and MDS-inspired biosensing platforms (e.g., CRISPR-Cas systems, nanopore-forming proteins, and gas vesicles), followed by a critical discussion of their functions and properties. We then discuss several transduction mechanisms (optical, acoustic, magnetic, and electrical) involved in MDS-inspired biosensing. We further detail the applications of the MDS-inspired biosensors to detect a variety of analytes (nucleic acids, peptides, proteins, pathogens, cells, small molecules, and metal ions). In the end, we propose the key challenges and future perspectives in seeking new and improved MDS tools that can potentially lead to breakthrough discoveries in developing a new generation of biosensors with a combination of low cost; high sensitivity, accuracy, and precision; and fast detection. Overall, this Review gives a historical review of MDS, elucidates the principles of emulating MDS to develop biosensors, and analyzes the recent advancements, current challenges, and future trends in this field. It provides a unique critical analysis of emulating MDS to develop robust biosensors and discusses the design of such biosensors using elements found in MDS, showing that emulating MDS is a promising approach to conceptually advancing the design of biosensors.

Methods to Improve Molecular Diagnosis in Genomic Cold Cases in Pediatric Neurology

During the last decade, genetic testing has emerged as an important etiological diagnostic tool for Mendelian diseases, including pediatric neurological conditions. A genetic diagnosis has a considerable impact on disease management and treatment; however, many cases remain undiagnosed after applying standard diagnostic sequencing techniques. This review discusses various methods to improve the molecular diagnostic rates in these genomic cold cases. We discuss extended analysis methods to consider, non-Mendelian inheritance models, mosaicism, dual/multiple diagnoses, periodic re-analysis, artificial intelligence tools, and deep phenotyping, in addition to integrating various omics methods to improve variant prioritization. Last, novel genomic technologies, including long-read sequencing, artificial long-read sequencing, and optical genome mapping are discussed. In conclusion, a more comprehensive molecular analysis and a timely re-analysis of unsolved cases are imperative to improve diagnostic rates. In addition, our current understanding of the human genome is still limited due to restrictions in technologies. Novel technologies are now available that improve upon some of these limitations and can capture all human genomic variation more accurately. Last, we recommend a more routine implementation of high molecular weight DNA extraction methods that is coherent with the ability to use and/or optimally benefit from these novel genomic methods.

Using synthetic chromosome controls to evaluate the sequencing of difficult regions within the human genome

BACKGROUND

Next-generation sequencing (NGS) can identify mutations in the human genome that cause disease and has been widely adopted in clinical diagnosis. However, the human genome contains many polymorphic, low-complexity, and repetitive regions that are difficult to sequence and analyze. Despite their difficulty, these regions include many clinically important sequences that can inform the treatment of human diseases and improve the diagnostic yield of NGS.

RESULTS

To evaluate the accuracy by which these difficult regions are analyzed with NGS, we built an in silico decoy chromosome, along with corresponding synthetic DNA reference controls, that encode difficult and clinically important human genome regions, including repeats, microsatellites, HLA genes, and immune receptors. These controls provide a known ground-truth reference against which to measure the performance of diverse sequencing technologies, reagents, and bioinformatic tools. Using this approach, we provide a comprehensive evaluation of short- and long-read sequencing instruments, library preparation methods, and software tools and identify the errors and systematic bias that confound our resolution of these remaining difficult regions.

CONCLUSIONS

This study provides an analytical validation of diagnosis using NGS in difficult regions of the human genome and highlights the challenges that remain to resolve these difficult regions.

Use of MALDI-TOF mass spectrometry for virus identification: a review

The possibilities of virus identification, including SARS-CoV-2, by MALDI-TOF mass spectrometry are discussed in this review.

Application of long-read sequencing to elucidate complex pharmacogenomic regions: a proof of principle

The use of pharmacogenomics in clinical practice is becoming standard of care. However, due to the complex genetic makeup of pharmacogenes, not all genetic variation is currently accounted for. Here, we show the utility of long-read sequencing to resolve complex pharmacogenes by analyzing a well-characterised sample. This data consists of long reads that were processed to resolve phased haploblocks. 73% of pharmacogenes were fully covered in one phased haploblock, including 9/15 genes that are 100% complex. Variant calling accuracy in the pharmacogenes was high, with 99.8% recall and 100% precision for SNVs and 98.7% precision and 98.0% recall for Indels. For the majority of gene-drug interactions in the DPWG and CPIC guidelines, the associated genes could be fully resolved (62% and 63% respectively). Together, these findings suggest that long-read sequencing data offers promising opportunities in elucidating complex pharmacogenes and haplotype phasing while maintaining accurate variant calling.

Using DNA testing for the precise, definite, and low-cost diagnosis of sickle cell disease and other Haemoglobinopathies: findings from Tanzania

Background

Sickle cell disease (SCD) is an important cause of under-five mortality. Tanzania is the 5th country in the world with the highest births prevalence of SCD individuals. Significant advances in the neonatal diagnosis of SCD using rapid point-of-care testing have been made. However genetic confirmation is still required for positive cases, in uncertain cases, in multiply transfused patients, to resolve compound heterozygosity (Hb S/ β⁰ Thal or Hb S/ β⁺ thal) not uncommon in the coastal regions of East Africa and increasingly also for pre-marital counselling and potentially for future curative approaches such as gene therapy. The currently available DNA tests are prohibitively expensive. Here, we describe an easy-to-use, affordable and accurate β-globin sequencing approach that can be easily integrated within existing NBS for SCD and other haemoglobinopathies especially in Low- and Middle-income Countries.

Aim

To evaluate an affordable DNA technology for the diagnosis of Sickle cell disease and other haemoglobinopathies in a resource-limited setting.

Methods

Laboratory-based validation study was conducted by Muhimbili University of Health and Allied Sciences and the University of Oxford involving sequencing of the entire β -haemoglobin locus using the Oxford Nanopore MinION platform. A total number of 36 Dried blood spots and whole blood samples were subjected to conventional protein-based methods (isoelectric focusing, HPLC), and/or sequenced by the Sanger method as comparators.

Results

Sequencing results for SCD using the MinION were 100% concordant with those from the Sanger method. In addition, the long-read DNA sequencing method enabled the resolution of cases with unusual phenotypes which make up 1% of all children in Tanzania. The cost is £11/ sample for consumables, which is cheaper compared to other sequencing platforms.

Conclusions

This is the first report of a comprehensive single DNA assay as a definitive diagnostic test for SCD and other haemoglobinopathies. The test is fast, precise, accurate and affordable.

Father-to-offspring transmission of extremely long NOTCH2NLC repeat expansions with contractions: genetic and epigenetic profiling with long-read sequencing

Background

GGC repeat expansions in NOTCH2NLC are associated with neuronal intranuclear inclusion disease. Very recently, asymptomatic carriers with NOTCH2NLC repeat expansions were reported. In these asymptomatic individuals, the CpG island in NOTCH2NLC is hypermethylated, suggesting that two factors repeat length and DNA methylation status should be considered to evaluate pathogenicity. Long-read sequencing can be used to simultaneously profile genomic and epigenomic alterations. We analyzed four sporadic cases with NOTCH2NLC repeat expansion and their phenotypically normal parents. The native genomic DNA that retains base modification was sequenced on a per-trio basis using both PacBio and Oxford Nanopore long-read sequencing technologies. A custom workflow was developed to evaluate DNA modifications. With these two technologies combined, long-range DNA methylation information was integrated with complete repeat DNA sequences to investigate the genetic origins of expanded GGC repeats in these sporadic cases.

Results

In all four families, asymptomatic fathers had longer expansions (median: 522, 390, 528 and 650 repeats) compared with their affected offspring (median: 93, 117, 162 and 140 repeats, respectively). These expansions are much longer than the disease-causing range previously reported (in general, 41–300 repeats). Repeat lengths were extremely variable in the father, suggesting somatic mosaicism. Instability is more frequent in alleles with uninterrupted pure GGCs. Single molecule epigenetic analysis revealed complex DNA methylation patterns and epigenetic heterogeneity. We identified an aberrant gain-of-methylation region (2.2 kb in size beyond the CpG island and GGC repeats) in asymptomatic fathers. This methylated region was unmethylated in the normal allele with bilateral transitional zones with both methylated and unmethylated CpG dinucleotides, which may be protected from methylation to ensure NOTCH2NLC expression.

Conclusions

We clearly demonstrate that the four sporadic NOTCH2NLC-related cases are derived from the paternal GGC repeat contraction associated with demethylation. The entire genetic and epigenetic landscape of the NOTCH2NLC region was uncovered using the custom workflow of long-read sequence data, demonstrating the utility of this method for revealing epigenetic/mutational changes in repetitive elements, which are difficult to characterize by conventional short-read/bisulfite sequencing methods. Our approach should be useful for biomedical research, aiding the discovery of DNA methylation abnormalities through the entire genome.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-021-01192-5.

SpecHap: a diploid phasing algorithm based on spectral graph theory

Haplotype phasing plays an important role in understanding the genetic data of diploid eukaryotic organisms. Different sequencing technologies (such as next-generation sequencing or third-generation sequencing) produce various genetic data that require haplotype assembly. Although multiple diploid haplotype phasing algorithms exist, only a few will work equally well across all sequencing technologies. In this work, we propose SpecHap, a novel haplotype assembly tool that leverages spectral graph theory. On both in silico and whole-genome sequencing datasets, SpecHap consumed less memory and required less CPU time, yet achieved comparable accuracy with state-of-art methods across all the test instances, which comprises sequencing data from next-generation sequencing, linked-reads, high-throughput chromosome conformation capture, PacBio single-molecule real-time, and Oxford Nanopore long-reads. Furthermore, SpecHap successfully phased an individual Ambystoma mexicanum, a species with gigantic diploid genomes, within 6 CPU hours and 945MB peak memory usage, while other tools failed to yield results either due to memory overflow (40GB) or time limit exceeded (5 days). Our results demonstrated that SpecHap is scalable, efficient, and accurate for diploid phasing across many sequencing platforms.