Contig Length Research Articles

Chromosome-level assembly of the common vetch (Vicia sativa) reference genome.

Vicia sativa L. (common vetch, n =6) is an annual, herbaceous, climbing legume, originating in the Fertile Crescent of the Middle East and now widespread in the Mediterranean basin, West, Central and Eastern Asia, North and South America. V.sativa is of economic importance as a forage legume in countries such as Australia, China, and the USA, and contributes valuable nitrogen to agricultural rotation cropping systems. To accelerate precision genome breeding and genomics-based selection of this legume, we present a chromosome-level reference genome sequence for V.sativa, constructed using a combination of long-read Oxford Nanopore sequencing, short-read Illumina sequencing, and high-throughput chromosome conformation data (CHiCAGO and Hi-C) analysis. The chromosome-level assembly of six pseudo-chromosomes has a total genome length of 1.65Gbp, with a median contig length of 684Kbp. BUSCO analysis of the assembly demonstrated very high completeness of 98% of the dicotyledonous orthologs. RNA-seq analysis and gene modelling enabled the annotation of 53,218 protein-coding genes. This V.sativa assembly will provide insights into vetch genome evolution and be a valuable resource for genomic breeding, genetic diversity and for understanding adaption to diverse arid environments.

Slag: A program for seeded local assembly of genes in complex genomes.

Although finished genomes have become more common, there is still a need for assemblies of individual genes or chromosomal regions when only unassembled reads are available. slag (Seeded Local Assembly of Genes) fulfils this need by performing iterative local assembly based on cycles of matching‐read retrieval with blast and assembly with cap3, phrap, spades, canu or unicycler. The target sequence can be nucleotide or protein. Read fragmentation allows slag to use phrap or cap3 to assemble long reads at lower coverage (e.g., 5×) than is possible with canu or unicycler. In simple, nonrepetitive genomes, a slag assembly can cover a whole chromosome, but in complex genomes the growth of target‐matching contigs is limited as additional reads are consumed by consensus contigs consisting of repetitive elements. Apart from genomic complexity, contig length and correctness depend on read length and accuracy. With pyrosequencing or Illumina reads, slag‐assembled contigs are accurate enough to allow design of PCR primers, while contigs assembled from Oxford Nanopore or pre‐HiFi Pacific Biosciences long reads are generally only accurate enough to design baiting sequences for further targeted sequencing. In an application with real reads, slag successfully extended sequences for four wheat genes, which were verified by cloning and Sanger sequencing of overlapping amplicons. slag is a robust alternative to atram 2 for local assemblies, especially for read sets with less than 20× coverage. slag is freely available at https://github.com/cfcrane/SLAG.

A Superior Contiguous Whole Genome Assembly for Shrimp (Penaeus indicus)

Penaeid shrimp fishery and culture is a commercial enterprise contributing to employment, nutritional security and foreign exchange of developing countries. The genetic improvement programs being operated in shrimp benefit hugely from genomic resources. We report here a high-quality genome assembly for a penaeid shrimp, Penaeus indicus, which is the only Crustacean assembly to meet the reference standards of 1 and 10 Mb N50 lengths for contigs and scaffolds, respectively, among genomes of >1.5 Gb assembly length. The assembly is 1.93 Gb length (34.4 Mb scaffold N50) with 28,720 protein-coding genes and 49.31% repeat elements. The P. indicus assembly has 31.99% of simple sequence repeats, the highest among sequenced animal genomes. In comparison to other shrimp genomes having short contig lengths, the P. indicus assembly has 346 un-gapped contigs of over 1 Mb length and betters other shrimp genomes on sequence contiguity. This contiguous genome revealed 15,563 coding single nucleotide polymorphisms (SNPs) of which 2,572 are non-synonymous. The assembly and the SNP data resources have applications to genetic improvement programs, evolutionary studies and stock management.

Combined assembly of long and short sequencing reads improve the efficiency of exploring the soil metagenome

