Sequencing error profiles of Illumina sequencing instruments

"Outlaw" mutations in quasispecies of SARS-CoV-2 inhibit replication

The evolution of SARS-CoV-2, the agent of COVID-19, has been remarkable for its high mutation potential, leading to the appearance of variants. Some mutations have never appeared in the published genomes, which represent consensus, or bona fide genomes. Here we tested the hypothesis that mutations that did not appear in consensus genomes were, in fact, as frequent as the mutations that appeared during the various epidemic episodes, but were not expressed because lethal. To identify these mutations, we analysed the genomes of 90 nasopharyngeal samples and the quasispecies determined by next-generation sequencing. Mutations observed in the quasispecies and not in the consensus genomes were considered to be lethal, what we called "outlaw" mutations. Among these mutations, we analysed the 21 most frequent. Eight of these "outlaws" were in the RNA polymerase and we were able to use a structural biology model and molecular dynamics simulations to demonstrate the functional incapacity of these mutated RNA polymerases. Three other mutations affected the spike, a major protein involved in the pathogenesis of COVID-19. Overall, by analysing the SARS-CoV-2 quasispecies obtained during sequencing, this method made it possible to identify "outlaws," showing areas that could potentially become the target of treatments.

Genome sequencing analysis and validation of infestation-related functional genes of Vibrio parahaemolyticus LG2206 isolated from the hepatopancreas of diseased mud crab (Scylla paramamosain) in South China

Vibrio parahaemolyticus (V. parahaemolyticus) is a major bacterial pathogen found in brackish environments, leading to disease outbreaks and great economic losses in the mud crab industry. This study investigated the molecular mechanism of V. parahaemolyticus infecting mud crabs through genome sequencing analysis, survival experiments, and the expression patterns of related functional genes. A strain of V. parahaemolyticus with high pathogenicity and lethality was isolated from diseased mud crab in South China. The genome sequencing results showed that the genome size of V. parahaemolyticus was a circular chromosome of 3,357,271 bp, with a GC content of 45 %, containing 2985 protein-coding genes, denoted as V. parahaemolyticus LG2206. Genome analysis data revealed that a total of 113 adherence coding genes were obtained, including 120 virulence factor coding genes, 37 type III secretion system (T3SS) coding genes, and 277 sequences of T3SS effectors. Survival experiments showed that the mortality was 20 % within 96 h in the 1 × 104 CFU/mL infection group, 90 % in the 3.2 × 105 CFU/mL treatment group, and 100 % in the 1 × 106 CFU/mL treatment group. The LD50 of V. parahaemolyticus LG2206 was determined as 4.6 × 104 CFU/mL. Six genes of znuA and fliD (flagellin encoding genes), yscE and yscR (T3SS encoding genes), and nfuA and htpX (virulence factor encoding genes) were selected and validated by quantitative real-time PCR analysis after infection with 4.6 × 104 CFU/mL of V. parahaemolyticus LG2206 for 96 h. The expression of the six genes exhibited a significant up-regulation trend at all tested time points. The results indicated that the infestation-related genes screened in the experiment play important roles in the infestation process. This study provides timely and effective information to further analyze the molecular mechanism of V. parahaemolyticus infection and develop comprehensive measures for disease prevention and control.

Somatic mutations in aging and disease

Time always leaves its mark, and our genome is no exception. Mutations in the genome of somatic cells were first hypothesized to be the cause of aging in the 1950s, shortly after the molecular structure of DNA had been described. Somatic mutation theories of aging are based on the fact that mutations in DNA as the ultimate template for all cellular functions are irreversible. However, it took until the 1990s to develop the methods to test if DNA mutations accumulate with age in different organs and tissues and estimate the severity of the problem. By now, numerous studies have documented the accumulation of somatic mutations with age in normal cells and tissues of mice, humans, and other animals, showing clock-like mutational signatures that provide information on the underlying causes of the mutations. In this review, we will first briefly discuss the recent advances in next-generation sequencing that now allow quantitative analysis of somatic mutations. Second, we will provide evidence that the mutation rate differs between cell types, with a focus on differences between germline and somatic mutation rate. Third, we will discuss somatic mutational signatures as measures of aging, environmental exposure, and activities of DNA repair processes. Fourth, we will explain the concept of clonally amplified somatic mutations, with a focus on clonal hematopoiesis. Fifth, we will briefly discuss somatic mutations in the transcriptome and in our other genome, i.e., the genome of mitochondria. We will end with a brief discussion of a possible causal contribution of somatic mutations to the aging process.

Development of tools for the genetic manipulation of Campylobacter and their application to the N-glycosylation system of Campylobacter hepaticus, an emerging pathogen of poultry

Various species of campylobacters cause significant disease problems in both humans and animals. The continuing development of tools and methods for genetic and molecular manipulation of campylobacters enables the detailed study of bacterial virulence and disease pathogenesis. Campylobacter hepaticus is an emerging pathogen that causes spotty liver disease (SLD) in poultry. SLD has a significant economic and animal welfare impact as the disease results in elevated mortalities and significant decreases in egg production. Although potential virulence genes of C. hepaticus have been identified, they have not been further studied and characterized, as appropriate genetic tools and methods to transform and perform mutagenesis studies in C. hepaticus have not been available. In this study, the genetic manipulation of C. hepaticus is reported, with the development of novel plasmid vectors, methods for transformation, site-specific mutagenesis, and mutant complementation. These tools were used to delete the pglB gene, an oligosaccharyltransferase, a central enzyme of the N-glycosylation pathway, by allelic exchange. In the mutant strain, N-glycosylation was completely abolished. The tools and methods developed in this study represent innovative approaches that can be applied to further explore important virulence factors of C. hepaticus and other closely related Campylobacter species.

IMPORTANCE

Spotty liver disease (SLD) of layer chickens, caused by infection with Campylobacter hepaticus, is a significant economic and animal welfare burden on an important food production industry. Currently, SLD is controlled using antibiotics; however, alternative intervention methods are needed due to increased concerns associated with environmental contamination with antibiotics, and the development of antimicrobial resistance in many bacterial pathogens of humans and animals. This study has developed methods that have enabled the genetic manipulation of C. hepaticus. To validate the methods, the pglB gene was inactivated by allelic exchange to produce a C. hepaticus strain that could no longer N-glycosylate proteins. Subsequently, the mutation was complemented by reintroduction of the gene in trans, on a plasmid vector, to demonstrate that the phenotypic changes noted were caused by the mutation of the targeted gene. The tools developed enable ongoing studies to understand other virulence mechanisms of this important emerging pathogen.

A Novel Next-Generation Sequencing Assay for the Identification of BCR::ABL1 Transcript Type and Accurate and Sensitive Detection of TKI-Resistant Mutations

BACKGROUND

The clinical management of chronic myeloid leukemia (CML) patients requires the identification of the type of BCR::ABL1 transcript at diagnosis and the monitoring of its expression and potential tyrosine kinase inhibitor (TKI) resistance mutations during treatment. Detection of resistant mutation requires transcript type-specific amplification of BCR::ABL1 from RNA.

METHODS

In this study, a custom RNA-based next-generation sequencing (NGS) assay (Dup-Seq BCR::ABL1) that enables (a) the identification of BCR::ABL1 transcript type and (b) the detection of resistance mutations from common and atypical BCR::ABL1 transcript types was developed and validated. The assay design covers BCR exon 1 to ABL1 exon 10 and employs duplicate PCR amplification for error correction. The custom data analysis pipeline enables breakpoint determination and overlapped mutation calling from duplicates, which minimizes the low-level mutation artifacts.

RESULTS

This study demonstrates that this novel assay achieves high accuracy (positive percent agreement (PPA) for fusion: 98.5%; PPA and negative percent agreement (NPA) for mutation at 97.8% and 100.0%, respectively) and sensitivity (limit of detection (LOD) for mutation detection at 3% from 10 000 copies of BCR::ABL1 input).

CONCLUSIONS

The Dup-Seq BCR::ABL1 assay not only allows for the identification of BCR::ABL1 typical and atypical transcript types and accurate and sensitive detection of TKI-resistant mutations but also simplifies molecular testing work flow for the clinical management of CML patients.

Sequencing by binding rivals SMOR error-corrected sequencing by synthesis technology for accurate detection and quantification of minor (< 0.1%) subpopulation variants

BACKGROUND

Detecting very minor (< 1%) subpopulations using next-generation sequencing is a critical need for multiple applications, including the detection of drug resistant pathogens and somatic variant detection in oncology. A recently available sequencing approach termed 'sequencing by binding (SBB)' claims to have higher base calling accuracy data "out of the box." This paper evaluates the utility of using SBB for the detection of ultra-rare drug resistant subpopulations in Mycobacterium tuberculosis (Mtb) using a targeted amplicon assay and compares the performance of SBB to single molecule overlapping reads (SMOR) error corrected sequencing by synthesis (SBS) data.

