Validation of Golden Gate Assemblies Using Highly Multiplexed Nanopore Amplicon Sequencing.

Golden Gate cloning has revolutionized synthetic biology. Its concept of modular, highly characterized libraries of parts that can be combined into higher order assemblies allows engineering principles to be applied to biological systems. The basic parts, typically stored in Level 0 plasmids, are sequence validated by the method of choice and can be combined into higher order assemblies on demand. Higher order assemblies are typically transcriptional units, and multiple transcriptional units can be assembled into multi-gene constructs. Higher order Golden Gate assembly based on defined and validated parts usually does not introduce sequence changes. Therefore, simple validation of the assemblies, e.g., by colony polymerase chain reaction (PCR) or restriction digest pattern analysis is sufficient. However, in many experimental setups, researchers do not use defined parts, but rather part libraries, resulting in assemblies of high combinatorial complexity where sequencing again becomes mandatory. Here, we present a detailed protocol for the use of a highly multiplexed dual barcode amplicon sequencing using the Nanopore sequencing platform for in-house sequence validation. The workflow, called DuBA.flow, is a start-to-finish procedure that provides all necessary steps from a single colony to the final easy-to-interpret sequencing report.

Promoter Deletion Leading to Allele Specific Expression in a Genetically Unsolved Case of Primary Ciliary Dyskinesia.

Variation in the non-coding genome represents an understudied mechanism of disease and it remains challenging to predict if single nucleotide variants, small insertions and deletions, or structural variants in non-coding genomic regions will be detrimental. Our approach using complementary RNA-seq and targeted long-read DNA sequencing can prioritize identification of non-coding variants that lead to disease via alteration of gene splicing or expression. We have identified a patient with primary ciliary dyskinesia with a pathogenic coding variant on one allele of the SPAG1 gene, while the second allele appears normal by whole exome sequencing despite an autosomal recessive inheritance pattern. RNA sequencing revealed reduced SPAG1 transcript levels and exclusive allele specific expression of the known pathogenic allele, suggesting the presence of a non-coding variant on the second allele that impacts transcription. Targeted long-read DNA sequencing identified a heterozygous 3 kilobase deletion of the 5' untranslated region of SPAG1, overlapping the promoter and first non-coding exon. This non-coding deletion was missed by whole exome sequencing and gene-specific deletion/duplication analysis, highlighting the importance of investigating the non-coding genome in patients with "missing" disease-causing variation. This paradigm demonstrates the utility of both RNA and long-read DNA sequencing in identifying pathogenic non-coding variants in patients with unexplained genetic disease.

Interchromosomal segmental duplication drives translocation and loss of P. falciparum histidine-rich protein 3.

Most malaria rapid diagnostic tests (RDTs) detect Plasmodium falciparum histidine-rich protein 2 (PfHRP2) and PfHRP3, but deletions of pfhrp2 and phfrp3 genes make parasites undetectable by RDTs. We analyzed 19,313 public whole-genome-sequenced P. falciparum field samples to understand these deletions better. Pfhrp2 deletion only occurred by chromosomal breakage with subsequent telomere healing. Pfhrp3 deletions involved loss from pfhrp3 to the telomere and showed three patterns: no other associated rearrangement with evidence of telomere healing at breakpoint (Asia; Pattern 13-TARE1); associated with duplication of a chromosome 5 segment containing multidrug-resistant-1 gene (Asia; Pattern 13-5++); and most commonly, associated with duplication of a chromosome 11 segment (Americas/Africa; Pattern 13-11++). We confirmed a 13-11 hybrid chromosome with long-read sequencing, consistent with a translocation product arising from recombination between large interchromosomal ribosome-containing segmental duplications. Within most 13-11++ parasites, the duplicated chromosome 11 segments were identical. Across parasites, multiple distinct haplotype groupings were consistent with emergence due to clonal expansion of progeny from intrastrain meiotic recombination. Together, these observations suggest negative selection normally removes 13-11++pfhrp3 deletions, and specific conditions are needed for their emergence and spread including low transmission, findings that can help refine surveillance strategies.

Prokaryotic communities profiling of Indonesian hot springs using long-read Oxford Nanopore sequencing.