BackgroundAdvances in DNA sequencing technologies have transformed our capacity to perform life science research, decipher the dynamics of complex soil microbial communities and exploit them for plant disease management. However, soil is a complex conglomerate, which makes functional metagenomics studies very challenging.ResultsMetagenomes were assembled by long-read (PacBio, PB), short-read (Illumina, IL), and mixture of PB and IL (PI) sequencing of soil DNA samples were compared. Ortholog analyses and functional annotation revealed that the PI approach significantly increased the contig length of the metagenomic sequences compared to IL and enlarged the gene pool compared to PB. The PI approach also offered comparable or higher species abundance than either PB or IL alone, and showed significant advantages for studying natural product biosynthetic genes in the soil microbiomes.ConclusionOur results provide an effective strategy for combining long and short-read DNA sequencing data to explore and distill the maximum information out of soil metagenomics.

Development of polymorphic SSR markers and their applicability in genetic diversity evaluation in Euptelea pleiosperma

Euptelea pleiosperma is a characteristic species of East Asian flora with both ornamental and scientific values. Based on the reduced-representation sequencing (RRS) technology of RAD-Seq, this study conducted high-throughput Illumina paired-end sequencing to find SSR marker information in the genome of E. pleiosperma, and to screen and verify polymorphism of SSR markers. We obtained 5.5G of high-quality data using RAD-Seq. The total number of contigs of the RAD tags was 299,376, with the maximum contig length of 2,062 bp and the average length of 445 bp. From these sequences, we identified 20,718 SSR loci, with a distribution density of one SSR per 6.45 kb (1/6.45 kb). Among all SSRs, dinucleotides (52.00%) were the most detected SSRs, followed by mononucleotides (21.63%). AG/CT was the dominant motif in the SSR loci, accounting for 34.8%. Primers were successfully designed for 14,593 loci, and 100 pairs of these primers were randomly selected for chemical synthesis and validated by SSR-PCR amplification in 20 individuals of E. pleiosperma. Seventy-nine primers were able to amplify the target bands. Cervus 3.0 software was used to analyze the selected 20 SSR loci with good polymorphism. For the 20 SSR markers, the number of alleles ranged from 4 to 9, and the observed heterozygosity and expected heterozygosity were from 0.35 to 0.75 and 0.541 to 0.875, respectively. The information content of polymorphic loci ranged from 0.463 to 0.848, with an average value of 0.638. Among them, there were 18 highly polymorphic loci, and 20 SSR loci did not deviate from the HardyWeinberg equilibrium. Furthermore, the 20 pairs of SSR primers were used to conduct PCoA analysis based on Nei’s genetic distance of 51 individuals from three populations. The results showed that these SSR markers could distinguish genetic differences based on different geographical locations.

Genome Sequence of Chrysotila roscoffensis, a Coccolithphore Contributed to Global Biogeochemical Cycles

Chrysotila is a genus of coccolithophores. Together with Emiliania, it is one of the representative genera in the Haptophyta which have been extensively studied. They are photosynthetic unicellular marine algae sharing the common characteristic of the production of CaCO3 platelets (coccoliths) on the surface of their cells and are crucial contributors to global biogeochemical cycles. Here, we report the genome assembly of Chrysotila roscoffensis. The assembled genome size was ~636 Mb distributed across 769 scaffolds with N50 of 1.63 Mb, and maximum contig length of ~2.6 Mb. Repetitive elements accounted for approximately 59% of the genome. A total of 23,341 genes were predicted from C. roscoffensis genome. The divergence time between C. roscoffensis and Emiliania huxleyi was estimated to be around 537.6 Mya. Gene families related to cytoskeleton, cellular motility and morphology, and ion transport were expanded. The genome of C. roscoffensis will provide a foundation for understanding the genetic and phenotypic diversification and calcification mechanisms of coccolithophores.

The Genome of the "Sea Vomit" Didemnum vexillum.

