General properties and phylogenetic utilities of nuclear ribosomal DNA and mitochondrial DNA commonly used in molecular systematics

A footworm in the door: revising Onchocerca phylogeny with previously unknown cryptic species in wild North American ungulates

Onchocerca is an important genus of vector-borne filarial nematodes that infect both humans and animals worldwide. Many Onchocerca spp., most of medical and veterinary health relevance, are the focus of a variety of diagnostic and molecular research. However, despite the importance of these parasites, there is growing evidence of previously unexplored genetic diversity of these nematodes, particularly among wild ungulate hosts in North America. These understudied parasites prevent us from comprehending the evolutionary history of the genus Onchocerca, monitoring potential One Health threats, and improving our filarioid diagnostic capabilities. In order to fill these knowledge gaps, we identified five uncharacterized Onchocerca lineages and compared them with other well-known filarioid species using single and concatenated gene regions (i.e., nd5, cox1, 12S, 18S, 28S, hsp70, MyoHC, rbp1). Phylogenetic analyses revealed that the novel Onchocerca lineages of wild North American ungulates segregate into two clades. One clade comprised Onchocerca lineages II, IV, and V and other species found mainly in domestic animals and humans, and the second comprised Onchocerca lineages I and III and other species from a variety of hosts including cervids, bovids, and equids. The formation of two clearly separate clades supports the idea of at least two independent expansion events of ancestral Onchocerca spp. into the North American continent via the Bering land bridge. Cophylogenetic analysis shows evidence of ancestral Onchocerca spp. of Bovidae host-switching to wild Cervidae and giving rise to the novel Onchocerca spp. Lastly, pairwise analysis confirms informative molecular markers of diagnostic relevance in both mitochondrial and nuclear gene regions of filarioid nematodes. The overall information provides greater context to the genus Onchocerca and emphasizes the need to discover, characterize, and monitor neglected parasites, especially those of wildlife origin.

The first complete mitochondrial genome of Loimia arborea (Polychaeta: Terebellidae) and phylogenetic analysis

In order to understand the molecular insights within the Terebellida, the complete mitochondrial genome of Loimia arborea was sequenced. The mitochondrial genome is 16,023 bp, with 13 protein-coding genes (PCGs), 23 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, and a non-coding region (D-loop). Notably, two adjacent copies of methionine tRNA genes (trnMs) were detected. The phylogeny of Terebellida was constructed based on 13 PCGs from 13 species, two main clades were strongly supported, i.e., Cirratuliformia (clade A) and Terebelliformia (clade B). And the monophyly of the Terebellidae was restored.

Genetic Variation among the Partial Gene Sequences of the Ribosomal Protein Large-Two, the Internal Transcribed Spacer, and the Small Ribosomal Subunit of Blastocystis sp. from Human Fecal Samples

In the present study, we compared the genetic variability of fragments from the internal transcribed spacer region (ITS) and the small subunit ribosomal DNA (SSUrDNA) as nuclear markers, in contrast with the ribosomal protein large two (rpl2) loci, placed in the mitochondrion-related organelles (MROs) within and among human fecal samples with Blastocystis. Samples were analyzed using polymerase chain reaction (PCR)-sequencing, phylogenies, and genetics of population structure analyses were performed. In total, 96 sequences were analyzed, i.e., 33 of SSUrDNA, 35 of rpl2, and 28 of ITS. Only three subtypes (STs) were identified, i.e., ST1 (11.4%), ST2 (28.6%), and ST3 (60%); in all cases, kappa indexes were 1, meaning a perfect agreement among ST assignations. The topologies of phylogenetic inferences were similar among them, clustering to each ST in its specific cluster; discrepancies between phylogeny and assignment of STs were not observed. The STRUCTURE v2.3.4 software assigned three subpopulations corresponding to the STs 1-3, respectively. The population indices were consistent with those previously reported by other groups. Our results suggest the potential use of the ITS and rpl2 genes as molecular markers for Blastocystis subtyping as an alternative approach for the study of the genetic diversity observed within and between human isolates of this microorganism.

Mitochondrial and ribosomal markers in the identification of nematodes of clinical and veterinary importance

BACKGROUND

Nematodes of the Ascarididae, Ancylostomatidae and Onchocercidae families are parasites of human and veterinary importance causing infections with high prevalence worldwide. Molecular tools have significantly improved the diagnosis of these helminthiases, but the selection of genetic markers for PCR or metabarcoding purposes is often challenging because of the resolution these may show.

METHODS

Nuclear 18S rRNA, internal transcribed spacers 1 (ITS-1) and 2 (ITS-2), mitochondrial gene cytochrome oxidase 1 (cox1) and mitochondrial rRNA genes 12S and 16S loci were studied for 30 species of the mentioned families. Accordingly, their phylogenetic interspecies resolution, pairwise nucleotide p-distances and sequence availability in GenBank were analyzed.

RESULTS

The 18S rRNA showed the least interspecies resolution since separate species of the Ascaris, Mansonella, Toxocara or Ancylostoma genus were intermixed in phylogenetic trees as opposed to the ITS-1, ITS-2, cox1, 12S and 16S loci. Moreover, pairwise nucleotide p-distances were significantly different in the 18S compared to the other loci, with an average of 99.1 ± 0.1%, 99.8 ± 0.1% and 98.8 ± 0.9% for the Ascarididae, Ancylostomatidae and Onchocercidae families, respectively. However, ITS-1 and ITS-2 average pairwise nucleotide p-distances in the three families ranged from 72.7% to 87.3%, and the cox1, 12S and 16S ranged from 86.4% to 90.4%. Additionally, 2491 cox1 sequences were retrieved from the 30 analyzed species in GenBank, whereas 212, 1082, 994, 428 and 143 sequences could be obtained from the 18S, ITS-1, ITS-2, 12S and 16S markers, respectively.

CONCLUSIONS

The use of the cox1 gene is recommended because of the high interspecies resolution and the large number of sequences available in databases. Importantly, confirmation of the identity of an unknown specimen should always be complemented with the careful morphological examination of worms and the analysis of other markers used for specific parasitic groups.

Genetic variation of the freshwater snail Indoplanorbis exustus (Gastropoda: Planorbidae) in Thailand, inferred from 18S and 28S rDNA sequences

Indoplanorbis exustus, a freshwater pulmonate snail, is widely distributed in tropical and subtropical zones and plays a significant role as an intermediate host for trematode parasites. Various genetic markers have been used for species identification and phylogenetic studies of this snail. However, there are limited studies about their molecular genetics based on nuclear ribosomal DNA (rDNA) genes. A genetic analysis of I. exustus in Thailand was conducted based on the nuclear 18S rDNA (339 bp) and 28S rDNA (1036 bp) genes. Indoplanorbis snails were collected from 29 localities in 21 provinces covering six regions of Thailand. Nucleotide sequences from 44 snails together with sequences from the GenBank database were examined for phylogenetic relationships and genetic diversity. All sequences of the selected nucleotide regions exhibited a high level of similarity (99%) to the sequences of I. exustus in the GenBank database. The maximum likelihood tree based on the 18S and 28S rDNA fragment sequences of I. exustus in Thailand revealed only one group with clear separation from another genus in the family Planorbidae. The I. exustus 28S rDNA sequences showed intraspecific genetic divergence ranging from 0 to 0.78% and were classified into 8 different haplotypes. Conversely, the 18S rDNA data showed lower variation than the 28S rDNA data and revealed a single haplotype and intraspecific distances of zero among all sampled individuals. The haplotype network of 28S rDNA sequences of I. exustus in Thailand revealed six unique haplotypes and two haplotypes shared by at least two regions. Overall, both markers were successful in the identification of I. exustus. However, these markers, particularly the 18S rDNA, may not be suitable for genetic analysis within the species, particularly for population genetic studies, due to their limited variation as seen in this study. In summary, this study not only enhances understanding of genetic variation in I. exustus but is also useful for the selection of molecular markers in future genetic research.

