Characterization of the two intra-individual sequence variants in the 18S rRNA gene in the plant parasitic nematode, Rotylenchulus reniformis

Evolution of myxozoan mitochondrial genomes: insights from myxobolids

BACKGROUND

Myxozoa is a class of cnidarian parasites that encompasses over 2,400 species. Phylogenetic relationships among myxozoans remain highly debated, owing to both a lack of informative morphological characters and a shortage of molecular markers. Mitochondrial (mt) genomes are a common marker in phylogeny and biogeography. However, only five complete myxozoan mt genomes have been sequenced: four belonging to two closely related genera, Enteromyxum and Kudoa, and one from the genus Myxobolus. Interestingly, while cytochrome oxidase genes could be identified in Enteromyxum and Kudoa, no such genes were found in Myxobolus squamalis, and another member of the Myxobolidae (Henneguya salminicola) was found to have lost its entire mt genome. To evaluate the utility of mt genomes to reconstruct myxozoan relationships and to understand if the loss of cytochrome oxidase genes is a characteristic of myxobolids, we sequenced the mt genome of five myxozoans (Myxobolus wulii, M. honghuensis, M. shantungensis, Thelohanellus kitauei and, Sphaeromyxa zaharoni) using Illumina and Oxford Nanopore platforms.

RESULTS

Unlike Enteromyxum, which possesses a partitioned mt genome, the five mt genomes were encoded on single circular chromosomes. An mt plasmid was found in M. wulii, as described previously in Kudoa iwatai. In all new myxozoan genomes, five protein-coding genes (cob, cox1, cox2, nad1, and nad5) and two rRNAs (rnl and rns) were recognized, but no tRNA. We found that Myxobolus and Thelohanellus species shared unidentified reading frames, supporting the view that these mt open reading frames are functional. Our phylogenetic reconstructions based on the five conserved mt genes agree with previously published trees based on the 18S rRNA gene.

CONCLUSIONS

Our results suggest that the loss of cytochrome oxidase genes is not a characteristic of all myxobolids, the ancestral myxozoan mt genome was likely encoded on a single circular chromosome, and mt plasmids exist in a few lineages. Our findings indicate that myxozoan mt sequences are poor markers for reconstructing myxozoan phylogenetic relationships because of their fast-evolutionary rates and the abundance of repeated elements, which complicates assembly.

Wild Caught Nematode Identification and Early Embryo Development: An accessible undergraduate research experience

Ample evidence suggests that participation in undergraduate research in community college is critical for stimulating interest and retention in STEM careers. Guided skill development and practice in a collaborative lab setting allows students to be competitive when applying to future research opportunities. The goals of this undergraduate research experience (URE) was for student-driven discovery with unknown outcomes including: introduction to primary literature, developmental biology, developing hypotheses, learning worm maintenance, microscopy, PCR, and sequencing analysis. The use of C. elegans and wild caught nematodes facilitated an exciting and affordable project that can be built on in future UREs.

Characterizing Repeats in Two Whole-Genome Amplification Methods in the Reniform Nematode Genome

One of the major problems in the U.S. and global cotton production is the damage caused by the reniform nematode, Rotylenchulus reniformis. Amplification of DNA from single nematodes for further molecular analysis can be challenging sometimes. In this research, two whole-genome amplification (WGA) methods were evaluated for their efficiencies in DNA amplification from a single reniform nematode. The WGA was carried out using both REPLI-g Mini and Midi kits, and the GenomePlex single cell whole-genome amplification kit. Sequence analysis produced 4 Mb and 12 Mb of genomic sequences for the reniform nematode using REPLI-g and SIGMA libraries. These sequences were assembled into 28,784 and 24,508 contigs, respectively, for REPLI-g and SIGMA libraries. The highest repeats in both libraries were of low complexity, and the lowest for the REPLI-g library were for satellites and for the SIGMA library, RTE/BOV-B. The same kind of repeats were observed for both libraries; however, the SIGMA library had four other repeat elements (Penelope (long interspersed nucleotide element (LINE)), RTE/BOV-B (LINE), PiggyBac, and Mirage/P-element/Transib), which were not seen in the REPLI-g library. DNA transposons were also found in both libraries. Both reniform nematode 18S rRNA variants (RN_VAR1 and RN_VAR2) could easily be identified in both libraries. This research has therefore demonstrated the ability of using both WGA methods, in amplification of gDNA isolated from single reniform nematodes.