Tunicates are the sister group of vertebrates and thus occupy a key position for investigations into vertebrate innovations as well as into the consequences of the vertebrate-specific genome duplications. Nevertheless, tunicate genomes have not been studied extensively in the past, and comparative studies of tunicate genomes have remained scarce. The carpet sea squirt Didemnum vexillum, commonly known as “sea vomit”, is a colonial tunicate considered an invasive species with substantial ecological and economical risk. We report the assembly of the D. vexillum genome using a hybrid approach that combines 28.5 Gb Illumina and 12.35 Gb of PacBio data. The new hybrid scaffolded assembly has a total size of 517.55 Mb that increases contig length about eightfold compared to previous, Illumina-only assembly. As a consequence of an unusually high genetic diversity of the colonies and the moderate length of the PacBio reads, presumably caused by the unusually acidic milieu of the tunic, the assembly is highly fragmented (L50 = 25,284, N50 = 6539). It is sufficient, however, for comprehensive annotations of both protein-coding genes and non-coding RNAs. Despite its shortcomings, the draft assembly of the “sea vomit” genome provides a valuable resource for comparative tunicate genomics and for the study of the specific properties of colonial ascidians.

Ancestral Flowering Plant Chromosomes and Gene Orders Based on Generalized Adjacencies and Chromosomal Gene Co-Occurrences.

Recurrent whole genome duplication and the ensuing loss of redundant genes-fractionation-complicate efforts to reconstruct the gene orders and chromosomes of the ancestors associated with the nodes of a phylogeny. Loss of genes disrupts the gene adjacencies key to current techniques. With our RACCROCHE pipeline, instead of starting with the inference of short ancestral segments, we suggest delaying the choice of gene adjacencies while we accumulate many more syntenically validated generalized (gapped) adjacencies. We obtain longer ancestral contigs using maximum weight matching (MWM). Similarly, we do not construct chromosomes by successively piecing together contigs into larger segments, but rather compile counts of pairwise contig co-occurrences on the set of extant genomes and use these to cluster the contigs. Chromosome-level contig assemblies for a monoploid genome emerge naturally at each node of the phylogeny and the contigs then can be ordered along the chromosome. Sampling alternative MWM solutions, visualizing heat maps, and applying gap statistics allow us to estimate the number of chromosomes in the reconstruction. We introduce several measures of quality: length of contigs, continuity of contig structure on successive ancestors, coverage of the extant genome by the reconstruction, and rearrangement relations among the inferred chromosomes. The reconstructed ancestors are visualized by painting the ancestral projections on the descendant genomes. We submit genomes drawn from a broad range of monocot orders to our pipeline, confirming the tetraploidization event "tau" in the stem lineage between the alismatids and the lilioids. We show additional applications to the Solanaceae and to four Brassica genomes, producing evidence about the monoploid ancestor in each case.

Genome-Wide Transcriptome Profiling, Characterization, and Functional Identification of NAC Transcription Factors in Sorghum under Salt Stress.

Salinity stress has become a significant concern to global food security. Revealing the mechanisms that enable plants to survive under salinity has immense significance. Sorghum has increasingly attracted researchers interested in understanding the survival and adaptation strategies to high salinity. However, systematic analysis of the DEGs (differentially expressed genes) and their relative expression has not been reported in sorghum under salt stress. The de novo transcriptomic analysis of sorghum under different salinity levels from 60 to 120 mM NaCl was generated using Illumina HiSeq. Approximately 323.49 million high-quality reads, with an average contig length of 1145 bp, were assembled de novo. On average, 62% of unigenes were functionally annotated to known proteins. These DEGs were mainly involved in several important metabolic processes, such as carbohydrate and lipid metabolism, cell wall biogenesis, photosynthesis, and hormone signaling. SSG 59-3 alleviated the adverse effects of salinity by suppressing oxidative stress (H2O2) and stimulating enzymatic and non-enzymatic antioxidant activities (SOD, APX, CAT, APX, POX, GR, GSH, ASC, proline, and GB), as well as protecting cell membrane integrity (MDA and electrolyte leakage). Significant up-regulation of transcripts encoding the NAC, MYB, and WRYK families, NHX transporters, the aquaporin protein family, photosynthetic genes, antioxidants, and compatible osmolyte proteins were observed. The tolerant line (SSG 59-3) engaged highly efficient machinery in response to elevated salinity, especially during the transport and influx of K+ ions, signal transduction, and osmotic homeostasis. Our data provide insights into the evolution of the NAC TFs gene family and further support the hypothesis that these genes are essential for plant responses to salinity. The findings may provide a molecular foundation for further exploring the potential functions of NAC TFs in developing salt-resistant sorghum lines.

