Meat-Borne-Parasite: A Nanopore-Based Meta-Barcoding Work-Flow for Parasitic Microbiodiversity Assessment in the Wild Fauna of French Guiana.

French Guiana, located in the Guiana Shield, is a natural reservoir for many zoonotic pathogens that are of considerable medical or veterinary importance. Until now, there has been limited data available on the description of parasites circulating in this area, especially on protozoan belonging to the phylum Apicomplexa; conversely, the neighbouring countries describe a high parasitic prevalence in animals and humans. Epidemiological surveillance is necessary, as new potentially virulent strains may emerge from these forest ecosystems, such as Amazonian toxoplasmosis. However, there is no standard tool for detecting protozoa in wildlife. In this study, we developed Meat-Borne-Parasite, a high-throughput meta-barcoding workflow for detecting Apicomplexa based on the Oxford Nanopore Technologies sequencing platform using the 18S gene of 14 Apicomplexa positive samples collected in French Guiana. Sequencing reads were then analysed with MetONTIIME pipeline. Thanks to a scoring rule, we were able to classify 10 samples out of 14 as Apicomplexa positive and reveal the presence of co-carriages. The same samples were also sequenced with the Illumina platform for validation purposes. For samples identified as Apicomplexa positive by both platforms, a strong positive correlation at up to the genus level was reported. Overall, the presented workflow represents a reliable method for Apicomplexa detection, which may pave the way for more comprehensive biomonitoring of zoonotic pathogens.

Benchmarking of computational methods for m6A profiling with Nanopore direct RNA sequencing.

N6-methyladenosine (m6A) is the most abundant internal eukaryotic mRNA modification, and is involved in the regulation of various biological processes. Direct Nanopore sequencing of native RNA (dRNA-seq) emerged as a leading approach for its identification. Several software were published for m6A detection and there is a strong need for independent studies benchmarking their performance on data from different species, and against various reference datasets. Moreover, a computational workflow is needed to streamline the execution of tools whose installation and execution remains complicated. We developed NanOlympicsMod, a Nextflow pipeline exploiting containerized technology for comparing 14 tools for m6A detection on dRNA-seq data. NanOlympicsMod was tested on dRNA-seq data generated from in vitro (un)modified synthetic oligos. The m6A hits returned by each tool were compared to the m6A position known by design of the oligos. In addition, NanOlympicsMod was used on dRNA-seq datasets from wild-type and m6A-depleted yeast, mouse and human, and each tool's hits were compared to reference m6A sets generated by leading orthogonal methods. The performance of the tools markedly differed across datasets, and methods adopting different approaches showed different preferences in terms of precision and recall. Changing the stringency cut-offs allowed for tuning the precision-recall trade-off towards user preferences. Finally, we determined that precision and recall of tools are markedly influenced by sequencing depth, and that additional sequencing would likely reveal additional m6A sites. Thanks to the possibility of including novel tools, NanOlympicsMod will streamline the benchmarking of m6A detection tools on dRNA-seq data, improving future RNA modification characterization.

The Oxford Nanopore MinION as a Versatile Technology for the Diagnosis and Characterization of Emerging Plant Viruses.

The emergence of novel viral epidemics that could affect major crops represents a serious threat to global food security. The early and accurate identification of the causative viral agent is the most important step for a rapid and effective response to disease outbreaks. Over the last years, the Oxford Nanopore Technologies (ONT) MinION sequencer has been proposed as an effective diagnostic tool for the early detection and identification of emerging viruses in plants, providing many advantages compared with different high-throughput sequencing (HTS) technologies. Here, we provide a step-by-step protocol that we optimized to obtain the virome of "Lamon bean" plants (Phaseolus vulgaris L.), an agricultural product with Protected Geographical Indication (PGI) in North-East of Italy, which is frequently subjected to multiple infections caused by different RNA viruses. The conversion of viral RNA in ds-cDNA enabled the use of Genomic DNA Ligation Sequencing Kit and Native Barcoding DNA Kit, which have been originally developed for DNA sequencing. This allowed the simultaneous diagnosis of both DNA- and RNA-based pathogens, providing a more versatile alternative to the use of direct RNA and/or direct cDNA sequencing kits.

