Sequencing by avidity enables high accuracy with low reagent consumption

Human de novo mutation rates from a four-generation pedigree reference

Understanding the human de novo mutation (DNM) rate requires complete sequence information1. Here using five complementary short-read and long-read sequencing technologies, we phased and assembled more than 95% of each diploid human genome in a four-generation, twenty-eight-member family (CEPH 1463). We estimate 98-206 DNMs per transmission, including 74.5 de novo single-nucleotide variants, 7.4 non-tandem repeat indels, 65.3 de novo indels or structural variants originating from tandem repeats, and 4.4 centromeric DNMs. Among male individuals, we find 12.4 de novo Y chromosome events per generation. Short tandem repeats and variable-number tandem repeats are the most mutable, with 32 loci exhibiting recurrent mutation through the generations. We accurately assemble 288 centromeres and six Y chromosomes across the generations and demonstrate that the DNM rate varies by an order of magnitude depending on repeat content, length and sequence identity. We show a strong paternal bias (75-81%) for all forms of germline DNM, yet we estimate that 16% of de novo single-nucleotide variants are postzygotic in origin with no paternal bias, including early germline mosaic mutations. We place all this variation in the context of a high-resolution recombination map (~3.4 kb breakpoint resolution) and find no correlation between meiotic crossover and de novo structural variants. These near-telomere-to-telomere familial genomes provide a truth set to understand the most fundamental processes underlying human genetic variation.

Nanopore-based consensus sequencing enables accurate multimodal tumor cell-free DNA profiling

Shallow genome-wide cell-free DNA sequencing holds great promise for noninvasive cancer monitoring by providing reliable copy number alteration (CNA) and fragmentomic profiles. Single-nucleotide variations (SNVs) are, however, much harder to identify with low sequencing depth due to sequencing errors. Here, we present Nanopore Rolling Circle Amplification (RCA)-enhanced Consensus Sequencing (NanoRCS), which leverages RCA and consensus calling based on genome-wide long-read nanopore sequencing to enable simultaneous multimodal tumor fraction (TF) estimation through SNVs, CNAs, and fragmentomics. The efficacy of NanoRCS is tested on 18 cancer patient samples and seven healthy controls, demonstrating its ability to reliably detect TFs as low as 0.24%. In vitro experiments confirm that SNV measurements are essential for detecting TFs below 3%. NanoRCS provides an opportunity for cost-effective and rapid sample processing, which aligns well with clinical needs, particularly in settings where quick and accurate cancer monitoring is essential for personalized treatment strategies.

Engineering Gizmos for Short Cancer Genetic Fragments Discrimination

Currently, mankind is fiercely struggling with cancer. Recently, we have been winning the battle against cancer through precision medicine and accompanying diagnostic methods, and we are raising many hopes with blockbuster drugs. It would be even better if we could read the cancer nucleotide sequence, identify them in advance, and suggest treatments simultaneously. However, this may be an impossible dream because it takes a lot of time and effort to diagnose and ensure all the long gene sequences of cancer at once. Thus, victory will be even closer if a rapid and accurate diagnosis of the cancer-specific gene biomarkers that will soon be imprinted can be made. With the advent of nanotechnology, a new short cancer diagnostic toolkit has been proposed to achieve the goal. This review presents a small diagnostic device that detects certain cancers' genetic fragments (simply 'Gizmo'). The development of numerous diagnostic methods has focused on (1) directly detecting pre-selectively targeted genes using novel diagnostic systems, and (2) indirectly detecting substantial improvements in diagnostic sensitivity only through detection signal amplification without existing gene amplification steps. Our fight against cancer is not a dream, but the result of success, and it is assumed that victory will accelerate as soon as possible.

Implementation of Multi-Criteria Decision-Making for Selecting Most Effective Genome Sequencing Technology

