MultiNanopolish: Refined grouping method for reducing redundant calculations in nanopolish

FPGA-based accelerator for adaptive banded event alignment in nanopore sequencing data analysis

BACKGROUND

Adaptive Banded Event Alignment (ABEA) stands as a critical algorithmic component in sequence polishing and DNA methylation detection, employing dynamic programming to align raw Nanopore signal with reference reads. Motivated by the observation that, compared to CPUs and GPUs, cutting-edge FPGAs demonstrate-in certain cases-superior performance at a reduced cost and energy consumption, this paper presents an efficient FPGA-based accelerator for ABEA, leveraging the inherent high parallelism and sequential access pattern within ABEA.

RESULT

Our proposed FPGA-based ABEA accelerator significantly enhances ABEA performance compared to the original CPU-based implementation in Nanopolish as well as the state-of-art acceleration on GPU and FPGA platforms. Specifically, targeting Xilinx VU9P, our accelerator achieves an average throughput speedup of 10.05 × over the CPU-only implementation, an average 1.81 × speedup over the state-of-art GPU acceleration with only 7.2% of the energy, and a speedup of 10.11 × compared to an existing FPGA accelerator.

CONCLUSION

Our work demonstrates that intensive genome analysis can benefit significantly from cutting-edge FPGAs, offering improvements in both performance and energy consumption.

Hierarchical architecture of neo-sex chromosomes and accelerated adaptive evolution in tortricid moths

Sex chromosomes can expand through fusion with autosomes, thereby acquiring unique evolutionary patterns. In butterflies and moths (Lepidoptera), these sex chromosome-autosome (SA) fusions occur relatively frequently, suggesting possible evolutionary advantages. Here, we investigated how SA fusion affects chromosome features and molecular evolution in leafroller moths (Lepidoptera: Tortricidae). Phylogenomic analysis showed that Tortricidae diverged ∼124 million years ago, accompanied by an SA fusion between the Merian elements M(20 + 17) and MZ. In contrast to partial autosomal fusions, the fused neo-Z Chromosome developed a hierarchical architecture, in which the three elements exhibit heterogeneous sequence features and evolutionary patterns. Specifically, the M17 part had a distinct base composition and chromatin domains. Unlike M20 and MZ, M17 was expressed at the same levels as autosomes in both sexes, compensating for the lost gene dosage in females. Concurrently, the SA fusion drove M17 as an evolutionary hotspot, accelerating the evolution of several genes related to ecological adaptation (e.g., ABCCs) and facilitating the divergence of closely related species, whereas the undercompensated M20 did not show such an effect. Thus, accelerated evolution under a novel pattern of dosage compensation may have favored the adaptive radiation of this group. This study demonstrates the association between a karyotype variant and adaptive evolution and explains the recurrent SA fusion in the Lepidoptera.

Simple, reference-independent assessment to empirically guide correction and polishing of hybrid microbial community metagenomic assembly

Hybrid metagenomic assembly of microbial communities, leveraging both long- and short-read sequencing technologies, is becoming an increasingly accessible approach, yet its widespread application faces several challenges. High-quality references may not be available for assembly accuracy comparisons common for benchmarking, and certain aspects of hybrid assembly may benefit from dataset-dependent, empiric guidance rather than the application of a uniform approach. In this study, several simple, reference-free characteristics-particularly coding gene content and read recruitment profiles-were hypothesized to be reliable indicators of assembly quality improvement during iterative error-fixing processes. These characteristics were compared to reference-dependent genome- and gene-centric analyses common for microbial community metagenomic studies. Two laboratory-scale bioreactors were sequenced with short- and long-read platforms, and assembled with commonly used software packages. Following long read assembly, long read correction and short read polishing were iterated up to ten times to resolve errors. These iterative processes were shown to have a substantial effect on gene- and genome-centric community compositions. Simple, reference-free assembly characteristics, specifically changes in gene fragmentation and short read recruitment, were robustly correlated with advanced analyses common in published comparative studies, and therefore are suitable proxies for hybrid metagenome assembly quality to simplify the identification of the optimal number of correction and polishing iterations. As hybrid metagenomic sequencing approaches will likely remain relevant due to the low added cost of short-read sequencing for differential coverage binning or the ability to access lower abundance community members, it is imperative that users are equipped to estimate assembly quality prior to downstream analyses.

Induction of phytoextraction, phytoprotection and growth promotion activities in Lupinus albus under mercury abiotic stress conditions by Peribacillus frigoritolerans subsp., mercuritolerans subsp. nov