Phylogeography of Two Enigmatic Sulphur Butterflies, Colias mongola Alphéraky, 1897 and Colias tamerlana Staudinger, 1897 (Lepidoptera, Pieridae), with Relations to Wolbachia Infection

The genus Colias Fabricius, 1807 includes numerous taxa and forms with uncertain status and taxonomic position. Among such taxa are Colias mongola Alphéraky, 1897 and Colias tamerlana Staudinger, 1897, interpreted in the literature either as conspecific forms, as subspecies of different but morphologically somewhat similar Colias species or as distinct species-level taxa. Based on mitochondrial and nuclear DNA markers, we reconstructed a phylogeographic pattern of the taxa in question. We recover and include in our analysis DNA barcodes of the century-old type specimens, the lectotype of C. tamerlana deposited in the Natural History Museum (Museum für Naturkunde), Berlin, Germany (ZMHU) and the paralectotype of C. tamerlana and the lectotype of C. mongola deposited in the Zoological Institute, Russian Academy of Sciences, St. Petersburg, Russia (ZISP). Our analysis grouped all specimens within four (HP_I-HP_IV) deeply divergent but geographically poorly structured clades which did not support nonconspecifity of C. mongola-C. tamerlana. We also show that all studied females of the widely distributed haplogroup HP_II were infected with a single Wolbachia strain belonging to the supergroup B, while the males of this haplogroup, as well as all other investigated specimens of both sexes, were not infected. Our data highlight the relevance of large-scale sampling dataset analysis and the need for testing for Wolbachia infection to avoid erroneous phylogenetic reconstructions and species misidentification.

Complete mitochondrial genome of a golden orb-web spider Trichonephila clavata (Chelicerata, Arachnida) from South Korea

The mitochondrial genome of a golden orb-web spider Trichonephila clavata (L. Koch, 1878) from South Korea is determined and characterized in detail, which is the second mitochondrial genome reported from this species: the first was published from the Chinese sample by Pan et al. (2016). It was 14,436 bp in length being composed of 13 protein-coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes, and one control region (CR). It has a base composition of 35.99% for 'A,' 14.88% for 'G,' 9.09% for 'C,' and 40.04% for 'T.' Comparing the South Korean and Chinese mitochondrial genomes, we observed 8% nucleotide sequence differences between their CRs, caused by the different numbers and sorts of possessed tandem repeats, suggesting a promising molecular marker to distinguish South Korean individuals from Chinese ones. The phylogenetic trees using the maximum likelihood (ML) method were reconstructed with nucleotides (without 3rd codon position) and amino acids from 13 PCGs, respectively, which consistently confirmed that T. clavata (Subfamily Nephilinae) from South Korea and China are clustered together, distinctly separated from the other subfamily Araneinae in the monophyletic family Araneidae.

Molecular Identification and Phylogenetic Analysis of Laelapidae Mites (Acari: Mesostigmata)

The family Laelapidae (Dermanyssoidea) is morphologically and ecologically the most diverse group of Mesostigmata mites. Although molecular genetic data are widely used in taxonomic identification and phylogenetic analysis, most classifications in Mesostigmata mites are based solely on morphological characteristics. In the present study, eight species of mites from the Laelapidae (Dermanyssoidea) family collected from different species of small rodents in Lithuania, Norway, Slovakia, and the Czech Republic were molecularly characterized using the nuclear (28S ribosomal RNA) and mitochondrial (cytochrome oxidase subunit I gene) markers. Obtained molecular data from 113 specimens of mites were used to discriminate between species and investigate the phylogenetic relationships and genetic diversity among Laelapidae mites from six genera. This study provides new molecular data on Laelaps agilis, Laelaps hilaris, Laelaps jettmari, Haemogamasus nidi, Eulaelaps stabularis, Hyperlaelaps microti, Myonyssus gigas, and Hirstionyssus sp. mites collected from different rodent hosts and geographical regions in Europe.

Comparative analysis of mitochondrial genomes reveals family-specific architectures and molecular features in scorpions (Arthropoda: Arachnida: Scorpiones)

Scorpions are a group of arachnids with great evolutionary success that comprise more than 2,000 described species. Mitochondrial genomes have been little studied in this clade. We describe and compare different scorpion mitochondrial genomes and analyze their architecture and molecular characteristics. We assembled eight new scorpion mitochondrial genomes from transcriptomic datasets, annotated them, predicted the secondary structures of tRNAs, and compared the nucleotide composition, codon usage, and relative synonymous codon usage of 16 complete scorpion mitochondrial genomes. Lastly, we provided a phylogeny based on all mitochondrial protein coding genes. We characterized the mitogenomes in detail and reported particularities such as dissimilar synteny in the family Buthidae compared to other scorpions, unusual tRNA secondary structures, and unconventional start and stop codons in all scorpions. Our comparative analysis revealed that scorpion mitochondrial genomes exhibit different architectures and features depending on taxonomic identity. We highlight the parvorder Buthida, particularly the family Buthidae, as it invariably exhibited different mitogenome features such as synteny, codon usage, and AT-skew compared to the parvorder Iurida that included the rest of the scorpion families we analyzed in this study. Our results provide a better understanding of the evolution of mitogenome features and phylogenetic relationships in scorpions.

Possible species discrimination of a blotched nerite Nerita albicilla with their distribution pattern and demographic history in the Indo-Pacific

The blotched nerite Nerita albicilla (Linnaeus 1758) is distributed in intertidal areas of the Indo-Pacific. In South Korea, it has been found only in the southernmost region of Jeju Island so far. Owing to its limited distribution, it can be a promising intertidal species helpful for monitoring global climate change effects in the Korean Peninsula. We performed population genetic analyses based on 393 COI haplotypes from 697 N. albicilla, including 167 from this study and 530 from public databases. The results showed that there are two distinct genetic lineages in N. albicilla: PAIO (Palearctic, Australasia, Indo-Malay, and Oceania) and Afrotropic lineages. DNA barcoding gap analyses indicated that the two lineages could be differentiated into two different species: N. albicilla (PAIO) and N. originalis sp. nov. (Afrotropic) (3.96%). Additionally, it was revealed that their divergence time was ca. 5.96 Ma and dramatic diversification of COI haplotypes occurred during the late Pliocene and Pleistocene. The results of MDA, BSP, and neutrality test implied recent population size expansion, which was estimated to be ca. 250 Ka. Finally, we discussed whether the observation of N. originalis sp. nov. in South Korea is due to the northward migration through ocean currents caused by global warming or due to artificial activity through marine transportation.

Molecular taxonomy of phlebotomine sand flies (Diptera, Psychodidae) with emphasis on DNA barcoding: A review