OBJECTIVES: Indonesia's location at the convergence of multiple tectonic plates results in a unique geomorphological feature with abundant hot springs. This study pioneers the metagenomic exploration of Indonesian hot springs, harbouring unique life forms despite high temperatures. The microbial community of hot springs is taxonomically versatile and biotechnologically valuable. 16s rRNA amplicon sequencing of the metagenome is a viable option for the microbiome investigation. This study utilized Oxford Nanopore's long-read 16 S rRNA sequencing for enhanced species identification, improved detection of rare members, and a more detailed community composition profile. DATA DESCRIPTION: Water samples were taken from three hot springs of the Bali, Indonesia (i) Angseri, 8.362503 S, 115.133452 E; (ii) Banjar, 8.210270 S, 114.967063 E; and (iii) Batur, 8.228806 S, 115.404829 E. BioLit Genomic DNA Extraction Kit (SRL, Mumbai, India) was used to isolate DNA from water samples. The quantity and quality of the DNA were determined using a NanoDrop™ spectrophotometer and a Qubit fluorometer (Thermo Fisher Scientific, USA). The library was created using Oxford Nanopore Technology kits, and the sequencing was done using Oxford Nanopore's GridION platform. All sequencing data was obtained in FASTQ files and filtered using NanoFilt software. This dataset is valuable for searching novel bacteria diversity and their existence.

Description of the complete rDNA repeat unit structure of Coturnixjaponica Temminck et Schlegel, 1849 (Aves).

Ribosomal RNA (18S, 5.8S, 28S) gene clusters in genomes form regions that consist of multiple tandem repeats. They are located on a single or several pairs of chromosomes and play an important role in the formation of the nucleolus responsible for the assembly of ribosome subunits. The rRNA gene cluster sequences are widely used for taxonomic studies, however at present, complete information on the avian rDNA repeat unit structure including intergenic spacer sequence is available only for the chicken (Gallusgallusdomesticus Linnaeus, 1758). The GC enrichment and high-order repeats peculiarities within the intergenic spacer described for the chicken rDNA cluster may be responsible for these failures. The karyotype of the Japanese quail (Coturnixjaponica Temminck et Schlegel, 1849) deserves close attention because, unlike most birds, it has three pairs of nucleolar organizer bearing chromosomes, two of which are microchromosomes enriched in repeating elements and heterochromatin that carry translocated terminal nucleolar organizers. Here we assembled and annotated the complete Japanese quail ribosomal gene cluster sequence of 21166 base pairs (GenBank under the registration tag BankIt2509210 CoturnixOK523374). This is the second deciphered avian rDNA cluster after the chicken. Despite the revealed high similarity with the chicken corresponding sequence, it has a number of specific features, which include a slightly lower degree of GC content and the presence of bendable elements in the content of both the transcribed spacer I and the non-transcribed intergenic spacer.

Characterization of the Rat Osteosarcoma Cell Line UMR-106 by Long-Read Technologies Identifies a Large Block of Amplified Genes Associated with Human Disease.

BACKGROUND/OBJECTIVES: The rat osteosarcoma cell line UMR-106 is widely used for the study of bone cancer biology but it has not been well characterized with modern genomic methods. METHODS: To better understand the biology of UMR-106 cells we used a combination of optical genome mapping (OGM), long-read sequencing nanopore sequencing and RNA sequencing.The UMR-106 genome was compared to a strain-matched Sprague-Dawley rat for variants associated with human osteosarcoma while expression data were contrasted with a public osteoblast dataset. RESULTS: Using the COSMIC database to identify the most affected genes in human osteosarcomas we found somatic mutations in Tp53 and H3f3a. OGM identified a relatively small number of differences between the cell line and a strain-matched control animal but did detect a ~45 Mb block of amplification that included Myc on chromosome 7 which was confirmed by long-read sequencing. The amplified region showed several blocks of non-contiguous rearranged sequence implying complex rearrangements during their formation and included 14 genes reported as biomarkers in human osteosarcoma, many of which also showed increased transcription. A comparison of 5mC methylation from the nanopore reads of tumor and control samples identified genes with distinct differences including the OS marker Cdkn2a. CONCLUSIONS: This dataset illustrates the value of long DNA methods for the characterization of cell lines and how inter-species analysis can inform us about the genetic nature underlying mutations that underpin specific tumor types. The data should be a valuable resource for investigators studying osteosarcoma, in general, and specifically the UMR-106 model.