Plant-parasitic nematodes associated with sugarcane in Kilimanjaro, Tanzania

Morphological and molecular analyses of plant-parasitic nematodes (PPN) from 12 sugarcane plantation sites of Tanganyika Planting Company (TPC) Limited in Kilimanjaro region of Tanzania revealed the presence of six PPN genera, i.e. Helicotylenchus, Hemicycliophora, Pratylenchus, Rotylenchulus, Scutellonema, and Tylenchorhynchus. The genera with the highest densities and present in virtually all samples were Pratylenchus and Rotylenchulus, and the most important species appeared to be R. parvus, P. zeae, T. crassicaudatus, and T. ventrosignatus. A total sequences of 11 partial ITS, 15 D2-D3 of 28S, and 6 partial 18S of rRNA gene, and 7 partial COI gene of mtDNA of these species were obtained in this study. Morphology and molecular data comparisons between the Tanzanian R. parvus and the South African R. parvus indicated that R. parvus is a cryptic species complex. Based on the results of morphological and molecular analyses of T. crassicaudatus and T. agri from China, Haiti, Indonesia, Iran, Niger and the USA, T. agri syn. n. is proposed as a junior synonym of T. crassicaudatus.

Widespread prevalence but contrasting patterns of intragenomic rRNA polymorphisms in nematodes: Implications for phylogeny, species delimitation and life history inference

Ribosomal RNA genes have long been a favoured locus in phylogenetic and metabarcoding studies. Within a genome, rRNA loci are organized as tandem repeated arrays and the copies are homogenized through the process of concerted evolution. However, some level of rRNA variation (intragenomic polymorphism) is known to persist and be maintained in the genomes of many species. In nematode worms, the extent of rRNA polymorphism (RP) across species and the evolutionary and life history factors that contribute to the maintenance of intragenomic RP is largely unknown. Here, we present an extensive analysis across 30 terrestrial nematode species representing a range of free-living and parasitic taxa isolated worldwide. Our results indicate that RP is common and widespread, ribosome function appears to be maintained despite mutational changes, and intragenomic variants are stable in the genome and neutrally evolving. However, levels of variation were varied widely across rRNA locus and species, with some taxa observed to lack RP entirely. Higher levels of RP were significantly correlated with shorter generation time and high reproductive rates, and population-level factors may play a role in the geographic and phylogenetic structuring of rRNA variants observed in genera such as Rotylenchulus and Pratylenchus. Although RP did not dramatically impact the clustering and recovery of taxa in mock metabarcoding analyses, the present study has significant implications for global biodiversity estimates of nematode species derived from environmental rRNA amplicon studies, as well as our understanding of the evolutionary and ecological factors shaping genetic diversity across the nematode Tree of Life.

Evidence of Intra-individual SSU Polymorphisms in Dark-spored Myxomycetes (Amoebozoa)

The nuclear small subunit rRNA gene (SSU or 18S) is a marker frequently used in phylogenetic and barcoding studies in Amoebozoa, including Myxomycetes. Despite its common usage and the confirmed existence of divergent copies of ribosomal genes in other protists, the potential presence of intra-individual SSU variability in Myxomycetes has never been studied before. Here we investigated the pattern of nucleotide polymorphism in the 5' end fragment of SSU by cloning and sequencing a total of 238 variants from eight specimens, each representing a species of the dark-spored orders Stemonitidales and Physarales. After excluding singletons, a relatively low SSU intra-individual variability was found but our data indicate that this might be a widely distributed phenomenon in Myxomycetes as all samples analyzed possessed various ribotypes. To determine if the occurrence of multiple SSU variants within a single specimen has a negative effect on the circumscription of species boundaries, we conducted phylogenetic analyses that revealed that clone variation may be detrimental for inferring phylogenetic relationships among some of the specimens analyzed. Despite that intra-individual variability should be assessed in additional taxa, our results indicate that special care should be taken for species identification when working with closely related species.

Genotypic Characterization of Torymus sinensis (Hymenoptera: Torymidae) After Its Introduction in Tuscany (Italy) for the Biological Control of Dryocosmus kuriphilus (Hymenoptera: Cynipidae)