The taxonomy and systematics of sand flies (Diptera, Psychodidae, Phlebotominae) are one of the pillars of research aimed to identifying vector populations and the agents transmitted by these insects. Traditionally, the use of morphological traits has been the main line of evidence for the definition of species, but the use of DNA sequences is useful as an integrative approach for their delimitation. Here, we discuss the current status of the molecular taxonomy of sand flies, including their most sequenced molecular markers and the main results. Only about 37% of all sand fly species have been processed for any molecular marker and are publicly available in the NCBI GenBank or BOLD Systems databases. The genera Phlebotomus, Nyssomyia, Psathyromyia and Psychodopygus are well-sampled, accounting for more than 56% of their sequenced species. However, less than 34% of the species of Sergentomyia, Lutzomyia, Trichopygomyia and Trichophoromyia have been sampled, representing a major gap in the knowledge of these groups. The most sequenced molecular markers are those within mtDNA, especially the DNA barcoding fragment of the cytochrome c oxidase subunit I (coi) gene, which has shown promising results in detecting cryptic diversity within species. Few sequences of conserved genes have been generated, which hampers higher-level phylogenetic inferences. We argue that sand fly species should be sequenced for at least the coi DNA barcoding marker, but multiple markers with different mutation rates should be assessed, whenever possible, to generate multilocus analysis.

Mitochondrial DNA of Sardinian and North-West Italian Populations Revealed a New Piece in the Mosaic of Phylogeography and Phylogeny of Salariopsis fluviatilis (Blenniidae)

The genus Salariopsis (Blenniidae) comprises freshwater blenny fish that inhabits Mediterranean Sea, Black Sea, and north-east Atlantic areas. Three species were formally described to date: Salariopsis fluviatilis. S. economidisi, and S. atlantica. In this study, 103 individuals were collected from different Italian regions (Sardinia, Liguria, Piedmont, Lombardy) and analyzed using the mtDNA Control Region and the ribosomal 16s gene. We aimed (i) to depict the phylogeographic patterns of S. fluviatilis in northern Italy and Sardinia and (ii) to compare the genetic structure of Italian samples with those from other Mediterranean regions. Results obtained showed the presence of a well-supported genetic structuring among Italian S. fluviatilis populations, shedding new light on the phylogeographic patterns of northern Italian populations of S. fluviatilis sensu stricto across the Ligurian Alpine ridge and the Sardinia Island-mainland dispersal patterns. Furthermore, our species delimitation analysis was consistent in supporting results of previous research about the presence of genetic differentiation among S. fluviatilis, evidencing: (i) a large group of S. fluviatilis sensu stricto that includes two sub-groups (Occidental and Oriental), (ii) one group comprising populations from the Middle East of a taxonomic entity corresponding to Salariopsis cf. fluviatilis, and (iii) one group of Iberian individuals from the Guadiana River.

The complete mitochondrial genome of the woodwasp Euxiphydria potanini (Hymenoptera, Xiphydrioidea) and phylogenetic implications for symphytans

The long-necked woodwasp superfamily Xiphydrioidea belongs to the suborder Symphyta (Hymenoptera). Here we newly characterize the complete mitochondrial genome of the South Korean Euxiphydria potanini (Xiphydriidae) using next-generation sequencing: 16,500 bp long with 84.27% A + T content and 37 typical mitochondrial genes including those encoding 13 PCGs, 2 rRNAs, 22 tRNAs, and one A + T rich region. We compare the patterns of symphytan mitochondrial gene arrangement with those of an ancestral insect form and found some synapomorphic rearrangements in phylogenetic context. We use a variety of nucleotide and amino acid sequence alignments (thirteen mtPCGs and/or eight nDNAs) alongside step-by-step exclusions of long-branched taxa to elucidate the phylogenetic position of Xiphydrioidea and phylogenetic relationships among the seven symphytan superfamilies, except for Anaxyeloidea of which no mtgenome was available. The monophyly of symphytan superfamilies (with weak support for Pamphilioidea), sister-group relationship of Xiphydrioidea and Cephoidea, and Symphyta being paraphyletic to Apocrita, etc. are consistently supported by maximum likelihood and Bayesian inference trees. We also discuss the problematic phylogenetic positions of Orussoidea and Siricoidea and propose a hypothetical scenario of morphological character transition during hymenopteran evolution based on morphological key characteristics, such as the cenchrus and the wasp-waist.

Complete mitochondrial genome of a malaria vector mosquito Anopheles sinensis from South Korea

We present the complete mitochondrial genome of a Malaria vector Mosquito Anopheles sinensis Wiedemann, 1828 from South Korea. The mitochondrial genome is about 15,421 bp long and contains 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes, and an A-T rich region. The present data were compared with those from China with respect to PCG sequence differences, tRNA structure, gene order, and control region (CR) structure. An. sinensis mitochondrial genomes from northeast Asia share identical gene composition and gene order. In contrast, they have distinct differences in the CR within the range 8.75% (51/583 bp) to 9.95% (58/583 bp). The phylogenetic analysis showed that An. sinensis from South Korea was clustered together with those from China, but there existed distinct genetic distance between the two. Likewise, mitochondrial genome sequences from other Anopheles species were employed to infer phylogenetic relationships among the members of the genus Anopheles. This study further promotes the enrichment of An. sinensis mitochondrial genome data, providing useful information for their mitochondrial genetic differences along with geographical distances in northeast Asia.

Partial Ribosomal Nontranscribed Spacer Sequences Distinguish Rhagoletis zephyria (Diptera: Tephritidae) From the Apple Maggot, R. pomonella

The apple maggot, Rhagoletis pomonella (Walsh), was introduced into the apple-growing regions of the Pacific Northwest in the U.S.A. during the past 60-100 yr. Apple maggot (larvae, puparia, and adults) is difficult to distinguish from its morphologically similar sister species, Rhagoletis zephyria Snow, which is native and abundant in the Pacific Northwest. While morphological identifications are common practice, a simple, inexpensive assay based on genetic differences would be very useful when morphological traits are unclear. Here we report nucleotide substitution and insertion-deletion mutations in the nontranscribed spacer (NTS) of the ribosomal RNA gene cistron of R. pomonella and R. zephyria that appear to be diagnostic for these two fly species. Insertion-deletion variation is substantial and results in a 49 base-pair difference in PCR amplicon size between R. zephyria and R. pomonella that can be scored using agarose gel electrophoresis. PCR amplification and DNA sequencing of 766 bp of the NTS region from 38 R. pomonella individuals and 35 R. zephyria individuals from across their geographic ranges led to the expected PCR fragments of approx. 840 bp and 790 bp, respectively, as did amplification and sequencing of a smaller set of 26 R. pomonella and 16 R. zephyria flies from a sympatric site in Washington State. Conversely, 633 bp mitochondrial COI barcode sequences from this set of flies were polyphyletic with respect to R. pomonella and R. zephyria. Thus, differences in NTS PCR products on agarose gels potentially provide a simple way to distinguish between R. pomonella and R. zephyria.

Dactylogyridae 2022: a meta-analysis of phylogenetic studies and generic diagnoses of parasitic flatworms using published genetic and morphological data

Dactylogyridae is one of the most studied families of parasitic flatworms with more than 1000 species and 166 genera described to date including ecto- and endoparasites. Dactylogyrid monogeneans were suggested as model organisms for host-parasite macroevolutionary and biogeographical studies due to the scientific and economic importance of some of their host lineages. Consequently, an array of phylogenetic research into different dactylogyrid lineages has been produced over the past years but the last family-wide study was published 16 years ago. Here, we provide a meta-analysis of the phylogenetic relationships of Dactylogyridae including representatives of all genera with available molecular data (n=67). First, we investigate the systematic informativeness of morphological characters widely used to diagnose dactylogyrid genera through a parsimony analysis of the characters, character mapping, and phylogenetic comparative methods. Second, we provide an overview of the current state of the systematics of the family and its subfamilies, and summarise potentially poly- and paraphyletic genera. Third, we elaborate on the implications of taxonomic, citation, and confirmation bias in past studies. Fourth, we discuss host range, biogeographical, and freshwater-marine patterns. We found two well-supported macroclades which we assigned to the subfamilies Dactylogyrinae and Ancyrocephalinae. These subfamilies further include 16 well-supported clades with only a few synapomorphies that could be deduced from generic diagnoses in the literature. Furthermore, few morphological characters considered systematically informative at the genus level display a strong phylogenetic signal. However, the parsimony analysis suggests that these characters provide little information on the relationships between genera. We conclude that a strong taxonomic bias and low coverage of DNA sequences and regions limit knowledge on morphological and biogeographical evolutionary patterns that can be inferred from these results. We propose addressing potential citation and confirmation biases through a 'level playing field' multiple sequence alignment as provided by this study.