Revolutionising High Resolution HLA Genotyping for Transplant Assessment: Validation, Implementation and Challenges of Oxford Nanopore Technologies' Q20+ Sequencing.

The advent of third-generation sequencing (TGS) represents a significant shift in the field of genetic sequencing, enabling single-molecule sequencing to overcome limitations of short-read NGS platforms. Several studies have assessed the utilisation of TGS in HLA genotyping, though many of these studies have described the high error rate as an obstacle to achieving a robust and highly accurate HLA typing assay. In 2021, Oxford Nanopore Technologies (ONT) released the high-accuracy sequencing Kit 14 and the MinION flow cell model R10.4.1, which were reported to achieve sequencing accuracies up to 99%. The aim of this study was to validate this novel high-accuracy sequencing kit for HLA genotyping coupled with a full-gene HLA PCR assay. Comparison with historical data obtained using legacy flow cell models such as R9.4, R10.3 and R10.4 was also done to assess progressive improvement in sequencing performance with each sequential release. The workflow was validated based on data throughput, sequence quality and accuracy, and HLA genotyping resolution. An initial validation was performed using an internal reference panel of 42 samples representing common HLA allele groups, followed by an analysis of data obtained from 111 sequencing batch runs since the implementation, to assess assay performance and define quality control metrics to assess inter-run variability and monitor quality. Furthermore, challenges arising from MinION flow cell stability and use, and assessment of barcode contamination are discussed. The findings of this study highlight advantages of ONT sequencing kit 14/R10.4.1 for HLA genotyping and the implementation considerations for the routine diagnostic HLA laboratory.

Microbiome-induced reprogramming in post-transcriptional landscape using nanopore direct RNA sequencing.

It has been widely recognized that the microbiota has the capacity to shape host gene expression and physiological functions. However, there remains a paucity of comprehensive study revealing the host transcriptional landscape regulated by the microbiota. Here, we comprehensively examined mRNA landscapes in mouse tissues (brain and cecum) from specific-pathogen-free and germ-free mice using nanopore direct RNA sequencing. Our results show that the microbiome has global influence on a host's RNA modifications (m6A, m5C, Ψ), isoform generation, poly(A) tail length, and transcript abundance in both brain and cecum tissues. Moreover, the microbiome exerts tissue-specific effects on various post-transcriptional regulatory processes. In addition, the microbiome impacts the coordination of multiple RNA modifications in host brain and cecum tissues. In conclusion, we establish the relationship between microbial regulation and gene expression. Our results help the understanding of the mechanisms by which the microbiome reprograms host gene expression.

Investigating the Quantification Capabilities of a Nanopore-Based Sequencing Platform for Food Safety Application via External Standards of Lambda DNA and Lambda Spiked Beef.

Six hundred million cases of disease and roughly 420,000 deaths occur globally each year due to foodborne pathogens. Current methods to screen and identify pathogens in swine, poultry, and cattle products include immuno-based techniques (e.g., immunoassay integrated biosensors), molecular methods (e.g., DNA hybridization and PCR assays), and traditional culturing. These methods are often used in tandem to screen, quantify, and characterize samples, prolonging real-time comprehensive analysis. Next-generation sequencing (NGS) is a relatively new technology that combines DNA-sequencing chemistry and bioinformatics to generate and analyze large amounts of short- or long-read DNA sequences and whole genomes. The goal of this project was to evaluate the quantitative capabilities of the real-time NGS Oxford Nanopore Technologies' MinION sequencer through a shotgun-based sequencing approach. This investigation explored the correlation between known amounts of the analyte (lambda DNA as a pathogenic bacterial surrogate) with data output, in both the presence and absence of a background matrix (Bos taurus DNA). A positive linear correlation was observed between the concentration of analyte and the amount of data produced, number of bases sequenced, and number of reads generated in both the presence and absence of a background matrix. In the presence of bovine DNA, the sequenced data were successfully mapped to the NCBI lambda reference genome. Furthermore, the workflow from pre-extracted DNA to target identification took less than 3 h, demonstrating the potential of long-read sequencing in food safety as a rapid method for screening, identification, and quantification.