Torymus sinensis Kamijo (Hymenoptera: Torymidae) is an alien parasitoid that is used in many areas of the world for biological control the Asian chestnut gall wasp, Dryocosmus kuriphilus Yasumatsu (Hymenoptera: Cynipidae). In Italy, this parasitoid was imported from Japan in 2003 and subsequently multiplied and released throughout the country. In this study, a phylogenetic investigation was carried out on insects from three different sites in northern Tuscany (Italy). Moreover, the possible hybridization between T. sinensis and some native Torymus species was evaluated. The conserved region 18S rRNA gene and the hypervariable ITS2 (Internal Transcribed Spacer 2) region of the ribosomal cistrone were selected as molecular markers. Sequencing the amplified products, after cloning, ruled out any hybridization between T. sinensis and the native Torymus species, and also confirmed the presence of two haplotypes for the Tuscan population of T. sinensis both for the region of the 18S rRNA gene as well as for the ITS2 region. These results confirm that the environmental impact of the alien parasitoid T. sinensis in the study site is acceptable, although an extensive and repeated monitoring would be desirable.

Characterization of a non-sexual population of Strongyloides stercoralis with hybrid 18S rDNA haplotypes in Guangxi, Southern China

Strongyloidiasis is a much-neglected but sometimes fatal soil born helminthiasis. The causing agent, the small intestinal parasitic nematode Strongyloides stercoralis can reproduce sexually through the indirect/heterogonic life cycle, or asexually through the auto-infective or the direct/homogonic life cycles. Usually, among the progeny of the parasitic females both, parthenogenetic parasitic (females only) and sexual free-living (females and males) individuals, are present simultaneously. We isolated S. stercoralis from people living in a village with a high incidence of parasitic helminths, in particular liver flukes (Clonorchis sinensis) and hookworms, in the southern Chinese province Guangxi. We determined nuclear and mitochondrial DNA sequences of individual S. stercoralis isolated from this village and from close by hospitals and we compared these S. stercoralis among themselves and with selected published S. stercoralis from other geographic locations. For comparison, we also analyzed the hookworms present in the same location. We found that, compared to earlier studies of S. stercoralis populations in South East Asia, all S. stercoralis sampled in our study area were very closely related, suggesting a recent common source of infection for all patients. In contrast, the hookworms from the same location, while all belonging to the species Necator americanus, showed rather extensive genetic diversity even within host individuals. Different from earlier studies conducted in other geographic locations, almost all S. stercoralis in this study appeared heterozygous for different sequence variants of the 18S rDNA hypervariable regions (HVR) I and IV. In contrast to earlier investigations, except for three males, all S. stercoralis we isolated in this study were infective larvae, suggesting that the sampled population reproduces predominantly, if not exclusively through the clonal life cycles. Consistently, whole genome sequencing of individual worms revealed higher heterozygosity than reported earlier for likely sexual populations of S. stercoralis. Elevated heterozygosity is frequently associated with asexual clonal reproduction.

A Rapid Diagnostic for Detection of Aphelenchoides besseyi and A. fujianensis Based on Real-Time PCR

Aphelenchoides besseyi and A. fujianensis have been frequently found in mixed populations associated with forage grass seed in Brazil. The morphological similarity between both species has previously led A. fujianensis to be erroneously identified as A. besseyi. A. besseyi is a quarantine pest in many countries that import Brazilian forage seed; however, there is no current evidence suggesting that A. fujianensis is a plant-parasitic species. Two real-time polymerase chain reaction (qPCR) diagnostics were developed to detect each species and an operational envelope was established. A set of primers and hydrolysis probes for each species was designed targeting the large subunit (LSU) region. To assess their specificity, primers and probes sets were tested with samples of nontarget Aphelenchoides and Paraphelenchus sp. also frequently associated with forage seed. Experiments using dilutions of purified plasmid standards underpinned the sensitivity of the qPCR assays, which detected as few as 10 copies of target nematode ribosomal DNA. Thus, the developed diagnostics were sufficiently sensitive to detect DNA extracted from a fragment of a single target nematode. There was a positive correlation between copy number of the target species and nematode abundance, suggesting the potential of this method for quantification. Evidence of intra-individual variability among cloned sequences of the LSU region in a single A. besseyi population is also reported.