Background/Objectives: In recent years, molecular diagnosis has become increasingly critical in identifying human pathogens with unknown genes. Methods: An innovative approach, the fuzzy-based preference ranking organization method for enrichment evaluation (PROMETHEE) technique, one of the most effective multi-criteria decision-making (MCDM) methods, was used to evaluate criteria, including portability, generation type, max read/run, max output data/run, processing time/run, read length, accuracy, diagnostic sensitivity, test minimum loading volume, test cost/run, instrument cost, error rate, throughput capability, ability to sequence the large whole genome, small whole genome, and exome and large panel, mutation detection ability, whole-genome sequencing with single-stranded sequencing, and single-stranded sequencing accuracy, to determine the most suitable sequencing technology. Results: Based on the analysis, the Avidite Base Chemistry (ABC), Nanopore, and Illumina sequencing platforms sequentially emerged as the most favorable options based on their net flows of 0.0346, 0.0041, and 0.0003, respectively. Conclusions: Our findings provide important data to facilitate the selection of genome detection technologies. Through the use of innovative approaches, complex evaluations can be analyzed and the right choices can be made. Importantly, the technique has a degree of subjectivity, so varying conditions may lead to different findings.

Scalable co-sequencing of RNA and DNA from individual nuclei

The ideal technology for directly investigating the relationship between genotype and phenotype would analyze both RNA and DNA genome-wide and with single-cell resolution; however, existing tools lack the throughput required for comprehensive analysis of complex tumors and tissues. We introduce a highly scalable method for jointly profiling DNA and expression following nucleosome depletion (DEFND-seq). In DEFND-seq, nuclei are nucleosome-depleted, tagmented and separated into individual droplets for messenger RNA and genomic DNA barcoding. Once nuclei have been depleted of nucleosomes, subsequent steps can be performed using the widely available 10x Genomics droplet microfluidic technology and commercial kits. We demonstrate the production of high-complexity mRNA and gDNA sequencing libraries from thousands of individual nuclei from cell lines, fresh and archived surgical specimens for associating gene expression with both copy number and single-nucleotide variants.

Comparison of Sequencing-by-Synthesis and Avidity Base Chemistry Next-Generation Sequencing Platforms in Identifying Somatic Variants of Hematological Malignancies

BACKGROUND

This study evaluates the performance of the Illumina NextSeq™ 550 and the Element Biosciences AVITI™ next-generation sequencing (NGS) system, in detecting single nucleotide variants (SNVs) and gene fusions.

METHODS

A set of 66 NGS libraries, consisting of 33 DNA, 24 cDNA, and triplicates of 3 control libraries, were prepared from bone marrow samples targeting 63 genes and related fusions, and initially sequenced using the NextSeq 550 in the Cleveland Clinic's molecular diagnostic laboratory. The same libraries were subsequently sequenced on the AVITI. The resulting data were analyzed using a combination of Cleveland Clinic developed pipelines and ArcherDx virtual machine software.

RESULTS

The study found that all 105 SNVs and 39 gene fusions identified by the NextSeq 550 were also detected in the AVITI, demonstrating a high degree of concordance between the platforms. The analyses revealed R2 values of 0.86 for read depth and 0.96 for VAF of the 105 DNA variants, and 0.95 for read depth and 0.97 for fusion percentage of the 39 fusion variants. In the reproducibility studies, the VAF and fusion percentage of all variants were within 2 standard deviations of the mean when the same positive controls were sequenced 3 times on the AVITI.

CONCLUSIONS

These results indicate that the NextSeq 550 and the AVITI provide comparable performance in terms of accuracy and sensitivity for variant detection. Notably, the AVITI chemistry requires substantially lower PhiX input than the NextSeq 550 needs for this application. This results in substantial cost and efficiency benefits.

Small variant benchmark from a complete assembly of X and Y chromosomes

The sex chromosomes contain complex, important genes impacting medical phenotypes, but differ from the autosomes in their ploidy and large repetitive regions. To enable technology developers along with research and clinical laboratories to evaluate variant detection on male sex chromosomes X and Y, we create a small variant benchmark set with 111,725 variants for the Genome in a Bottle HG002 reference material. We develop an active evaluation approach to demonstrate the benchmark set reliably identifies errors in challenging genomic regions and across short and long read callsets. We show how complete assemblies can expand benchmarks to difficult regions, but highlight remaining challenges benchmarking variants in long homopolymers and tandem repeats, complex gene conversions, copy number variable gene arrays, and human satellites.

Innovations in Short-Read Sequencing Technologies and Their Applications to Clinical Genomics

BACKGROUND