Liolophura species discrimination with geographical distribution patterns and their divergence and expansion history on the northwestern Pacific coast

The chiton Liolophura japonica (Lischke 1873) is distributed in intertidal areas of the northwestern Pacific. Using COI and 16S rRNA, we found three genetic lineages, suggesting separation into three different species. Population genetic analyses, the two distinct COI barcoding gaps albeit one barcoding gap in the 16S rRNA, and phylogenetic relationships with a congeneric species supported this finding. We described L. koreana, sp. nov. over ca. 33°24′ N (JJ), and L. sinensis, sp. nov. around ca. 27°02′–28°00′ N (ZJ). We confirmed that these can be morphologically distinguished by lateral and dorsal black spots on the tegmentum and the shape of spicules on the perinotum. We also discuss species divergence during the Plio-Pleistocene, demographic expansions following the last interglacial age in the Pleistocene, and augmentation of COI haplotype diversity during the Pleistocene. Our study sheds light on the potential for COI in examining marine invertebrate species discrimination and distribution in the northwestern Pacific.

Assessing the suitability of mitochondrial and nuclear DNA genetic markers for molecular systematics and species identification of helminths

BACKGROUND

Genetic markers are employed widely in molecular studies, and their utility depends on the degree of sequence variation, which dictates the type of application for which they are suited. Consequently, the suitability of a genetic marker for any specific application is complicated by its properties and usage across studies. To provide a yardstick for future users, in this study we assess the suitability of genetic markers for molecular systematics and species identification in helminths and provide an estimate of the cut-off genetic distances per taxonomic level.

METHODS

We assessed four classes of genetic markers, namely nuclear ribosomal internal transcribed spacers, nuclear rRNA, mitochondrial rRNA and mitochondrial protein-coding genes, based on certain properties that are important for species identification and molecular systematics. For molecular identification, these properties are inter-species sequence variation; length of reference sequences; easy alignment of sequences; and easy to design universal primers. For molecular systematics, the properties are: average genetic distance from order/suborder to species level; the number of monophyletic clades at the order/suborder level; length of reference sequences; easy alignment of sequences; easy to design universal primers; and absence of nucleotide substitution saturation. Estimation of the cut-off genetic distances was performed using the 'K-means' clustering algorithm.

RESULTS

The nuclear rRNA genes exhibited the lowest sequence variation, whereas the mitochondrial genes exhibited relatively higher variation across the three groups of helminths. Also, the nuclear and mitochondrial rRNA genes were the best possible genetic markers for helminth molecular systematics, whereas the mitochondrial protein-coding and rRNA genes were suitable for molecular identification. We also revealed that a general gauge of genetic distances might not be adequate, using evidence from the wide range of genetic distances among nematodes.

CONCLUSION

This study assessed the suitability of DNA genetic markers for application in molecular systematics and molecular identification of helminths. We provide a novel way of analyzing genetic distances to generate suitable cut-off values for each taxonomic level using the 'K-means' clustering algorithm. The estimated cut-off genetic distance values, together with the summary of the utility and limitations of each class of genetic markers, are useful information that can benefit researchers conducting molecular studies on helminths.

Molecular and morphological description of the first Hepatozoon (Apicomplexa: Hepatozoidae) species infecting a neotropical turtle, with an approach to its phylogenetic relationships

Haemogregarines (Adeleorina) have a high prevalence in turtles. Nevertheless, there is only one Hepatozoon species described that infects Testudines so far; it is Hepatozoon fitzsimonsi which infects the African tortoise Kinixys belliana. Colombia harbours a great diversity of chelonians; however, most of them are threatened. It is important to identify and characterize chelonian haemoparasite infections to improve the clinical assessments, treatments and the conservation and reintroduction programs of these animals. To evaluate such infections for the Colombian wood turtle Rhinoclemmys melanosterna, we analysed blood from 70 individuals. By using the morphological characteristics of blood stages as well as molecular information (18S rRNA sequences), here we report a new Hepatozoon species that represents the first report of a hepatozoid species infecting a semi-aquatic continental turtle in the world. Although the isolated lineage clusters within the phylogenetic clades that have morphological species of parasites already determined, their low nodal support makes their position within each group inconclusive. It is important to identify new molecular markers to improve parasite species identification. In-depth research on blood parasites infecting turtles is essential for increasing knowledge that could assess this potential unknown threat, to inform the conservation of turtles and for increasing the state of knowledge on parasites.

The mitochondrial genome of a giant water bug Lethocerus deyrollei (Hemiptera: Belostomatidae) from South Korea

A giant water bug Lethocerus deyrollei (Hemiptera: Belostomatidae) is a large, predatory, and nocturnal hemipteran insect, which has been considered threatened and thus enrolled as an endangered species in South Korea and Japan. Here, we characterized the complete mitochondrial genome of L. deyrollei, which has a circular form with 19,295 bp in length, which is the longest when compared to those of the 111 hemipteran species reported so far. Its longest genome size is due to the extremely extended CR (4686 bp), which is much longer than those of China and Japan. It consisted of a total of 37 genes (13 PCGs, 22 tRNA genes, and two rRNA genes) and one control region (CR). The genome composition and gene order were identical to those previously reported from the same species of China and Japan with over 99.7% sequence similarities except for CR and trnI. The nucleotide composition was highly A + T biased, accounting for 71% of the whole mitochondrial genome, as in other species of Nepoidea. Based on the aa sequences of 13 PCGs, we reconstructed a maximum likelihood tree, which indicated that the three mitochondrial genomes of L. deyrollei from South Korea, China, and Japan are grouped, and also Lethocerus, Belostomatidae, Nepoidea, Nepomorpha, Heteroptra are strong monophyletic groups, respectively.

The complete mitochondrial genome of an endangered triton snail Charonia lampas (Littorinimorpha: Charoniidae) from South Korea

The complete mitochondrial genome of an endangered triton snail Charonia lampas (Littorinimorpha: Charoniidae) from South Korea was determined with the size of 15,330 bp, which contained 13 protein-coding genes, 22 tRNAs, and two rRNAs. The overall features of the mitochondrial genome are similar to those of typical gastropod mitochondrial genomes, except for the absence of a control region (CR): a conserved gene order, a high A + T content (67.45% for C. lampas), preference for AT-rich codons, etc. It is the second mitochondrial genome reported from this species after Cho et al. (2017), and there is ca. 98% sequence similarity between those of the species. According to the maximum-likelihood tree, C. lampas placed within the monophyletic Tonnoidea, forming a monoclade with Cypraeoidea and Stromboidea within Littorinimorpha. Interestingly, we confirmed that Vermetoidea is placed at the basal position of the monophyletic Littorinimorpha.

Characterization of Euplotes lynni nov. spec., E. indica nov. spec. and description of E. aediculatus and E. woodruffi (Ciliophora, Euplotidae) using an integrative approach