RESULTS

SBS displayed an elevated error rate when compared to SMOR error-corrected SBS and SBB techniques. SMOR error-corrected SBS and SBB technologies performed similarly within the linear range studies and error rate studies.

CONCLUSIONS

With lower sequencing error rates within SBB sequencing, this technique looks promising for both targeted and unbiased whole genome sequencing, leading to the identification of minor (< 1%) subpopulations without the need for error correction methods.

CLEMENT: genomic decomposition and reconstruction of non-tumor subclones

Genome-level clonal decomposition of a single specimen has been widely studied; however, it is mostly limited to cancer research. In this study, we developed a new algorithm CLEMENT, which conducts accurate decomposition and reconstruction of multiple subclones in genome sequencing of non-tumor (normal) samples. CLEMENT employs the Expectation-Maximization (EM) algorithm with optimization strategies specific to non-tumor subclones, including false variant call identification, non-disparate clone fuzzy clustering, and clonal allele fraction confinement. In the simulation and in vitro cell line mixture data, CLEMENT outperformed current cancer decomposition algorithms in estimating the number of clones (root-mean-square-error = 0.58-0.78 versus 1.43-3.34) and in the variant-clone membership agreement (∼85.5% versus 70.1-76.7%). Additional testing on human multi-clonal normal tissue sequencing confirmed the accurate identification of subclones that originated from different cell types. Clone-level analysis, including mutational burden and signatures, provided a new understanding of normal-tissue composition. We expect that CLEMENT will serve as a crucial tool in the currently emerging field of non-tumor genome analysis.

Utility analyses of AVITI sequencing chemistry

BACKGROUND

DNA sequencing is a critical tool in modern biology. Over the last two decades, it has been revolutionized by the advent of massively parallel sequencing, leading to significant advances in the genome and transcriptome sequencing of various organisms. Nevertheless, challenges with accuracy, lack of competitive options and prohibitive costs associated with high throughput parallel short-read sequencing persist.

RESULTS

Here, we conduct a comparative analysis using matched DNA and RNA short-reads assays between Element Biosciences' AVITI and Illumina's NextSeq 550 chemistries. Similar comparisons were evaluated for synthetic long-read sequencing for RNA and targeted single-cell transcripts between the AVITI and Illumina's NovaSeq 6000. For both DNA and RNA short-read applications, the study found that the AVITI produced significantly higher per sequence quality scores. For PCR-free DNA libraries, we observed an average 89.7% lower experimentally determined error rate when using the AVITI chemistry, compared to the NextSeq 550. For short-read RNA quantification, AVITI platform had an average of 32.5% lower error rate than that for NextSeq 550. With regards to synthetic long-read mRNA and targeted synthetic long read single cell mRNA sequencing, both platforms' respective chemistries performed comparably in quantification of genes and isoforms. The AVITI displayed a marginally lower error rate for long reads, with fewer chemistry-specific errors and a higher mutation detection rate.

CONCLUSION

These results point to the potential of the AVITI platform as a competitive candidate in high-throughput short read sequencing analyses when juxtaposed with the Illumina NextSeq 550.

Limitations in next-generation sequencing-based genotyping of breast cancer polygenic risk score loci

Considering polygenic risk scores (PRSs) in individual risk prediction is increasingly implemented in genetic testing for hereditary breast cancer (BC) based on next-generation sequencing (NGS). To calculate individual BC risks, the Breast and Ovarian Analysis of Disease Incidence and Carrier Estimation Algorithm (BOADICEA) with the inclusion of the BCAC 313 or the BRIDGES 306 BC PRS is commonly used. The PRS calculation depends on accurately reproducing the variant allele frequencies (AFs) and, consequently, the distribution of PRS values anticipated by the algorithm. Here, the 324 loci of the BCAC 313 and the BRIDGES 306 BC PRS were examined in population-specific database gnomAD and in real-world data sets of five centers of the German Consortium for Hereditary Breast and Ovarian Cancer (GC-HBOC), to determine whether these expected AFs can be reproduced by NGS-based genotyping. Four PRS loci were non-existent in gnomAD v3.1.2 non-Finnish Europeans, further 24 loci showed noticeably deviating AFs. In real-world data, between 11 and 23 loci were reported with noticeably deviating AFs, and were shown to have effects on final risk prediction. Deviations depended on the sequencing approach, variant caller and calling mode (forced versus unforced) employed. Therefore, this study demonstrates the necessity to apply quality assurance not only in terms of sequencing coverage but also observed AFs in a sufficiently large cohort, when implementing PRSs in a routine diagnostic setting. Furthermore, future PRS design should be guided by the technical reproducibility of expected AFs across commonly used genotyping methods, especially NGS, in addition to the observed effect sizes.

Scalable assessment of genome editing off-targets associated with genetic variants

Genome editing with RNA-guided DNA binding factors carries risk of off-target editing at homologous sequences. Genetic variants may introduce sequence changes that increase homology to a genome editing target, thereby increasing risk of off-target editing. Conventional methods to verify candidate off-targets rely on access to cells with genomic DNA carrying these sequences. However, for candidate off-targets associated with genetic variants, appropriate cells for experimental verification may not be available. Here we develop a method, Assessment By Stand-in Off-target LentiViral Ensemble with sequencing (ABSOLVE-seq), to integrate a set of candidate off-target sequences along with unique molecular identifiers (UMIs) in genomes of primary cells followed by clinically relevant gene editor delivery. Gene editing of dozens of candidate off-target sequences may be evaluated in a single experiment with high sensitivity, precision, and power. We provide an open-source pipeline to analyze sequencing data. This approach enables experimental assessment of the influence of human genetic diversity on specificity evaluation during gene editing therapy development.

The Advantage of Targeted Next-Generation Sequencing over qPCR in Testing for Druggable EGFR Variants in Non-Small-Cell Lung Cancer