Transcriptome profile of Haemaphysalis longicornis (Acari: Ixodidae) exposed to Cymbopogon citratus essential oil and citronellal suggest a cytotoxic mode of action involving mitochondrial Ca2+ overload and depolarization

Haemaphysalis longicornis is an ixodid tick species of medical and veterinary importance. Investigation into the acaricidal activities of botanicals have increased recently but information about their molecular mechanism of action is scarce. Here, RNA-seq analysis of the ticks exposed to Cymbopogon citratus essential oil and citronellal was performed and the responsive genes were identified. More than 6.39 G clean reads with Q20 ≥ 94.88% were obtained for each H. longicornis sample, with an average GC content of 50.94%. Using the Trinity method, 166,710 unigenes with a mean length of 869 bp and a maximum contig length of 29,156 bp were obtained. The upregulation of genes was concentration-dependent in most of the treated groups. Many genes responsive to C. citratus oil and citronellal were stress-related and they include genes associated with adrenergic signaling/calcium channels, cGMP-PKG signaling, apoptosis, focal adhesion, ECM-receptor interaction, ubiquitin-mediated proteolysis, mTOR signaling pathway, and longevity regulating pathway. The upregulation of genes (CACNAID, ADCY9, TPM1, and MYH6) associated with calcium channels suggests a neurotoxic mode of action, whereas, the upregulation of apoptosis-associated genes (CYC, DRONC, CASP7, CASP9, BCL2L1, bcl-xL, etc.) suggests a cytotoxic mode of action. The metabolism of C. citratus essential oil generates oxidative stress which increases the intra-mitochondrial free Ca2+ and triggers the formation of reactive oxygen species (ROS) that culminates to mitochondrial depolarization, ATP depletion, and either necrotic or apoptotic death. The neurotoxic and cytotoxic effects exhibited by the monoterpenes in H. longicornis is vital and could be exploited for the advancement of acaricide development and eco-friendly tick control.

Long-read metagenomics of multiple displacement amplified DNA of low-biomass human gut phageomes by SACRA pre-processing chimeric reads.

The human gut bacteriophage community (phageome) plays an important role in the host's health and disease; however, the entire structure is poorly understood, partly owing to the generation of many incomplete genomes in conventional short-read metagenomics. Here, we show long-read metagenomics of amplified DNA of low-biomass phageomes with multiple displacement amplification (MDA), involving the development of a novel bioinformatics tool, split amplified chimeric read algorithm (SACRA), that efficiently pre-processed numerous chimeric reads generated through MDA. Using five samples, SACRA markedly reduced the average chimera ratio from 72% to 1.5% in PacBio reads with an average length of 1.8 kb. De novo assembly of chimera-less PacBio long reads reconstructed contigs of ≥5 kb with an average proportion of 27%, which was 1% in contigs from MiSeq short reads, thereby dramatically improving contig length and genome completeness. Comparison of PacBio and MiSeq contigs found MiSeq contig fragmentations frequently near local repeats and hypervariable regions in the phage genomes, and those caused by multiple homologous phage genomes coexisting in the community. We also developed a reference-independent method to assess the completeness of the linear phage genomes. Overall, we established a SACRA-coupled long-read metagenomics robust to highly diverse gut phageomes, identifying high-quality circular and linear phage genomes with adequate sequence quantity.

Application of shotgun metagenomics sequencing and targeted sequence capture to detect circulating porcine viruses in the Dutch-German border region.