Four species belonging to the genus Euplotes have been investigated, namely: E. lynni nov. spec., E. indica nov. spec., E. aediculatus, and E. woodruffi. All populations are from India and were investigated using morphological and molecular markers. The phylogenetic relationships were inferred from small subunit ribosomal rRNA gene (SSU rRNA), internal transcribed spacer (ITS) region, and mitochondrial cytochrome c oxidase subunit I (COI) gene. Predicted secondary structure models for two new species using the hypervariable region of the SSU rRNA gene and ITS2 region support the distinctness of both species. Morphological characters were subjected to principal component analysis (PCA) and genetic variations were studied in-depth to analyze the relatedness of the two new species with their congeners. An integrative approach combining morphological features, molecular analysis, and ecological characteristics was carried out to understand the phylogenetic position of the reported species within the different clades of the genus Euplotes.

Molecular forms of Anopheles subpictus and Anopheles sundaicus in the Indian subcontinent

Background

Anopheles subpictus and Anopheles sundaicus are closely related species, each comprising several sibling species. Ambiguities exist in the classification of these two nominal species and the specific status of members of these species complexes. Identifying fixed molecular forms and mapping their spatial distribution will help in resolving the taxonomic ambiguities and understanding their relative epidemiological significance.

Methods

DNA sequencing of Internal Transcribed Spacer-2 (ITS2), 28S-rDNA (D1-to-D3 domains) and cytochrome oxidase-II (COII) of morphologically identified specimens of two nominal species, An. subpictus sensu lato (s.l.) and An. sundaicus s.l., collected from the Indian subcontinent, was performed and subjected to genetic distance and molecular phylogenetic analyses.

Results

Molecular characterization of mosquitoes for rDNA revealed the presence of two molecular forms of An. sundaicus s.l. and three molecular forms of An. subpictus s.l. (provisionally designated as Form A, B and C) in the Indian subcontinent. Phylogenetic analyses revealed two distinct clades: (i) subpictus clade, with a single molecular form of An. subpictus (Form A) prevalent in mainland India and Sri Lanka, and (ii) sundaicus clade, comprising of members of Sundaicus Complex, two molecular forms of An. subpictus s.l. (Form B and C), prevalent in coastal areas or islands in Indian subcontinent, and molecular forms of An. subpictus s.l. reported from Thailand and Indonesia. Based on the number of float-ridges on eggs, all An. subpictus molecular Form B were classified as Species B whereas majority (80%) of the molecular Form A were classified as sibling species C. Fixed intragenomic sequence variation in ITS2 with the presence of two haplotypes was found in molecular Form A throughout its distribution.

Conclusion

A total of three molecular forms of An. subpictus s.l. and two molecular forms of An. sundaicus s.l. were recorded in the Indian subcontinent. Phylogenetically, two forms of An. subpictus s.l. (Form B and C) prevalent in coastal areas or islands in the Indian subcontinent and molecular forms reported from Southeast Asia are members of Sundaicus Complex. Molecular Form A of An. subpictus is distantly related to all other forms and deserve a distinct specific status.

Two Distinct Genotypes of Spissistilus festinus (Say, 1830) (Hemiptera, Membracidae) in the United States Revealed by Phylogenetic and Morphological Analyses

Spissistilus festinus (Say, 1830) (Hemiptera: Membracidae) is a frequent pest of leguminous crops in the Southern United States, and a vector of grapevine red blotch virus. There is currently no information on the genetic diversity of S. festinus. In this study, populations of S. festinus were collected in 2015-2017 from various crops and geographic locations in the United States, and fragments of the mitochondrial cytochrome C oxidase 1 (mt-COI) gene and the nuclear internal transcribed spacer 2 (ITS2) region were characterized by polymerase chain reaction and sequencing. Maximum-likelihood and Bayesian analyses of the mt-COI and ITS2 sequences yielded similar phylogenetic tree topologies, revealing two distinct genetic S. festinus lineages with all of the specimens from California comprising one phylogenetic clade, alongside a single GenBank entry from Arizona, and all specimens from the Southeastern United States comprising a statistically-supported distinct clade, regardless of host and year of collection. The mt-COI gene fragment showed up to 10.8% genetic distance between the two phylogenetic clades. These results suggest the existence of two genotypes within S. festinus in the United States. The only distinct morphological trait between the two genotypes was a less elevated pronotum in the representative specimens from California, compared to the representative specimens from the Southeastern United States. Since this phenotypic feature is inconspicuous, a diagnostic polymerase chain reaction targeting a variable region of the mt-COI fragment was developed to reliably distinguish between the specimens of the two genotypes of S. festinus and to facilitate their specific identification.

Population genetic structure of the thick-tailed bushbaby (Otolemur crassicaudatus) from the Soutpansberg Mountain range, Northern South Africa, based on four mitochondrial DNA regions

Greater bushbabies, strepsirrhine primates, that are distributed across central, eastern and southern Africa, with northern and eastern South Africa representing the species' most southerly distribution. Greater bushbabies are habitat specialists whose naturally fragmented habitats are getting even more fragmented due to anthropogenic activities. Currently, there is no population genetic data or study published on the species. The aim of our study was to investigate the genetic variation in a thick-tailed bushbaby, Otolemur crassicaudatus, population in the Soutpansberg mountain range, Limpopo Province, South Africa. Four mitochondrial regions, ranging from highly conserved to highly variable, were sequenced from 47 individuals. The sequences were aligned and genetic diversity, structure, as well as demographic analyses were performed. Low genetic diversity (π = 0.0007-0.0038 in coding regions and π = 0.0127 in non-coding region; Hd = 0.166-0.569 in coding regions and Hd = 0.584 in non-coding region) and sub-structuring (H = 2-3 in coding regions and H = 4 in non-coding region) was observed with two divergent haplogroups (haplotype pairwise distance = 3-5 in coding region and 6-10 in non-coding region) being identified. This suggests the population may have experienced fixation of mitochondrial haplotypes due to limited female immigration, which is consistent with philopatric species, that alternative haplotypes are not native to this population, and that there may be male mobility from adjacent populations. This study provides the first detailed insights into the mitochondrial genetic diversity of a continental African strepsirrhine primate and demonstrates the utility of mitochondrial DNA in intraspecific genetic population analyses of these primates.

Molecular and Histopathological Data on Levisunguis subaequalis Curran, Overstreet, Collins & Benz, 2014 (Pentastomida: Eupentastomida: Porocephalida: Porocephaloidea: Sebekidae: Sebekinae) from Gambusia affinis in Alabama, USA

Levisunguis subaequalis Curran, Overstreet, Collins & Benz, 2014 , was recently described from the lungs of the definitive hosts, softshell turtles, Apalone ferox (Schneider, 1783), and Apalone spinifera aspera (Agassiz, 1857) as well as the viscera of an intermediate host, the western mosquitofish, Gambusia affinis (Baird and Girard, 1853). However, the original account lacked molecular data. Furthermore, histological examination of infected host tissues in the original account of L. subaequalis did not reveal any pathological changes in the intermediate host. The present work provides a robust morphological description of the nymph and novel molecular data from the 18S and 28S ribosomal gene regions and the cytochrome c oxidase subunit 1 (COI) mitochondrial gene. Phylogenetic analyses using Bayesian inference and maximum likelihood analysis with concatenated sequence data from these 3 regions, as well as each region individually, placed the turtle pentastomid L. subaequalis as a sister clade to the crocodilian pentastomids of the genus Sebekia Sambon, 1922. While only concatenated phylogenetic analyses agreed with the currently accepted classification of the Eupentastomida and phylogenetic signal assessment indicated that the concatenated data set yielded the most phylogenetic signal, data from more taxa are still needed for robust phylogenetic inferences to be made. The intensity of infection ranged from 2 to 171 nymphs per fish, compared with the highest previously reported intensity of 6. These high-intensity infections with L. subaequalis were characterized by the nymphs occupying 5-50% of the coelomic cavity of G. affinis. However, despite this heavy parasite infection, fish exhibited minimal pathology. Observed pathology was characterized by compression or effacement of organs adjacent to the nymphs, particularly liver, swim bladder, and intestines, as well as the formation of granulomas around shed pentastomid cuticles. Nonetheless, the morphological and molecular data provided in the present work will bolster future efforts to identify this pentastomid in other hosts where pathology may be present in addition to aiding in the advancement of the field of molecular pentastomid systematics.