Principal Component Analysis and Molecular Characterization of Reniform Nematode Populations in Alabama

U.S. cotton production is suffering from the yield loss caused by the reniform nematode (RN), Rotylenchulus reniformis. Management of this devastating pest is of utmost importance because, no upland cotton cultivar exhibits adequate resistance to RN. Nine populations of RN from distinct regions in Alabama and one population from Mississippi were studied and thirteen morphometric features were measured on 20 male and 20 female nematodes from each population. Highly correlated variables (positive) in female and male RN morphometric parameters were observed for body length (L) and distance of vulva from the lip region (V) (r = 0.7) and tail length (TL) and c' (r = 0.8), respectively. The first and second principal components for the female and male populations showed distinct clustering into three groups. These results show pattern of sub-groups within the RN populations in Alabama. A one-way ANOVA on female and male RN populations showed significant differences (p ≤ 0.05) among the variables. Multiple sequence alignment (MSA) of 18S rRNA sequences (421) showed lengths of 653 bp. Sites within the aligned sequences were conserved (53%), parsimony-informative (17%), singletons (28%), and indels (2%), respectively. Neighbor-Joining analysis showed intra and inter-nematodal variations within the populations as clone sequences from different nematodes irrespective of the sex of nematode isolate clustered together. Morphologically, the three groups (I, II and III) could not be distinctly associated with the molecular data from the 18S rRNA sequences. The three groups may be identified as being non-geographically contiguous.

Contrasting evolutionary patterns of 28S and ITS rRNA genes reveal high intragenomic variation in Cephalenchus (Nematoda): Implications for species delimitation

Concerted evolution is often assumed to be the evolutionary force driving multi-family genes, including those from ribosomal DNA (rDNA) repeat, to complete homogenization within a species, although cases of non-concerted evolution have been also documented. In this study, sequence variation of 28S and ITS ribosomal RNA (rRNA) genes in the genus Cephalenchus is assessed at three different levels, intragenomic, intraspecific, and interspecific. The findings suggest that not all Cephalenchus species undergo concerted evolution. High levels of intraspecific polymorphism, mostly due to intragenomic variation, are found in Cephalenchus sp1 (BRA-01). Secondary structure analyses of both rRNA genes and across different species show a similar substitution pattern, including mostly compensatory (CBC) and semi-compensatory (SBC) base changes, thus suggesting the functionality of these rRNA copies despite the variation found in some species. This view is also supported by low sequence variation in the 5.8S gene in relation to the flanking ITS-1 and ITS-2 as well as by the existence of conserved motifs in the former gene. It is suggested that potential cross-fertilization in some Cephalenchus species, based on inspection of female reproductive system, might contribute to both intragenomic and intraspecific polymorphism of their rRNA genes. These results reinforce the potential implications of intragenomic and intraspecific genetic diversity on species delimitation, especially in biodiversity studies based solely on metagenetic approaches. Knowledge of sequence variation will be crucial for accurate species diversity estimation using molecular methods.

Morphological and molecular characterisation of one new and several known species of the reniform nematode, Rotylenchulus Linford & Oliveira, 1940 (Hoplolaimidae: Rotylenchulinae), and a phylogeny of the genus

The reniform nematodes of the genusRotylenchulusare semi-endoparasites of numerous herbaceous and woody plant roots and are mainly distributed in tropical and subtropical regions. In this study, we provide morphological and molecular characterisation of six out of ten presently known valid species ofRotylenchulus:R. clavicaudatus,R. leptus,R. macrodoratus,R. macrosoma,R. reniformisandR. saccharifrom South Africa, USA, Italy and Spain.Rotylenchulus parvuswas only studied morphologically. A new species,R. macrosomoidessp. n., isolated from soil and roots of sugarcane in South Africa, is described. The phylogeny ofRotylenchulus, as inferred from the analyses of D2-D3 of 28S rRNA, ITS rRNA,coxImtDNA andhsp90gene sequences, is presented. The study revealed thatR. reniformisandR. macrosomahave a sister relationship, but that relationships between otherRotylenchulusspecies remain unresolved. The phylogenetic analysis also confirmed the hypothesis that this genus originated from the Afrotropical zoogeographical region. Our study revealed thatR. reniformisandR. macrosomoidessp. n. have two distinct rRNA gene types andR. macrosomahave three rRNA gene types in their genomes. PCR with species-specific primers was developed for rapid diagnostics ofR. reniformis.