Porcine viruses have been emerging in recent decades, threatening animal and human health, as well as economic stability for pig farmers worldwide. Next‐generation sequencing (NGS) can detect and characterize known and unknown viruses but has limited sensitivity when an unbiased approach, such as shotgun metagenomics sequencing, is used. To increase the sensitivity of NGS for the detection of viruses, we applied and evaluated a broad viral targeted sequence capture (TSC) panel and compared it to an unbiased shotgun metagenomic approach. A cohort of 36 pooled porcine nasal swab and blood serum samples collected from both sides of the Dutch–German border region were evaluated. Overall, we detected 46 different viral species using TSC, compared to 40 viral species with a shotgun metagenomics approach. Furthermore, we performed phylogenetic analysis on recovered influenza A virus (FLUAV) genomes from Germany and revealed a close similarity to a zoonotic influenza strain previously detected in the Netherlands. Although TSC introduced coverage bias within the detected viruses, it improved sensitivity, genome sequence depth and contig length. In‐depth characterization of the swine virome, coupled with developing new enrichment techniques, can play a crucial role in the surveillance of circulating porcine viruses and emerging zoonotic pathogens.

Unraveling the Genome of a High Yielding Colombian Sugarcane Hybrid.

Recent developments in High Throughput Sequencing (HTS) technologies and bioinformatics, including improved read lengths and genome assemblers allow the reconstruction of complex genomes with unprecedented quality and contiguity. Sugarcane has one of the most complicated genomes among grassess with a haploid length of 1Gbp and a ploidies between 8 and 12. In this work, we present a genome assembly of the Colombian sugarcane hybrid CC 01-1940. Three types of sequencing technologies were combined for this assembly: PacBio long reads, Illumina paired short reads, and Hi-C reads. We achieved a median contig length of 34.94 Mbp and a total genome assembly of 903.2 Mbp. We annotated a total of 63,724 protein coding genes and performed a reconstruction and comparative analysis of the sucrose metabolism pathway. Nucleotide evolution measurements between orthologs with close species suggest that divergence between Saccharum officinarum and Saccharum spontaneum occurred <2 million years ago. Synteny analysis between CC 01-1940 and the S. spontaneum genome confirms the presence of translocation events between the species and a random contribution throughout the entire genome in current sugarcane hybrids. Analysis of RNA-Seq data from leaf and root tissue of contrasting sugarcane genotypes subjected to water stress treatments revealed 17,490 differentially expressed genes, from which 3,633 correspond to genes expressed exclusively in tolerant genotypes. We expect the resources presented here to serve as a source of information to improve the selection processes of new varieties of the breeding programs of sugarcane.

Full-length transcriptome assembly of andrias davidianus (amphibia: caudata) skin via hybrid sequencing.

The Chinese giant salamander, Andrias davidianus, is the largest amphibian species in the world; it is thus an economically and ecologically important species. The skin of A. davidianus exhibits complex adaptive structural and functional adaptations to facilitate survival in aquatic and terrestrial ecosystems. Here, we report the first full-length amphibian transcriptome from the dorsal skin of A. davidianus, which was assembled using hybrid sequencing and the PacBio and Illumina platforms. A total of 153,038 transcripts were hybrid assembled (mean length of 2039 bp and N50 of 2172 bp), and 133,794 were annotated in at least one database (nr, Swiss-Prot, KEGG, KOGs, GO, and nt). A total of 58,732, 68,742, and 115,876 transcripts were classified into 24 KOG categories, 1903 GO term categories, and 46 KEGG pathways (level 2), respectively. A total of 207,627 protein-coding regions, 785 transcription factors, 27,237 potential long non-coding RNAs, and 8299 simple sequence repeats were also identified. The hybrid-assembled transcriptome recovered more full-length transcripts, had a higher N50 contig length, and a higher annotation rate of unique genes compared with that assembled in previous studies using next-generation sequencing. The high-quality full-length reference gene set generated in this study will help elucidate the genetic characteristics of A. davidianus skin and aid the identification of functional skin proteins.