First molecular assessment of the interrelationships of cladorchiid digeneans (Digenea: Paramphistomoidea), parasites of Neotropical fishes, including descriptions of three new species and new host and geographical records

The first molecular assessment of phylogenetic relationships of cladorchiid digeneans (superfamily Paramphistomoidea Fischoeder, 1901) from freshwater fishes based on 28S rDNA, ITS2 and cox1 sequences reveals the subfamilies Dadayiinae Fukui, 1929 and Kalitrematinae Travassos, 1933 as non-monophyletic, whereas Dadaytrema Travassos, 1931 represented by three species is monophyletic. Fourteen species of cladorchiids were found in characiform, perciform and siluriform fishes in the Neotropical Region (Brazil and Peru), with numerous new host and geographical records. The first scanning electron micrographs of seven species are presented. Two new species of dadayiine and one new species of kalitrematine paramphistomes are described. Microrchis macrovarium sp. n. from Pimelodella cristata (Müller et Troschel) (type host), Tetranematichthys quadrifilis (Kner) and Pterodoras granulosus (Valenciennes) in Brazil and Peru differs from all three congeners in that the testes are directly tandem, not oblique, and that the ovary, between the caecal ends, is widely separated from the testes. Pronamphistoma philippei sp. n. from Heros sp. in Brazil is distinguished from the type and only species, Pronamphistoma cichlasomae Thatcher, 1992, by the absence of the anterior collar-like expansion present in the type species, the presence of extramural rather than intramural pharyngeal sacs, and the unusual development of the dorsal and ventral exterior circular muscle fibre series in the acetabulum. Pseudocladorchis romani sp. n. from P. granulosus (type host), Brachyplatystoma vaillantii (Valenciennes), Calophysus macropterus (Lichtenstein), Megalodoras uranoscopus (Eigenmann et Eigenmann) and Oxydoras niger (Valenciennes) in Brazil and Peru, is most similar to Pseudocladorchis nephrodorchis Daday, 1907 but differs in the shape of the testes (irregular, versus reniform in the latter species) and the size of the ovary (as large as, or larger than, the testes in the new species). The generic diagnosis of Pronamphistoma Thatcher, 1992 is amended. Dadaytremoides parauchenipteri Lunaschi, 1989 is transferred to Doradamphistoma Thatcher, 1979 as D. parauchenipteri (Lunaschi, 1989) comb. n. based on morphological and molecular evidence.

Identification and Phylogenetic Classification of Fasciola species Isolated from Sheep and Cattle by PCR-RFLP in Zabol, in Sistan and Baluchistan Province, Southeast Iran

BACKGROUND

The detection of Fasciola species in various geographical regions is essential for health policymaking. Here, we aimed to identify livestock (cattle and sheep) related Fasciola genotypes by restriction fragment length polymorphism PCR.

METHODS

Seventy adult Fasciola flukes were collected from 70 infected livers of 35 cattle and 35 sheep slaughtered in Zabol abattoir, south-east Iran (Jan-Jul 2017). Fasciola species were determined based on molecular features. For molecular detection, Fasciola ITS1 region was amplified and sequenced. A 700 bp fragment was amplified. These were digested with RasΙ enzyme. F. hepatica specific fragments were 47, 59, 68, 104, and 370, while those related to F. gigantica had 45, 55, 170, 370.

RESULTS

The two main species of F. hepatica and F. gigantica are responsible for fasciolosis in sheep and cattle in our region. From 35 Fasciola isolated from cattle, 3 and 32 were F. hepatica and F. giagantica respectively. From 35 Fasciola isolated from sheep, 4 were F. hepatica and 31 were F. gigantica.

CONCLUSION

All Seventy Fasciola samples from two different hosts (cattle and sheep) were identified as either F. hepatica or F. gigantica by PCR-RFLP. Genotypic variability of Fasciola species was high in our region. It is recommended to assess molecular variation of Fasciola isolates in other host livestock.

Prevalence of Animal Fasciolosis and Specification of Fasciola spp. Isolated from Sheep, Goats and Cattle by Molecular Method: Hamadan Province, West of Iran

BACKGROUND

Fascioliasis is a common disease among humans and animals. Having global distribution, disease is developed by hepatic trematodes, Fasciola hepatica and F. gigantica. The main objective of this research was determining the prevalence of Fasciola species in Hamadan livestock and identifying those using PCR-RFLP.

METHODS

Overall, 13607 livestock livers in the slaughterhouse of Hamadan, west of Iran including 10846 sheep, 995 cattle and 1766 goats were examined in 2015. In addition, 75 Fasciola (41 worms from sheep, 22 worms from goats and 12 worms from cattle) were examined by PCR-RFLP method.

RESULTS

Totally, 100 livers were infected to Fasciola species (total prevalence: 0.74%; sheep 0.5%, goats 1.4%, cattle 1.5%). In the molecular results, prevalence of F. hepatica was higher than (92.5% in the sheep) F. gigantic, and in the cattle 91.5% and in the goats was 54.5% F. hepatica. Genotyping identified species confirmed intermediated types as F. hepatica.

CONCLUSION

Prevalence of Fasciola in this province is not so high. Intermediate types identified by PCR-RFLP method determined as F. hepatica by nucleotide sequencing. Because of morphological differences and interspecies variety, the accurate identification of Fasciola species needs using nucleotide sequencing.

DNA Barcoding of Metazoan Zooplankton Copepods from South Korea

Copepods, small aquatic crustaceans, are the most abundant metazoan zooplankton and outnumber every other group of multicellular animals on earth. In spite of ecological and biological importance in aquatic environment, their morphological plasticity, originated from their various lifestyles and their incomparable capacity to adapt to a variety of environments, has made the identification of species challenging, even for expert taxonomists. Molecular approaches to species identification have allowed rapid detection, discrimination, and identification of cryptic or sibling species based on DNA sequence data. We examined sequence variation of a partial mitochondrial cytochrome C oxidase I gene (COI) from 133 copepod individuals collected from the Korean Peninsula, in order to identify and discriminate 94 copepod species covering six copepod orders of Calanoida, Cyclopoida, Harpacticoida, Monstrilloida, Poecilostomatoida and Siphonostomatoida. The results showed that there exists a clear gap with ca. 20 fold difference between the averages of within-specific sequence divergence (2.42%) and that of between-specific sequence divergence (42.79%) in COI, suggesting the plausible utility of this gene in delimitating copepod species. The results showed, with the COI barcoding data among 94 copepod species, that a copepod species could be distinguished from the others very clearly, only with four exceptions as followings: Mesocyclops dissimilis–Mesocyclops pehpeiensis (0.26% K2P distance in percent) and Oithona davisae–Oithona similis (1.1%) in Cyclopoida, Ostrincola japonica–Pseudomyicola spinosus (1.5%) in Poecilostomatoida, and Hatschekia japonica–Caligus quadratus (5.2%) in Siphonostomatoida. Thus, it strongly indicated that COI may be a useful tool in identifying various copepod species and make an initial progress toward the construction of a comprehensive DNA barcode database for copepods inhabiting the Korean Peninsula.