The emergence of targeted therapies in non-small-cell lung cancer (NSCLC), including inhibitors of epidermal growth factor receptor (EGFR) tyrosine kinase, has increased the need for robust companion diagnostic tests. Nowadays, detection of actionable variants in exons 18-21 of the EGFR gene by qPCR and direct DNA sequencing is often replaced by next-generation sequencing (NGS). In this study, we evaluated the diagnostic usefulness of targeted NGS for druggable EGFR variants testing in clinical NSCLC material previously analyzed by the IVD-certified qPCR test with respect to DNA reference material. We tested 59 NSCLC tissue and cytology specimens for EGFR variants using the NGS 'TruSight Tumor 15' assay (Illumina) and the qPCR 'cobas EGFR mutation test v2' (Roche Diagnostics). The sensitivity and specificity of targeted NGS assay were evaluated using the biosynthetic and biological DNA reference material with known allelic frequencies (VAF) of EGFR variants. NGS demonstrated a sufficient lower detection limit for diagnostic applications (VAF < 5%) in DNA reference material; all EGFR variants were correctly identified. NGS showed high repeatability of VAF assessment between runs (CV% from 0.02 to 3.98). In clinical material, the overall concordance between NGS and qPCR was 76.14% (Cohen's Kappa = 0.5933). The majority of discordant results concerned false-positive detection of EGFR exon 20 insertions by qPCR. A total of 9 out of 59 (15%) clinical samples showed discordant results for one or more EGFR variants in both assays. Additionally, we observed TP53 to be a frequently co-mutated gene in EGFR-positive NSCLC patients. In conclusion, targeted NGS showed a number of superior features over qPCR in EGFR variant detection (exact identification of variants, calculation of allelic frequency, high analytical sensitivity), which might enhance the basic diagnostic report.

SARS-CoV-2 population dynamics in immunocompetent individuals in a closed transmission chain shows genomic diversity over the course of infection

BACKGROUND

SARS-CoV-2 remains rapidly evolving, and many biologically important genomic substitutions/indels have characterised novel SARS-CoV-2 lineages, which have emerged during successive global waves of the pandemic. Worldwide genomic sequencing has been able to monitor these waves, track transmission clusters, and examine viral evolution in real time to help inform healthcare policy. One school of thought is that an apparent greater than average divergence in an emerging lineage from contemporary variants may require persistent infection, for example in an immunocompromised host. Due to the nature of the COVID-19 pandemic and sampling, there were few studies that examined the evolutionary trajectory of SARS-CoV-2 in healthy individuals.

METHODS

We investigated viral evolutionary trends and participant symptomatology within a cluster of 16 SARS-CoV-2 infected, immunocompetent individuals with no co-morbidities in a closed transmission chain. Longitudinal nasopharyngeal swab sampling allowed characterisation of SARS-CoV-2 intra-host variation over time at both the dominant and minor genomic variant levels through Nimagen-Illumina sequencing.

RESULTS

A change in viral lineage assignment was observed in individual infections; however, there was only one indel and no evidence of recombination over the period of an acute infection. Minor and dominant genomic modifications varied between participants, with some minor genomic modifications increasing in abundance to become the dominant viral sequence during infection.

CONCLUSIONS

Data from this cohort of SARS-CoV-2-infected participants demonstrated that long-term persistent infection in an immunocompromised host was not necessarily a prerequisite for generating a greater than average frequency of amino acid substitutions. Amino acid substitutions at both the dominant and minor genomic sequence level were observed in immunocompetent individuals during infection showing that viral lineage changes can occur generating viral diversity.

Benchmarking short and long read polishing tools for nanopore assemblies: achieving near-perfect genomes for outbreak isolates

BACKGROUND

Oxford Nanopore provides high throughput sequencing platforms able to reconstruct complete bacterial genomes with 99.95% accuracy. However, even small levels of error can obscure the phylogenetic relationships between closely related isolates. Polishing tools have been developed to correct these errors, but it is uncertain if they obtain the accuracy needed for the high-resolution source tracking of foodborne illness outbreaks.

RESULTS

We tested 132 combinations of assembly and short- and long-read polishing tools to assess their accuracy for reconstructing the genome sequences of 15 highly similar Salmonella enterica serovar Newport isolates from a 2020 onion outbreak. While long-read polishing alone improved accuracy, near perfect accuracy (99.9999% accuracy or ~ 5 nucleotide errors across the 4.8 Mbp genome, excluding low confidence regions) was only obtained by pipelines that combined both long- and short-read polishing tools. Notably, medaka was a more accurate and efficient long-read polisher than Racon. Among short-read polishers, NextPolish showed the highest accuracy, but Pilon, Polypolish, and POLCA performed similarly. Among the 5 best performing pipelines, polishing with medaka followed by NextPolish was the most common combination. Importantly, the order of polishing tools mattered i.e., using less accurate tools after more accurate ones introduced errors. Indels in homopolymers and repetitive regions, where the short reads could not be uniquely mapped, remained the most challenging errors to correct.

CONCLUSIONS

Short reads are still needed to correct errors in nanopore sequenced assemblies to obtain the accuracy required for source tracking investigations. Our granular assessment of the performance of the polishing pipelines allowed us to suggest best practices for tool users and areas for improvement for tool developers.

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.

Tracking-seq reveals the heterogeneity of off-target effects in CRISPR-Cas9-mediated genome editing

The continued development of novel genome editors calls for a universal method to analyze their off-target effects. Here we describe a versatile method, called Tracking-seq, for in situ identification of off-target effects that is broadly applicable to common genome-editing tools, including Cas9, base editors and prime editors. Through tracking replication protein A (RPA)-bound single-stranded DNA followed by strand-specific library construction, Tracking-seq requires a low cell input and is suitable for in vitro, ex vivo and in vivo genome editing, providing a sensitive and practical genome-wide approach for off-target detection in various scenarios. We show, using the same guide RNA, that Tracking-seq detects heterogeneity in off-target effects between different editor modalities and between different cell types, underscoring the necessity of direct measurement in the original system.

Sequencing of Kaposi's Sarcoma Herpesvirus (KSHV) genomes from persons of diverse ethnicities and provenances with KSHV-associated diseases demonstrate multiple infections, novel polymorphisms, and low intra-host variance

Recently published near full-length KSHV genomes from a Cameroon Kaposi sarcoma case-control study showed strong evidence of viral recombination and mixed infections, but no sequence variations associated with disease. Using the same methodology, an additional 102 KSHV genomes from 76 individuals with KSHV-associated diseases have been sequenced. Diagnoses comprise all KSHV-associated diseases (KAD): Kaposi sarcoma (KS), primary effusion lymphoma (PEL), KSHV-associated large cell lymphoma (KSHV-LCL), a type of multicentric Castleman disease (KSHV-MCD), and KSHV inflammatory cytokine syndrome (KICS). Participants originated from 22 different countries, providing the opportunity to obtain new near full-length sequences of a wide diversity of KSHV genomes. These include near full-length sequence of genomes with KSHV K1 subtypes A, B, C, and F as well as subtype E, for which no full sequence was previously available. High levels of recombination were observed. Fourteen individuals (18%) showed evidence of infection with multiple KSHV variants (from two to four unique genomes). Twenty-six comparisons of sequences, obtained from various sampling sites including PBMC, tissue biopsies, oral fluids, and effusions in the same participants, identified near complete genome conservation between different biological compartments. Polymorphisms were identified in coding and non-coding regions, including indels in the K3 and K15 genes and sequence inversions here reported for the first time. One such polymorphism in KSHV ORF46, specific to the KSHV K1 subtype E2, encoded a mutation in the leucine loop extension of the uracil DNA glycosylase that results in alteration of biochemical functions of this protein. This confirms that KSHV sequence variations can have functional consequences warranting further investigation. This study represents the largest and most diverse analysis of KSHV genome sequences to date among individuals with KAD and provides important new information on global KSHV genomics.

nf-core/airrflow: An adaptive immune receptor repertoire analysis workflow employing the Immcantation framework

Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) is a valuable experimental tool to study the immune state in health and following immune challenges such as infectious diseases, (auto)immune diseases, and cancer. Several tools have been developed to reconstruct B cell and T cell receptor sequences from AIRR-seq data and infer B and T cell clonal relationships. However, currently available tools offer limited parallelization across samples, scalability or portability to high-performance computing infrastructures. To address this need, we developed nf-core/airrflow, an end-to-end bulk and single-cell AIRR-seq processing workflow which integrates the Immcantation Framework following BCR and TCR sequencing data analysis best practices. The Immcantation Framework is a comprehensive toolset, which allows the processing of bulk and single-cell AIRR-seq data from raw read processing to clonal inference. nf-core/airrflow is written in Nextflow and is part of the nf-core project, which collects community contributed and curated Nextflow workflows for a wide variety of analysis tasks. We assessed the performance of nf-core/airrflow on simulated sequencing data with sequencing errors and show example results with real datasets. To demonstrate the applicability of nf-core/airrflow to the high-throughput processing of large AIRR-seq datasets, we validated and extended previously reported findings of convergent antibody responses to SARS-CoV-2 by analyzing 97 COVID-19 infected individuals and 99 healthy controls, including a mixture of bulk and single-cell sequencing datasets. Using this dataset, we extended the convergence findings to 20 additional subjects, highlighting the applicability of nf-core/airrflow to validate findings in small in-house cohorts with reanalysis of large publicly available AIRR datasets.

Comparative analysis of bottom trawl and nanopore sequencing in fish biodiversity assessment: The sylt outer reef example

The assessment of fish diversity is crucial for effective conservation and management strategies, especially in ecologically sensitive regions such as marine protected areas. This study contrasts the effectiveness of environmental DNA (eDNA) metabarcoding analysis employing Nanopore technology with compare beam trawl surveys at the Sylt Outer Reef, a Natura 2000 site in the North Sea, Germany. Out of the 17 fish species caught in a bottom trawl (using a 3m beam trawl), 14 were also identified through eDNA extracted from water samples. The three species not detected in the eDNA results were absent because they lacked representation in public DNA databases. The eDNA method detected twice as many fish species as the beam trawl, totalling 36 species, of which 14 were also detected by the trawl. Additionally, the selection of primers (Mifish) facilitated the identification of one marine mammal species, the harbour porpoise. In conclusion, the findings underscore the potential of eDNA coupled with MinION sequencing (Long read technology) as a robust tool for biodiversity assessment, surpassing traditional methods in detecting species richness.

Next-generation data filtering in the genomics era

Genomic data are ubiquitous across disciplines, from agriculture to biodiversity, ecology, evolution and human health. However, these datasets often contain noise or errors and are missing information that can affect the accuracy and reliability of subsequent computational analyses and conclusions. A key step in genomic data analysis is filtering - removing sequencing bases, reads, genetic variants and/or individuals from a dataset - to improve data quality for downstream analyses. Researchers are confronted with a multitude of choices when filtering genomic data; they must choose which filters to apply and select appropriate thresholds. To help usher in the next generation of genomic data filtering, we review and suggest best practices to improve the implementation, reproducibility and reporting standards for filter types and thresholds commonly applied to genomic datasets. We focus mainly on filters for minor allele frequency, missing data per individual or per locus, linkage disequilibrium and Hardy-Weinberg deviations. Using simulated and empirical datasets, we illustrate the large effects of different filtering thresholds on common population genetics statistics, such as Tajima's D value, population differentiation (FST), nucleotide diversity (π) and effective population size (Ne).