Divergence thresholds and divergent biodiversity estimates: can metabarcoding reliably describe zooplankton communities?

DNA metabarcoding is a promising method for describing communities and estimating biodiversity. This approach uses high-throughput sequencing of targeted markers to identify species in a complex sample. By convention, sequences are clustered at a predefined sequence divergence threshold (often 3%) into operational taxonomic units (OTUs) that serve as a proxy for species. However, variable levels of interspecific marker variation across taxonomic groups make clustering sequences from a phylogenetically diverse dataset into OTUs at a uniform threshold problematic. In this study, we use mock zooplankton communities to evaluate the accuracy of species richness estimates when following conventional protocols to cluster hypervariable sequences of the V4 region of the small subunit ribosomal RNA gene (18S) into OTUs. By including individually tagged single specimens and “populations” of various species in our communities, we examine the impact of intra- and interspecific diversity on OTU clustering. Communities consisting of single individuals per species generated a correspondence of 59–84% between OTU number and species richness at a 3% divergence threshold. However, when multiple individuals per species were included, the correspondence between OTU number and species richness dropped to 31–63%. Our results suggest that intraspecific variation in this marker can often exceed 3%, such that a single species does not always correspond to one OTU. We advocate the need to apply group-specific divergence thresholds when analyzing complex and taxonomically diverse communities, but also encourage the development of additional filtering steps that allow identification of artifactual rRNA gene sequences or pseudogenes that may generate spurious OTUs.

Characterization of the reniform nematode genome by shotgun sequencing

The reniform nematode (RN), a major agricultural pest particularly on cotton in the United States, is among the major plant-parasitic nematodes for which limited genomic information exists. In this study, over 380 Mb of sequence data were generated from pooled DNA of four adult female RNs and assembled into 67,317 contigs, including 25,904 (38.5%) predicted coding contigs and 41,413 (61.5%) noncoding contigs. Most of the characterized repeats were of low complexity (88.9%), and 0.9% of the contigs matched with 53.2% of GenBank ESTs. The most frequent Gene Ontology (GO) terms for molecular function and biological process were protein binding (32%) and embryonic development (20%). Further analysis showed that 741 (1.1%), 94 (0.1%), and 169 (0.25%) RN genomic contigs matched with 1328 (13.9%), 1480 (5.4%), and 1330 (7.4%) supercontigs of Meloidogyne incognita, Brugia malayi, and Pristionchus pacificus, respectively. Chromosome 5 of Caenorhabditis elegans had the highest number of hits to the RN contigs. Seven putative detoxification genes and three carbohydrate-active enzymes (CAZymes) involved in cell wall degradation were studied in more detail. Additionally, kinases, G protein-coupled receptors, and neuropeptides functioning in physiological, developmental, and regulatory processes were identified in the RN genome.

Phylogenetic analysis of entomoparasitic nematodes, potential control agents of flea populations in natural foci of plague

Entomoparasitic nematodes are natural control agents for many insect pests, including fleas that transmit Yersinia pestis, a causative agent of plague, in the natural foci of this extremely dangerous zoonosis. We examined the flea samples from the Volga-Ural natural focus of plague for their infestation with nematodes. Among the six flea species feeding on different rodent hosts (Citellus pygmaeus, Microtus socialis, and Allactaga major), the rate of infestation varied from 0 to 21%. The propagation rate of parasitic nematodes in the haemocoel of infected fleas was very high; in some cases, we observed up to 1,000 juveniles per flea specimen. Our study of morphology, life cycle, and rDNA sequences of these parasites revealed that they belong to three distinct species differing in the host specificity. On SSU and LSU rRNA phylogenies, these species representing three genera (Rubzovinema, Psyllotylenchus, and Spilotylenchus), constitute a monophyletic group close to Allantonema and Parasitylenchus, the type genera of the families Allantonematidae and Parasitylenchidae (Nematoda: Tylenchida). We discuss the SSU-ITS1-5.8S-LSU rDNA phylogeny of the Tylenchida with a special emphasis on the suborder Hexatylina.