Massively parallel sequencing (MPS) of nucleic acids has been a transformative technology for basic and applied genomic science, increasing efficiencies and decreasing costs to enable studies of unprecedented scope and impact. In clinical settings, these technological and scientific advances have led to the development of tests that are increasingly fast, comprehensive, and more frequently employed. Practitioners of genomic medicine have applied these tools across clinical settings, including diagnosis of inherited disorders and cancers and infectious disease detection and surveillance. In recent years, the commercial marketplace for MPS sequencers and reagents has been dominated by a few companies. The growing demand for sequencing has led to the recent emergence of several new sequencing platforms with techniques that may provide alternatives or improvements to existing workflows or allow the adoption of sequencing workflows in new settings. Clinical genomics laboratories will evaluate these platforms from a unique perspective, focusing on how technological advancements can improve patient care.

CONTENT

This review describes short-read sequencing platforms provided by Illumina, Element Biosciences, MGI, PacBio, Singular Genomics, Thermo Fisher Scientific, and Ultima Genomics. This review discusses their innovative approaches, principles, workflows, and applications.

SUMMARY

This review aims to inform laboratory geneticists, clinicians, and researchers about emerging short-read technologies and their applications in clinical genomics. By highlighting their principles and potential contributions, we aim to assist laboratories in selecting suitable solutions for their sequencing needs considering key factors such as applications, throughput, and integration with existing laboratory workflows.

Sequencing and Optical Genome Mapping for the Adventurous Chemist

This review provides a comprehensive overview of the chemistries and workflows of the sequencing methods that have been or are currently commercially available, providing a very brief historical introduction to each method. The main optical genome mapping approaches are introduced in the same manner, although only a subset of these are or have ever been commercially available. The review comes with a deck of slides containing all of the figures for ease of access and consultation.

Improved characterization of 3' single-cell RNA-seq libraries with paired-end avidity sequencing

Prevailing poly(dT)-primed 3' single-cell RNA-seq protocols generate barcoded cDNA fragments containing the reverse transcriptase priming site or in principle the polyadenylation site. Direct sequencing across this site was historically difficult because of DNA sequencing errors induced by the homopolymeric primer at the 'barcode' end. Here, we evaluate the capability of 'avidity base chemistry' DNA sequencing from Element Biosciences to sequence through the primer and enable accurate paired-end read alignment and precise quantification of polyadenylation sites. We find that the Element Aviti instrument sequences through the thymine homopolymer into the subsequent cDNA sequence without detectable loss of accuracy. The additional sequence enables direct and independent assignment of reads to polyadenylation sites, which bypasses the complexities and limitations of conventional approaches but does not consistently improve read mapping rates compared to single-end alignment. We also characterize low-level artifacts and demonstrate necessary adjustments to adapter trimming and sequence alignment regardless of platform, particularly in the context of extended read lengths. Our analyses confirm that Element avidity sequencing is an effective alternative to Illumina sequencing for standard single-cell RNA-seq, particularly for polyadenylation site measurement but do not rule out the potential for similar performance from other emerging platforms.

Efficient small fragment sequencing of human, cattle, and bison miRNA, small RNA, or csRNA-seq libraries using AVITI

BACKGROUND

Next-Generation Sequencing (NGS) catalyzed breakthroughs across various scientific domains. Illumina's sequencing by synthesis method has long been central to NGS, but new sequencing methods like Element Biosciences' AVITI technology are emerging. AVITI is reported to offer improved signal-to-noise ratios and cost reductions. However, its reliance on rolling circle amplification, which can be affected by polymer size, raises questions about its effectiveness in sequencing small RNAs (sRNAs) such as microRNAs (miRNAs), small nucleolar RNAs (snoRNAs), and many others. These sRNAs are crucial regulators of gene expression and involved in various biological processes. Additionally, capturing capped small RNAs (csRNA-seq) is a powerful method for mapping active or "nascent" RNA polymerase II transcription initiation in tissues and clinical samples.

RESULTS

Here, we report a new protocol for seamlessly sequencing short fragments on the AVITI and demonstrate that AVITI and Illumina sequencing technologies equivalently capture human, cattle (Bos taurus), and bison (Bison bison) sRNA or csRNA sequencing libraries, increasing confidence in both sequencing approaches. Additionally, analysis of generated nascent transcription start site (TSS) data for cattle and bison revealed inaccuracies in their current genome annotations, underscoring the potential and necessity to translate small and nascent RNA sequencing methodologies to livestock.

CONCLUSIONS

Our accelerated and optimized protocol bridges the advantages of AVITI sequencing with critical methods that rely on sequencing short fragments. This advance bolsters the utility of AVITI technology alongside traditional Illumina platforms, offering new opportunities for NGS applications.