mtDNA analysis using Mitopore

Mitochondrial DNA (mtDNA) analysis is crucial for the diagnosis of mitochondrial disorders, forensic investigations, and basic research. Existing pipelines are complex, expensive, and require specialized personnel. In many cases, including the diagnosis of detrimental single nucleotide variants (SNVs), mtDNA analysis is still carried out using Sanger sequencing. Here, we developed a simple workflow and a publicly available webserver named Mitopore that allows the detection of mtDNA SNVs, indels, and haplogroups. To simplify mtDNA analysis, we tailored our workflow to process noisy long-read sequencing data for mtDNA analysis, focusing on sequence alignment and parameter optimization. We implemented Mitopore with eliBQ (eliminate bad quality reads), an innovative quality enhancement that permits the increase of per-base quality of over 20% for low-quality data. The whole Mitopore workflow and webserver were validated using patient-derived and induced pluripotent stem cells harboring mtDNA mutations. Mitopore streamlines mtDNA analysis as an easy-to-use fast, reliable, and cost-effective analysis method for both long- and short-read sequencing data. This significantly enhances the accessibility of mtDNA analysis and reduces the cost per sample, contributing to the progress of mtDNA-related research and diagnosis.

High heteroplasmy is associated with low mitochondrial copy number and selection against non-synonymous mutations in the snail Cepaea nemoralis

Molluscan mitochondrial genomes are unusual because they show wide variation in size, radical genome rearrangements and frequently show high variation (> 10%) within species. As progress in understanding this variation has been limited, we used whole genome sequencing of a six-generation matriline of the terrestrial snail Cepaea nemoralis, as well as whole genome sequences from wild-collected C. nemoralis, the sister species C. hortensis, and multiple other snail species to explore the origins of mitochondrial DNA (mtDNA) variation. The main finding is that a high rate of SNP heteroplasmy in somatic tissue was negatively correlated with mtDNA copy number in both Cepaea species. In individuals with under ten mtDNA copies per nuclear genome, more than 10% of all positions were heteroplasmic, with evidence for transmission of this heteroplasmy through the germline. Further analyses showed evidence for purifying selection acting on non-synonymous mutations, even at low frequency of the rare allele, especially in cytochrome oxidase subunit 1 and cytochrome b. The mtDNA of some individuals of Cepaea nemoralis contained a length heteroplasmy, including up to 12 direct repeat copies of tRNA-Val, with 24 copies in another snail, Candidula rugosiuscula, and repeats of tRNA-Thr in C. hortensis. These repeats likely arise due to error prone replication but are not correlated with mitochondrial copy number in C. nemoralis. Overall, the findings provide key insights into mechanisms of replication, mutation and evolution in molluscan mtDNA, and so will inform wider studies on the biology and evolution of mtDNA across animal phyla.

Manufacturing DNA in E. coli yields higher-fidelity DNA than in vitro enzymatic synthesis

Biotechnologies such as gene therapy have brought DNA vectors to the forefront of pharmaceuticals. The quality of starting material plays a pivotal role in determining final product quality. Here, we examined the fidelity of DNA replication using enzymatic methods (in vitro) compared to plasmid DNA produced in vivo in E. coli. Next-generation sequencing approaches rely on in vitro polymerases, which have inherent limitations in sensitivity. To address this challenge, we introduce a novel assay based on loss-of-function (LOF) mutations in the conditionally toxic sacB gene. Our findings show that DNA production in E. coli results in significantly fewer LOF mutations (80- to 3,000-fold less) compared to enzymatic DNA replication methods such as polymerase chain reaction (PCR) and rolling circle amplification (RCA). These results suggest that using DNA produced by PCR or RCA may introduce a substantial number of mutation impurities, potentially affecting the quality and yield of final pharmaceutical products. Our study underscores that DNA synthesized in vitro has a significantly higher mutation rate than DNA produced traditionally in E. coli. Therefore, utilizing in vitro enzymatically produced DNA in biotechnology and biomanufacturing may entail considerable fidelity-related risks, while using DNA starting material derived from E. coli substantially mitigates this risk.

OliTag-seq enhances in cellulo detection of CRISPR-Cas9 off-targets

The potential for off-target mutations is a critical concern for the therapeutic application of CRISPR-Cas9 gene editing. Current detection methodologies, such as GUIDE-seq, exhibit limitations in oligonucleotide integration efficiency and sensitivity, which could hinder their utility in clinical settings. To address these issues, we introduce OliTag-seq, an in-cellulo assay specifically engineered to enhance the detection of off-target events. OliTag-seq employs a stable oligonucleotide for precise break tagging and an innovative triple-priming amplification strategy, significantly improving the scope and accuracy of off-target site identification. This method surpasses traditional assays by providing comprehensive coverage across various sgRNAs and genomic targets. Our research particularly highlights the superior sensitivity of induced pluripotent stem cells (iPSCs) in detecting off-target mutations, advocating for using patient-derived iPSCs for refined off-target analysis in therapeutic gene editing. Furthermore, we provide evidence that prolonged Cas9 expression and transient HDAC inhibitor treatments enhance the assay's ability to uncover off-target events. OliTag-seq merges the high sensitivity typical of in vitro assays with the practical application of cellular contexts. This approach significantly improves the safety and efficacy profiles of CRISPR-Cas9 interventions in research and clinical environments, positioning it as an essential tool for the precise assessment and refinement of genome editing applications.

Benchmarking short-, long- and hybrid-read assemblers for metagenome sequencing of complex microbial communities

Metagenome community analyses, driven by the continued development in sequencing technology, is rapidly providing insights in many aspects of microbiology and becoming a cornerstone tool. Illumina, Oxford Nanopore Technologies (ONT) and Pacific Biosciences (PacBio) are the leading technologies, each with their own advantages and drawbacks. Illumina provides accurate reads at a low cost, but their length is too short to close bacterial genomes. Long reads overcome this limitation, but these technologies produce reads with lower accuracy (ONT) or with lower throughput (PacBio high-fidelity reads). In a critical first analysis step, reads are assembled to reconstruct genomes or individual genes within the community. However, to date, the performance of existing assemblers has never been challenged with a complex mock metagenome. Here, we evaluate the performance of current assemblers that use short, long or both read types on a complex mock metagenome consisting of 227 bacterial strains with varying degrees of relatedness. We show that many of the current assemblers are not suited to handle such a complex metagenome. In addition, hybrid assemblies do not fulfil their potential. We conclude that ONT reads assembled with CANU and Illumina reads assembled with SPAdes offer the best value for reconstructing genomes and individual genes of complex metagenomes, respectively.

Overriding Mendelian inheritance in Arabidopsis with a CRISPR toxin-antidote gene drive that impairs pollen germination

Synthetic gene drives, inspired by natural selfish genetic elements and transmitted to progeny at super-Mendelian (>50%) frequencies, present transformative potential for disseminating traits that benefit humans throughout wild populations, even facing potential fitness costs. Here we constructed a gene drive system in plants called CRISPR-Assisted Inheritance utilizing NPG1 (CAIN), which uses a toxin-antidote mechanism in the male germline to override Mendelian inheritance. Specifically, a guide RNA-Cas9 cassette targets the essential No Pollen Germination 1 (NPG1) gene, serving as the toxin to block pollen germination. A recoded, CRISPR-resistant copy of NPG1 serves as the antidote, providing rescue only in pollen cells that carry the drive. To limit potential consequences of inadvertent release, we used self-pollinating Arabidopsis thaliana as a model. The drive demonstrated a robust 88-99% transmission rate over two successive generations, producing minimal resistance alleles that are unlikely to inhibit drive spread. Our study provides a strong basis for rapid genetic modification or suppression of outcrossing plant populations.

MHConstructor: A high-throughput, haplotype-informed solution to the MHC assembly challenge

The extremely high levels of genetic polymorphism within the human major histocompatibility complex (MHC) limit the usefulness of reference-based alignment methods for sequence assembly. We incorporate a short read de novo assembly algorithm into a workflow for novel application to the MHC. MHConstructor is a containerized pipeline designed for high-throughput, haplotype-informed, reproducible assembly of both whole genome sequencing and target-capture short read data in large, population cohorts. To-date, no other self-contained tool exists for the generation of de novo MHC assemblies from short read data. MHConstructor facilitates wide-spread access to high quality, alignment-free MHC sequence analysis.