CRISPR/Cas9-Mediated Enrichment Coupled to Nanopore Sequencing Provides a Valuable Tool for the Precise Reconstruction of Large Genomic Target Regions.

Complete and accurate identification of genetic variants associated with specific phenotypes can be challenging when there is a high level of genomic divergence between individuals in a study and the corresponding reference genome. We have applied the Cas9-mediated enrichment coupled to nanopore sequencing to perform a targeted de novo assembly and accurately reconstruct a genomic region of interest. This approach was used to reconstruct a 250-kbp target region on chromosome 5 of the common bean genome (Phaseolus vulgaris) associated with the shattering phenotype. Comparing a non-shattering cultivar (Midas) with the reference genome revealed many single-nucleotide variants and structural variants in this region. We cut five 50-kbp tiled sub-regions of Midas genomic DNA using Cas9, followed by sequencing on a MinION device and de novo assembly, generating a single contig spanning the whole 250-kbp region. This assembly increased the number of Illumina reads mapping to genes in the region, improving their genotypability for downstream analysis. The Cas9 tiling approach for target enrichment and sequencing is a valuable alternative to whole-genome sequencing for the assembly of ultra-long regions of interest, improving the accuracy of downstream genotype-phenotype association analysis.

STArS (STrain-Amplicon-Seq), a targeted nanopore sequencing workflow for SARS-CoV-2 diagnostics and genotyping.

Diagnostic tests based on reverse transcription-quantitative polymerase chain reaction (RT-qPCR) are the gold standard approach to detect severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection from clinical specimens. However, unless specifically optimized, this method is usually unable to recognize the specific viral strain responsible of coronavirus disease 2019, a crucial information that is proving increasingly important in relation to virus spread and treatment effectiveness. Even if some RT-qPCR commercial assays are currently being developed for the detection of viral strains, they focus only on single/few genetic variants that may not be sufficient to uniquely identify a specific strain. Therefore, genome sequencing approaches remain the most comprehensive solution for virus genotyping and to recognize viral strains, but their application is much less widespread due to higher costs. Starting from the well-established ARTIC protocol coupled to nanopore sequencing, in this work, we developed STArS (STrain-Amplicon-Seq), a cost/time-effective sequencing-based workflow for both SARS-CoV-2 diagnostics and genotyping. A set of 10 amplicons was initially selected from the ARTIC tiling panel, to cover: (i) all the main biologically relevant genetic variants located on the Spike gene; (ii) a minimal set of variants to uniquely identify the currently circulating strains; (iii) genomic sites usually amplified by RT-qPCR method to identify SARS-CoV-2 presence. PCR-amplified clinical samples (both positive and negative for SARS-CoV-2 presence) were pooled together with a serially diluted exogenous amplicon at known concentration and sequenced on a MinION device. Thanks to a scoring rule, STArS had the capability to accurately classify positive samples in agreement with RT-qPCR results, both at the qualitative and quantitative level. Moreover, the method allowed to effectively genotype strain-specific variants and thus also return the phylogenetic classification of SARS-CoV-2-postive samples. Thanks to the reduced turnaround time and costs, the proposed approach represents a step towards simplifying the clinical application of sequencing for viral genotyping, hopefully aiding in combatting the global pandemic.

Characterization of full-length CNBP expanded alleles in myotonic dystrophy type 2 patients by Cas9-mediated enrichment and nanopore sequencing.