Detection and analysis of complex structural variation in human genomes across populations and in brains of donors with psychiatric disorders

Complex structural variations (cxSVs) are often overlooked in genome analyses due to detection challenges. We developed ARC-SV, a probabilistic and machine-learning-based method that enables accurate detection and reconstruction of cxSVs from standard datasets. By applying ARC-SV across 4,262 genomes representing all continental populations, we identified cxSVs as a significant source of natural human genetic variation. Rare cxSVs have a propensity to occur in neural genes and loci that underwent rapid human-specific evolution, including those regulating corticogenesis. By performing single-nucleus multiomics in postmortem brains, we discovered cxSVs associated with differential gene expression and chromatin accessibility across various brain regions and cell types. Additionally, cxSVs detected in brains of psychiatric cases are enriched for linkage with psychiatric GWAS risk alleles detected in the same brains. Furthermore, our analysis revealed significantly decreased brain-region- and cell-type-specific expression of cxSV genes, specifically for psychiatric cases, implicating cxSVs in the molecular etiology of major neuropsychiatric disorders.

Hypothesis-based investigation of known AD risk variants reveals the genetic underpinnings of neuropathological lesions observed in Alzheimer's-type dementia

Alzheimer's disease (AD) is the leading cause of dementia worldwide. Besides neurofibrillary tangles and amyloid beta (Aβ) plaques, a wide range of co-morbid neuropathological features can be observed in AD brains. Since AD has a very strong genetic background and displays a wide phenotypic heterogeneity, this study aims at investigating the genetic underpinnings of co-morbid and hallmark neuropathological lesions. This was realized by obtaining the genotypes for 75 AD risk variants from low-coverage whole-genome sequencing data for 325 individuals from the Leuven Brain Collection. Association testing with deeply characterized neuropathological lesions revealed a strong and likely direct effect of rs117618017, a SNP in exon 1 of APH1B, with tau-related pathology. Second, a relation between APOE and granulovacuolar degeneration, a proxy for necroptosis, was also discovered in addition to replication of the well-known association of APOE with AD hallmark neuropathological lesions. Additionally, several nominal associations with AD risk genes were detected for pTDP pathology, α-synuclein lesions and pTau-related pathology. These findings were confirmed in a meta-analysis with three independent cohorts. For example, we replicated a prior association between TPCN1 (rs6489896) and LATE-NC risk. Furthermore, we identified new putative LATE-NC-linked SNPs, including rs7068231, located upstream of ANK3. We found association between BIN1 (rs6733839) and α-synuclein pathology, and replicated a prior association between USP6NL (rs7912495) and Lewy body pathology. Additionally, we also found that UMAD1 (rs6943429) was nominally associated with Lewy body pathology. Overall, these results contribute to a broader general understanding of how AD risk variants discovered in large-scale clinical genome-wide association studies are involved in the pathological mechanisms of AD and indicate the importance of downstream elimination of phenotypic heterogeneity introduced in these studies.

Parallel molecular data storage by printing epigenetic bits on DNA

DNA storage has shown potential to transcend current silicon-based data storage technologies in storage density, longevity and energy consumption1-3. However, writing large-scale data directly into DNA sequences by de novo synthesis remains uneconomical in time and cost4. We present an alternative, parallel strategy that enables the writing of arbitrary data on DNA using premade nucleic acids. Through self-assembly guided enzymatic methylation, epigenetic modifications, as information bits, can be introduced precisely onto universal DNA templates to enact molecular movable-type printing. By programming with a finite set of 700 DNA movable types and five templates, we achieved the synthesis-free writing of approximately 275,000 bits on an automated platform with 350 bits written per reaction. The data encoded in complex epigenetic patterns were retrieved high-throughput by nanopore sequencing, and algorithms were developed to finely resolve 240 modification patterns per sequencing reaction. With the epigenetic information bits framework, distributed and bespoke DNA storage was implemented by 60 volunteers lacking professional biolab experience. Our framework presents a new modality of DNA data storage that is parallel, programmable, stable and scalable. Such an unconventional modality opens up avenues towards practical data storage and dual-mode data functions in biomolecular systems.

Exploring the unmapped cysteine redox proteoform landscape