High-Quality Genome Resource of Clonostachys rosea Strain CanS41 by Oxford Nanopore Long-Read Sequencing.

Clonostachys rosea is a necrotrophic mycoparasitic fungus with excellent biological control ability against numerous fungal plant pathogens. Here, we performed genomic sequencing of C. rosea strain CanS41 using Oxford Nanopore sequencing technology. We generated a high-quality genome assembly (>99.99% accuracy), which comprised 26 contigs containing 60.68 Mb sequences with a GC content of 48.55% and a repeat content of 8.38%. The N50 contig length is 3.02 Mb. In total, 20,818 protein-coding genes were identified and functionally annotated. Genes encoding secreted proteins and carbohydrate-active enzymes as well as secondary metabolic gene clusters were also identified and analyzed. In summary, the high-quality genome assembly and gene annotation provided here will allow further exploration of biological functions and enhance biological control ability of C. rosea.

The karyotype, genome survey, and assembly of Mud artemisia (Artemisia selengensis).

Artemisia selengensis is traditional Chinese medicine and phytochemical analysis indicated that A. selengensis contains essential oils, fatty acids and phenolic acids. The lack of reference genomic information may lead to tardiness in molecular biology research of A. selengensis. Karyotype analysis, genome survey, and genome assembly was employed to acquire information on the genome structure of A. selengensis. The chromosome number is 2n = 2x = 36, karyotype formula is 28m + 8Sm, karyotype asymmetry coefficient is 58.8%, and karyotypes were symmetric to Stebbins' type 2A. Besides, the flow cytometry findings reported that the mean peak value of fluorescent intensity is 1,170,677, 2C DNA content is 12pg and the genome size was estimated to be approximately 5.87Gb. Furthermore, the genome survey generates 341,478,078 clean reads, unfortunately, after K-mer analysis, no significant peak can be observed, the heterozygosity, repetitive rate and genome size was unable to estimated. It is speculated that this phenomenon might be due to the complexity of genome structure. 37,266 contigs are preliminary assembled with Oxford Nanopore Technology (ONT) sequencing, totaling 804Mb and GC content was 34.08%. The total length is 804,475,881bp, N50 is 29,624bp, and the largest contig length is 239,792bp. This study reveals the preliminary information of genome size of A. selengensis. These findings may provide supportive information for sequencing and assembly of whole-genome sequencing and encourage the progress of functional gene discovery, genetic improvement, evolutionary study, and structural studies of A. selengensis.

De novo genome assembly of the marine teleost, bluefin trevally (Caranx melampygus).

The bluefin trevally, Caranx melampygus, also known as the bluefin kingfish or bluefin jack, is known for its remarkable, bright-blue fins. This marine teleost is a widely prized sportfish, but few resources have been devoted to the genomics and conservation of this species because it is not targeted by large-scale commercial fisheries. Population declines from recreational and artisanal overfishing have been observed in Hawai‘i, USA, resulting in both an interest in aquaculture and concerns about the long-term conservation of this species. Most research to-date has been performed in Hawai‘i, raising questions about the status of bluefin trevally populations across its Indo-Pacific range. Genomic resources allow for expanded research on stock status, genetic diversity, and population demography. We present a high quality, 711 Mb nuclear genome assembly of a Hawaiian bluefin trevally from noisy long-reads with a contig NG50 of 1.2 Mb and longest contig length of 8.9 Mb. As measured by single-copy orthologs, the assembly was 95% complete, and the genome is comprised of 16.9% repetitive elements. The assembly was annotated with 33.1 K protein-coding genes, 71.4% of which were assigned putative functions, using RNA-seq data from eight tissues from the same individual. This is the first whole-genome assembly published for the carangoid genus Caranx. Using this assembled genome, a multiple sequentially Markovian coalescent model was implemented to assess population demography. Estimates of effective population size suggest population expansion has occurred since the Late Pleistocene. This genome will be a valuable resource for comparative phylogenomic studies of carangoid fishes and will help elucidate demographic history and delineate stock structure for bluefin trevally populations throughout the Indo-Pacific.