Variation and plasticity in life-history traits and fitness of wild Arabidopsis thaliana populations are not related to their genotypic and ecological diversity

BACKGROUND

Despite its implications for population dynamics and evolution, the relationship between genetic and phenotypic variation in wild populations remains unclear. Here, we estimated variation and plasticity in life-history traits and fitness of the annual plant Arabidopsis thaliana in two common garden experiments that differed in environmental conditions. We used up to 306 maternal inbred lines from six Iberian populations characterized by low and high genotypic (based on whole-genome sequences) and ecological (vegetation type) diversity.

RESULTS

Low and high genotypic and ecological diversity was found in edge and core Iberian environments, respectively. Given that selection is expected to be stronger in edge environments and that ecological diversity may enhance both phenotypic variation and plasticity, we expected genotypic diversity to be positively associated with phenotypic variation and plasticity. However, maternal lines, irrespective of the genotypic and ecological diversity of their population of origin, exhibited a substantial amount of phenotypic variation and plasticity for all traits. Furthermore, all populations harbored maternal lines with canalization (robustness) or sensitivity in response to harsher environmental conditions in one of the two experiments.

CONCLUSIONS

Overall, we conclude that the environmental attributes of each population probably determine their genotypic diversity, but all populations maintain substantial phenotypic variation and plasticity for all traits, which represents an asset to endure in changing environments.

Metabarcoding read abundances of orchid mycorrhizal fungi are correlated to copy numbers estimated using ddPCR

Quantifying the abundances of fungi is key to understanding natural variation in mycorrhizal communities in relation to plant ecophysiology and environmental heterogeneity. High-throughput metabarcoding approaches have transformed our ability to characterize and compare complex mycorrhizal communities. However, it remains unclear how well metabarcoding read counts correlate with actual read abundances in the sample, potentially limiting their use as a proxy for species abundances. Here, we use droplet digital PCR (ddPCR) to evaluate the reliability of ITS2 metabarcoding data for quantitative assessments of mycorrhizal communities in the orchid species Neottia ovata sampled at multiple sites. We performed specific ddPCR assays for eight families of orchid mycorrhizal fungi and compared the results with read counts obtained from metabarcoding. Our results demonstrate a significant correlation between DNA copy numbers measured by ddPCR assays and metabarcoding read counts of major mycorrhizal partners of N. ovata, highlighting the usefulness of metabarcoding for quantifying the abundance of orchid mycorrhizal fungi. Yet, the levels of correlation between the two methods and the numbers of false zero values varied across fungal families, which warrants cautious evaluation of the reliability of low-abundance families. This study underscores the potential of metabarcoding data for more quantitative analyses of mycorrhizal communities and presents practical workflows for metabarcoding and ddPCR to achieve a more comprehensive understanding of orchid mycorrhizal communities.

Template-specific optimization of NGS genotyping pipelines reveals allele-specific variation in MHC gene expression

Using high-throughput sequencing for precise genotyping of multi-locus gene families, such as the major histocompatibility complex (MHC), remains challenging, due to the complexity of the data and difficulties in distinguishing genuine from erroneous variants. Several dedicated genotyping pipelines for data from high-throughput sequencing, such as next-generation sequencing (NGS), have been developed to tackle the ensuing risk of artificially inflated diversity. Here, we thoroughly assess three such multi-locus genotyping pipelines for NGS data, the DOC method, AmpliSAS and ACACIA, using MHC class IIβ data sets of three-spined stickleback gDNA, cDNA and "artificial" plasmid samples with known allelic diversity. We show that genotyping of gDNA and plasmid samples at optimal pipeline parameters was highly accurate and reproducible across methods. However, for cDNA data, the gDNA-optimal parameter configuration yielded decreased overall genotyping precision and consistency between pipelines. Further adjustments of key clustering parameters were required tο account for higher error rates and larger variation in sequencing depth per allele, highlighting the importance of template-specific pipeline optimization for reliable genotyping of multi-locus gene families. Through accurate paired gDNA-cDNA typing and MHC-II haplotype inference, we show that MHC-II allele-specific expression levels correlate negatively with allele number across haplotypes. Lastly, sibship-assisted cDNA-typing of MHC-I revealed novel variants linked in haplotype blocks, and a higher-than-previously-reported individual MHC-I allelic diversity. In conclusion, we provide novel genotyping protocols for the three-spined stickleback MHC-I and -II genes, and evaluate the performance of popular NGS-genotyping pipelines. We also show that fine-tuned genotyping of paired gDNA-cDNA samples facilitates amplification bias-corrected MHC allele expression analysis.

Scuphr: A probabilistic framework for cell lineage tree reconstruction

Cell lineage tree reconstruction methods are developed for various tasks, such as investigating the development, differentiation, and cancer progression. Single-cell sequencing technologies enable more thorough analysis with higher resolution. We present Scuphr, a distance-based cell lineage tree reconstruction method using bulk and single-cell DNA sequencing data from healthy tissues. Common challenges of single-cell DNA sequencing, such as allelic dropouts and amplification errors, are included in Scuphr. Scuphr computes the distance between cell pairs and reconstructs the lineage tree using the neighbor-joining algorithm. With its embarrassingly parallel design, Scuphr can do faster analysis than the state-of-the-art methods while obtaining better accuracy. The method's robustness is investigated using various synthetic datasets and a biological dataset of 18 cells.

SNP-SVant: A Computational Workflow to Predict and Annotate Genomic Variants in Organisms Lacking Benchmarked Variants

Whole-genome sequencing is widely used to investigate population genomic variation in organisms of interest. Assorted tools have been independently developed to call variants from short-read sequencing data aligned to a reference genome, including single nucleotide polymorphisms (SNPs) and structural variations (SVs). We developed SNP-SVant, an integrated, flexible, and computationally efficient bioinformatic workflow that predicts high-confidence SNPs and SVs in organisms without benchmarked variants, which are traditionally used for distinguishing sequencing errors from real variants. In the absence of these benchmarked datasets, we leverage multiple rounds of statistical recalibration to increase the precision of variant prediction. The SNP-SVant workflow is flexible, with user options to tradeoff accuracy for sensitivity. The workflow predicts SNPs and small insertions and deletions using the Genome Analysis ToolKit (GATK) and predicts SVs using the Genome Rearrangement IDentification Software Suite (GRIDSS), and it culminates in variant annotation using custom scripts. A key utility of SNP-SVant is its scalability. Variant calling is a computationally expensive procedure, and thus, SNP-SVant uses a workflow management system with intermediary checkpoint steps to ensure efficient use of resources by minimizing redundant computations and omitting steps where dependent files are available. SNP-SVant also provides metrics to assess the quality of called variants and converts between VCF and aligned FASTA format outputs to ensure compatibility with downstream tools to calculate selection statistics, which are commonplace in population genomics studies. By accounting for both small and large structural variants, users of this workflow can obtain a wide-ranging view of genomic alterations in an organism of interest. Overall, this workflow advances our capabilities in assessing the functional consequences of different types of genomic alterations, ultimately improving our ability to associate genotypes with phenotypes. © 2024 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol: Predicting single nucleotide polymorphisms and structural variations Support Protocol 1: Downloading publicly available sequencing data Support Protocol 2: Visualizing variant loci using Integrated Genome Viewer Support Protocol 3: Converting between VCF and aligned FASTA formats.

Advancements in long-read genome sequencing technologies and algorithms

The recent advent of long read sequencing technologies, such as Pacific Biosciences (PacBio) and Oxford Nanopore technology (ONT), have led to substantial improvements in accuracy and computational cost in sequencing genomes. However, de novo whole-genome assembly still presents significant challenges related to the quality of the results. Pursuing de novo whole-genome assembly remains a formidable challenge, underscored by intricate considerations surrounding computational demands and result quality. As sequencing accuracy and throughput steadily advance, a continuous stream of innovative assembly tools floods the field. Navigating this dynamic landscape necessitates a reasonable choice of sequencing platform, depth, and assembly tools to orchestrate high-quality genome reconstructions. This comprehensive review delves into the intricate interplay between cutting-edge long read sequencing technologies, assembly methodologies, and the ever-evolving field of genomics. With a focus on addressing the pivotal challenges and harnessing the opportunities presented by these advancements, we provide an in-depth exploration of the crucial factors influencing the selection of optimal strategies for achieving robust and insightful genome assemblies.

Systematic mapping of TF-mediated cell fate changes by a pooled induction coupled with scRNA-seq and multi-omics approaches