Strain SAICEUPBMT was isolated from soils of Almadén (Ciudad Real, Spain), subjected to a high mercury concentration. SAICEUPBMT significantly increased aerial plant weight, aerial plant length and the development of secondary roots under mercury stress; increased twice the absorption of mercury by the plant, while favoring its development in terms of biomass. Similarly, plants inoculated with SAICEUPBMT and grown in soils contaminated with mercury, express a lower activity of antioxidant enzymes; catalase enzymes (CAT), superoxide dismutase (SOD), ascorbate peroxidase (APX), glutathione reductase (GR) for defense against ROS (reactive oxygen species). Whole genome analysis showed that ANI (95. 96 %), dDDH (72.9 %), AAI (93.3 %) and TETRA (0.99) values were on the thresholds established for differentiation a subspecies. The fatty acids analysis related the strain with the Peribacillus frigoritolerans species. And the synapomorphic analysis reveals a common ancestor with analysis related the strain with the Peribacillus frigoritolerans species. Results from genomic analysis together with differences in phenotypic features and chemotaxonomic analysis support the proposal of strain SAICEUPBMT as the type strain of a novel subspecies for which the name Peribacillus frigoritolerans subps. mercuritolerans sp. nov is proposed. The absence of virulence genes and transmissible resistance mechanisms reveals its safety for agronomic uses, under mercury stress conditions. The ability of Peribacillus frigoritolerans subsp. mercuritolerans subsp. nov to improve plant development was tested in a Lupinus albus model, demonstrating a great potential for plant phytoprotection against mercury stress.

Cryptic transmission and novel introduction of Dengue 1 and 2 genotypes in Colombia

Dengue fever remains as a public health challenge in Colombia, standing as the most prevalent infectious disease in the country. The cyclic nature of dengue epidemics, occurring approximately every 3 years, is intricately linked to meteorological events like El Niño Southern Oscillation (ENSO). Therefore, the Colombian system faces challenges in genomic surveillance. This study aimed to evaluate local dengue virus (DENV) transmission and genetic diversity in four Colombian departments with heterogeneous incidence patterns (department is first-level territorial units in Colombia). For this study, we processed 266 serum samples to identify DENV. Subsequently, we obtained 118 genome sequences by sequencing DENV genomes from serum samples of 134 patients infected with DENV-1 and DENV-2 serotypes. The predominant serotype was DENV-2 (108/143), with the Asian-American (AA) genotype (91/118) being the most prevalent one. Phylogenetic analysis revealed concurrent circulation of two lineages of both DENV-2 AA and DENV-1 V, suggesting ongoing genetic exchange with sequences from Venezuela and Cuba. The continuous migration of Venezuelan citizens into Colombia can contribute to this exchange, emphasizing the need for strengthened prevention measures in border areas. Notably, the time to most recent common ancestor analysis identified cryptic transmission of DENV-2 AA since approximately 2015, leading to the recent epidemic. This challenges the notion that major outbreaks are solely triggered by recent virus introductions, emphasizing the importance of active genomic surveillance. The study also highlighted the contrasting selection pressures on DENV-1 V and DENV-2 AA, with the latter experiencing positive selection, possibly influencing its transmissibility. The presence of a cosmopolitan genotype in Colombia, previously reported in Brazil and Peru, raises concerns about transmission routes, emphasizing the necessity for thorough DENV evolution studies. Despite limitations, the study underscores genomic epidemiology's crucial role in early detection and comprehension of DENV genotypes, recommending the use of advanced sequencing techniques as an early warning system to help prevent and control dengue outbreaks in Colombia and worldwide.

Salpa genome and developmental transcriptome analyses reveal molecular flexibility enabling reproductive success in a rapidly changing environment

Ocean warming favors pelagic tunicates, such as salps, that exhibit increasingly frequent and rapid population blooms, impacting trophic dynamics and composition and human marine-dependent activities. Salp blooms are a result of their successful reproductive life history, alternating seasonally between asexual and sexual protogynous (i.e. sequential) hermaphroditic stages. While predicting future salp bloom frequency and intensity relies on an understanding of the transitions during the sexual stage from female through parturition and subsequent sex change to male, these transitions have not been explored at the molecular level. Here we report the development of the first complete genome of S. thompsoni and the North Atlantic sister species S. aspera. Genome and comparative analyses reveal an abundance of repeats and G-quadruplex (G4) motifs, a highly stable secondary structure, distributed throughout both salp genomes, a feature shared with other tunicates that perform alternating sexual-asexual reproductive strategies. Transcriptional analyses across sexual reproductive stages for S. thompsoni revealed genes associated with male sex differentiation and spermatogenesis are expressed as early as birth and before parturition, inconsistent with previous descriptions of sequential sexual differentiation in salps. Our findings suggest salp are poised for reproductive success at birth, increasing the potential for bloom formation as ocean temperatures rise.

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.

GPU-accelerated and pipelined methylation calling

Motivation

The third-generation DNA sequencing technologies, such as Nanopore Sequencing, can operate at very high speeds and produce longer reads, which in turn results in a challenge for the computational analysis of such massive data. Nanopolish is a software package for signal-level analysis of Oxford Nanopore sequencing data. Call-methylation module of Nanopolish can detect methylation based on Hidden Markov Model (HMM). However, Nanopolish is limited by the long running time of some serial and computationally expensive processes. Among these, Adaptive Banded Event Alignment (ABEA) is the most time-consuming step, and the prior work, f5c, has already parallelized and optimized ABEA on GPU. As a result, the remaining methylation score calculation part, which uses HMM to identify if a given base is methylated or not, has become the new performance bottleneck.