Cysteine redox proteoforms define the diverse molecular states that proteins with cysteine residues can adopt. A protein with one cysteine residue must adopt one of two binary proteoforms: reduced or oxidized. Their numbers scale: a protein with 10 cysteine residues must assume one of 1,024 proteoforms. Although they play pivotal biological roles, the vast cysteine redox proteoform landscape comprising vast numbers of theoretical proteoforms remains largely uncharted. Progress is hampered by a general underappreciation of cysteine redox proteoforms, their intricate complexity, and the formidable challenges that they pose to existing methods. The present review advances cysteine redox proteoform theory, scrutinizes methodological barriers, and elaborates innovative technologies for detecting unique residue-defined cysteine redox proteoforms. For example, chemistry-enabled hybrid approaches combining the strengths of top-down mass spectrometry (TD-MS) and bottom-up mass spectrometry (BU-MS) for systematically cataloguing cysteine redox proteoforms are delineated. These methods provide the technological means to map uncharted redox terrain. To unravel hidden redox regulatory mechanisms, discover new biomarkers, and pinpoint therapeutic targets by mining the theoretical cysteine redox proteoform space, a community-wide initiative termed the "Human Cysteine Redox Proteoform Project" is proposed. Exploring the cysteine redox proteoform landscape could transform current understanding of redox biology.

A comprehensive adsorption and desorption study on the interaction of DNA oligonucleotides with TiO2 nanolayers

The utilization of TiO2 nanolayers that possess excellent biocompatibility and physical properties in DNA sensing and sequencing remains largely to be explored. To examine their applicability in gene sequencing, a comprehensive study on the interaction of DNA oligonucleotides with TiO2 nanolayers was performed through adsorption and desorption experiments. TiO2 nanolayers with 10 nm thickness were fabricated via magnetron sputtering onto a 6-inch silicon wafer. A simple chip block method, validated via quartz crystal microbalance experiments with dissipation monitoring (QCM-D), was proposed to study the adsorption behaviors and interaction mechanisms under a variety of critical influencing factors, including DNA concentration, length, and type, adsorption time, pH, and metal ions. It is determined that the adsorption takes 2 h to reach saturation in the MES solution and the adsorption capacity is significantly enhanced by lowering the pH due to the isoelectric point being pH = 6 for TiO2. The adsorption percentages of nucleobases are largely similar in the MES solution while following 5T = 5G > 5C > 5A in HEPES buffer for an adsorption duration of 2.5 h. Through pre-adsorption experiments, it is deduced that DNA oligonucleotides are horizontally adsorbed on the nanolayer. This further demonstrates that mono-, di-, and tri-valent metal ions promote the adsorption, whereas Zn2+ has strong adsorption by inducing DNA condensation. Based on the desorption experiments, it is revealed that electrostatic force dominates the adsorption over van der Waals force and hydrogen bonds. The phosphate group is the main functional group for adsorption, and the adsorption strength increases with the length of the oligonucleotide. This study provides comprehensive data on the adsorption of DNA oligonucleotides onto TiO2 nanolayers and clarifies the interaction mechanisms therein, which will be valuable for applications of TiO2 in DNA-related applications.

A familial, telomere-to-telomere reference for human de novo mutation and recombination from a four-generation pedigree

Using five complementary short- and long-read sequencing technologies, we phased and assembled >95% of each diploid human genome in a four-generation, 28-member family (CEPH 1463) allowing us to systematically assess de novo mutations (DNMs) and recombination. From this family, we estimate an average of 192 DNMs per generation, including 75.5 de novo single-nucleotide variants (SNVs), 7.4 non-tandem repeat indels, 79.6 de novo indels or structural variants (SVs) originating from tandem repeats, 7.7 centromeric de novo SVs and SNVs, and 12.4 de novo Y chromosome events per generation. STRs and VNTRs are the most mutable with 32 loci exhibiting recurrent mutation through the generations. We accurately assemble 288 centromeres and six Y chromosomes across the generations, documenting de novo SVs, and demonstrate that the DNM rate varies by an order of magnitude depending on repeat content, length, and sequence identity. We show a strong paternal bias (75-81%) for all forms of germline DNM, yet we estimate that 17% of de novo SNVs are postzygotic in origin with no paternal bias. We place all this variation in the context of a high-resolution recombination map (~3.5 kbp breakpoint resolution). We observe a strong maternal recombination bias (1.36 maternal:paternal ratio) with a consistent reduction in the number of crossovers with increasing paternal (r=0.85) and maternal (r=0.65) age. However, we observe no correlation between meiotic crossover locations and de novo SVs, arguing against non-allelic homologous recombination as a predominant mechanism. The use of multiple orthogonal technologies, near-telomere-to-telomere phased genome assemblies, and a multi-generation family to assess transmission has created the most comprehensive, publicly available "truth set" of all classes of genomic variants. The resource can be used to test and benchmark new algorithms and technologies to understand the most fundamental processes underlying human genetic variation.