Transcriptional regulation controls cellular functions through interactions between transcription factors (TFs) and their chromosomal targets. However, understanding the fate conversion potential of multiple TFs in an inducible manner remains limited. Here, we introduce iTF-seq as a method for identifying individual TFs that can alter cell fate toward specific lineages at a single-cell level. iTF-seq enables time course monitoring of transcriptome changes, and with biotinylated individual TFs, it provides a multi-omics approach to understanding the mechanisms behind TF-mediated cell fate changes. Our iTF-seq study in mouse embryonic stem cells identified multiple TFs that trigger rapid transcriptome changes indicative of differentiation within a day of induction. Moreover, cells expressing these potent TFs often show a slower cell cycle and increased cell death. Further analysis using bioChIP-seq revealed that GCM1 and OTX2 act as pioneer factors and activators by increasing gene accessibility and activating the expression of lineage specification genes during cell fate conversion. iTF-seq has utility in both mapping cell fate conversion and understanding cell fate conversion mechanisms.

Effects of parental age and polymer composition on short tandem repeat de novo mutation rates

Short tandem repeats (STRs) are hotspots of genomic variability in the human germline because of their high mutation rates, which have long been attributed largely to polymerase slippage during DNA replication. This model suggests that STR mutation rates should scale linearly with a father's age, as progenitor cells continually divide after puberty. In contrast, it suggests that STR mutation rates should not scale with a mother's age at her child's conception, since oocytes spend a mother's reproductive years arrested in meiosis II and undergo a fixed number of cell divisions that are independent of the age at ovulation. Yet, mirroring recent findings, we find that STR mutation rates covary with paternal and maternal age, implying that some STR mutations are caused by DNA damage in quiescent cells rather than polymerase slippage in replicating progenitor cells. These results echo the recent finding that DNA damage in oocytes is a significant source of de novo single nucleotide variants and corroborate evidence of STR expansion in postmitotic cells. However, we find that the maternal age effect is not confined to known hotspots of oocyte mutagenesis, nor are postzygotic mutations likely to contribute significantly. STR nucleotide composition demonstrates divergent effects on de novo mutation (DNM) rates between sexes. Unlike the paternal lineage, maternally derived DNMs at A/T STRs display a significantly greater association with maternal age than DNMs at G/C-containing STRs. These observations may suggest the mechanism and developmental timing of certain STR mutations and contradict prior attribution of replication slippage as the primary mechanism of STR mutagenesis.

Leveraging new methods for comprehensive characterization of mitochondrial DNA in esophageal squamous cell carcinoma

BACKGROUND

Mitochondria play essential roles in tumorigenesis; however, little is known about the contribution of mitochondrial DNA (mtDNA) to esophageal squamous cell carcinoma (ESCC). Whole-genome sequencing (WGS) is by far the most efficient technology to fully characterize the molecular features of mtDNA; however, due to the high redundancy and heterogeneity of mtDNA in regular WGS data, methods for mtDNA analysis are far from satisfactory.

METHODS

Here, we developed a likelihood-based method dMTLV to identify low-heteroplasmic mtDNA variants. In addition, we described fNUMT, which can simultaneously detect non-reference nuclear sequences of mitochondrial origin (non-ref NUMTs) and their derived artifacts. Using these new methods, we explored the contribution of mtDNA to ESCC utilizing the multi-omics data of 663 paired tumor-normal samples.

RESULTS

dMTLV outperformed the existing methods in sensitivity without sacrificing specificity. The verification using Nanopore long-read sequencing data showed that fNUMT has superior specificity and more accurate breakpoint identification than the current methods. Leveraging the new method, we identified a significant association between the ESCC overall survival and the ratio of mtDNA copy number of paired tumor-normal samples, which could be potentially explained by the differential expression of genes enriched in pathways related to metabolism, DNA damage repair, and cell cycle checkpoint. Additionally, we observed that the expression of CBWD1 was downregulated by the non-ref NUMTs inserted into its intron region, which might provide precursor conditions for the tumor cells to adapt to a hypoxic environment. Moreover, we identified a strong positive relationship between the number of mtDNA truncating mutations and the contribution of signatures linked to tumorigenesis and treatment response.

CONCLUSIONS

Our new frameworks promote the characterization of mtDNA features, which enables the elucidation of the landscapes and roles of mtDNA in ESCC essential for extending the current understanding of ESCC etiology. dMTLV and fNUMT are freely available from https://github.com/sunnyzxh/dMTLV and https://github.com/sunnyzxh/fNUMT , respectively.

Subclonal Cancer Driver Mutations Are Prevalent in the Unresected Peritumoral Edema of Adult Diffuse Gliomas

Adult diffuse gliomas commonly recur regardless of therapy. As recurrence typically arises from the peritumoral edema adjacent to the resected bulk tumor, the profiling of somatic mutations from infiltrative malignant cells within this critical, unresected region could provide important insights into residual disease. A key obstacle has been the inability to distinguish between next-generation sequencing (NGS) noise and the true but weak signal from tumor cells hidden among the noncancerous brain tissue of the peritumoral edema. Here, we developed and validated True2 sequencing to reduce NGS-associated errors to <1 false positive/100 kb panel positions while detecting 97.6% of somatic mutations with an allele frequency ≥0.1%. True2 was then used to study the tumor and peritumoral edema of 22 adult diffuse gliomas including glioblastoma, astrocytoma, oligodendroglioma, and NF1-related low-grade neuroglioma. The tumor and peritumoral edema displayed a similar mutation burden, indicating that surgery debulks these cancers physically but not molecularly. Moreover, variants in the peritumoral edema included unique cancer driver mutations absent in the bulk tumor. Finally, analysis of multiple samples from each patient revealed multiple subclones with unique mutations in the same gene in 17 of 22 patients, supporting the occurrence of convergent evolution in response to patient-specific selective pressures in the tumor microenvironment that may form the molecular foundation of recurrent disease. Collectively, True2 enables the detection of ultralow frequency mutations during molecular analyses of adult diffuse gliomas, which is necessary to understand cancer evolution, recurrence, and individual response to therapy.

SIGNIFICANCE

True2 is a next-generation sequencing workflow that facilitates unbiased discovery of somatic mutations across the full range of variant allele frequencies, which could help identify residual disease vulnerabilities for targeted adjuvant therapies.

Sexual transmission of urogenital bacteria: whole metagenome sequencing evidence from a sexual network study

Sexual transmission of the urogenital microbiota may contribute to adverse sexual and reproductive health outcomes. The extent of sexual transmission of the urogenital microbiota is unclear as prior studies largely investigated specific pathogens. We used epidemiologic data and whole metagenome sequencing to characterize urogenital microbiota strain concordance between participants of a sexual network study. Individuals who screened positive for genital Chlamydia trachomatis were enrolled and referred their sexual contacts from the prior 60-180 days. Snowball recruitment of sexual contacts continued for up to four waves. Vaginal swabs and penile urethral swabs were collected for whole metagenome sequencing. We evaluated bacterial strain concordance using inStrain and network analysis. We defined concordance as ≥99.99% average nucleotide identity over ≥50% shared coverage; we defined putative sexual transmission as concordance between sexual contacts with <5 single-nucleotide polymorphisms per megabase. Of 138 participants, 74 (54%) were female; 120 (87%) had genital chlamydia; and 43 (31%) were recruited contacts. We identified 115 strain-concordance events among 54 participants representing 25 bacterial species. Seven events (6%) were between sexual contacts including putative heterosexual transmission of Fannyhessea vaginae, Gardnerella leopoldii, Prevotella amnii, Sneathia sanguinegens, and Sneathia vaginalis (one strain each), and putative sexual transmission of Lactobacillus iners between female contacts. Most concordance events (108, 94%) were between non-contacts, including eight female participants connected through 18 Lactobacillus crispatus and 3 Lactobacillus jensenii concordant strains, and 14 female and 2 male participants densely interconnected through 52 Gardnerella swidsinskii concordance events.IMPORTANCEEpidemiologic evidence consistently indicates bacterial vaginosis (BV) is sexually associated and may be sexually transmitted, though sexual transmission remains subject to debate. This study is not capable of demonstrating BV sexual transmission; however, we do provide strain-level metagenomic evidence that strongly supports heterosexual transmission of BV-associated species. These findings strengthen the evidence base that supports ongoing investigations of concurrent male partner treatment for reducing BV recurrence. Our data suggest that measuring the impact of male partner treatment on F. vaginae, G. leopoldii, P. amnii, S. sanguinegens, and S. vaginalis may provide insight into why a regimen does or does not perform well. We also observed a high degree of strain concordance between non-sexual-contact female participants. We posit that this may reflect limited dispersal capacity of vaginal bacteria coupled with individuals' comembership in regional transmission networks where transmission may occur between parent and child at birth, cohabiting individuals, or sexual contacts.