Results

This article focuses on the call-methylation module that resides in the Nanopolish package. We propose Galaxy-methyl, which parallelizes and optimizes the methylation score calculation step on GPU and then pipelines the four steps of the call-methylation module. Galaxy-methyl increases the execution concurrency across CPUs and GPUs as well as hardware resource utilization for both. The experimental results collected indicate that Galaxy-methyl can achieve 3×-5× speedup compared with Nanopolish, and reduce the total execution time by 35% compared with f5c, on average.

Availability and implementation

The source code of Galaxy-methyl is available at https://github.com/fengyilin118/.

Comparing Long-Read Assemblers to Explore the Potential of a Sustainable Low-Cost, Low-Infrastructure Approach to Sequence Antimicrobial Resistant Bacteria With Oxford Nanopore Sequencing

Long-read sequencing (LRS) can resolve repetitive regions, a limitation of short read (SR) data. Reduced cost and instrument size has led to a steady increase in LRS across diagnostics and research. Here, we re-basecalled FAST5 data sequenced between 2018 and 2021 and analyzed the data in relation to gDNA across a large dataset (n = 200) spanning a wide GC content (25-67%). We examined whether re-basecalled data would improve the hybrid assembly, and, for a smaller cohort, compared long read (LR) assemblies in the context of antimicrobial resistance (AMR) genes and mobile genetic elements. We included a cost analysis when comparing SR and LR instruments. We compared the R9 and R10 chemistries and reported not only a larger yield but increased read quality with R9 flow cells. There were often discrepancies with ARG presence/absence and/or variant detection in LR assemblies. Flye-based assemblies were generally efficient at detecting the presence of ARG on both the chromosome and plasmids. Raven performed more quickly but inconsistently recovered small plasmids, notably a ∼15-kb Col-like plasmid harboring bla KPC . Canu assemblies were the most fragmented, with genome sizes larger than expected. LR assemblies failed to consistently determine multiple copies of the same ARG as identified by the Unicycler reference. Even with improvements to ONT chemistry and basecalling, long-read assemblies can lead to misinterpretation of data. If LR data are currently being relied upon, it is necessary to perform multiple assemblies, although this is resource (computing) intensive and not yet readily available/useable.

Long-Reads-Based Metagenomics in Clinical Diagnosis With a Special Focus on Fungal Infections

Identification of the causative infectious agent is essential in the management of infectious diseases, with the ideal diagnostic method being rapid, accurate, and informative, while remaining cost-effective. Traditional diagnostic techniques rely on culturing and cell propagation to isolate and identify the causative pathogen. These techniques are limited by the ability and the time required to grow or propagate an agent in vitro and the facts that identification based on morphological traits are non-specific, insensitive, and reliant on technical expertise. The evolution of next-generation sequencing has revolutionized genomic studies to generate more data at a cheaper cost. These are divided into short- and long-read sequencing technologies, depending on the length of reads generated during sequencing runs. Long-read sequencing also called third-generation sequencing emerged commercially through the instruments released by Pacific Biosciences and Oxford Nanopore Technologies, although relying on different sequencing chemistries, with the first one being more accurate both platforms can generate ultra-long sequence reads. Long-read sequencing is capable of entirely spanning previously established genomic identification regions or potentially small whole genomes, drastically improving the accuracy of the identification of pathogens directly from clinical samples. Long-read sequencing may also provide additional important clinical information, such as antimicrobial resistance profiles and epidemiological data from a single sequencing run. While initial applications of long-read sequencing in clinical diagnosis showed that it could be a promising diagnostic technique, it also has highlighted the need for further optimization. In this review, we show the potential long-read sequencing has in clinical diagnosis of fungal infections and discuss the pros and cons of its implementation.

Comparative evaluation of Nanopore polishing tools for microbial genome assembly and polishing strategies for downstream analysis

Assembling high-quality microbial genomes using only cost-effective Nanopore long-read systems such as Flongle is important to accelerate research on the microbial genome and the most critical point for this is the polishing process. In this study, we performed an evaluation based on BUSCO and Prokka gene prediction in terms of microbial genome assembly for eight state-of-the-art Nanopore polishing tools and combinations available. In the evaluation of individual tools, Homopolish, PEPPER, and Medaka demonstrated better results than others. In combination polishing, the second round Homopolish, and the PEPPER × medaka combination also showed better results than others. However, individual tools and combinations have specific limitations on usage and results. Depending on the target organism and the purpose of the downstream research, it is confirmed that there remain some difficulties in perfectly replacing the hybrid polishing carried out by the addition of a short-read. Nevertheless, through continuous improvement of the protein pores, related base-calling algorithms, and polishing tools based on improved error models, a high-quality microbial genome can be achieved using only Nanopore reads without the production of additional short-read data. The polishing strategy proposed in this study is expected to provide useful information for assembling the microbial genome using only Nanopore reads depending on the target microorganism and the purpose of the research.