Myotonic dystrophy type 2 (DM2) is caused by CCTG repeat expansions in the CNBP gene, comprising 75 to >11,000 units and featuring extensive mosaicism, making it challenging to sequence fully expanded alleles. To overcome these limitations, we used PCR-free Cas9-mediated nanopore sequencing to characterize CNBP repeat expansions at the single-nucleotide level in nine DM2 patients. The length of normal and expanded alleles can be assessed precisely using this strategy, agreeing with traditional methods, and revealing the degree of mosaicism. We also sequenced an entire ~50 kbp expansion, which has not been achieved previously for DM2 or any other repeat-expansion disorders. Our approach precisely counted the repeats and identified the repeat pattern for both short interrupted and uninterrupted alleles. Interestingly, in the expanded alleles, only two DM2 samples featured the expected pure CCTG repeat pattern, while the other seven presented also TCTG blocks at the 3' end, which have not been reported before in DM2 patients, but confirmed hereby with orthogonal methods. The demonstrated approach simultaneously determines repeat length, structure/motif, and the extent of somatic mosaicism, promising to improve the molecular diagnosis of DM2 and achieve more accurate genotype-phenotype correlations for the better stratification of DM2 patients in clinical trials.

The Virome of 'Lamon Bean': Application of MinION Sequencing to Investigate the Virus Population Associated with Symptomatic Beans in the Lamon Area, Italy.

'Lamon bean' is a protected geographical indication (PGI) for a product of four varieties of bean (Phaseolus vulgaris L.) grown in a specific area of production, which is located in the Belluno district, Veneto region (N.E. of Italy). In the last decade, the 'Lamon bean' has been threatened by severe virus epidemics that have compromised its profitability. In this work, the full virome of seven bean samples showing different foliar symptoms was obtained by MinION sequencing. Evidence that emerged from sequencing was validated through RT-PCR and ELISA in a large number of plants, including different ecotypes of Lamon bean and wild herbaceous hosts that may represent a virus reservoir in the field. Results revealed the presence of bean common mosaic virus (BCMV), cucumber mosaic virus (CMV), peanut stunt virus (PSV), and bean yellow mosaic virus (BYMV), which often occurred as mixed infections. Moreover, both CMV and PSV were reported in association with strain-specific satellite RNAs (satRNAs). In conclusion, this work sheds light on the cause of the severe diseases affecting the 'Lamon bean' by exploitation of MinION sequencing.

'Nebbiolo' genome assembly allows surveying the occurrence and functional implications of genomic structural variations in grapevines (Vitis vinifera L.).

BACKGROUND: 'Nebbiolo' is a grapevine cultivar typical of north-western Italy, appreciated for producing high-quality red wines. Grapevine cultivars are characterized by possessing highly heterozygous genomes, including a great incidence of genomic rearrangements larger than 50 bp, so called structural variations (SVs). Even though abundant, SVs are an under-explored source of genetic variation mainly due to methodological limitations at their detection. RESULTS: We employed a multiple platform approach to produce long-range genomic data for two different 'Nebbiolo' clones, namely: optical mapping, long-reads and linked-reads. We performed a haplotype-resolved de novo assembly for cultivar 'Nebbiolo' (clone CVT 71) and used an ab-initio strategy to annotate it. The annotated assembly enhanced our ability to detect SVs, enabling the study of genomic regions not present in the grapevines' reference genome and accounting for their functional implications. We performed variant calling analyses at three different organizational levels: i) between haplotypes of clone CVT 71 (primary assembly vs haplotigs), ii) between 'Nebbiolo' and 'Cabernet Sauvignon' assemblies and iii) between clones CVT 71 and CVT 185, representing different 'Nebbiolo' biotypes. The cumulative size of non-redundant merged SVs indicated a total of 79.6 Mbp for the first comparison and 136.1 Mbp for the second one, while no SVs were detected for the third comparison. Interestingly, SVs differentiating cultivars and haplotypes affected similar numbers of coding genes. CONCLUSIONS: Our results suggest that SVs accumulation rate and their functional implications in 'Nebbiolo' genome are highly-dependent on the organizational level under study. SVs are abundant when comparing 'Nebbiolo' to a different cultivar or the two haplotypes of the same individual, while they turned absent between the two analysed clones.