Genome-wide detection of somatic mosaicism at short tandem repeats

MOTIVATION

Somatic mosaicism has been implicated in several developmental disorders, cancers, and other diseases. Short tandem repeats (STRs) consist of repeated sequences of 1-6 bp and comprise >1 million loci in the human genome. Somatic mosaicism at STRs is known to play a key role in the pathogenicity of loci implicated in repeat expansion disorders and is highly prevalent in cancers exhibiting microsatellite instability. While a variety of tools have been developed to genotype germline variation at STRs, a method for systematically identifying mosaic STRs is lacking.

RESULTS

We introduce prancSTR, a novel method for detecting mosaic STRs from individual high-throughput sequencing datasets. prancSTR is designed to detect loci characterized by a single high-frequency mosaic allele, but can also detect loci with multiple mosaic alleles. Unlike many existing mosaicism detection methods for other variant types, prancSTR does not require a matched control sample as input. We show that prancSTR accurately identifies mosaic STRs in simulated data, demonstrate its feasibility by identifying candidate mosaic STRs in Illumina whole genome sequencing data derived from lymphoblastoid cell lines for individuals sequenced by the 1000 Genomes Project, and evaluate the use of prancSTR on Element and PacBio data. In addition to prancSTR, we present simTR, a novel simulation framework which simulates raw sequencing reads with realistic error profiles at STRs.

AVAILABILITY AND IMPLEMENTATION

prancSTR and simTR are freely available at https://github.com/gymrek-lab/trtools. Detailed documentation is available at https://trtools.readthedocs.io/.

Improved characterization of single-cell RNA-seq libraries with paired-end avidity sequencing

Prevailing poly(dT)-primed 3' single-cell RNA-seq protocols generate barcoded cDNA fragments containing the reverse transcriptase priming site, which is expected to be the poly(A) tail or a genomic adenine homopolymer. Direct sequencing across this priming site was historically difficult because of DNA sequencing errors induced by the homopolymeric primer at the 'barcode' end. Here, we evaluate the capability of "avidity base chemistry" DNA sequencing from Element Biosciences to sequence through this homopolymer accurately, and the impact of the additional cDNA sequence on read alignment and precise quantification of polyadenylation site usage. We find that the Element Aviti instrument sequences through the thymine homopolymer into the subsequent cDNA sequence without detectable loss of accuracy. The resulting paired-end alignments enable direct and independent assignment of reads to polyadenylation sites, which bypasses complexities and limitations of conventional approaches but does not consistently improve read mapping rates compared to single-end alignment. We also characterize low-level artifacts and arrive at an adjusted adapter trimming and alignment workflow that significantly improves the alignment of sequence data from Element and Illumina, particularly in the context of extended read lengths. Our analyses confirm that Element avidity sequencing is an effective alternative to Illumina sequencing for standard single-cell RNA-seq, particularly for polyadenylation site analyses but do not rule out the potential for similar performance from other emerging platforms.

Fine-scale characterization of the soybean rhizosphere microbiome via synthetic long reads and avidity sequencing

BACKGROUND

The rhizosphere microbiome displays structural and functional dynamism driven by plant, microbial, and environmental factors. While such plasticity is a well-evidenced determinant of host health, individual and community-level microbial activity within the rhizosphere remain poorly understood, due in part to the insufficient taxonomic resolution achieved through traditional marker gene amplicon sequencing. This limitation necessitates more advanced approaches (e.g., long-read sequencing) to derive ecological inferences with practical application. To this end, the present study coupled synthetic long-read technology with avidity sequencing to investigate eukaryotic and prokaryotic microbiome dynamics within the soybean (Glycine max) rhizosphere under field conditions.

RESULTS