Genetic variants of Kudoa septempunctata (Myxozoa: Multivalvulida), a flounder parasite causing foodborne disease

Foodborne disease outbreaks caused by raw olive flounders (Paralichthys olivaceus) parasitized with Kudoa septempunctata have been reported in Japan. Origins of olive flounders consumed in Japan vary, being either domestic or imported, and aquaculture-raised or natural. Although it is unknown whether different sources are associated with different outcomes, it is desirable to identify whether this is the case by determining whether unique K. septempunctata strains occur and if so, whether some are associated with foodborne illness. We here developed an intraspecific genotyping method, using the sequence variation of mitochondrial genes. We collected olive flounder samples from foodborne disease outbreaks, domestic fish farms or quarantine offices and investigated whether K. septempunctata genotype is associated with pathogenicity or geographic origin. The 104 samples were classified into three genotypes, ST1, ST2 and ST3. Frequency of symptomatic cases differed by genotypes, but the association was not statistically significant. Whereas K. septempunctata detected from aquaculture-raised and natural fish from Japan were either ST1 or ST2, those from fish inspected at quarantine from Korea to Japan were ST3. Our method can be applied to phylogeographic analysis of K. septempunctata and contribute to containing the foodborne disease. The genotype database is hosted in the PubMLST website (http://pubmlst.org/kseptempunctata/).

Molecular characterization ofOpisthorchis noverca(Digenea: Opisthorchiidae) based on nuclear ribosomal ITS2 and mitochondrial COI genes

Opisthorchiasis is a public health problem in South-East Asian countries and Eastern Europe. The infection implicates mainly two species ofOpisthorchis, namelyO. viverriniandO. felineus,that occur mostly in fish-eating mammals and humans, although there are rare reports of human cases involving two other species,O. novercaandO. guayaquilensis.Opisthorchis novercahas been reported frequently in dogs and pigs from the Indian subcontinent, with rare reports from cattle and human subjects. With a view to supplementing morphology-based identification of this species, the present study aimed to provide molecular characterization ofO. noverca, using rDNA internal transcribed spacer 2 (ITS2) and mitochondrial cytochrome oxidase I (mt COI) markers so as to determine its genetic correlation with other species of Opisthorchiidae, and also to generate a taxon-specific molecular marker based on the ITS2 region. The phylogenetic relationship betweenO. novercaand other species of the genus was determined using molecular sequence data. To strengthen the result, secondary structure sequence analyses of ITS2 with hemi-compensatory base changes (hCBCs), and amino acid sequence analyses, were also evaluated. Our results confirm thatO. novercais a distinct and valid species.

Speciation Progress: A Case Study on the Bushcricket Poecilimon veluchianus

Different mechanisms such as selection or genetic drift permitted e.g. by geographical isolation can lead to differentiation of populations and could cause subsequent speciation. The two subspecies of Poecilimon veluchianus, a bushcricket endemic to central Greece, show a parapatric distribution and are partially reproductively isolated. Therefore, P. veluchianus is suitable to investigate an ongoing speciation process. We based our analysis on sequences of the internal transcribed spacer (ITS) and the mitochondrial control region (CR). The population genetic analysis based on the nuclear marker ITS revealed a barrier to gene flow within the range of Poecilimon veluchianus, which corresponds well to the described subspecies. In contrast to the results based on the nuclear ITS marker, the mitochondrial CR marker does not clearly support the separation into two subspecies with restricted gene flow and a clear contact zone. Furthermore, we could identify isolation by distance (IBD) as one important mechanism responsible for the observed genetic structure (based on the ITS marker). The population genetic analysis based on the nuclear marker ITS also suggests the existence of hybrids in the wild. Furthermore, the simultaneous lack of strong prezygotic barriers and the presence of postzygotic mating barriers, observed in previous laboratory experiments, suggest that a secondary contact after an allopatric phase is more likely than parapatric speciation.

The Mitochondrial Genomes of a Myxozoan Genus Kudoa Are Extremely Divergent in Metazoa

The Myxozoa are oligo-cellular parasites with alternate hosts--fish and annelid worms--and some myxozoan species harm farmed fish. The phylum Myxozoa, comprising 2,100 species, was difficult to position in the tree of life, due to its fast evolutionary rate. Recent phylogenomic studies utilizing an extensive number of nuclear-encoded genes have confirmed that Myxozoans belong to Cnidaria. Nevertheless, the evolution of parasitism and extreme body simplification in Myxozoa is not well understood, and no myxozoan mitochondrial DNA sequence has been reported to date. To further elucidate the evolution of Myxozoa, we sequenced the mitochondrial genomes of the myxozoan species Kudoa septempunctata, K. hexapunctata and K. iwatai and compared them with those of other metazoans. The Kudoa mitochondrial genomes code for ribosomal RNAs, transfer RNAs, eight proteins for oxidative phosphorylation and three proteins of unknown function, and they are among the metazoan mitochondrial genomes coding the fewest proteins. The mitochondrial-encoded proteins were extremely divergent, exhibiting the fastest evolutionary rate in Metazoa. Nevertheless, the dN/dS ratios of the protein genes in genus Kudoa were approximately 0.1 and similar to other cnidarians, indicating that the genes are under negative selection. Despite the divergent genetic content, active oxidative phosphorylation was indicated by the transcriptome, metabolism and structure of mitochondria in K. septempunctata. As possible causes, we attributed the divergence to the population genetic characteristics shared between the two most divergent clades, Ctenophora and Myxozoa, and to the parasitic lifestyle of Myxozoa. The fast-evolving, functional mitochondria of the genus Kudoa expanded our understanding of metazoan mitochondrial evolution.

Molecular phylogeny of commercially important lobster species from Indian coast inferred from mitochondrial and nuclear DNA sequences

Lobsters constitute low-volume high-value crustacean fishery resource along Indian coast. For the conservation and management of this declining resource, accurate identification of species and larvae is essential. The objectives of this work were to generate species-specific molecular signatures of 11 commercially important species of lobsters of families Palinuridae and Scyllaridae and to reconstruct a phylogeny to clarify the evolutionary relationships among genera and species included in this study. Partial sequences were generated for all the candidate species from sampling sites along the Indian coast using markers like Cytochrome oxidase I (COI), 16SrRNA, 12SrRNA, and 18SrRNA genes, and analyzed. The genetic identities of widely distributed Thenus species along the Indian coast to be Thenus unimaculatus and the sub-species of Panulirus homarus to be P. homarus homarus were confirmed. Phylogeny reconstruction using the individual gene and concatenated mtDNA data set were carried out. The overall results suggested independent monophyly of Scyllaridae and Stridentes of Palinuridae. The interspecific divergence was found to be highest for the 12SrRNA compared with other genes. Significant incongruence between mtDNA and nuclear 18SrRNA gene tree topologies was observed. The results hinted an earlier origin for Palinuridae compared with Scyllaridae. The DNA sequence data generated from this study will aid in the correct identification of lobster larvae and will find application in research related to larval transport and distribution.

Phylogenetic information content of Copepoda ribosomal DNA repeat units: ITS1 and ITS2 impact

The utility of various regions of the ribosomal repeat unit for phylogenetic analysis was examined in 16 species representing four families, nine genera, and two orders of the subclass Copepoda (Crustacea). Fragments approximately 2000 bp in length containing the ribosomal DNA (rDNA) 18S and 28S gene fragments, the 5.8S gene, and the internal transcribed spacer regions I and II (ITS1 and ITS2) were amplified and analyzed. The DAMBE (Data Analysis in Molecular Biology and Evolution) software was used to analyze the saturation of nucleotide substitutions; this test revealed the suitability of both the 28S gene fragment and the ITS1/ITS2 rDNA regions for the reconstruction of phylogenetic trees. Distance (minimum evolution) and probabilistic (maximum likelihood, Bayesian) analyses of the data revealed that the 28S rDNA and the ITS1 and ITS2 regions are informative markers for inferring phylogenetic relationships among families of copepods and within the Cyclopidae family and associated genera. Split-graph analysis of concatenated ITS1/ITS2 rDNA regions of cyclopoid copepods suggested that the Mesocyclops, Thermocyclops, and Macrocyclops genera share complex evolutionary relationships. This study revealed that the ITS1 and ITS2 regions potentially represent different phylogenetic signals.

A Consistent PCR-RFLP Assay Based on ITS-2 Ribosomal DNA for Differentiation of Fasciola Species

OBJECTIVE(S)

Fascioliasis is a zoonotic parasitic disease caused by liver fluke species of Fasciola hepatica and Fasciola gigantica. Differentiation of these two species, based on their morphological characteristics, is difficult. The current study aimed to use PCR-RFLP assay to distinguish between F. hepatica and F. gigantica, based on profiles of RFLP, produced by effect of endonucleases on ITS2 of the ribosomal DNA genes from these two species.

MATERIALS AND METHODS

Adult Fasciola spp. were isolated from bile duct of naturally infected animals. The species of Fasciola were confirmed by sequencing the 505 bp region of the ITS2 gene in the isolates. By running the sequences of the samples in NEBcutter, suitable restriction enzymes (MspI and KpnI) were selected. Eight F. gigantica and eighteen F. hepatica samples were evaluated.

RESULTS

While RFLP pattern with MspI produced a profile by which it was difficult to differentiate these two species, KpnI along with MspI, produced a consistent pattern of a 231, 212 and 93 bp fragments in F. hepatica. This pattern was not seen in F. gigantica.

CONCLUSION

Findings of this study demonstrated that RFLP with KpnI and MspI produce a suitable pattern which simply differentiates F. hepatica from F. gigantica.

Utility of divergent domains of 28S ribosomal RNA in species discrimination of paramphistomes (Trematoda: Digenea: Paramphistomoidea)

Among the digenetic trematodes, paramphistomes are known to be the causative agent of "amphistomiasis" or the stomach fluke disease of domestic and wild animals, mainly ruminants. The use of 28S (divergent domains) and 18S rRNA for phylogenetic inference is significantly warranted for these flukes since it is as yet limited to merely the exploration of the second internal transcribed spacer (ITS2) region. The present study intended to explore the divergent domains (D1-D3) of 28S rRNA and simultaneously equate the phylogenetic information with 18S rRNA in paramphistomes. Divergence of the 28S rRNA domains was evident amongst the divergent (D) domains, where D1 domain emerged as the most variable and D2, the most robust domain, since the latter could provide a higher resolution of the species. D2 was the only domain that comprised compensatory mutations in the helices of its structural constraints; this domain is thus well suited for species distinction and may be considered a potential DNA barcode complementary to mitochondrial DNA. 28S (D1 + D2 + D3) rRNA provided a significant resolution of the taxa corroborating with the taxonomy of these flukes and thus proved to be more robust as a phylogenetic marker for lower levels than 18S rRNA. Phylogenetic inferences of paramphitomes are still scarcely explored; additional data from other taxa belonging to this family may estimate better the biodiversity of these flukes.

Low genetic diversity in wide-spread Eurasian liver fluke Opisthorchis felineus suggests special demographic history of this trematode species

Opisthorchis felineus or Siberian liver fluke is a trematode parasite (Opisthorchiidae) that infects the hepato-biliary system of humans and other mammals. Despite its public health significance, this wide-spread Eurasian species is one of the most poorly studied human liver flukes and nothing is known about its population genetic structure and demographic history. In this paper, we attempt to fill this gap for the first time and to explore the genetic diversity in O. felineus populations from Eastern Europe (Ukraine, European part of Russia), Northern Asia (Siberia) and Central Asia (Northern Kazakhstan). Analysis of marker DNA fragments from O. felineus mitochondrial cytochrome c oxidase subunit 1 and 3 (cox1, cox3) and nuclear rDNA internal transcribed spacer 1 (ITS1) sequences revealed that genetic diversity is very low across the large geographic range of this species. Microevolutionary processes in populations of trematodes may well be influenced by their peculiar biology. Nevertheless, we suggest that lack of population genetics structure observed in O. felineus can be primarily explained by the Pleistocene glacial events and subsequent sudden population growth from a very limited group of founders. Rapid range expansion of O. felineus through Asian and European territories after severe bottleneck points to a high dispersal potential of this trematode species.

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.

DNA barcoding methods for invertebrates

Invertebrates comprise approximately 34 phyla, while vertebrates represent one subphylum and insects a (very large) class. Thus, the clades excepting vertebrates and insects encompass almost all of animal diversity. Consequently, the barcoding challenge in invertebrates is that of barcoding animals in general. While standard extraction, cleaning, PCR methods, and universal primers work for many taxa, taxon-specific challenges arise because of the shear genetic and biochemical diversity present across the kingdom, and because problems arising as a result of this diversity, and solutions to them, are still poorly characterized for many metazoan clades. The objective of this chapter is to emphasize general approaches, and give practical advice for overcoming the diverse challenges that may be encountered across animal taxa, but we stop short of providing an exhaustive inventory. Rather, we encourage researchers, especially those working on poorly studied taxa, to carefully consider methodological issues presented below, when standard approaches perform poorly.

Complete mitochondrial genome of Acheilognathus signifer (Cypriniformes, Cyprinidae): comparison of light-strand replication origins from six Acheilognathus species

Acheilognathus signifer (Cypriniformes, Cyprinidae) is an endemic, endangered fresh water fish species in the Korean peninsula. In this study, the mitochondrial genome of A. signifer was completely sequenced, which is 16,566 bp in length. The characteristics of the complete mitochondrial genome were described in detail.

Quick detection of Dryocosmus kuriphilus Yasumatsu (Hymenoptera: Cynipidae) in chestnut dormant buds by nested PCR

Dryocosmus kuriphilus Yasumatsu (Hymenoptera: Cynipidae) develops in chestnut buds that remain asymptomatic from oviposition (June-July) until budburst; it is, thus, easily spread by plant material used in propagation. Therefore, it is particularly interesting to identify infested plant batches before their movement. Unfortunately, a non-destructive method for checking buds has not yet been developed, and the only technique available is the screening of a bud sample. The visual investigation is long and requires highly skilled and trained staff. The purpose of this work was to set up an effective and fast method able to identify the presence of first instar larvae of D. kuriphilus in a large number of chestnut buds by PCR. Four primer pairs were designed on nuclear and mitochondrial sequences of a set of seven gall wasp taxa and tested on five different cynipid's DNA. Nested diagnostic PCR was carried out on DNA extracted from samples of 2 g buds simulating four levels of infestation (larvae were added to uninfested buds); 320 bp amplicon of 28S sequence was chosen as a marker to detect one larva out of 2 g buds. The method showed a potential efficiency of 5000 to 15,000 buds per week, depending on bud size.