Length variation in 18S rRNA expansion segment 43/e4 of Daphnia obtusa: ancient or recurring polymorphism?

Copy number and sequence variation in rDNA of Daphnia pulex from natural populations: insights from whole-genome sequencing

Ribosomal DNA (rDNA) has a vital role in ribosome biogenesis as it contains the genes that encode ribosomal RNA (rRNA) separated by intergenic spacers (IGSs). The rRNA genes occur in hundreds to tens of thousands of copies per haploid genome in eukaryotes and are generally highly conserved with low variation within species. Due to the repetitive nature and large size of rDNA arrays, detecting intraindividual variation can be difficult. In this study, we use whole-genome sequences of 169 Daphnia pulex individuals from 10 natural populations to measure the copy number and sequence variation in rDNA. This revealed that variation in rDNA copy number between individuals spans an order of magnitude. We further observed a substantial level of sequence variation within individual genomes. As expected, single-nucleotide polymorphisms occurred in regions of lower functional constraint such as the IGS and expansion segments of the rRNA genes. The presence of strong linkage disequilibrium among variants facilitated identification of haplotypes within each population. Although there was evidence of recombination among haplotypes from different populations, it is insufficient to eliminate linkage disequilibrium within populations. Estimating copy number and haplotype diversity within individuals revealed that the level of intraindividual sequence variation is not strongly correlated with copy number. The observed patterns of variation highlight a complex evolutionary history of rDNA in D. pulex. Future research should explore the functional implications of rDNA copy number and sequence variation on organismal phenotypes.

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.

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 two intra-individual sequence variants in the 18S rRNA gene in the plant parasitic nematode, Rotylenchulus reniformis

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.