Synthetic long-read sequencing permitted de novo reconstruction of the entire 18S-ITS1-ITS2 region of the eukaryotic rRNA operon as well as all nine hypervariable regions of the 16S rRNA gene. All full-length, mapped eukaryotic amplicon sequence variants displayed genus-level classification, and 44.77% achieved species-level classification. The resultant eukaryotic microbiome encompassed five kingdoms (19 genera) of protists in addition to fungi - a depth unattainable with conventional short-read methods. In the prokaryotic fraction, every full-length, mapped amplicon sequence variant was resolved at the species level, and 23.13% at the strain level. Thirteen species of Bradyrhizobium were thereby distinguished in the prokaryotic microbiome, with strain-level identification of the two Bradyrhizobium species most reported to nodulate soybean. Moreover, the applied methodology delineated structural and compositional dynamism in response to experimental parameters (i.e., growth stage, cultivar, and biostimulant application), unveiled a saprotroph-rich core microbiome, provided empirical evidence for host selection of mutualistic taxa, and identified key microbial co-occurrence network members likely associated with edaphic and agronomic properties.

CONCLUSIONS

This study is the first to combine synthetic long-read technology and avidity sequencing to profile both eukaryotic and prokaryotic fractions of a plant-associated microbiome. Findings herein provide an unparalleled taxonomic resolution of the soybean rhizosphere microbiota and represent significant biological and technological advancements in crop microbiome research.

ProcaryaSV: structural variation detection pipeline for bacterial genomes using short-read sequencing

BACKGROUND

Structural variations play an important role in bacterial genomes. They can mediate genome adaptation quickly in response to the external environment and thus can also play a role in antibiotic resistance. The detection of structural variations in bacteria is challenging, and the recognition of even small rearrangements can be important. Even though most detection tools are aimed at and benchmarked on eukaryotic genomes, they can also be used on prokaryotic genomes. The key features of detection are the ability to detect small rearrangements and support haploid genomes. Because of the limiting performance of a single detection tool, combining the detection abilities of multiple tools can lead to more robust results. There are already available workflows for structural variation detection for long-reads technologies and for the detection of single-nucleotide variation and indels, both aimed at bacteria. Yet we are unaware of structural variations detection workflows for the short-reads sequencing platform. Motivated by this gap we created our workflow. Further, we were interested in increasing the detection performance and providing more robust results.

RESULTS

We developed an open-source bioinformatics pipeline, ProcaryaSV, for the detection of structural variations in bacterial isolates from paired-end short sequencing reads. Multiple tools, starting with quality control and trimming of sequencing data, alignment to the reference genome, and multiple structural variation detection tools, are integrated. All the partial results are then processed and merged with an in-house merging algorithm. Compared with a single detection approach, ProcaryaSV has improved detection performance and is a reproducible easy-to-use tool.

CONCLUSIONS

The ProcaryaSV pipeline provides an integrative approach to structural variation detection from paired-end next-generation sequencing of bacterial samples. It can be easily installed and used on Linux machines. It is publicly available on GitHub at https://github.com/robinjugas/ProcaryaSV .

Chasing Sequencing Perfection: Marching Toward Higher Accuracy and Lower Costs

Next-generation sequencing (NGS), represented by Illumina platforms, has been an essential cornerstone of basic and applied research. However, the sequencing error rate of 1 per 1000 bp (10-3) represents a serious hurdle for research areas focusing on rare mutations, such as somatic mosaicism or microbe heterogeneity. By examining the high-fidelity sequencing methods developed in the past decade, we summarized three major factors underlying errors and the corresponding 12 strategies mitigating these errors. We then proposed a novel framework to classify 11 preexisting representative methods according to the corresponding combinatory strategies and identified three trends that emerged during methodological developments. We further extended this analysis to eight long-read sequencing methods, emphasizing error reduction strategies. Finally, we suggest two promising future directions that could achieve comparable or even higher accuracy with lower costs in both NGS and long-read sequencing.

Genetic analysis of seven patients with inherited ichthyosis and Nagashima‑type palmoplantar keratoderma

Inherited ichthyosis comprises a series of heterogeneous dermal conditions; it mainly manifests as widespread hyperkeratosis, xerosis and scaling of the skin. At times, overlapping symptoms require differential diagnosis between ichthyosis and several other similar disorders. The present study reports seven patients with confirmed or suspected to be associated with ichthyosis by conducting a thorough clinical and genetic investigation. Genetic testing was conducted using whole‑exome sequencing, with Sanger sequencing as the validation method. The MEGA7 program was used to analyze the conservation of amino acid residues affected by the detected missense variants. The enrolled patients exhibited ichthyosis‑like but distinct clinical manifestations. Genetic analysis identified diagnostic variations in the FLG, STS, KRT10 and SERPINB7 genes and clarified the carrying status of each variant in the respective family members. The two residues affected by the detected missense variants remained conserved across multiple species. Of note, the two variants, namely STS: c.452C>T(p.P151L) and c.647_650del(p.L216fs) are novel. In conclusion, a clear genetic differential diagnosis was made for the enrolled ichthyosis‑associated patients; the study findings also extended the mutation spectrum of ichthyosis and provided solid evidence for the counseling of the affected families.

