<abstract xmlns="http://www.w3.org/1999/xhtml">

<p><italic>Anguina tritici</italic> is the first plant-parasitic nematode described in literature, dating back to the year 1743. It is responsible for causing earcockle (seed gall) and tundu diseases in wheat and rye. Notably, this nematode has been observed to survive in an anhydrobiotic state for up to 32 years within wheat seed galls. These exceptional characteristics have inspired the sequencing of the <italic>A. tritici</italic> genome. In this study, we present the initial draft genome of <italic>A. tritici</italic>, obtained using the Illumina MiSeq platform with coverage of 60-fold. The genome is estimated to have a size of 164 Mb and comprises 39,965 protein-coding genes, exhibiting a GC content of 39.1%. The availability of this genome data will serve as a foundation for future functional biological investigations, particularly for genes whose functions remain unknown to this day.</p>
</abstract>

Anguina tritici is the first plant-parasitic nematode described in literature, dating back to the year 1743. It is responsible for causing earcockle (seed gall) and tundu diseases in wheat and rye. Notably, this nematode has been observed to survive in an anhydrobiotic state for up to 32 years within wheat seed galls. These exceptional characteristics have inspired the sequencing of the A. tritici genome. In this study, we present the initial draft genome of A. tritici, obtained using the Illumina MiSeq platform with coverage of 60-fold. The genome is estimated to have a size of 164 Mb and comprises 39,965 protein-coding genes, exhibiting a GC content of 39.1%. The availability of this genome data will serve as a foundation for future functional biological investigations, particularly for genes whose functions remain unknown to this day.