Resolving Sulfation Posttranslational Modifications on a Peptide Hormone using Nanopores.

Peptide hormones are decorated with post-translational modifications (PTMs) that are crucial for receptor recognition. Tyrosine sulfation on plant peptide hormones is, for example, essential for plant growth and development. Measuring the occurrence and position of sulfotyrosine is, however, compromised by major technical challenges during isolation and detection. Nanopores can sensitively detect protein PTMs at the single-molecule level. By translocating PTM variants of the plant pentapeptide hormone phytosulfokine (PSK) through a nanopore, we here demonstrate the accurate identification of sulfation and phosphorylation on the two tyrosine residues of PSK. Sulfation can be clearly detected and distinguished (>90%) from phosphorylation on the same residue. Moreover, the presence or absence of PTMs on the two close-by tyrosine residues can be accurately determined (>96% accuracy). Our findings demonstrate the extraordinary sensitivity of nanopore protein measurements, providing a powerful tool for identifying position-specific sulfation on peptide hormones and promising wider applications to identify protein PTMs.

Efficient, High-Quality Engineering of Therapeutic Extracellular Vesicles on an Integrated Nanoplatform.

Engineered extracellular vesicles (EVs) have been recognized as important therapeutics for gene and cell therapy. To achieve clinically desired therapy, technologies for EV engineering have high demands on the efficacy in producing EVs and their qualities, which, however, remain challenging to conventional routes due to their limited control on therapeutic payload delivery, EV secretion, and extracellular microenvironments. Here, we report a nanoplatform (denoted as PURE) that enables efficient electro-transfection while stimulating cells to produce high-quality EVs carrying functional RNAs. PURE further employs an ammonium removal zone to maintain the physiological conditions of the extracellular microenvironment and an EV uptake zone that efficiently (87.1%) captures EVs in situ with porous hydrogels. The platform achieved about a 12-fold higher yield of engineered EVs and a 146-fold abundance of desired therapeutics compared to those naturally secreted from cells. The PURE-engineered miR-130a-EVs were validated for effectively upregulating the mTOR signaling pathway in both in vitro and in vivo. Their therapeutic capability was then verified by enhancing the in vitro activation of primordial follicles. In vivo applications further highlighted the therapeutic effects of miR-130a-EVs in restoring ovary function in aged mice. The PURE platform represents a strategy for the clinical translation of EV-mediated therapy.

Anion Selectivities in Zwitterion Grafted Nanopores: Effect of Zwitterion Architecture.

The separation of ions of similar charge is a crucial challenge in many applications, from water treatment to precious metal recovery. Membranes with cross-linked zwitterionic amphiphilic copolymer (ZAC-X) selective layers, which feature self-assembled, zwitterion-lined nanodomains for permeation, offer unique permselectivity between monovalent anions (e.g., Cl-/F-). This has motivated studies on the mechanisms of transport and selectivity in this family of materials. In this study, we conducted molecular dynamics simulations of aqueous salt solutions within zwitterion-functionalized nanopores to elucidate the influence of dipole orientation of the ZI ligands on anion diffusivities, partitioning, and permeabilities. Our model compares systems with contrasting ZI organization: surface-cation-anion (S-ZI+-ZI-, Motif A) and surface-anion-cation (S-ZI--ZI+, Motif B). Our results reveal that Motif A exhibits less pronounced ion pairing due to a spatial separation in the radial profiles of cations and anions. Motif B demonstrates prominent ion pairing for smaller anions owing to their overlap with cation distributions. Further, our potential of mean force profiles reveals that anion partitioning increases with anion size in both ligand motifs, whereas Motif B exhibits significantly higher partitioning selectivity toward larger anions compared to Motif A. Our results for ion diffusivities show that the self-diffusivities of both anions and cations are lower for Motif B compared to Motif A. Such trends in anion partitioning and diffusivities can be explained by differences in the interactions and steric hindrance experienced by the anionic species in Motifs A and B. Finally, our results for anion permselectivity, obtained by combining partitioning and diffusivity, indicate that partitioning trends dominate over diffusivity trends. Consequently, anion permeability increases with anion size, and ligand Motif B yields much higher permselectivity toward larger anions compared to ligand Motif A.