Efficient small fragment sequencing of human, cow, and bison miRNA, small RNA or csRNA-seq libraries using AVITI

Next-Generation Sequencing (NGS) catalyzed breakthroughs across various scientific domains. Illumina's sequencing by synthesis method has long been essential for NGS but emerging technologies like Element Biosciences' sequencing by avidity (AVITI) represent a novel approach. It has been reported that AVITI offers improved signal-to-noise ratios and cost reductions. However, the method relies on rolling circle amplification which can be impacted by polymer size, raising questions about its efficacy sequencing small RNAs (sRNA) molecules including microRNAs (miRNAs), piwi-interacting RNAs (piRNAs), and others that are crucial regulators of gene expression and involved in various biological processes. In addition, capturing capped small RNAs (csRNA-seq) has emerged as a powerful method to map active or "nascent" RNA polymerase II transcription initiation in tissues and clinical samples. Here, we report a new protocol for seamlessly sequencing short DNA fragments on the AVITI and demonstrate that AVITI and Illumina sequencing technologies equivalently capture human, cattle (Bos taurus) and the bison (Bison bison) sRNA or csRNA sequencing libraries, augmenting the confidence in both approaches. Additionally, analysis of generated nascent transcription start sites (TSSs) data for cattle and bison revealed inaccuracies in their current genome annotations and highlighted the possibility and need to translate small RNA sequencing methodologies to livestock. Our accelerated and optimized protocol therefore bridges the advantages of AVITI sequencing and critical methods that rely on sequencing short DNA fragments.

ArCH: improving the performance of clonal hematopoiesis variant calling and interpretation

MOTIVATION

The acquisition of somatic mutations in hematopoietic stem and progenitor stem cells with resultant clonal expansion, termed clonal hematopoiesis (CH), is associated with increased risk of hematologic malignancies and other adverse outcomes. CH is generally present at low allelic fractions, but clonal expansion and acquisition of additional mutations leads to hematologic cancers in a small proportion of individuals. With high depth and high sensitivity sequencing, CH can be detected in most adults and its clonal trajectory mapped over time. However, accurate CH variant calling is challenging due to the difficulty in distinguishing low frequency CH mutations from sequencing artifacts. The lack of well-validated bioinformatic pipelines for CH calling may contribute to lack of reproducibility in studies of CH.

RESULTS

Here, we developed ArCH, an Artifact filtering Clonal Hematopoiesis variant calling pipeline for detecting single nucleotide variants and short insertions/deletions by combining the output of four variant calling tools and filtering based on variant characteristics and sequencing error rate estimation. ArCH is an end-to-end cloud-based pipeline optimized to accept a variety of inputs with customizable parameters adaptable to multiple sequencing technologies, research questions, and datasets. Using deep targeted sequencing data generated from six acute myeloid leukemia patient tumor: normal dilutions, 31 blood samples with orthogonal validation, and 26 blood samples with technical replicates, we show that ArCH improves the sensitivity and positive predictive value of CH variant detection at low allele frequencies compared to standard application of commonly used variant calling approaches.

AVAILABILITY AND IMPLEMENTATION

The code for this workflow is available at: https://github.com/kbolton-lab/ArCH.

Nanopore DNA sequencing technologies and their applications towards single-molecule proteomics

Sequencing of nucleic acids with nanopores has emerged as a powerful tool offering rapid readout, high accuracy, low cost and portability. This label-free method for sequencing at the single-molecule level is an achievement on its own. However, nanopores also show promise for the technologically even more challenging sequencing of polypeptides, something that could considerably benefit biological discovery, clinical diagnostics and homeland security, as current techniques lack portability and speed. Here we survey the biochemical innovations underpinning commercial and academic nanopore DNA/RNA sequencing techniques, and explore how these advances can fuel developments in future protein sequencing with nanopores.