<div xmlns="http://www.w3.org/1999/xhtml">
<p>At least four distinct major lineages of plant parasites have emerged within the species-rich and trophically diverse phylum Nematoda (<a ref-type="bibr" href="#j_jofnem-2023-0031_ref_015">Quist et al., 2015</a>). <italic>A. tritici</italic>, commonly known as the seed gall nematode, is a plant parasitic nematode (PPN) that infects the seeds of wheat and barley, leading to the development of the earcockle disease. This nematode also acts as a vector for the bacteria <italic>Clavibacter tritici</italic>, causing tundu disease in their association (<a ref-type="bibr" href="#j_jofnem-2023-0031_ref_001">Swarup and Gupta, 1971</a>). According to the phylogenetic analysis conducted by <a ref-type="bibr" href="#j_jofnem-2023-0031_ref_022">van Megen et al. (2009)</a>, <italic>Anguina tritici</italic> is classified within the Anguinidae family, along with its closely related species <italic>Ditylenchus dipsaci</italic> and <italic>D. destructor</italic>. In the past decade, the genomes of several parasitic nematodes have been sequenced, enabling comparative genomic analyses that shed light on the gain or loss of novel genes among nematode species (<a ref-type="bibr" href="#j_jofnem-2023-0031_ref_002">Wasmuth et al., 2008</a>). As more genome sequences become available, researchers will be able to identify species-specific peculiarities in biological characteristics.</p>
<p>The biology of <italic>A. tritici</italic> differs from other PPNs, which usually feed on or in the roots of the crop plants but do not infest the seeds. <italic>A. tritici</italic> is an obligate parasite specifically targeting the seeds of its host plant wheat and barley, unlike the rice white tip nematode, <italic>Aphelenchoides besseyi</italic>, which also survives in an anhydrobiotic state but in dry tissues, under hulls of rice grains. <italic>A. tritici</italic> exhibits the remarkable ability to survive for several decades in an anhydrobiotic state without losing its parasitic capability (<a ref-type="bibr" href="#j_jofnem-2023-0031_ref_003">Limber, 1973</a>) but <italic>A. besseyi</italic> is reported to survive for a maximum of three years without losing its parasitic ability (<a ref-type="bibr" href="#j_jofnem-2023-0031_ref_023">Cralley, 1949</a>; <a ref-type="bibr" href="#j_jofnem-2023-0031_ref_024">Yoshi and Yamamoto, 1950</a>) Among plant parasitic nematodes, <italic>A. tritici</italic> exhibits notable variability in certain biological traits, distinguishing it from the other PPN species and making it an attractive model system for investigating fundamental questions related to survival, ageing, seed parasitism, bacterial associations, and more. Currently, limited genomic information is available regarding plant-parasitic nematodes that specifically target crop seed.</p>
<p>Recently, the genome of <italic>Aphelenchoides besseyi</italic>, a nematode that feeds on the growing tips of rice seedlings, causing the white-tip disease, has been published (<a ref-type="bibr" href="#j_jofnem-2023-0031_ref_020">Lai et al., 2023</a>; <a ref-type="bibr" href="#j_jofnem-2023-0031_ref_021">Ji et al., 2023</a>). However, unlike <italic>A. tritici</italic>, it does not enter the grain but remains localized in the husk part of the paddy seed. The genome sequences of two closely related species of <italic>Anguina tritici</italic>, namely <italic>Ditylenchus dipsaci</italic> (<a ref-type="bibr" href="#j_jofnem-2023-0031_ref_016">Mimee et al., 2019</a>) and <italic>D. destructor</italic> (<a ref-type="bibr" href="#j_jofnem-2023-0031_ref_019">Zheng et al., 2016</a>) are also available. These genome sequences provide valuable resources for comparative analysis and serve as important references for understanding the genomic features of <italic>A. tritici</italic> in a broader context. The objective of this study was to present the draft genome of the seed gall nematode, <italic>A. tritici</italic>. To propagate the <italic>A. tritici</italic>, cockled seeds were crushed and placed in a modified Baermann funnel for the extraction of second-stage juveniles (J2s). The extracted J2s were then inoculated onto greenhouse-grown wheat cultivar ‘PBW 343’ for multiple generations to generate an inbred population of the nematode. Before DNA extraction, J2s were surface sterilized using M9 buffer and nematode sterilization buffer to remove any contaminants (<a ref-type="bibr" href="#j_jofnem-2023-0031_ref_004">Joshi et al., 2019</a>). High-quality genomic DNA was extracted using the modified fast CTAB method (<a ref-type="bibr" href="#j_jofnem-2023-0031_ref_005">Thomas et al., 1997</a>). The isolated high-quality DNA was utilized for the preparation of a genomic DNA sequencing library using the TruSeq DNA sample preparation kit from Illumina technologies, USA by following the manufacturer’s protocol. The quality of DNA libraries was assessed using the Bioanalyzer 2100 from Agilent technologies, USA. Paired-end sequencing of the genomic library was performed on the Illumina MiSeq (2X300) platform. The genome size was estimated using the K-mer value of the Abyss Version 1.3.7 (<a ref-type="bibr" href="#j_jofnem-2023-0031_ref_008">Simpson et al., 2009</a>). Paired-end files were merged using PandaSeq 2.5 (<a ref-type="bibr" href="#j_jofnem-2023-0031_ref_006">Masella et al., 2012</a>), and the reads were subjected to quality checking using FastQC software v3.0 (Andrew, 2012). High-quality reads with a Phred score of ≥20, after trimming adapters using trimmomatic version 0.39 (<a ref-type="bibr" href="#j_jofnem-2023-0031_ref_007">Bolger et al., 2014</a>), were retained. <italic>De novo</italic> assembly was carried out using Abyss Version 1.3.7 (<a ref-type="bibr" href="#j_jofnem-2023-0031_ref_008">Simpson et al., 2009</a>) at different K-mers (45, 50, 55, and 61). The quality assessment of the assembled genome was performed using BUSCO version 3.0.2 (<a ref-type="bibr" href="#j_jofnem-2023-0031_ref_009">Robert et al., 2017</a>). Protein-coding genes were predicted using the <italic>ab initio</italic> approach of AUGUSTUS version 3.0.3 (<a ref-type="bibr" href="#j_jofnem-2023-0031_ref_010">Stanke et al., 2004</a>) with default parameters. Predicted genes with a length of less than 400 bp were excluded, and the remaining genes were annotated against the NCBI redundant nucleotide database using RAPsearch2 and an in-house perl script (<a ref-type="bibr" href="#j_jofnem-2023-0031_ref_011">Yuzhen et al., 2011</a>).</p>
<p>The Illumina MiSeq sequencing strategy generated a total of 11,065,381,294 bases, with a read count of 40,210,848. The GC percentage, Q20 percentage, and Q30 percentage were calculated as 39.1%, 92.73%, and 84.08%, respectively. By merging the paired-end reads using the Pandaseq program, a total of 7,938,027,993 bases and 18,913,374 read counts were obtained, with a GC percentage of 39.1%. The draft genome assembly of <italic>A. tritici</italic> resulted in a size of up to 164 Mb with a GC content of 39.1%. The genome sequence data has been submitted to GenBank with the following identifiers: BioProject ID: PRJNA802691, BioSample: SAMN25565353, and accession ID: JAKMXJ000000000. The genome achieved a total coverage level of 60-fold. Quantitative assessment of assembly completeness was performed using BUSCO, revealing C: 36.6% [S: 35.5%, D: 1.1%], F: 14.8%, M: 48.6%, and n: 982 (<a ref-type="table" href="#j_jofnem-2023-0031_tab_001">Table 1</a>). A total of 44,954 protein coding genes were predicted using the <italic>ab initio</italic> Augustus version 3.0.3. After filtering genes with a length below 400 base pairs, a total of 39,965 predicted genes were identified, resulting in an overall gene density of 26.5784. The average gene length was 1,347 bp, while the average exon and intron lengths were 168 bp and 280 bp, respectively. Of the 41,746 predicted transcripts, 62% were successfully annotated with functional annotations against the NCBI database using RAP search and an in-house perl script. However, 14,856 (37.18%) genes remained unannotated in the NCBI database. In comparison with closely related species, the genome size of <italic>Anguina tritici</italic> (164 Mb) falls between that of <italic>Ditylenchus dipsaci</italic> (227.2 Mb) (<a ref-type="bibr" href="#j_jofnem-2023-0031_ref_016">Mimee et al., 2019</a>) and <italic>Ditylenchus destructor</italic> (112 Mb) (<a ref-type="bibr" href="#j_jofnem-2023-0031_ref_019">Zheng et al., 2016</a>) indicating that <italic>Anguina tritici</italic> has a moderately sized genome in comparison to the other two species. The GC content, shows slight variations among the three species. <italic>Anguina tritici</italic> has the highest GC content, indicating a relatively higher abundance of G and C nucleotides in its genome compared to the <italic>Ditylenchus dipsaci</italic> (37.5) and <italic>Ditylenchus destructor</italic> (36.6). <italic>Anguina tritici</italic> has the highest number of coding genes (39,965) among the three species, indicating a more complex and diverse genetic repertoire. <italic>Ditylenchus dipsaci</italic> has an intermediate number of genes (26,428), while <italic>Ditylenchus destructor</italic> has the fewest genes (13,938). These variations may reflect evolutionary adaptations and functional differences between these nematode species.</p>
<table-wrap id="j_jofnem-2023-0031_tab_001" position="float"><label>Table 1.</label>
<caption><p>Characteristics of <italic>A. tritici</italic> genome.</p></caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th align="left" valign="middle"><bold>Number of contigs</bold></th>
<th align="left" valign="middle"><bold>843604</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td align="left" valign="top"><bold>Assembled genome size</bold></td>
<td align="left" valign="top">164398910 (164 Mb)</td>
</tr>
<tr>
<td align="left" valign="top"><bold>N50</bold></td>
<td align="left" valign="top">3232</td>
</tr>
<tr>
<td align="left" valign="top"><bold>Number of gene models</bold></td>
<td align="left" valign="top">39965</td>
</tr>
<tr>
<td align="left" valign="top"><bold>Median intron, bp</bold></td>
<td align="left" valign="top">65</td>
</tr>
<tr>
<td align="left" valign="top"><bold>G+C %</bold></td>
<td align="left" valign="top">39.1%</td>
</tr>
<tr>
<td align="left" valign="top"><bold>Gene density</bold></td>
<td align="left" valign="top">26.5784</td>
</tr>
<tr>
<td align="left" valign="top"><bold>Average gene size (bp)</bold></td>
<td align="left" valign="top">1347</td>
</tr>
<tr>
<td align="left" valign="top"><bold>Average CDS length (bp)</bold></td>
<td align="left" valign="top">168.88</td>
</tr>
<tr>
<td align="left" valign="top"><bold>Average number of exon/genes</bold></td>
<td align="left" valign="top">3.18</td>
</tr>
<tr>
<td align="left" valign="top"><bold>Average intron length (bp)</bold></td>
<td align="left" valign="top">280</td>
</tr>
<tr>
<td align="left" valign="top"><bold>BUSCO</bold></td>
<td align="left" valign="top">C: 36.6% [S: 35.5%, D: 1.1%], F: 14.8%, M: 48.6%, n: 982</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>Broadly, the genome size of <italic>A. tritici</italic> was determined to be 164 Mb, which is larger than most of the published genomes of PPNs, except for <italic>Meloidogyne arenaria</italic>, <italic>Rotylenchulus reniformis</italic>, and which are 284 Mb and 380 Mb, respectively (<a ref-type="bibr" href="#j_jofnem-2023-0031_ref_013">Nyaku et al., 2014</a>; <a ref-type="bibr" href="#j_jofnem-2023-0031_ref_012">Sato et al., 2018</a>). The GC content was measured to be 39.1%, which is higher than most PPNs, except for <italic>Bursaphelenchus xylophilus</italic> (<a ref-type="bibr" href="#j_jofnem-2023-0031_ref_014">Kikuchi et al., 2011</a>). From the <italic>A. tritici</italic> genome, a total of 39,965 protein-coding genes were predicted. This number is higher than most of the PPN genomes sequenced, except for <italic>Meloidogyne javanica</italic>, <italic>Meloidogyne floridensis</italic>, and <italic>Meloidogyne arenaria</italic> (<a ref-type="bibr" href="#j_jofnem-2023-0031_ref_017">Blanc-Mathieu et al., 2017</a>). On an average, three exons per gene were observed, compared to two exons per gene in <italic>Pratylenchus coffeae</italic> (<a ref-type="bibr" href="#j_jofnem-2023-0031_ref_018">Burke et al., 2015</a>) and four in most PPNs. The <italic>A. tritici</italic> genome exhibits characteristics of a compact parasitic genome, similar to other nematodes, with a low number of small introns and a high number of repeats. Out of the 41,746 predicted transcripts, 62% were successfully annotated with functional annotations. These unannotated genes may possess unique functions in the <italic>A. tritici</italic> genome since they do not show significant matches in the NCBI database. Their validation could shed light on their species-specific roles. The <italic>A. tritici</italic> genome will be valuable in unravelling the mechanisms involved in anhydrobiosis and the genes responsible for parasitizing the seed parts of plants, an area of nematology that has been less explored. Integrating RNAseq data in the annotation process could enhance gene prediction and annotation accuracy, thereby improving the overall quality and completeness of the final genome assembly. Future studies should consider incorporating RNAseq data to further enhance the understanding of the <italic>A. tritici</italic> genome.</p>
</div>