TRMT1L-catalyzed m2 2G27 on tyrosine tRNA is required for efficient mRNA translation and cell survival under oxidative stress.

tRNA modifications are critical for several aspects of their functions, including decoding, folding, and stability. Using a multifaceted approach encompassing eCLIP-seq and Nanopore tRNA-seq, we show that the human tRNA methyltransferase TRMT1L interacts with component of the Rix1 ribosome biogenesis complex and binds to the 28S rRNA, as well as to a subset of tRNAs. Mechanistically, we demonstrate that TRMT1L is responsible for catalyzing m2 2G solely at position 27 of tRNA-Tyr-GUA. Surprisingly, TRMT1L depletion also impaired the deposition of acp3U and dihydrouridine on tRNA-Tyr-GUA, Cys-GCA, and Ala-CGC. TRMT1L knockout cells have a marked decrease in tRNA-Tyr-GUA levels, coinciding with a reduction in global translation rates and hypersensitivity to oxidative stress. Our results establish TRMT1L as the elusive methyltransferase catalyzing the m2 2G27 modification on tRNA Tyr, resolving a long-standing gap of knowledge and highlighting its potential role in a tRNA modification circuit crucial for translation regulation and stress response.

3-hour genome sequencing and targeted analysis to rapidly assess genetic risk.

Purpose: Rapid genetic testing in the critical care setting may guide diagnostic evaluation, direct therapies, and help families and care providers make informed decisions about goals of care. We tested whether a simplified DNA extraction and library preparation process would enable us to perform ultra-rapid assessment of genetic risk for a Mendelian condition, based on information from an affected sibling, using long-read genome sequencing and targeted analysis. Methods: Following extraction of DNA from cord blood and rapid library preparation, genome sequencing was performed on an Oxford Nanopore PromethION. FASTQ files were generated from original sequencing data in near real-time and aligned to a reference genome. Variant calling and analysis were performed at timed intervals. Results: We optimized the DNA extraction and library preparation methods to create sufficient library for sequencing from 500 µL of blood. Real-time, targeted analysis was performed to determine that the newborn was neither affected nor a heterozygote for variants underlying a Mendelian condition. Phasing of the target region and prior knowledge of the affected haplotypes supported our interpretation despite a low level of coverage at 3 hours of life. Conclusion: This proof-of-concept experiment demonstrates how prior knowledge of haplotype structure or familial variants can be used to rapidly evaluate an individual at risk for a genetic disease. While ultra-rapid sequencing remains both complex and cost prohibitive, our method is more easily automated than prior approaches and uses smaller volumes of blood, thus may be more easily adopted for future studies of ultra-rapid genome sequencing in the clinical setting.

Empirical assessment of the enrichment-based metagenomic methods in identifying diverse respiratory pathogens.

Probe-based nucleic acid enrichment represents an effective route to enhance the detection capacity of next-generation sequencing (NGS) in a set of clinically diverse and relevant microbial species. In this study, we assessed the effect of the enrichment-based sequencing on identifying respiratory infections using tiling RNA probes targeting 76 respiratory pathogens and sequenced using both Illumina and Oxford Nanopore platforms. Forty respiratory swab samples pre-tested for a panel of respiratory pathogens by qPCR were used to benchmark the sequencing data. We observed a general improvement in sensitivity after enrichment. The overall detection rate increased from 73 to 85% after probe capture detected by Illumina. Moreover, enrichment with probe sets boosted the frequency of unique pathogen reads by 34.6 and 37.8-fold for Illumina DNA and cDNA sequencing, respectively. This also resulted in significant improvements on genome coverage especially in viruses. Despite these advantages, we found that library pooling may cause reads mis-assignment, probably due to crosstalk issues arise from post-capture PCR and from pooled sequencing, thus increasing the risk of bleed-through signal. Taken together, an overall improvement in the breadth and depth of pathogen coverage is achieved using enrichment-based sequencing method. For future applications, automated library processing and pooling-free sequencing could enhance the precision and timeliness of probe enrichment-based clinical metagenomics.

