Literature DB >> 23593343

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

Seloame T Nyaku1, Venkateswara R Sripathi, Ramesh V Kantety, Yong Q Gu, Kathy Lawrence, Govind C Sharma.   

Abstract

The 18S rRNA gene is fundamental to cellular and organismal protein synthesis and because of its stable persistence through generations it is also used in phylogenetic analysis among taxa. Sequence variation in this gene within a single species is rare, but it has been observed in few metazoan organisms. More frequently it has mostly been reported in the non-transcribed spacer region. Here, we have identified two sequence variants within the near full coding region of 18S rRNA gene from a single reniform nematode (RN) Rotylenchulus reniformis labeled as reniform nematode variant 1 (RN_VAR1) and variant 2 (RN_VAR2). All sequences from three of the four isolates had both RN variants in their sequences; however, isolate 13B had only RN variant 2 sequence. Specific variable base sites (96 or 5.5%) were found within the 18S rRNA gene that can clearly distinguish the two 18S rDNA variants of RN, in 11 (25.0%) and 33 (75.0%) of the 44 RN clones, for RN_VAR1 and RN_VAR2, respectively. Neighbor-joining trees show that the RN_VAR1 is very similar to the previously existing R. reniformis sequence in GenBank, while the RN_VAR2 sequence is more divergent. This is the first report of the identification of two major variants of the 18S rRNA gene in the same single RN, and documents the specific base variation between the two variants, and hypothesizes on simultaneous co-existence of these two variants for this gene.

Entities:  

Mesh:

Substances:

Year:  2013        PMID: 23593343      PMCID: PMC3623918          DOI: 10.1371/journal.pone.0060891

Source DB:  PubMed          Journal:  PLoS One        ISSN: 1932-6203            Impact factor:   3.240


Introduction

The ribosomal DNA (rDNA) plays a pivotal role in protein synthesis in eukaryotes and changes in these genes can profoundly affect ecological interactions, host range, trophic production, and thus the overall growth and resource requirement of the organism [1]. Therefore, rRNA variation has effects far beyond cell biology alone. The rDNA region is composed of numerous copies of tandemly repeated transcription units within the genome. It consists of the 18S, 5.8S, and 28S genes as well as internal and external transcribed spacers (ITS and ETS) [2]. Concerted evolution is implicated in the conservative nature of the repeats [3], [4], through molecular mechanisms such as gene conversion [4], [5] and unequal crossing-over [6]. Since various regions of this gene evolve at different rates [2], the sequence homology provides comparative phylogenetic differences between or within organisms. Intra and inter variations within the 18S rRNA gene sequences of species have been observed in genomes of some organisms [7], [8], [9] often attributed to frequent recombination events involving unequal crossovers and gene conversions producing rDNA units having similar sequences [9]. The RN is endemic to the southern U.S. and in a number of tropical and subtropical regions of the world [10], including Africa, Asia and Australia [11]. The absence or infrequent usage of non-hosts, such as corn (Zea mays) in crop rotations with cotton (Gossypium hirsutum), aggravates the incidence of RN. This organism also has a diverse host range that includes nearly 300 plant species [10], [12] and possesses the ability to live in prolonged dry and harsh environments [12]. This anhydrobiosis feature facilitates its dispersal to distant areas by dust storms [13]. Considerable damage is caused by the RN and this reduces cotton production in Alabama [14], Louisiana [15], and Mississippi [16] where in 2006 approximately 328,073 bales of cotton or $128 million in revenue were lost. Monoculture is the main reason why RN populations keep increasing [17], [18]. The rDNA serves as a key molecular marker for establishing phylogenetic relationships. Within a number of organisms, the rRNA genes can provide the evolutionary history [19]. This gene has variable regions which have been used to elucidate molecular phylogenies in organisms including those of free-living, animal, and plant parasitic nematodes [20]–[22]. In choosing a molecular marker or region for phylogenetic analysis, considerable care is required to accurately reflect the true evolutionary relationships among taxa. The 18S rDNA region, because of its conservative nature, has been a staple in such phylogenic analyses within and between various taxa [23]. The presence of more than one variant therefore, is an infrequent rarity and demands further investigation. The objectives of this study were to: 1. Elucidate the full-length 18S rRNA gene nucleotide sequence of a single RN, and 2. Characterize the 18S rRNA gene for occurrence of any potential variation in this gene.

Results