Real-Time On-Site Diagnosis of Quarantine Pathogens in Plant Tissues by Nanopore-Based Sequencing.

Rapid and sensitive assays for the identification of plant pathogens are necessary for the effective management of crop diseases. The main limitation of current diagnostic testing is the inability to combine broad and sensitive pathogen detection with the identification of key strains, pathovars, and subspecies. Such discrimination is necessary for quarantine pathogens, whose management is strictly dependent on genotype identification. To address these needs, we have established and evaluated a novel all-in-one diagnostic assay based on nanopore sequencing for the detection and simultaneous characterization of quarantine pathogens, using Xylella fastidiosa as a case study. The assay proved to be at least as sensitive as standard diagnostic tests and the quantitative results agreed closely with qPCR-based analysis. The same sequencing results also allowed discrimination between subspecies when present either individually or in combination. Pathogen detection and typing were achieved within 13 min of sequencing owing to the use of an internal control that allowed to stop sequencing when sufficient data had accumulated. These advantages, combined with the use of portable equipment, will facilitate the development of next-generation diagnostic assays for the efficient monitoring of other plant pathogens.

A new giant keelback slug of the genus Limax from the Balkans, described by citizen scientists.

Background: Despite their large size, striking colouration and genital extravagance, the taxonomy of the European giant keelback slugs of the genus Limax is still poorly understood. Preliminary morphological and molecular data suggest that many unnamed or unrecognised species exist, especially in the Alps, the Mediterranean and the Balkans. New information: We organised a citizen science expedition to Durmitor National Park in Montenegro and discovered a new species, genetically distinct, but morphologically similar to the sympatric L.cinereoniger Wolf 1803 and describe it as L.pseudocinereoniger.

Characterization of FMR1 Repeat Expansion and Intragenic Variants by Indirect Sequence Capture.

Traditional methods for the analysis of repeat expansions, which underlie genetic disorders, such as fragile X syndrome (FXS), lack single-nucleotide resolution in repeat analysis and the ability to characterize causative variants outside the repeat array. These drawbacks can be overcome by long-read and short-read sequencing, respectively. However, the routine application of next-generation sequencing in the clinic requires target enrichment, and none of the available methods allows parallel analysis of long-DNA fragments using both sequencing technologies. In this study, we investigated the use of indirect sequence capture (Xdrop technology) coupled to Nanopore and Illumina sequencing to characterize FMR1, the gene responsible of FXS. We achieved the efficient enrichment (> 200×) of large target DNA fragments (~60-80 kbp) encompassing the entire FMR1 gene. The analysis of Xdrop-enriched samples by Nanopore long-read sequencing allowed the complete characterization of repeat lengths in samples with normal, pre-mutation, and full mutation status (> 1 kbp), and correctly identified repeat interruptions relevant for disease prognosis and transmission. Single-nucleotide variants (SNVs) and small insertions/deletions (indels) could be detected in the same samples by Illumina short-read sequencing, completing the mutational testing through the identification of pathogenic variants within the FMR1 gene, when no typical CGG repeat expansion is detected. The study successfully demonstrated the parallel analysis of repeat expansions and SNVs/indels in the FMR1 gene at single-nucleotide resolution by combining Xdrop enrichment with two next-generation sequencing approaches. With the appropriate optimization necessary for the clinical settings, the system could facilitate both the study of genotype-phenotype correlation in FXS and enable a more efficient diagnosis and genetic counseling for patients and their relatives.

ACoRE: Accurate SARS-CoV-2 genome reconstruction for the characterization of intra-host and inter-host viral diversity in clinical samples and for the evaluation of re-infections.