Sequencing an F1 hybrid of Silurus asotus and S. meridionalis enabled the assembly of high-quality parental genomes

Genome complexity such as heterozygosity may heavily influence its de novo assembly. Sequencing somatic cells of the F1 hybrids harboring two sets of genetic materials from both of the paternal and maternal species may avoid alleles discrimination during assembly. However, the feasibility of this strategy needs further assessments. We sequenced and assembled the genome of an F1 hybrid between Silurus asotus and S. meridionalis using the SequelII platform and Hi-C scaffolding technologies. More than 300 Gb raw data were generated, and the final assembly obtained 2344 scaffolds composed of 3017 contigs. The N50 length of scaffolds and contigs was 28.55 Mb and 7.49 Mb, respectively. Based on the mapping results of short reads generated for the paternal and maternal species, each of the 29 chromosomes originating from S. asotus and S. meridionalis was recognized. We recovered nearly 94% and 96% of the total length of S. asotus and S. meridionalis. BUSCO assessments and mapping analyses suggested that both genomes had high completeness and accuracy. Further analyses demonstrated the high collinearity between S. asotus, S. meridionalis, and the related Pelteobagrus fulvidraco. Comparison of the two genomes with that assembled only using the short reads from non-hybrid parental species detected a small portion of sequences that may be incorrectly assigned to the different species. We supposed that at least part of these situations may have resulted from mitotic recombination. The strategy of sequencing the F1 hybrid genome can recover the vast majority of the parental genomes and may improve the assembly of complex genomes.

Draft Genome Assembly and Transcriptome Dataset for European Turnip (Brassica rapa L. ssp. rapifera), ECD4 Carrying Clubroot Resistance.

Draft Genome Assembly and Transcriptome Dataset for European Turnip (Brassica rapa L. ssp. rapifera), ECD4 Carrying Clubroot Resistance.

Simulation study and comparative evaluation of viral contiguous sequence identification tools

BackgroundViruses, including bacteriophages, are important components of environmental and human associated microbial communities. Viruses can act as extracellular reservoirs of bacterial genes, can mediate microbiome dynamics, and can influence the virulence of clinical pathogens. Various targeted metagenomic analysis techniques detect viral sequences, but these methods often exclude large and genome integrated viruses. In this study, we evaluate and compare the ability of nine state-of-the-art bioinformatic tools, including Vibrant, VirSorter, VirSorter2, VirFinder, DeepVirFinder, MetaPhinder, Kraken 2, Phybrid, and a BLAST search using identified proteins from the Earth Virome Pipeline to identify viral contiguous sequences (contigs) across simulated metagenomes with different read distributions, taxonomic compositions, and complexities.ResultsOf the tools tested in this study, VirSorter achieved the best F1 score while Vibrant had the highest average F1 score at predicting integrated prophages. Though less balanced in its precision and recall, Kraken2 had the highest average precision by a substantial margin. We introduced the machine learning tool, Phybrid, which demonstrated an improvement in average F1 score over tools such as MetaPhinder. The tool utilizes machine learning with both gene content and nucleotide features. The addition of nucleotide features improves the precision and recall compared to the gene content features alone.Viral identification by all tools was not impacted by underlying read distribution but did improve with contig length. Tool performance was inversely related to taxonomic complexity and varied by the phage host. For instance, Rhizobium and Enterococcus phages were identified consistently by the tools; whereas, Neisseria prophage sequences were commonly missed in this study.ConclusionThis study benchmarked the performance of nine state-of-the-art bioinformatic tools to identify viral contigs across different simulation conditions. This study explored the ability of the tools to identify integrated prophage elements traditionally excluded from targeted sequencing approaches. Our comprehensive analysis of viral identification tools to assess their performance in a variety of situations provides valuable insights to viral researchers looking to mine viral elements from publicly available metagenomic data.