Amplification of RN DNA using primer pairs SSUF05 and SSUR81 gave a band of approximately 1800 bp and 49 clones were sequenced from four female RNs. Sequencing of the 1800 bp band was accomplished by using a combination of six primers ( ). Multiple sequence alignment (MSA) was performed using ClustalW (http://align.genome.jp/ ) with all the clone sequences generated from each nematode using default parameters and viewed using Bioedit software [24]. The MSA table is provided as an additional file (). Two types of 18S rDNA sequences were observed from MSA, these were labeled as RN_VAR1 and RN_VAR2 referred to as the first and second variants of the 18S rDNA sequence of the RN, respectively. Five of the clones (SSU1B2, SSU1B3, SSU1B7, SSU1B9, and SSU1B10) were identified as putative chimeras because these had their preference scores above 1.0. These clones were therefore not used for further analysis. The 44 clone sequences (100%) were made up of 5 (11.3%), 12 (27.3%), 19 (43.2%), and 8 (18.2%) clones from individual nematodes 1, 2, 3, and 4 designated as SSU1B, SSU12A, SSU13B, and SSU25A, respectively. Reniform nematode variant 1 (RN_VAR1) consisted of 11 (25.0%) clones out of the 44 clones sequenced and these were very similar in terms of variable sites, and consisted of 1, 4, 0, and 6 clones for SSU1B, SSU12A, SSU13B, and SSU25A, respectively. A majority of clones sequenced (33 or 75%) were classified as reniform variant 2 (RN_VAR2) sequences. These 33 clone sequences (100%) consisted of 4, 8, 19, and 2 clones for SSU1B, SSU12A, SSU13B, and SSU25A, respectively. Comparisons of the consensus sequences from the variant clones of the RN to GenBank gave varying hits to nematode sequences. Variant 1 (RN_VAR1) sequence had the top-most hit to Rotylenchulus reniformis isolate wb2 18S ribosomal RNA gene (GenBank accession number EU306342) with a 98% similarity, and a bit score of 3,125. Variant 2 (RN_VAR2) sequence had a top-most hit to RN_VAR1, followed by Heterodera schachtii isolate wb15 (GenBank accession number EU306355), with a 95% similarity, bit score was 2,791. The MSA analysis identified 96 (5.5%) nucleotide sites that distinguish the two 18S rRNA variants ( and ). These positions were dispersed widely between 95 to 1673 bp within the 18S rRNA gene. Seven distinct sites (0.4%) had indels within the consensus sequence at base positions 155, 156, 769, 1,330, 1,331, 1,633 and 1,634. The number of conserved regions among the two variants of the RNs was determined using default parameters [24]. Twenty-six conserved regions were present among both RN variants ( ). Regions 774–980 bp and 1,084–1283 bp had the highest stretch of conserved bases of 207 and 200, respectively.
Table 1

Oligonucleotide primers for sequencing of the RN 18S rRNA gene.

Oligo namePrimer sequence (5′-3′)Direction
SSUF05GAAACTGCGAA(CT)GGCTCATTForward
SSUR81TGATCC(AT) (GT)C(CT)GCAGGTTCACReverse
SSUF22 TCCAAGGAAGGCAGCAGGC Forward
SSUR13 GGGCATCACAGACCTGTTA Reverse
T3 ATTAACCCTCACTAAAGGGA Reverse
T7 TAATACGACTCACTATAGGG Forward
Table 2

Variable sites identified among reniform nematode (RN) 18S rRNA variants.

Variations at specific base locations in 18S clones (15) of rRNA variant 1 (RN_VAR1)
96105106110119123152154155156165253265294313316331334379408409410417418419
SSU1B8AGCTATCCATTGCCGCGCTGCTGAT
SSU12A4AGCTATCCATTGCCGCGCTAGTGAT
SSU12A9AGCTATCCATTGCCGCGCTAGTGAT
SSU12A10AGCTATCCATTGCCGCGCTAGTGAT
SSU12A12AGCTATCCATTGCCGCGCTAGTGAT
SSU25A1AGCTATCCATTGCCGCGCTAGTGAT
SSU25A4AGCTATCCATTGCCGCGCTAGTGAT
SSU25A5AGCTATCCATTGCCGCGCTAGTGAT
SSU25A6AGCTATCCATTGCCGCGCTGCTGAT
SSU25A8AGCTATCCATTGCCGCGCTAGTGAT
SSU25A10AGCTATCCA-TGCCGCGCTACTGAT
CONSENSUS A G C T A T C C A T T G C C G C G C T A G T G A T
Variations at specific base locations in 18S clones (34) of rRNA variant 2 (RN_VAR2)
SSU1B1CAAAGCTAAATAAGCTGACCATC
SSU1B4CAACGCTAAATAAGCTGGCCATC
SSU1B5CAACGCTAAATAAGCTGACCATC
SSU1B6CAACGCTAAATAAGCTGACCATC
SSU12A1CAACGCTAAATAAGCTGACCATC
SSU12A2CAACGCTAAATAAGCTGACCATC
SSU12A3CAACGCTAAATAAGCTGACCATC
SSU12A5CAACGCTAAATAAGCTGACCATC
SSU12A6CAACGCTAAATAAGCTGACCATC
SSU12A7CAACGCTAAATAAGCTGACCATC
SSU12A8CAACGCTAAATAAGCTGACCATC
SSU12A13CAACGCTAAATAAGCTGACCATC
SSU25A11CAACGCTAAATAAGCTGACCATC
SSU25A12CAACGCTAAATAAGCTGACCATC
SSU13B1CAACGCTAAATAAGCTGACCATC
SSU13B2CAACGCTAAATAAGCTGACCATC
SSU13B3CAACGCTAAATAAGCTGACCATC
SSU13B4CAACGCTAAATAAGCTGACCATC
SSU13B5CAACGCTAAATAAGCTGACCATC
SSU13B6CAACGCTAAATAAGCTGACCATC
SSU13B7CAACGCTAAATAAGCTGACCATC
SSU13B9CAACGCTAAATAAGCTGACCATC
SSU13B10CAACGCTAAATAAGCTGACCATC
SSU13B11CAACGCTAAATAAGCTGACCATC
SSU13B12CAACGCTAAATAAGCTGACCATC
SSU13B13CAACGCTAAATAAGCTGACCATC
SSU13B14CAACGCTAAATAAGCTGACCATC
SSU13B15CAACGCTAAATAAGCTGACCATC
SSU13B16CAACGCTAAATAAGCTGACCATC
SSU13B17CAACGCTAAATAAGCTGACCATC
SSU13B18CAACGCTAAATAAGCTGACCATC
SSU13B19CAACGCTAAATAAGCTGACCATC
SSU13B20CAACGCTAAATAAGCTGACCATC
CONSENSUS C A A A G C T A A A T A A G C T G A C C A T C
Table 3

Variable sites identified among reniform nematode (RN) 18S rRNA variants.

Variations at specific base locations in 18S clones (15) of rRNA variant 1 (RN_VAR1)
671672673769981983984987992996128412861300133013311335141114121415142114281439
SSU1B8AGCTATCCTAGCTT-CGACGA
SSU12A4AGCTATCCTAGCTTACGACGA
SSU12A9AGCTATCCTAGCTTACGACGA
SSU12A10AGCTATCCTAGCTTACGACGA
SSU12A12AGCTATCCTAGCTTACGACGA
SSU25A1AGCTATCCTAGCTTACGACGA
SSU25A4AGCTATCCTAGCTTACGACGA
SSU25A5AGCTATCCTAGCTTACGACGA
SSU25A6AGCTATCCTAGCTTACGACGA
SSU25A8AGCTATCCTAGCTTACGACGA
SSU25A10AGCTATCCTAGCTTACGATGA
CONSENSUS A G C T A T C C T A G C T T A C G A C G A
Variations at specific base locations in 18S clones (34) of rRNA variant 2 (RN_VAR2)
SSU1B1GAG-GCGTTCGATGATGTCG
SSU1B4GAG-GCGTTCGATGATGTCG
SSU1B5GAG-GCGTTCGATGATGTCG
SSU1B6GAG-GCGTTCGATGATGTCG
SSU12A1GAG-GCGTTCGATGATGTCG
SSU12A2GAG-GCGTTCGATGATGCCG
SSU12A3GAG-GCGTTCGATGATGCCG
SSU12A5GAG-GCGTTCGATGATGTCG
SSU12A6GAG-GCGTTCGATGATGTCG
SSU12A7GAG-GCGTTCGATGATGTCG
SSU12A8GAG-GCGTTCGATGATGTCG
SSU12A13GAG-GCGTTCGATGATGTCG
SSU25A11GAG-GCGTTCGATGATGTCG
SSU25A12GAG-GCGTTCGATGATGTCG
SSU13B1GAG-GCGTTCGATGATGTCG
SSU13B2GAG-GCGTTCGATGATGTCG
SSU13B3GAG-GCGTTCGATGATGTCG
SSU13B4GAG-GCGTTCGATGATGTCG
SSU13B5GAG-GCGTTCGATGATGTCG
SSU13B6GAG-GCGTTCGATGATGTCG
SSU13B7GAG-GCGTTCGATGATGTCG
SSU13B9GAG-GCGTTCGATGATGTCG
SSU13B10GAG-GCGTTCGATGATGTCG
SSU13B11GAG-GCGTTCGATGATGTCG
SSU13B12GAG-GCGTTCGATGATGTCG
SSU13B13GAG-GCGTTCGATGATGTCG
SSU13B14GAG-GCGTTCGATGATGTCG
SSU13B15GAG-GCGTTCGATGATGTCG
SSU13B16GAG-GCGTTCGGTGATGTCG
SSU13B17GAG-GCGTTCGATGATGTCG
SSU13B18GAG-GCGTTCGATGATGTCG
SSU13B19GAG-GCGTTCGATGATGTCG
SSU13B20GAGGGCGTTCGATGATGTCG
CONSENSUS G A G - G C G T T C G A T G A T G T C G
Figure 1

Twenty-six conserved regions within the 18S rRNA gene of four juvenile female RNs.

Initial phylogenetic analysis was based on 44 RN clone sequences, RN_VAR1, RN_VAR2 consensus sequences, two sequences each from a taxa in the super-family Heteroderidae (Afenestrata koreana isolate wb17, and Globodera achilleae isolate 1008), Hoplolaimidae (Scutellonema bradys 1G and Rotylenchus goodeyi isolate RotyGoo), and one sequence each from Rotylenchulidae (Rotylenchulus reniformis isolate wb2), and Meloidogynidae super-family Meloidogyne incognita isolate MeloInc1 ( ). A neighbor-joining (NJ) tree was generated with bootstrap test of phylogeny with 10,000 replications, random seed of 64,238 model used was nucleotide: maximum composite likelihood, substitutions (d:transitions + transversions), pattern among lineages was homogeneous, rates among sites were uniform, Gaps/missing data- complete deletion. The NJ tree showed a major split between clones consisting of RN_VAR1, and RN_VAR2. Other nematodes belonging to the various super-families clustered together. R. reniformis isolate wb2 grouped with RN_VAR1 and M. incognita isolate MeloInc1 did not group with any other taxa ( ). Bootstrapping using the NJ clearly separated the clones (SSU12A, SSU25A, SSU1B, and SSU13B). The inclusion of the consensus sequences from these clones, thus (RN_VAR1 and RN_VAR2) in the NJ tree resulted in these clustering with their resultant clone sequences. Interestingly, clones from nematode isolate 13B were of RN_VAR2, however, clones from the other three nematodes were either of RN_VAR1 or RN_VAR2 and clustered as such. A second phylogram was generated using a NJ tree, with the same parameters as the first. Nematodes belonging to the super-families Heteroderidae, Hoplolaimidae, Rotylenchulidae and Meloidogynidae were used in the analysis ( ). The NJ tree showed taxa belonging to the same super-families being clustered together. The RN_VAR1 grouped with Rotylenchulus reniformis isolate wb2 with a bootstrap percentage of 100, and Meloidogyne incognita isolate MeloInc1, the out-group was not closely related to any other taxa ( ).
Figure 2

Neighbor-Joining (NJ) analyses for 44 RN 18S rRNA clones and nematodes in the super-families, Hoplolaimidae, Heteroderidae, Rotylenchulidae, and Meloidogynidae.

Table 4

Nematode species belonging to Tylenchoidea, GenBank accession numbers, and length of 18S rRNA sequence.

Super-familyNematodesGenBank No.Sequence length (bp)
Heteroderidae Heterodera avenae isolate 848FJ0404031,725
Heterodera betae isolate 973FJ0404041,725
Heterodera goettingiana strain HeDeGoe1EU6699151,725
Heterodera hordecalis isolate 1063FJ0404051,725
Heterodera mani strain HeDeMan1EU6699161,725
Heterodera schachtii isolate HeDeSchAY2846171,744
Heterodera schachtii isolate wb15 EU3063551,772
Heterodera trifolii isolate 846FJ0404021,724
Afenestrata koreana isolate wb17EU3063571,788
Globodera achilleae isolate 1008FJ0403991,726
Globodera artemisiae isolate 1022FJ0404001,726
Globodera pallida isolate GlobPal1AY2846181,746
Globodera rostochiensis strain GlobRos4AY5938801,728
Globodera tabacum isolate GlobTab4FJ0404011,726
Punctodera stonei EU6823911,727
Hoplolaimidae Rotylenchus goodeyi isolate RotyGooAY2846091,634
Rotylenchus robustus AJ9665031,796
Rotylenchus uniformis isolate wb16EU3063561,796
Scutellonema bradys 1GAY2717231,811
Scutellonema bradys 2GAJ9665041,803
Helicotylenchus Canadensis isolate HeliCanAY2846051,635
Helicotylenchus dihystera AJ9664861,790
Helicotylenchus pseudorobustus isolate HeliPseAY2846061,755
Helicotylenchus varicaudatus isolate wb14EU3063541,802
Helicotylenchus vulgaris isolate HeliVulAY2846071,751
Rotylenchulidae Rotylenchulus reniformis isolate wb2EU3063421,873
Meloidogynidae Meloidogyne incognita isolate MeloInc1AY2846211,730
Figure 3

Neighbor-Joining (NJ) analysis for RN variants (RN_VAR1 and RN_VAR2) and nematodes in the super-families, Hoplolaimidae, Heteroderidae, Rotylenchulidae, and Meloidogynidae.

BLASTN results of RN_VAR1 and RN_VAR2 to REPLIg genomic library gave hits of 83 (26%) and 236 (74%), respectively. The average alignment length and bit scores for both variants were 400.43 bp and 749.64 for RN_VAR1 and 441.40 bp and 832.79 for RN_VAR2, respectively. Comparison of our RN 18S rDNA variant sequences to our genomic sequences to determine the presence of these variants in the RN genome, showed RN_VAR2 with high unique hits, average % identity, and low mis-matches within the REPLIg genomic sequences ( ). Similarly, both RN_VAR1 and RN_VAR2 had 250 hits to our RN ESTs, respectively. Among these hits, 40(16%), 42(17%), 23(9%), and 145 (58%) had 0, 1, 2, and 3–21 mis-matches to the RN ESTs for RN_VAR1 respectively. For RN_VAR2, 111(44%), 88(35%), 35(14%), and 16(7%) had 0, 1, 2, and 3–5 mis-matches, respectively. The lesser number of mis-matches were observed for RN_VAR2. The secondary structures generated for RN_VAR1 and RN_VAR2 show 12 common sub-structures among these two variants (). These variants also show two and three sub-structures that were unrelated among RN_VAR1 and RN_VAR2 respectively. Sixty-three sub-structures were used in the generation of this folding pattern from 1,027 matched nucleotides. Multiple alignments of the RNA molecules showed structurally conserved areas among the RN variants and 110 single nucleotide polymorphisms (SNPs) (). The consensus structure obtained for RN_VAR1 and RN_VAR2 showed the most-likely structure both variants form together (). Multiple sequence alignment of 18S rDNA between the RN variants (VAR 1 and VAR 2) with 9 other nematodes provided evidence of the presence of hairpin 17 in the RN sequences on positions 391–455 ( ). Folding patterns of hairpin 17 in these nematodes reflects the conservative nature of this region ( ). Base pairs at positions 17a, 17b, and 17c were conserved among the RN_VAR1 and RN_VAR2 and in the other nine nematodes, and substitutions were noted at positions 17d and 17e.
Table 5

Blast hits among RN 18S rRNA variants, and REPLIg 454 genomic library.

Query and SubjectTotal number of unique hitsAverage % IDAverage alignment lengthAverage mis-matchesAverage number of gapsAverage bit score
RN_VAR1 vs REPLIg genomic sequences83.0098.97400.432.231.55749.64
RN_VAR2 vs REPLIg genomic sequences236.0099.25441.401.072.08832.79
Figure 4

Multiple sequence alignment of the primary structure of the 18S rRNA gene in RN showing hairpin 17 corresponding to positions 391 - 455 in RN_VAR1 and RN_VAR2.

Figure 5

The secondary structure of hairpin 17 of the 18S rDNA for RN_VAR1, RN_VAR2 of RN and 9 other nematode species.

The pairs of nucleotides peculiar to certain taxa are boxed in red at positions 17a, 17b, and 17c. Base pair for Daptonema, CA at position 17a is boxed in violet color.

The secondary structure of hairpin 17 of the 18S rDNA for RN_VAR1, RN_VAR2 of RN and 9 other nematode species.

The pairs of nucleotides peculiar to certain taxa are boxed in red at positions 17a, 17b, and 17c. Base pair for Daptonema, CA at position 17a is boxed in violet color.

Discussion

Ribosomal DNA (rDNA) codes for RNA molecules and these are separated by intergenic spacers that control transcription of the rRNA genes. Both the intergenic spacers and the transcribed regions of the rRNA gene show evidence of greater divergence between and considerably less within a species, a phenomenon known as the concerted evolution. The 18S rRNA gene regions of the rDNA evolve differentially and this attribute has been successfully used in phylogenetic analysis in kingdoms, phyla, classes, and orders [23], [25]. This high copy number gene of eukaryotic nuclear rDNA is generally tandemly organized on a chromosome. The gene region is divided into domains that are differentially conserved. Expansion segments evolve at high rates and are relatively more variable, and differences have been observed within the 18S rRNA gene, due to the presence of variable regions which expand within this region [26]. The comparative analysis of 44 newly generated near full-length 18S rDNA sequences from four female RNs distinctly show the presence of two sequence variants (RN_VAR1 and RN_VAR2). The presence of intra-individual variation within the 18S rDNA of animal genomes is rare. Variation occurring within the 18S rDNA of a metazoan species was first reported for Dugesia (Schmidtea) mediterranea which is a free-living platyhelminth. Two types of rDNA (type I and II) were observed in the genome of this organism, with 8% of nucleotides distinguishing the variants [2]. Variation has also been observed in the 18S rDNA region of Strongyloides where four hyper-variable regions (HVR-I to IV) were identified within the 18S rDNA among 34 isolates from 15 species [27]. Hyper-variable regions I to III differed from IV because of nucleotide lengths which were unequal, e.g., HVR-IV region ranged from 23 to 39 nucleotides. Similar to our study, these variable regions could serve as markers for identification of Strongyloides species. A study conducted using the rRNA gene loci of 12 Drosophila species showed the presence of variants within 7.8–8.2 kb region [9]. A range of 3 to 18 variants were present in more than 3% of the rDNA region in 11 species. Unlike our study, which focuses exclusively on the coding region, in Drosophila the non-coding regions had more variants than in the coding regions, and these variants were 3 to 8 times more prevalent in the expansion region. Our ongoing work also suggests that the presence of coding region variation seems to occur in individually-sequenced single nematodes from locations across Alabama and Mississippi. Variants have also been observed in the Plasmodium within its life cycle stages [28] and suggestion have been made that such sequence variation may impart some advantage to the host in its ability to escape susceptibility based on such genetic change. Additionally, in freshwater planarians variation has been reported in 2.5% to 3.9% within the 18S rDNA region of various species and genera [29]. A phenomenon of the 18S rRNA gene is the presence of slower evolving areas ‘core regions’ and ‘expansion regions’. The former is involved in substrate binding, while the expansion region is a variable region [30], [31]. Thus, the nucleotide sites distinguishing our RN 18S rRNA variants could possibly occur within both the core and expansion regions. The specific nature of these sites may have functional roles within the 18S rRNA genes of this nematode. An interesting question is why does this sequence variation exist? Does this variation provide an evolutionary advantage to this organism? Some of the possible explanations can only be proposed. Models used in clarifying this concern relating to concerted evolution occurring within the rDNA locus have focused on unequal crossovers and gene conversions. Unequal crossover among sister chromatids has been known to occur more frequently as compared to exchanges between homologous chromosomes [32], [33]and this may have contributed to the emergence rDNA variants observed in this study. Phylogenetic analysis provided additional clarity to the distinct grouping of clones that clustered together with their respective variant when Meloidogyne incognita isolate MeloInc1 was included as an outgroup branched away from the other RN clones and other taxa The second phylogram generated had inclusion of 26 other taxa belonging to the superfamily Tylenchoidea, in an attempt to place these RN variants within this superfamily. The NJ tree showed that RN_VAR1 clustered with Rotylenchulus reniformis isolate wb2 with a bootstrap value of 100 while the RN_VAR2 was less similar to wb2. Presence of long expansion segments in freshwater crustaceans such as Daphnia have been utilized to study rDNA evolution and they have sequences that maintain their function and they evolve in unison. These organisms have regions or expansion segments where sequence variations are most prevalent [34]. Compared to Daphnia, our study focused on the whole coding region of 18S gene and not only on the expansion helices. The base variations were observed mainly at 96 nucleotide positions and these variations were uniformly distributed throughout the entire 18S rDNA coding region (95 to 1673 bp). This study and our ongoing research suggest that the two variants observed here are very likely ancient polymorphisms that have not escaped concerted evolution of this crucial rDNA multi-gene family. The natural question that arises with the presence of rDNA variants is whether the secondary structure is conserved and functional in a compensated manner, and does a helical structure exist in both variants of RN, or does it cause the secondary structure in the variants to become uncompensated and destabilized? The hairpin 17 region though highly conserved is still prone to substitutions which are retained by the 18S rDNA. Secondary structure analysis among Strongylida and Rhabditida distinctly revealed substitutions at positions 17a, 17b and 17c [35] The hairpin 17 secondary structures from nematodes studied showed conservative pairs UG in position 17a, GC in position 17b, and GC in position 17c, and the loss of substitutions in position17e. These substitutions were categorized as non-compensated because there were base pair losses at these positions. In our study however, the substitutions observed in positions 17d and 17e of RN_VAR1 and RN_VAR2 show compensated changes in nucleotide pairs in the conservative region of the hairpin 17 region. Our secondary structure RNA analysis showed that the SNPs among the two RN variants are compensatory. A possible explanation of two main types of sequence variants within an individual 18S gene is the presence of two locations of rDNA or nucleolar organizing regions, with occasional crossing-over between these two regions. This may be the reason why small gene conversion tracts have existed for the RN 18S rDNA. The conserved nature of nucleotides in helices or loops demonstrates a link between rRNA folding and the ribosome function [35]. As long as both sequences fold into stable secondary structures, and do not interfere with ribosome function, one of the two variants may not be eliminated. Studies involving diverging rDNA copies together with their expression potential have been reported [2], [36]. In our study it is clear that the two variants reside in the same nematode, and both 18S variants are expressed because their sequences were present in our RN ESTs. This gene will also be useful for genetic diversity studies in the various genera of nematodes. In summary, we have identified two distinct variants within the 18S rRNA gene within three individual RNs, with one variant sequence being the most similar to the Rotylenchulus reniformis isolate wb2 (RN_VAR1). This variation may be as a result of recombination events during unequal crossover occurring between sister chromatids, or between X and Y rDNA leading to the co-evolution of the rDNA sequences on these chromosomes. Our findings can be placed in the context of earlier molecular mechanisms and paralogous or ancestral (plesiomorphic) explanations proposed by other investigators. The role of gene conversion and unequal crossover mechanisms in RN needs to be further explored to fully understand the existence of these 18S variants. Short oligonucleotide probes designed from RN_VAR1 would yield specific markers that will differentiate this species from closely related taxa belonging to super families such as Heteroderidae, Hoplolaimidae, Rotylenchulidae, and Hoplolaimidae. In diagnostic research and environmental monitoring, probes to be developed will provide opportunities for the rapid diagnostic detection and quantification of RN compositions in environmental samples, through real-time PCR detections.. Although this study did not investigate location-specific differences in the presence of the different sequence variants, it would be interesting in the future to determine the geographic distribution of RN variants 1 and 2. This could be achieved using the full length 18S rDNA analysis of RNs collected from Africa, and other centers of RN diversity (e.g., Japan, Hawaii, and Southern half of U.S.).

Materials and Methods

Materials Statement

No specific permits were required for the described field studies. This is because soil samples were obtained from a single location within the state at Alabama Agricultural Research Station, BelleMina. This field is not privately owned or protected in any way, and the field studies did not involve endangered or protected species as determined by our soil diagnostic. It is regularly monitored by a Nematologist who insures that this area is extensively infested by RN.

Nematode Extraction from the Soil and DNA Extractions

Each soil sample was thoroughly mixed and kept at 4°C; a 150 cm3 subsample was used for the extraction of the nematodes as described [37]. These populations were further maintained on cotton cultivar ‘Delta and Pineland 425 BG/RR’ (DPL 425) in the greenhouse to grow and multiply. Extracted RN suspension was placed into a petri dish using a pipette, and placed on an Olympus SZH-ILLD dissection microscope (Olympus optical Co. Ltd. Japan) using a low magnification of 10. Individual nematodes were picked using a sterilized hook and then placed in a drop of water on a glass slide for morphological identification, under an Olympus IMT-2 compound microscope (Olympus optical Co. Ltd. Japan). The identified individual nematodes were then hand-picked and immediately placed into DNA lysis buffer. DNA was extracted separately from each of the four individual female RNs using a DNeasy Blood and Tissue Kit (Qiagen, Inc., Valencia, CA) according to the manufacturer’s protocol.

Polymerase Chain Reaction (PCR)

Extracted DNA (2.0 µl) of approximately 1ng/µl from single female RNs was transferred into PCR tubes containing 2.5 µl 10x High Fidelity PCR buffer, 1.0 µl MgCl2 (50 mM), 0.5 µl dNTPs (10 mM), 0.5 µl of forward and reverse primers each (10 µm) (synthesized by MWG-Biotech AG), 0.2 µl of high fidelity platinum taq (Invitrogen, Carlsbad, CA) and dNASE-free water added to a final volume of 25 µl. Primer pairs SSUR81 and SSUF05 ( ) were used in amplifying an approximately 1800 bp region of the 18S rRNA gene of four individual female RNs separately. Both primers SSUF05 and SSUR81 had some degeneracy in their bases. PCR was performed in a Peltier Thermal Cycler (PTC) tetrad 2 DNA engine (Bio-Rad, Hercules, CA): 94°C for 2 min, then 39 cycles of: 94°C for 15 sec, 55°C for 30 sec, and 68°C for 1 min and a final extension phase of 68°C for 7 min. The quality of PCR products was checked by electrophoresis of 6 µl of PCR reaction in 1% agarose gel with ethidium bromide staining. The bands were visualized and photographed under uv light Bio-Rad gel imager (Bio-Rad, CA). The PCR product size was determined by comparing with a 1 kb DNA marker.

Cloning of PCR Products

Polymerase Chain Reaction (PCR) products from four individual female nematodes were purified using a QIAquick PCR Purification Kit (Qiagen, Inc., Valencia, CA) according to the manufactures protocol. The fragments were then cloned into a plasmid vector using TOPO TA Cloning Kit (Invitrogen, Carlsbad, CA). The ligation reaction was made up of 4 µl of PCR product, 1 µl of salt solution (1.2 M NaCl and 0.06 M MgCl2), and 1 µl of TOPO vector. Several clones were picked for verification of inserts from PCR amplifications for each nematode. This was performed by amplification of the clonal DNA using M13 forward (5′-TGTAAAACGACGGCCAGT-3′) and reverse (5′-AGCGGATAACAATTTCACAC-3′) primers. PCR conditions were as follows: 94°C for 5 min, then 40 cycles of the following: 94°C for 30 sec, 55°C for 1 min, and 72°C for 1 min. The final extension phase was 72°C for 10 min. Individual bacterial colonies with inserts were picked and placed in separate 1.5 ml centrifuge tubes with 3 ml of liquid Luria- Bertani (LB ) media containing 100 µg/ml ampicillin and shaken at 37 ○C for 24 hours at 100 rpm in an Innova 4300 rotary incubator shaker (New Brunswick Scientific, Edison, NJ). Tubes containing bacterial cells were centrifuged for 3 min at 13,000 rpm in a Hermle MR-2 (National Labnet Company, Woodbridge, NJ) tabletop centrifuge to obtain a cell pellet. Plasmid DNA was isolated using a QIAprep Mini-prep kit (Qiagen, Inc., Valencia, CA) according to the manufacturer’s protocol.

Sequencing

The 18S rRNA gene of the RN was sequenced using six primers ( ). Plasmid inserts from at least ten colonies originating from each of the four nematodes were sequenced using the ABI PRISM Big Dye Terminator cycle sequencing ready reaction kit in an ABI 3700 nucleotide sequencer (Applied Biosystems, Foster City, CA) and screened for homology to nematoda sequences through BLASTN on the NCBI website (http://www.ncbi.nlm.nih.gov/Blast.cgi).

Alignment and Phylogenetic Analysis

The SeqMan Pro within the DNASTAR Lasergene v8.0 software (DNASTAR Inc., Madison, WI) was used in generating consensus sequences and in trimming extraneous sequences outside the respective amplification fragments. Individual clone sequences from each nematode were used for multiple sequence alignment (http://align.genome.jp/ ) using default parameters and viewed using BioEdit software [24]. All assembled 18S rRNA sequences have been deposited in GenBank under accession numbers JX406335–JX406383. Phylogenetic analyses were conducted using Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0 [38].

Sequence Examination for Chimera Presence

Potential production of chimeric DNA in our SSU rDNA clones was examined for all 49 clones. Analysis was carried out using Ballerophon (http://comp-bio.anu.edu.au/bellerophon/bellerophon.pl) with default setting (200 bp window) [39].

Secondary Structure Analysis

The ExpaRNA tool was used in comparative structural RNA analysis to compute the best arrangement of sequence-structure motifs common among the RN variants. The LocARNA tool was also used for multiple sequence alignment of RNA molecules from the RN variants and for the generation of a consensus structure [40], [41]. The common eukaryotic structure based on the known secondary structures of 9 other nematodes species were obtained from GenBank: Paracanthonchus caecus (AF047888), Daptonema procerus (AF047889), Chromadoropsis vivipara (AF047891), Enoplus brevis (U88336), Trichinella spiralis (U60231), Plectus sp., (U61761), Ascaris sp. (M58348, etc.); Brugia malayi (H30951, etc.), and Meloidogyne arenaria (U42342) and aligned together with RN_VAR1 and RN_VAR2 sequences. The secondary structures of hairpin 17 were then analyzed using the mfold program with set parameters.

Comparison of RN 18S rRNA variants to RN Genomic and EST Sequences

The RN variants were also used in performing blasts searches (BLASTN) against unassembled reniform genomic sequences ( ). Whole Genome Amplification (WGA) was performed using a REPLI-g midi kits (Qiagen, MD), for generating a RN genomic library used in 454 sequencing. Large-scale 454 sequencing of RN genome gave a total of 732 MB for the REPLIg genomic library prepared from DNA of four individual adult female RNs pooled together. Detailed blast hits are presented in the supplemental tables between REPLIg genomic sequences andRN_VAR1 (), and RN_VAR2 (. These RN raw reads have been deposited under the sequence read achieve (SRA) section of GenBank with accession numbers SRX099033. Furthermore, comparative analysis of the RN variant sequences to our RN ESTs generated from pooled eggs and vermiform stage cDNA libraries was undertaken. Over 50,000 sequence reads were generated from the RN transcriptome, resulting in more than 4 Mb (4,781,676 bases) of sequence data. Detailed blast hits are presented in the supplemental tables between RN ESTs and RN_VAR1 (), and RN_VAR2 (. The ESTs have been submitted to GenBank under the accession numbers SRX098224 and SRX098225. Predicted secondary structure of 18S rDNA for RN_VAR1. Numbers indicate similar sub-structures among RN_VAR1 and RN_VAR2, and arrows indicate dissimilar sub-structures among RN_VAR1 and RN_VAR2. (TIF) Click here for additional data file. Predicted secondary structure of 18S rDNA for RN_VAR2. Numbers indicate similar sub-structures among RN_VAR2 and RN_VAR1, and arrows indicate dissimilar sub-structures among RN_VAR2 and RN_VAR1. (TIF) Click here for additional data file. Multiple alignments of the 18S rDNA for RN_VAR1 and RN_VAR2. The bases highlighted in red indicate highly conserved areas between RN_VAR1 and RN_VAR2. The grooves in the grey bar represent dissimilarity among the nucleotides. (TIF) Click here for additional data file. The consensus structure formed both by RN_VAR1 and RN_VAR2 together. (TIF) Click here for additional data file. Multiple sequence alignment of 44 18S rDNA sequences from four female reniform nematodes. (XLSX) Click here for additional data file. Blast hits among RN_VAR1 and REPLIg genomic sequences. (XLSX) Click here for additional data file. Blast hits among RN_VAR2 and REPLIg genomic sequences. (XLSX) Click here for additional data file. Blast hits among RN_VAR1 and RN ESTs (XLSX) Click here for additional data file. Blast hits among RN_VAR2 and RN ESTs. (XLSX) Click here for additional data file.
  25 in total

Review 1.  Molecular analysis of nematode diversity and the evolution of parasitism.

Authors:  M Dorris; P De Ley; M L Blaxter
Journal:  Parasitol Today       Date:  1999-05

2.  18S rRNA gene variation among common airborne fungi, and development of specific oligonucleotide probes for the detection of fungal isolates.

Authors:  Zhihong Wu; Yoshihiko Tsumura; Göran Blomquist; Xiao-Ru Wang
Journal:  Appl Environ Microbiol       Date:  2003-09       Impact factor: 4.792

3.  LocARNA-P: accurate boundary prediction and improved detection of structural RNAs.

Authors:  Sebastian Will; Tejal Joshi; Ivo L Hofacker; Peter F Stadler; Rolf Backofen
Journal:  RNA       Date:  2012-03-26       Impact factor: 4.942

Review 4.  Ribosomal DNA: molecular evolution and phylogenetic inference.

Authors:  D M Hillis; M T Dixon
Journal:  Q Rev Biol       Date:  1991-12       Impact factor: 4.875

5.  Evolution of Hypervariable Regions, V4 and V7, of Insect 18S rRNA and Their Phylogenetic Implications.

Authors:  U W Hwang; H I Ree; W Kim
Journal:  Zoolog Sci       Date:  2000-01-01       Impact factor: 0.931

Review 6.  Molecular phylogeny of the animal kingdom.

Authors:  K G Field; G J Olsen; D J Lane; S J Giovannoni; M T Ghiselin; E C Raff; N R Pace; R A Raff
Journal:  Science       Date:  1988-02-12       Impact factor: 47.728

7.  Xenopus laevis 28S ribosomal RNA: a secondary structure model and its evolutionary and functional implications.

Authors:  C G Clark; B W Tague; V C Ware; S A Gerbi
Journal:  Nucleic Acids Res       Date:  1984-08-10       Impact factor: 16.971

8.  Hyper-variable regions in 18S rDNA of Strongyloides spp. as markers for species-specific diagnosis.

Authors:  Hideo Hasegawa; Shotaro Hayashida; Yatsukaho Ikeda; Hiroshi Sato
Journal:  Parasitol Res       Date:  2008-12-03       Impact factor: 2.289

9.  Rotylenchulus reniformis Management in Cotton with Crop Rotation.

Authors:  R F Davis; S R Koenning; R C Kemerait; T D Cummings; W D Shurley
Journal:  J Nematol       Date:  2003-03       Impact factor: 1.402

10.  Phylogenetic Relationships Among Xiphinema and Xiphidorus Nematode Species from Brazil Inferred from 18S rDNA Sequences.

Authors:  Claudio M G Oliveira; Judith Hübschen; Derek J F Brown; Luiz C C B Ferraz; Frank Wright; Roy Neilson
Journal:  J Nematol       Date:  2004-06       Impact factor: 1.402

View more
  5 in total

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

Authors:  Emily A Brown; Frédéric J J Chain; Teresa J Crease; Hugh J MacIsaac; Melania E Cristescu
Journal:  Ecol Evol       Date:  2015-05-13       Impact factor: 2.912

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

Authors:  E I Koshel; V V Aleshin; G A Eroshenko; V V Kutyrev
Journal:  Biomed Res Int       Date:  2014-04-03       Impact factor: 3.411

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

Authors:  Ambra Viviani; Rodolfo Bernardi; Andrea Cavallini; Elisabetta Rossi
Journal:  J Insect Sci       Date:  2019-07-01       Impact factor: 1.857

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

Authors:  Seloame T Nyaku; Ramesh V Kantety; Ernst Cebert; Kathy S Lawrence; Joseph O Honger; Govind C Sharma
Journal:  Plant Pathol J       Date:  2016-04-01       Impact factor: 1.795

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

Authors:  Siyu Zhou; Xiaoyin Fu; Pei Pei; Marek Kucka; Jing Liu; Lili Tang; Tingzheng Zhan; Shanshan He; Yingguang Frank Chan; Christian Rödelsperger; Dengyu Liu; Adrian Streit
Journal:  PLoS Negl Trop Dis       Date:  2019-05-06
  5 in total

北京卡尤迪生物科技股份有限公司 © 2022-2023.