Sequencing the SARS-CoV-2 genome from clinical samples can be challenging, especially in specimens with low viral titer. Here we report Accurate SARS-CoV-2 genome Reconstruction (ACoRE), an amplicon-based viral genome sequencing workflow for the complete and accurate reconstruction of SARS-CoV-2 sequences from clinical samples, including suboptimal ones that would usually be excluded even if unique and irreplaceable. The protocol was optimized to improve flexibility and the combination of technical replicates was established as the central strategy to achieve accurate analysis of low-titer/suboptimal samples. We demonstrated the utility of the approach by achieving complete genome reconstruction and the identification of false-positive variants in >170 clinical samples, thus avoiding the generation of inaccurate and/or incomplete sequences. Most importantly, ACoRE was crucial to identify the correct viral strain responsible of a relapse case, that would be otherwise mis-classified as a re-infection due to missing or incorrect variant identification by a standard workflow.

Exosomes from Plasma of Neuroblastoma Patients Contain Doublestranded DNA Reflecting the Mutational Status of Parental Tumor Cells.

Neuroblastoma (NB) is an aggressive infancy tumor, leading cause of death among preschool age diseases. Here we focused on characterization of exosomal DNA (exo-DNA) isolated from plasma cell-derived exosomes of neuroblastoma patients, and its potential use for detection of somatic mutations present in the parental tumor cells. Exosomes are small extracellular membrane vesicles secreted by most cells, playing an important role in intercellular communications. Using an enzymatic method, we provided evidence for the presence of double-stranded DNA in the NB exosomes. Moreover, by whole exome sequencing, we demonstrated that NB exo-DNA represents the entire exome and that it carries tumor-specific genetic mutations, including those occurring on known oncogenes and tumor suppressor genes in neuroblastoma (ALK, CHD5, SHANK2, PHOX2B, TERT, FGFR1, and BRAF). NB exo-DNA can be useful to identify variants responsible for acquired resistance, such as mutations of ALK, TP53, and RAS/MAPK genes that appear in relapsed patients. The possibility to isolate and to enrich NB derived exosomes from plasma using surface markers, and the quick and easy extraction of exo-DNA, gives this methodology a translational potential in the clinic. Exo-DNA can be an attractive non-invasive biomarker for NB molecular diagnostic, especially when tissue biopsy cannot be easily available.

Hydraena (s.str.) dinarica, new species (Coleoptera: Hydraenidae) along with further records of Hydraena spp. from Durmitor National Park, Montenegro and comments on the DNA barcoding problem with the genus.

BACKGROUND: Long-palped Water Beetles were collected during a taxon expedition in Montenegro which involved citizen scientists, students and taxonomists. The material was collected from springs, brooks, fens and the Tara River, at altitudes between 600 m and 1450 m above sea level, using fine-meshed hand-nets and by manual checking of submerged substrates. The morphological species delimitation was supplemented and congruent with mtDNA sequences mainly obtained in the field using the newly-developed MinION-based ONTrack pipeline. NEW INFORMATION: The new species Hydraena dinarica Freitag & de Vries, sp. n. from Durmitor Mt. is described, illustrated and compared in detail to closely-related congeners of the H. saga d'Orchymont, 1930/H. emarginata Rey, 1885 species complex. Five additional species and female specimens of two unidentified morphospecies of the genus were also recorded in the vicinity of Durmitor National Park. New records and the first DNA barcodes for Hydraena biltoni Jäch & Díaz, 2012 (endemic to Montenegro) and H. morio Kiesenwetter, 1849 are provided. Further records of H. nigrita Germar, 1824, H. minutissima Stephens, 1829, H. subintegra Ganglbauer, 1901 and females of two unidentified morphospecies are commented upon. The resulting inter- and intraspecific genetic distances and some observations of low or zero sequence divergence between recently-diverged species of Hydraena Kugelann, 1794 are briefly discussed.

A Long-Read Sequencing Approach for Direct Haplotype Phasing in Clinical Settings.