CRISPR-powered quantitative keyword search engine in DNA data storage

Despite the growing interest of archiving information in synthetic DNA to confront data explosion, quantitatively querying the data stored in DNA is still a challenge. Herein, we present Search Enabled by Enzymatic Keyword Recognition (SEEKER), which utilizes CRISPR-Cas12a to rapidly generate visible fluorescence when a DNA target corresponding to the keyword of interest is present. SEEKER achieves quantitative text searching since the growth rate of fluorescence intensity is proportional to keyword frequency. Compatible with SEEKER, we develop non-collision grouping coding, which reduces the size of dictionary and enables lossless compression without disrupting the original order of texts. Using four queries, we correctly identify keywords in 40 files with a background of ~8000 irrelevant terms. Parallel searching with SEEKER can be performed on a 3D-printed microfluidic chip. Overall, SEEKER provides a quantitative approach to conducting parallel searching over the complete content stored in DNA with simple implementation and rapid result generation.

Resistance-associated mutations to the anti-SARS-CoV-2 agent nirmatrelvir: Selection not induction

Mutations associated with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) resistance to antiprotease nirmatrelvir were reported. We aimed to detect them in SARS-CoV-2 genomes and quasispecies retrieved in our institute before drug availability in January 2022 and to analyze the impact of mutations on protease (3CLpro) structure. We sought for 38 3CLpro nirmatrelvir resistance mutations in a set of 62 673 SARS-CoV-2 genomes obtained in our institute from respiratory samples collected between 2020 and 2023 and for these mutations in SARS-CoV-2 quasispecies for 90 samples collected in 2020, using Python. SARS-CoV-2 protease with major mutation E166V was generated with Swiss Pdb Viewer and Molegro Molecular Viewer. We detected 22 (58%) of the resistance-associated mutations in 417 (0.67%) of the genomes analyzed; 325 (78%) of these genomes had been obtained from samples collected in 2020-2021. APOBEC signatures were found for 12/22 mutations. We also detected among viral quasispecies from 90 samples some minority reads harboring any of 15 nirmatrelvir resistance mutations, including E166V. Also, we predicted that E166V has a very limited effect on 3CLpro structure but may prevent drug attachment. Thus, we evidenced that mutations associated with nirmatrelvir resistance pre-existed in SARS-CoV-2 before drug availability. These findings further warrant SARS-CoV-2 genomic surveillance and SARS-CoV-2 quasispecies characterization.

CLAW: An automated Snakemake workflow for the assembly of chloroplast genomes from long-read data