Nematode Genome Announcement: A Draft Genome of Seed Gall Nematode, <italic>Anguina tritici</italic>

Journal of Nematology

This work is licensed under the Creative Commons Attribution 4.0 International License.

Nematode Genome Announcement: A Draft Genome of Seed Gall Nematode, 

Anguina tritici is the first plant-parasitic nematode described in literature, dating back to the year 1743. It is responsible for causing earcockle (seed...

Nematode Genome Announcement: A Draft Genome of Seed Gall Nematode, Anguina tritici

Article Category: Research Paper

Data publikacji: 15 lis 2023

Zakres stron: -

Otrzymano: 08 gru 2022

DOI: https://doi.org/10.2478/jofnem-2023-0031

Słowa kluczowe
earcockle, tundu disease, genome sequence, Illumina MiSeq

© 2023 Ashish Kumar Singh et al., published by Sciendo

This work is licensed under the Creative Commons Attribution 4.0 International License.

Characteristics of A. tritici genome.

Number of contigs	843604
Assembled genome size	164398910 (164 Mb)
N50	3232
Number of gene models	39965
Median intron, bp	65
G+C %	39.1%
Gene density	26.5784
Average gene size (bp)	1347
Average CDS length (bp)	168.88
Average number of exon/genes	3.18
Average intron length (bp)	280
BUSCO	C: 36.6% [S: 35.5%, D: 1.1%], F: 14.8%, M: 48.6%, n: 982

Nematode Genome Announcement: A Draft Genome of Seed Gall Nematode, Anguina tritici

Article Category: Research Paper

Data publikacji: 15 lis 2023

Zakres stron: -

Otrzymano: 08 gru 2022

DOI: https://doi.org/10.2478/jofnem-2023-0031

Słowa kluczoweearcockle, tundu disease, genome sequence, Illumina MiSeq

© 2023 Ashish Kumar Singh et al., published by Sciendo

This work is licensed under the Creative Commons Attribution 4.0 International License.

Characteristics of A. tritici genome.

Słowa kluczowe
earcockle, tundu disease, genome sequence, Illumina MiSeq