The reconstruction of individual haplotypes can facilitate the interpretation of disease risks; however, high costs and technical challenges still hinder their assessment in clinical settings. Second-generation sequencing is the gold standard for variant discovery but, due to the production of short reads covering small genomic regions, allows only indirect haplotyping based on statistical methods. In contrast, third-generation methods such as the nanopore sequencing platform developed by Oxford Nanopore Technologies (ONT) generate long reads that can be used for direct haplotyping, with fewer drawbacks. However, robust standards for variant phasing in ONT-based target resequencing efforts are not yet available. In this study, we presented a streamlined proof-of-concept workflow for variant calling and phasing based on ONT data in a clinically relevant 12-kb region of the APOE locus, a hotspot for variants and haplotypes associated with aging-related diseases and longevity. Starting with sequencing data from simple amplicons of the target locus, we demonstrated that ONT data allow for reliable single-nucleotide variant (SNV) calling and phasing from as little as 60 reads, although the recognition of indels is less efficient. Even so, we identified the best combination of ONT read sets (600) and software (BWA/Minimap2 and HapCUT2) that enables full haplotype reconstruction when both SNVs and indels have been identified previously using a highly-accurate sequencing platform. In conclusion, we established a rapid and inexpensive workflow for variant phasing based on ONT long reads. This allowed for the analysis of multiple samples in parallel and can easily be implemented in routine clinical practice, including diagnostic testing.

A Rapid and Accurate MinION-Based Workflow for Tracking Species Biodiversity in the Field.

Genetic markers (DNA barcodes) are often used to support and confirm species identification. Barcode sequences can be generated in the field using portable systems based on the Oxford Nanopore Technologies (ONT) MinION sequencer. However, to achieve a broader application, current proof-of-principle workflows for on-site barcoding analysis must be standardized to ensure a reliable and robust performance under suboptimal field conditions without increasing costs. Here, we demonstrate the implementation of a new on-site workflow for DNA extraction, PCR-based barcoding, and the generation of consensus sequences. The portable laboratory features inexpensive instruments that can be carried as hand luggage and uses standard molecular biology protocols and reagents that tolerate adverse environmental conditions. Barcodes are sequenced using MinION technology and analyzed with ONTrack, an original de novo assembly pipeline that requires as few as 1000 reads per sample. ONTrack-derived consensus barcodes have a high accuracy, ranging from 99.8 to 100%, despite the presence of homopolymer runs. The ONTrack pipeline has a user-friendly interface and returns consensus sequences in minutes. The remarkable accuracy and low computational demand of the ONTrack pipeline, together with the inexpensive equipment and simple protocols, make the proposed workflow particularly suitable for tracking species under field conditions.

New molecular insights on the response of the green alga Tetraselmis suecica to nitrogen starvation.

Microalgae are currently considered one of the most promising resources for biofuel production, aquaculture feedstock and new pharmaceuticals. Among them, green algae of the genus Tetraselmis are extensively studied for their lipid accumulation in nutrient-starvation conditions. In this paper, we present the full-transcriptome of Tetraselmis suecica and differential expression analysis between nitrogen-starved and -repleted conditions (at stationary phase) focusing not only on lipid metabolism but giving new insights on nutrient starvation responses. Transcripts involved in signal transduction pathways, stress and antioxidant responses and solute transport were strongly up-regulated when T. suecica was cultured under nitrogen starvation. On the contrary, transcripts involved in amino acid synthesis, degradation of sugars, secondary metabolite synthesis, as well as photosynthetic activity were down-regulated under the same conditions. Among differentially expressed transcripts, a polyketide synthase and three lipoxygenases (involved in the synthesis of secondary metabolites with antipredator, anticancer and anti-infective activities) were identified, suggesting the potential synthesis of bioactive compounds by this microalga. In addition, the transcript for a putative nitrilase, enzyme used in nitrile bioremediation, is here reported for the first time for T. suecica. These findings give new insights on T. suecica responses to nutrient starvation and on possible biotechnological applications for green algae.