Potential role of heteroplasmic mitochondrial DNA mutations in modulating the subtype-specific adaptation of oral squamous cell carcinoma to cisplatin therapy.

Cancer cells are constantly evolving to adapt to environmental changes, particularly during exposure to drug treatment. In this work, we aimed to characterize genetic and epigenetic changes in mitochondrial DNA (mtDNA) that may increase the resistance of oral squamous cell carcinoma (OSCC) to cisplatin. We first derived drug-resistant cells from two human OSCC cell lines, namely SAS and H103, by continual cisplatin treatments for about 4 months. To determine mtDNA changes induced by cisplatin, we performed nanopore sequencing and quantitative polymerase chain reaction analysis of mtDNA extracted from the cells pre- and post-treatment. We also assessed the mitochondrial functions of the cells and their capacity to generate intracellular reactive oxygen species (ROS). We found that in the cisplatin-resistant cells derived from SAS, there was a reduction in mtDNA content and significant enrichment of a m.3910G > C mutation in the MT-ND1 gene. However, such changes were not detected in cisplatin-resistant H103 cells. The cisplatin treatment also altered methylation patterns in both SAS and H103 cells and decreased their sensitivity to ROS-induced cytotoxicity. We suggest that the sequence alterations and epigenetic changes in mtDNA and the reduction in mtDNA content could be key drivers of cisplatin resistance in OSCC. These mtDNA alterations may participate in cellular adaptation that serves as a response to adverse changes in the environment, particularly exposure to cytotoxic agents. Importantly, the observed mtDNA changes may be influenced by the distinct genetic landscapes of various cancer subtypes. Overall, this study reveals significant insights into cisplatin resistance driven by complex mtDNA dynamics, particularly in OSCC. This underscores the need for targeted therapies tailored to the genetic profiles of individual OSCC patients to improve disease prognosis.

A Clickase-Mediated Immunoassay Based on Nanopore and Bionic Signal Labels for Ultrasensitive, Portable, and On-Site Detection of Ricin.

It is of particular importance to develop an effective method that possesses several merits simultaneously of rapid, ultrasensitive, portable, and on-site detection potential for food safety detection. Herein, we propose a clickase-mediated immunoassay based on nanopore and bionic signal labels for the detection of ricin. The introduction of Cu/Cys clickase and nanopore simultaneously effectively addressed the inherent limitations of natural enzymes and colorimetric signal output, respectively. Using this method, bionic signal labels can be easily formed through DNA and Gram-positive bacterial cell wall terminal peptide fragments (labeled by alkynyl and azide, respectively) and vancomycin. Translocation of the D-P@vancomycin through the nanopore generated highly specific oscillation current traces. This method showed a great on-site detection potential and superior analytical performance owing to the combination of the specificity of antibodies, high CuAAC click reaction catalytic efficiency of clickase, ultrasensitivity of the nanopore, and high signal resolution of D-P@vancomycin. Moreover, the practical applicability of the established method was also verified, achieving a limit of detection (LOD) down to 200.9 ag/mL with a wide linear relationship under the optimized conditions. In conclusion, this method is promising for rapid, portable, ultrasensitive, and on-site food safety detection.

What Quality Suffices for Nanopore Metabarcoding? Reconsidering Methodology and Ectomycorrhizae in Decaying Fagus sylvatica Bark as Case Study.