Chloroplasts are photosynthetic organelles in algal and plant cells that contain their own genome. Chloroplast genomes are commonly used in evolutionary studies and taxonomic identification and are increasingly becoming a target for crop improvement studies. As DNA sequencing becomes more affordable, researchers are collecting vast swathes of high-quality whole-genome sequence data from laboratory and field settings alike. Whole tissue read libraries sequenced with the primary goal of understanding the nuclear genome will inadvertently contain many reads derived from the chloroplast genome. These whole-genome, whole-tissue read libraries can additionally be used to assemble chloroplast genomes with little to no extra cost. While several tools exist that make use of short-read second generation and third-generation long-read sequencing data for chloroplast genome assembly, these tools may have complex installation steps, inadequate error reporting, poor expandability, and/or lack scalability. Here, we present CLAW (Chloroplast Long-read Assembly Workflow), an easy to install, customise, and use Snakemake tool to assemble chloroplast genomes from chloroplast long-reads found in whole-genome read libraries (https://github.com/aaronphillips7493/CLAW). Using 19 publicly available reference chloroplast genome assemblies and long-read libraries from algal, monocot and eudicot species, we show that CLAW can rapidly produce chloroplast genome assemblies with high similarity to the reference assemblies. CLAW was designed such that users have complete control over parameterisation, allowing individuals to optimise CLAW to their specific use cases. We expect that CLAW will provide researchers (with varying levels of bioinformatics expertise) with an additional resource useful for contributing to the growing number of publicly available chloroplast genome assemblies.

nf-core/airrflow: an adaptive immune receptor repertoire analysis workflow employing the Immcantation framework

Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) is a valuable experimental tool to study the immune state in health and following immune challenges such as infectious diseases, (auto)immune diseases, and cancer. Several tools have been developed to reconstruct B cell and T cell receptor sequences from AIRR-seq data and infer B and T cell clonal relationships. However, currently available tools offer limited parallelization across samples, scalability or portability to high-performance computing infrastructures. To address this need, we developed nf-core/airrflow, an end-to-end bulk and single-cell AIRR-seq processing workflow which integrates the Immcantation Framework following BCR and TCR sequencing data analysis best practices. The Immcantation Framework is a comprehensive toolset, which allows the processing of bulk and single-cell AIRR-seq data from raw read processing to clonal inference. nf-core/airrflow is written in Nextflow and is part of the nf-core project, which collects community contributed and curated Nextflow workflows for a wide variety of analysis tasks. We assessed the performance of nf-core/airrflow on simulated sequencing data with sequencing errors and show example results with real datasets. To demonstrate the applicability of nf-core/airrflow to the high-throughput processing of large AIRR-seq datasets, we validated and extended previously reported findings of convergent antibody responses to SARS-CoV-2 by analyzing 97 COVID-19 infected individuals and 99 healthy controls, including a mixture of bulk and single-cell sequencing datasets. Using this dataset, we extended the convergence findings to 20 additional subjects, highlighting the applicability of nf-core/airrflow to validate findings in small in-house cohorts with reanalysis of large publicly available AIRR datasets. nf-core/airrflow is available free of charge, under the MIT license on GitHub (https://github.com/nf-core/airrflow). Detailed documentation and example results are available on the nf-core website at (https://nf-co.re/airrflow).

Diversity of antibiotic resistance gene variants at subsequent stages of the wastewater treatment process revealed by a metagenomic analysis of PCR amplicons

Wastewater treatment plants have been recognised as point sources of various antibiotic-resistant bacteria (ARB) and antibiotic resistance genes (ARG) which are considered recently emerging biological contaminants. So far, culture-based and molecular-based methods have been successfully applied to monitor antimicrobial resistance (AMR) in WWTPs. However, the methods applied do not permit the comprehensive identification of the true diversity of ARGs. In this study we applied next-generation sequencing for a metagenomic analysis of PCR amplicons of ARGs from the subsequent stages of the analysed WWTP. The presence of 14 genes conferring resistance to different antibiotic families was screened by PCR. In the next step, three genes were selected for detailed analysis of changes of the profile of ARG variants along the process. A relative abundance of 79 variants was analysed. The highest diversity was revealed in the ermF gene, with 52 variants. The relative abundance of some variants changed along the purification process, and some ARG variants might be present in novel hosts for which they were currently unassigned. Additionally, we identified a pool of novel ARG variants present in the studied WWTP. Overall, the results obtained indicated that the applied method is sufficient for analysing ARG variant diversity.

An overlooked phenomenon: complex interactions of potential error sources on the quality of bacterial de novo genome assemblies

BACKGROUND

Parameters adversely affecting the contiguity and accuracy of the assemblies from Illumina next-generation sequencing (NGS) are well described. However, past studies generally focused on their additive effects, overlooking their potential interactions possibly exacerbating one another's effects in a multiplicative manner. To investigate whether or not they act interactively on de novo genome assembly quality, we simulated sequencing data for 13 bacterial reference genomes, with varying levels of error rate, sequencing depth, PCR and optical duplicate ratios.

RESULTS

We assessed the quality of assemblies from the simulated sequencing data with a number of contiguity and accuracy metrics, which we used to quantify both additive and multiplicative effects of the four parameters. We found that the tested parameters are engaged in complex interactions, exerting multiplicative, rather than additive, effects on assembly quality. Also, the ratio of non-repeated regions and GC% of the original genomes can shape how the four parameters affect assembly quality.

CONCLUSIONS

We provide a framework for consideration in future studies using de novo genome assembly of bacterial genomes, e.g. in choosing the optimal sequencing depth, balancing between its positive effect on contiguity and negative effect on accuracy due to its interaction with error rate. Furthermore, the properties of the genomes to be sequenced also should be taken into account, as they might influence the effects of error sources themselves.

Pangenome databases improve host removal and mycobacteria classification from clinical metagenomic data

BACKGROUND

Culture-free real-time sequencing of clinical metagenomic samples promises both rapid pathogen detection and antimicrobial resistance profiling. However, this approach introduces the risk of patient DNA leakage. To mitigate this risk, we need near-comprehensive removal of human DNA sequences at the point of sequencing, typically involving the use of resource-constrained devices. Existing benchmarks have largely focused on the use of standardized databases and largely ignored the computational requirements of depletion pipelines as well as the impact of human genome diversity.

RESULTS

We benchmarked host removal pipelines on simulated and artificial real Illumina and Nanopore metagenomic samples. We found that construction of a custom kraken database containing diverse human genomes results in the best balance of accuracy and computational resource usage. In addition, we benchmarked pipelines using kraken and minimap2 for taxonomic classification of Mycobacterium reads using standard and custom databases. With a database representative of the Mycobacterium genus, both tools obtained improved specificity and sensitivity, compared to the standard databases for classification of Mycobacterium tuberculosis. Computational efficiency of these custom databases was superior to most standard approaches, allowing them to be executed on a laptop device.

CONCLUSIONS

Customized pangenome databases provide the best balance of accuracy and computational efficiency when compared to standard databases for the task of human read removal and M. tuberculosis read classification from metagenomic samples. Such databases allow for execution on a laptop, without sacrificing accuracy, an especially important consideration in low-resource settings. We make all customized databases and pipelines freely available.

HAT: de novo variant calling for highly accurate short-read and long-read sequencing data

MOTIVATION

de novo variants (DNVs) are variants that are present in offspring but not in their parents. DNVs are both important for examining mutation rates as well as in the identification of disease-related variation. While efforts have been made to call DNVs, calling of DNVs is still challenging from parent-child sequenced trio data. We developed Hare And Tortoise (HAT) as an automated DNV detection workflow for highly accurate short-read and long-read sequencing data. Reliable detection of DNVs is important for human genomics and HAT addresses this need.

RESULTS

HAT is a computational workflow that begins with aligned read data (i.e. CRAM or BAM) from a parent-child sequenced trio and outputs DNVs. HAT detects high-quality DNVs from Illumina short-read whole-exome sequencing, Illumina short-read whole-genome sequencing, and highly accurate PacBio HiFi long-read whole-genome sequencing data. The quality of these DNVs is high based on a series of quality metrics including number of DNVs per individual, percent of DNVs at CpG sites, and percent of DNVs phased to the paternal chromosome of origin.

AVAILABILITY AND IMPLEMENTATION

https://github.com/TNTurnerLab/HAT.

The ALT pathway generates telomere fusions that can be detected in the blood of cancer patients

Telomere fusions (TFs) can trigger the accumulation of oncogenic alterations leading to malignant transformation and drug resistance. Despite their relevance in tumour evolution, our understanding of the patterns and consequences of TFs in human cancers remains limited. Here, we characterize the rates and spectrum of somatic TFs across >30 cancer types using whole-genome sequencing data. TFs are pervasive in human tumours with rates varying markedly across and within cancer types. In addition to end-to-end fusions, we find patterns of TFs that we mechanistically link to the activity of the alternative lengthening of telomeres (ALT) pathway. We show that TFs can be detected in the blood of cancer patients, which enables cancer detection with high specificity and sensitivity even for early-stage tumours and cancers of high unmet clinical need. Overall, we report a genomic footprint that enables characterization of the telomere maintenance mechanism of tumours and liquid biopsy analysis.

Sequencing by avidity enables high accuracy with low reagent consumption

We present avidity sequencing, a sequencing chemistry that separately optimizes the processes of stepping along a DNA template and that of identifying each nucleotide within the template. Nucleotide identification uses multivalent nucleotide ligands on dye-labeled cores to form polymerase-polymer-nucleotide complexes bound to clonal copies of DNA targets. These polymer-nucleotide substrates, termed avidites, decrease the required concentration of reporting nucleotides from micromolar to nanomolar and yield negligible dissociation rates. Avidity sequencing achieves high accuracy, with 96.2% and 85.4% of base calls having an average of one error per 1,000 and 10,000 base pairs, respectively. We show that the average error rate of avidity sequencing remained stable following a long homopolymer.

Effects of parental age and polymer composition on short tandem repeat de novo mutation rates

Short tandem repeats (STRs) are hotspots of genomic variability in the human germline because of their high mutation rates, which have long been attributed largely to polymerase slippage during DNA replication. This model suggests that STR mutation rates should scale linearly with a father's age, as progenitor cells continually divide after puberty. In contrast, it suggests that STR mutation rates should not scale with a mother's age at her child's conception, since oocytes spend a mother's reproductive years arrested in meiosis II and undergo a fixed number of cell divisions that are independent of the age at ovulation. Yet, mirroring recent findings, we find that STR mutation rates covary with paternal and maternal age, implying that some STR mutations are caused by DNA damage in quiescent cells rather than the classical mechanism of polymerase slippage in replicating progenitor cells. These results also echo the recent finding that DNA damage in quiescent oocytes is a significant source of de novo SNVs and corroborate evidence of STR expansion in postmitotic cells. However, we find that the maternal age effect is not confined to previously discovered hotspots of oocyte mutagenesis, nor are post-zygotic mutations likely to contribute significantly. STR nucleotide composition demonstrates divergent effects on DNM rates between sexes. Unlike the paternal lineage, maternally derived DNMs at A/T STRs display a significantly greater association with maternal age than DNMs at GC-containing STRs. These observations may suggest the mechanism and developmental timing of certain STR mutations and are especially surprising considering the prior belief in replication slippage as the dominant mechanism of STR mutagenesis.

Evaluation of Long-Read Sequencing Simulators to Assess Real-World Applications for Food Safety

Shiga toxin-producing Escherichia coli (STEC) and Listeria monocytogenes are routinely responsible for severe foodborne illnesses in the United States. Current identification methods utilized by the U.S. Food Safety Inspection Service require at least four days to identify STEC and six days for L. monocytogenes. Adoption of long-read, whole genome sequencing for food safety testing could significantly reduce the time needed for identification, but method development costs are high. Therefore, the goal of this project was to use NanoSim-H software to simulate Oxford Nanopore sequencing reads to assess the feasibility of sequencing-based foodborne pathogen detection and guide experimental design. Sequencing reads were simulated for STEC, L. monocytogenes, and a 1:1 combination of STEC and Bos taurus genomes using NanoSim-H. At least 2500 simulated reads were needed to identify the seven genes of interest targeted in STEC, and at least 500 reads were needed to detect the gene targeted in L. monocytogenes. Genome coverage of 30x was estimated at 21,521, and 11,802 reads for STEC and L. monocytogenes, respectively. Approximately 5-6% of reads simulated from both bacteria did not align with their respective reference genomes due to the introduction of errors. For the STEC and B. taurus 1:1 genome mixture, all genes of interest were detected with 1,000,000 reads, but less than 1x coverage was obtained. The results suggested sample enrichment would be necessary to detect foodborne pathogens with long-read sequencing, but this would still decrease the time needed from current methods. Additionally, simulation data will be useful for reducing the time and expense associated with laboratory experimentation.

Emergence of resistance-associated variants during sofosbuvir treatment in chronically infected hepatitis E patients

BACKGROUND AND AIMS

Chronic HEV infections remain a serious problem in immunocompromised patients, as specifically approved antiviral drugs are unavailable. In 2020, a 24-week multicenter phase II pilot trial was carried out, evaluating the nucleotide analog sofosbuvir by treating nine chronically HEV-infected patients with sofosbuvir (Trial Number NCT03282474). During the study, antiviral therapy reduced virus RNA levels initially but did not lead to a sustained virologic response. Here, we characterize the changes in HEV intrahost populations during sofosbuvir treatment to identify the emergence of treatment-associated variants.

APPROACH AND RESULTS

We performed high-throughput sequencing on RNA-dependent RNA polymerase sequences to characterize viral population dynamics in study participants. Subsequently, we used an HEV-based reporter replicon system to investigate sofosbuvir sensitivity in high-frequency variants. Most patients had heterogenous HEV populations, suggesting high adaptability to treatment-related selection pressures. We identified numerous amino acid alterations emerging during treatment and found that the EC 50 of patient-derived replicon constructs was up to ~12-fold higher than the wild-type control, suggesting that variants associated with lower drug sensitivity were selected during sofosbuvir treatment. In particular, a single amino acid substitution (A1343V) in the finger domain of ORF1 could reduce susceptibility to sofosbuvir significantly in 8 of 9 patients.

CONCLUSIONS

In conclusion, viral population dynamics played a critical role during antiviral treatment. High population diversity during sofosbuvir treatment led to the selection of variants (especially A1343V) with lower sensitivity to the drug, uncovering a novel mechanism of resistance-associated variants during sofosbuvir treatment.