Nanopore raw read accuracy has improved to over 99%, making it a potential tool for metabarcoding. For broad adoption, guidelines on quality filtering are needed to ensure reliable taxonomic unit recovery. This study aims to provide those guidelines for a fungal metabarcoding context and to apply them to a case study of ectomycorrhizae in the decaying bark of Fagus sylvatica. We introduce the eNano pipeline to test two standard metabarcoding approaches: (1) Reference-based mapping leveraging UNITE's species hypothesis system (SH approach); (2) Constructing 98% OTUs (OTU approach). Our results demonstrate that both approaches are effective with Nanopore data. When using a reference database, we recommend strict mapping criteria rather than Phred-based filtering. Leveraging the SH-system further enhances reproducibility and facilitates cross-study communication. For the 98% OTUs, filtering reads at ≥Q25 is recommended. Our case study reveals that the decay gradient is a primary determinant of community composition and that specific mycorrhizal fungi colonize decaying bark. Complementing our metabarcoding results with root tip morphotypification, we identify Laccaria amethystina and Tomentella sublilacina as key ectomycorrhizae of saplings on decaying logs. These findings demonstrate that Nanopore sequencing can provide valuable ecological insights and support its broader use in fungal metabarcoding as read quality continues to improve.

Genomic Analysis of Aeromonas salmonicida ssp. salmonicida Isolates Collected During Multiple Clinical Outbreaks Supports Association with a Single Epidemiological Unit.

Outbreaks of furunculosis cause significant losses in salmonid aquaculture worldwide. With a recent rise in antimicrobial resistance, regulatory measures to minimize the use of antibiotics in animal husbandry, including aquaculture, have increased scrutiny and availability of veterinary medical products to control this disease in production facilities. In such a regulatory environment, the utility of autogenous vaccines to assist with disease prevention and control as a veterinary-guided prophylactic measure is of high interest to the producers and veterinary services alike. However, evolving concepts of epidemiological units and epidemiological links need to be considered during approval and acceptance procedures for the application of autogenous vaccines in multiple aquaculture facilities. Here, we present the results of solid-state nanopore sequencing (Oxford Nanopore Technologies, ONT) performed on 54 isolates of Aeromonas salmonicida ssp. salmonicida sampled during clinical outbreaks of furunculosis in different aquaculture facilities from Bavaria, Germany, from 2017 to 2020. All of the performed analyses (phylogeny, single nucleotide polymorphism and 3D protein modeling for major immunogenic proteins) support a high probability that all studied isolates belong to the same epidemiological unit. Simultaneously, we describe a cost/effective method of whole genome analysis with the usage of ONT as a viable strategy to study outbreaks of other pathogens in the field of aquatic veterinary medicine for the purpose of developing the best autogenous vaccine candidates applicable to multiple aquaculture establishments.

Comparison of direct cDNA and PCR-cDNA Nanopore sequencing of RNA from Escherichia coli isolates.

Whole-transcriptome (long-read) RNA sequencing (Oxford Nanopore Technologies, ONT) holds promise for reference-agnostic analysis of differential gene expression in pathogenic bacteria, including for antimicrobial resistance genes (ARGs). However, direct cDNA ONT sequencing requires large concentrations of polyadenylated mRNA, and amplification protocols may introduce technical bias. Here we evaluated the impact of direct cDNA- and cDNA PCR-based ONT sequencing on transcriptomic analysis of clinical Escherichia coli. Four E. coli bloodstream infection-associated isolates (n=2 biological replicates per isolate) were sequenced using the ONT Direct cDNA Sequencing SQK-DCS109 and PCR-cDNA Barcoding SQK-PCB111.24 kits. Biological and technical replicates were distributed over eight flow cells using 16 barcodes to minimize batch/barcoding bias. Reads were mapped to a transcript reference and transcript abundance was quantified after in silico depletion of low-abundance and rRNA genes. We found there were strong correlations between read counts using both kits and when restricting the analysis to include only ARGs. We highlighted that correlations were weaker for genes with a higher GC content. Read lengths were longer for the direct cDNA kit compared to the PCR-cDNA kit whereas total yield was higher for the PCR-cDNA kit. In this small but methodologically rigorous evaluation of biological and technical replicates of isolates sequenced with the direct cDNA and PCR-cDNA ONT sequencing kits, we demonstrated that PCR-based amplification substantially improves yield with largely unbiased assessment of core gene and ARG expression. However, users of PCR-based kits should be aware of a small risk of technical bias which appears greater for genes with an unusually high (>52%)/low (<44%) GC content.