Secondary structure model for the ITS-2 precursor rRNA of strongyloid nematodes of equids: implications for phylogenetic inference

Genetic variability, cryptic species and phylogenetic relationship of six cyathostomin species based on mitochondrial and nuclear sequences

Cyathostomins are important intestinal nematode parasites of equines and include 50 accepted species. Their taxonomy has been frequently revised and the presence of cryptic species suggested. Furthermore, usually molecular- and morphology-based phylogenetic analyses give divergent results. In this study, the nucleotide sequences of the nuclear second internal transcribed spacer (ITS-2) and the mitochondrial partial cytochrome c oxidase subunit I (COI) were determined for adults of six cyathostomin species (Coronocyclus coronatus, Coronocyclus labiatus, Cylicocyclus nassatus, Cylicostephanus calicatus, Cylicostephanus longibursatus, Cylicostephanus minutus) collected from different equine species within two geographic regions. Maximum likelihood trees were calculated for ITS-2, COI, and concatenated data. No obvious differentiation was observed between geographic regions or equine host species. As previously reported, Coronocyclus coronatus and Cylicostephanus calicatus revealed a close relationship. Cryptic species were detected in Cylicostephanus minutus and Cylicostephanus calicatus. Cylicocyclus nassatus and Coronocyclus labiatus showed diverse mitochondrial and nuclear haplotypes occurring in different combinations, while Cylicostephanus longibursatus was comparatively homogenous. In conclusion, a combined analysis of nuclear and mitochondrial haplotypes improved resolution of the phylogeny and should be applied to the remaining cyathostomin species and across additional equine host species and geographic regions.

Nuclear and mitochondrial marker sequences reveal close relationship between Coronocyclus coronatus and a potential Cylicostephanus calicatus cryptic species complex

The Cyathostominae (Nematoda, Strongyloidea) parasitising equines represent a diverse group currently including 50 species. However, their taxonomy has been repeatedly revised and occasionally the presence of cryptic genospecies was suggested. Moreover, molecular- and morphology-based phylogenetic analyses give divergent results. For instance, molecular data have suggested close relationship between Coronocyclus coronatus and Cylicostephanus calicatus, although morphology-based taxonomy places them in different genera. Here, nuclear (internal transcribed spacer 2, ITS-2) and mitochondrial (cytochrome oxidase I, COI) sequences were obtained from the same individual, morphologically identified worms. In both morphospecies, two ITS-2 sequences types were observed: In Cor. coronatus, a small PCR product of 278 bp (nuclear haplotype group nHGBco) was always present but often in combination with a larger 369-370 bp fragment (nHGAco). In Cys. calicatus, either a large 370 bp product (nHGAca) or a short 281 bp amplicon (nHGBca) were found, but never both. Sequence identity between morphospecies was up to 100%. The smaller differed from the larger fragments by deletion of the region 110-198 in Cor. coronatus and 112-203 in Cys. calicatus. In COI, three and five mitochondrial haplotype groups (HGs), mtHG1co-mtHG3co and mtHG1ca-mtHG5ca were identified for Cor. coronatus and Cys. calicatus, respectively. In Cor. coronatus, there was no particular association of mtHG with nuclear genotypes (only nHGBco vs. both nHGBco plus nHGAco). In Cys. calicatus the nHGAca was always associated with the mtHG1ca, mtHG2ca or mtHG5ca whereas nHGBca was exclusively associated with mtHG3ca or mtHG4ca. Despite up to 100% identity in the nHGs, no mixing of mtHGs was observed between both species. Clear separation of certain nHGs with particular mtHGs in Cys. calicatus, despite the fact that the same host individuals were infected with both groups simultaneously, suggests presence of two non-interbreeding genospecies within Cys. calicatus, which needs further confirmation using additional samples from diverse geographical origins.

Systematic review of gastrointestinal nematodes of horses from Australia

Background

Equine gastrointestinal nematodes (GINs) have been the subject of intermittent studies in Australia over the past few decades. However, comprehensive information on the epidemiology of equine GINs, the efficacy of available anthelmintic drugs and the prevalence of anthelmintic resistance (AR) in Australasia is lacking. Herein, we have systematically reviewed existing knowledge on the horse GINs recorded in Australia, and main aspects of their pathogeneses, epidemiology, diagnoses, treatment and control.

Methods

Six electronic databases were searched for publications on GINs of Australian horses that met our inclusion criteria for the systematic review. Subsets of publications were subjected to review epidemiology, diagnoses, pathogeneses, treatment and control of GINs of horses from Australia.

Results

A total of 51 articles published between 1950 to 2018 were included. The main GINs reported in Australian horses were cyathostomins (at least 28 species), Draschia megastoma, Habronema muscae, H. majus, Oxyuris equi, Parascaris equorum, Strongyloides westeri and Trichostrongylus axei across different climatic regions of Queensland, New South Wales, Victoria, and Western Australia. Nematodes are diagnosed based on the traditional McMaster egg counting technique, though molecular markers to characterise common GINs of equines were characterised in 1990s. The use of anthelmintic drugs remains the most widely-used strategy for controlling equine GIN parasites in Australia; however, the threshold of faecal egg count that should trigger treatment in horses, remains controversial. Furthermore, anthelmintic resistance within GIN population of horses is becoming a common problem in Australia.

Conclusions

Although GINs infecting Australian horses have been the subject of occasional studies over the past few decades, the effective control of GIN infections is hampered by a generalised lack of knowledge in various disciplines of equine parasitology. Therefore, coordinated and focused research is required to fill our knowledge gaps in these areas to maximise equine health and minimise economic losses associated with the parasitic infections in Australia.

Electronic supplementary material

The online version of this article (10.1186/s13071-019-3445-4) contains supplementary material, which is available to authorized users.

Analysis of genetic variation in Globocephaloides populations from macropodid marsupials using a mutation scanning-based approach

Three species of Globocephaloides, parasitic nematodes occurring in macropodid marsupials in different areas of Australia, were characterized by the sequences of the first and second internal transcribed spacers (ITS-1 and ITS-2) of nuclear ribosomal DNA. Samples were subjected to PCR-coupled SSCP analysis and targeted sequencing, in order to assess genetic variation within and among individuals from different host species. Both SSCP and sequence data supported the current classification of morphospecies. Contrary to a previous hypothesis that cryptic species exist within the Globocephaloides trifidospicularis "complex", no or minor (0.2%) variation was detected among individuals from different hosts or geographical origins. Within G. macropodis populations, there was a consistent difference in both the ITS-1 (5.2%) and ITS-2 (7.1%) sequences between individuals derived from Macropus agilis and those from Macropus dorsalis. Although the results suggested that G. macropodis from each host species represented sibling species, future morphological study of individuals representing each G. macropodis genotype is warranted to provide further support for this hypothesis. (Nucleotide sequence data have been deposited in the GenBank database under accession nos. GQ131400-GQ131409.).

Comparative analysis of sequences and secondary structures of the rRNA internal transcribed spacer 2 (ITS2) in pollen beetles of the subfamily Meligethinae (Coleoptera, Nitidulidae): potential use of slippage-derived sequences in molecular systematics

A comparative analysis of ITS2 sequences and secondary structures in 89 species of pollen beetles of the subfamily Meligethinae (Coleoptera, Nitidulidae) was performed. The ITS2 folding pattern was highly conserved and comparable with the general model proposed for eukaryotes. Simple sequence repeats (SSRs) were responsible for most of the observed nucleotide variability (approximately 1-3%) and length variation (359-459bp). When plotted on secondary structures, SSRs mapped in expansion segments positioned at the apices of three ITS2 helices ('A', 'B' and 'D1') and appeared to have evolved under mechanisms of compensatory slippage. Homologies among SSRs nucleotides could not be unambiguously assigned, and thus were not useful to resolve phylogeny. However, slippage-derived motifs provided some preliminary genetic support for newly proposed taxonomic arrangements of several genera and subgenera of Meligethinae, corroborating existing morphological and ecological datasets.

Pan-eukaryote ITS2 homologies revealed by RNA secondary structure

For evolutionary comparisons, phylogenetics and evaluation of potential interbreeding taxa of a species, various loci have served for animals and plants and protistans. One [second internal transcribed spacer (ITS2) of the nuclear ribosomal DNA] is highly suitable for all. Its sequence is species specific. It has already been used extensively and very successfully for plants and some protistans, and a few animals (where historically, the mitochondrial genes have dominated species studies). Despite initial impressions that ITS2 is too variable, it has proven to provide useful biological information at higher taxonomic levels, even across all eukaryotes, thanks to the conserved aspects of its transcript secondary structure. The review of all eukaryote groups reveals that ITS2 is expandable, but always retains in its RNA transcript a common core structure of two helices with hallmark characteristics important for ribosomal RNA processing. This aspect of its RNA transcript secondary structure can rescue difficult alignment problems, making the ITS2 a more powerful tool for phylogenetics. Equally important, the recognition of eukaryote-wide homology regions provides extensive and detailed information to test experimental studies of ribosomal rRNA processing.

Nucleotide alterations in the D3 domain of the large subunit of ribosomal DNA among 21 species of equine strongyle

The expansion segments or divergent (D) domains in the large subunit (LSU) of the ribosomal DNA have been suggested as genetic markers for taxonomic and/or phylogenetic studies of parasites. In the present study, we assessed the degree of sequence variation in the D3 domain and flanking core regions of the LSU for 21 species of equine strongyles (Strongylida: Strongylidae) and determined which positions in the secondary structure of the LSU were associated with the nucleotide alterations. No intraspecific sequence variation was detected in 17 species, for which multiple individual worms were available. Mutations in sequence among species were detected at 19 nucleotide positions, most of which were located in the D3 domain. Fifteen alterations were transitions, three were transversions and one represented a site of multiple mutations. In relation to the secondary structure element of D3, 26% of these mutations were located in unpaired regions (i.e., end of loops, or in bulges of helices) and thus did not appear to alter the pairing arrangement in the helices of the secondary structure. Many of the other mutations represented partial or complete compensatory base pair changes. The magnitude of interspecific nucleotide variation in the D3 domain (0-4%) was considerably less than that recorded for some other nematode groups (enoplids and thelastomatoids), indicating that this region alone is of limited value for taxonomic and phylogenetic studies for strongyles of equids but is interesting in relation to the evolution of ribosomal DNA.

Analysis of the secondary structure of ITS1 in Pectinidae: implications for phylogenetic reconstruction and structural evolution

It is at present difficult to accurately position gaps in sequence alignment and to determine substructural homology in structure alignment when reconstructing phylogenies based on highly divergent sequences. Therefore, we have developed a new strategy for inferring phylogenies based on highly divergent sequences. In this new strategy, the whole secondary structure presented as a string in bracket notation is used as phylogenetic characters to infer phylogenetic relationships. It is no longer necessary to decompose the secondary structure into homologous substructural components. In this study, reliable phylogenetic relationships of eight species in Pectinidae were inferred from the structure alignment, but not from sequence alignment, even with the aid of structural information. The results suggest that this new strategy should be useful for inferring phylogenetic relationships based on highly divergent sequences. Moreover, the structural evolution of ITS1 in Pectinidae was also investigated. The whole ITS1 structure could be divided into four structural domains. Compensatory changes were found in all four structural domains. Structural motifs in these domains were identified further. These motifs, especially those in D2 and D3, may have important functions in the maturation of rRNAs.

Parasite histories and novel phylogenetic tools: Alternative approaches to inferring parasite evolution from molecular markers

Parasitological research is often contingent on the knowledge of the phylogeny/genealogy of the studied group. Although molecular phylogenetics has proved to be a powerful tool in such investigations, its application in the traditional fashion, based on a tree inference from the primary nucleotide sequences may, in many cases, be insufficient or even improper. These limitations are due to a number of factors, such as a scarcity/ambiguity of phylogenetic information in the sequences, an intricacy of gene relationships at low phylogenetic levels, or a lack of criteria when deciding among several competing coevolutionary scenarios. With respect to the importance of a precise and reliable phylogenetic background in many biological studies, attempts are being made to extend molecular phylogenetics with a variety of new data sources and methodologies. In this review, selected approaches potentially applicable to parasitological research are presented and their advantages as well as drawbacks are discussed. These issues include the usage of idiosyncratic markers (unique features with presumably low probability of homoplasy), such as insertion of mobile elements, gene rearrangements and secondary structure features; the problem of ancestral polymorphism and reticulate relationships at low phylogenetic levels; and the utility of a molecular clock to facilitate discrimination among alternative scenarios in host-parasite coevolution.

The use of nuclear ribosomal DNA markers for the identification of bursate nematodes (order Strongylida) and for the diagnosis of infections

Many bursate nematodes are of major importance to animal health. Animals are often parasitized by multiple species that differ in their prevalence, relative abundance and/or pathogenicity. Implementation of effective management strategies for these parasites requires reliable methods for their detection in hosts, identification to the species level and measurement of intensity of infection. One major problem is the difficulty of accurately identifying and distinguishing many species of bursate nematode because of the remarkable morphological similarity of their eggs and larvae. The inability to identify, with confidence, individual nematodes (irrespective of their life-cycle stage) to the species level by morphological methods has often led to a search for species-specific genetic markers. Studies over the past 15 years have shown that sequences of the internal transcribed spacers of ribosomal DNA provide useful genetic markers, providing the basis for the development of PCR-based diagnostic tools. Such molecular methods represent powerful tools for studying the systematics, epidemiology and ecology of bursate nematodes and, importantly, for the specific diagnosis of infections in animals and humans, thus contributing to improved control and prevention strategies for these parasites.

Advances in developing molecular-diagnostic tools for strongyloid nematodes of equids: fundamental and applied implications

Infections of equids with parasitic nematodes of the order Strongylida (subfamilies Strongylinae and Cyathostominae) are of major veterinary importance. In last decades, the widespread use of drugs against these parasites has led to problems of resistance within the Cyathostominae, and to an increase in their prevalence and intensity of infection. Novel control strategies, based on improved knowledge of parasite biology and epidemiology, have thus become important. However, there are substantial limitations in the understanding of fundamental biological and systematic aspects of these parasites, which have been due largely to limitations in their specific identification and diagnosis using traditional, morphological approaches. Recently, there has been progress in the development of DNA-based approaches for the specific identification of strongyloids of equids for systematic studies and disease diagnosis. The present article briefly reviews information on the classification, biology, pathogenesis, epidemiology of equine strongyloids and the diagnosis of infections, highlights knowledge gaps in these areas, describes recent advances in the use of molecular techniques for the genetic characterisation, specific identification and differentiation of strongyloids of equids as a basis for fundamental investigations of the systematics, population biology and ecology.

Phylogenetic analysis of the Metastrongyloidea (Nematoda: Strongylida) inferred from ribosomal RNA gene sequences

Phylogenetic relationships among nematodes of the strongylid superfamily Metastrongyloidea were analyzed using partial sequences from the large-subunit ribosomal RNA (LSU rRNA) and small-subunit ribosomal RNA (SSU rRNA) genes. Regions of nuclear ribosomal DNA (rDNA) were amplified by polymerase chain reaction, directly sequenced, aligned, and phylogenies inferred using maximum parsimony. Phylogenetic hypotheses inferred from the SSU rRNA gene supported the monophyly of representative taxa from each of the 7 currently accepted metastrongyloid families. Metastrongyloid taxa formed the sister group to representative trichostrongyloid sequences based on SSU data. Sequences from either the SSU or LSU RNA regions alone provided poor resolution for relationships within the Metastrongyloidea. However, a combined analysis using sequences from all rDNA regions yielded 3 equally parsimonious trees that represented the abursate Filaroididae as polyphyletic, Parafilaroides decorus as the sister species to the monophyletic Pseudaliidae, and a sister group relationship between Oslerus osleri and Metastrongylus salmi. Relationships among 3 members of the Crenosomatidae, and 1 representative of the Skrjabingylidae (Skrjabingylus chitwoodorum) were not resolved by these combined data. However, members of both these groups were consistently resolved as the sister group to the other metastrongyloid families. These relationships are inconsistent with traditional classifications of the Metastrongyloidea and existing hypotheses for their evolution.

Redescription of Zoniolaimus mawsonae Beveridge, 1983 (Nematoda: Strongyloidea) and the description of Z. latebrosus n. sp. from the red kangaroo Macropus rufus (Marsupialia: Macropodidae) based on morphological and molecular data

This study shows that the description of Z. mawsonae given by Beveridge (1983) represented a composite of two species. Hence, Z. mawsonae Beveridge, 1983 is re-described and a new species, Z. latebrosus, is described. Males of the two species differ in the lengths of their spicules (0.94-1.23 mm in Z. mawsonae compared with 1.53-1.95 mm in Z. latebrosus) and in several characteristics of the bursa and genital cone. Females of the two species can be identified based on the shape of the posterior end of the body, that of Z. mawsonae being markedly swollen. Molecular data, obtained from single strand conformation polymorphism (SSCP) analysis and DNA sequencing of the second internal transcribed spacer of ribosomal DNA, provided additional support for the separation of the two species.

Review of Papillostrongylus Johnston & Mawson, 1939 (Nematoda: Strongyloidea) from wallabies and kangaroos (Marsupialia: Macropodidae) using morphological and molecular techniques, with the description of P. barbatus n. sp

The strongyloid nematode genus Papillostrongylus Johnston & Mawson, 1939, from kangaroos and wallabies, is reviewed using morphological and molecular methods. P. labiatus Johnston & Mawson, 1939 is re-described from material from the type-host, the black-striped wallaby Macropus dorsalis, from eastern Queensland, Australia, in which it is a relatively common parasite. Additional records from M. parryi and Thylogale thetis are confirmed and considered to represent examples of host-switching. A geographically disjunct population of the nematode species occurs in M. bernardus and Petrogale brachyotis in Arnhem Land, Northern Territory, but assessment of its status requires additional material. Nematodes from M. rufus, M. giganteus, M. fuliginosus and M. robustus from inland regions of Australia, formerly attributed to P. labiatus, are here assigned to a new species, P. barbatus, distinguished by the presence of an external leaf-crown, larger size, by greater spicule length in the male and by a sinuous vagina in the female. Additional hosts of P. barbatus n. sp. are Petrogale assimilis and Pet. lateralis purpureicollis. Sequence analyses of the second internal transcribed spacer of ribosomal DNA (ITS-2) also showed that P. barbatus n. sp. differed at 40 (16.7%) of the 240 alignment positions when compared with P. labiatus. Most of these interspecific sequence differences occurred in loops or bulges of the predicted precursor rRNA secondary structure, or represented partial or total compenstory base pair changes in stems.

Applications of single-strand conformation polymorphism (SSCP) to taxonomy, diagnosis, population genetics and molecular evolution of parasitic nematodes

The analysis of genetic variation in parasitic nematodes has important implications for studying aspects of taxonomy, diagnosis, population genetics, drug resistance and molecular evolution. This article highlights some applications of PCR-based single-strand conformation polymorphism (SSCP) for the analysis of sequence variation in individual parasites (and their populations) to address some of these areas. It also describes the principles and advantages of SSCP, and provides some examples for future applications in parasitology.

Phylogeny of the Trichostrongylina (Nematoda) inferred from 28S rDNA sequences

We produced a molecular phylogeny of species within the order Strongylida (bursate nematodes) using the D1 and D2 domains of 28S rDNA, with 23 new sequences for each domain. A first analysis using Caenorhabditis elegans as an outgroup produced a tree with low resolution in which three taxa (Dictyocaulus filaria, Dictyocaulus noerneri, and Metastrongylus pudendotectus) showed highly divergent sequences. In a second analysis, these three species and C. elegans were removed and an Ancylostomatina, Bunostomum trigonocephalum, was chosen (on the basis of previous morphological analyses) as the outgroup for an analysis of the phylogenetic relationships between and within the Strongylina (strongyles) and Trichostrongylina (trichostrongyles). A very robust tree was obtained. The Trichostrongylina were monophyletic, but the Strongylina were paraphyletic, though this requires confirmation. Within the Trichostrongylina, the three superfamilies defined from morphological characters are confirmed, with the Trichostrongyloidea sister group to a clade including the Molineoidea and Heligmosomoidea. Within the Trichostrongyloidea, the Cooperiidae, Trichostrongylidae, and Haemonchidae were polytomous, the Haemonchinae were monophyletic, but the Ostertagiinae were paraphyletic. The sister-group relationships between Molineoidea and Heligmosomoidea were unsuspected from previous morphological analysis. No unequivocal morphological synapomorphy could be found for the grouping Molineoidea + Heligmosomoidea, but none was found which contradicted it.

Single-strand conformation polymorphism analysis of genetic variation in Labiostrongylus longispicularis from kangaroos

Single-strand conformation polymorphism (SSCP) analysis was employed to screen for sequence heterogeneity in the second internal transcribed spacer (ITS-2) of ribosomal (r) DNA of Labiostrongylus longispicularis, a parasitic strongylid nematode occuring in some species of kangaroo in different geographical regions of Australia. The results showed that most of the nematodes screened had different SSCP profiles, which were subsequently shown to correspond to polymorphisms and/or an indel in the ITS-2 sequence. These variable sites related mainly to unpaired regions of the predicted secondary structure of the precursor rRNA molecule. SSCP profiles could be used to distinguish L. longispicularis in Macropus robustus robustus (New South Wales) from L. longispicularis in Macropus robustus erubescens and Macropus rufus (South Australia). This difference corresponded to a transversional change in the ITS-2 sequence at alignment position 82. The study demonstrated clearly the effectiveness of SSCP analysis for future large-scale population genetic studies of L. longispicularis in order to test the hypothesis that L. longispicularis from different geographical regions represents multiple sibling species.

Evolutionary relationships of trichostrongyloid nematodes (Strongylida) inferred from ribosomal DNA sequence data

The evolutionary relationships of 21 species of trichostrongyloid nematodes were determined by use of sequence data of the second internal transcribed spacer of the ribosomal DNA aligned according to secondary structure information. Irrespective of the method of analysis used, the topologies of the phylogenetic trees derived from the molecular data differed with respect to all four hypotheses proposed previously for the evolutionary relationships of the different subfamilies within the Trichostrongylidae based on morphological data. Thus, the molecular data set did not resolve the conflict between the four previous proposals for the subfamilial relationships. Nonetheless, all trees derived from the molecular data showed strong support for the exclusion of the genera Filarinema and Amidostomum from the clade containing the species within the family Trichostrongylidae. This represents a major difference from the most recent proposal of the systematics of the Trichostrongyloidea in which these two genera were included within the Trichostrongylidae. Therefore, the molecular data support an earlier systematic framework in which Filarinema and Amidostomum were considered to be sister groups of the Trichostrongyloidea.

Mutation scanning analysis of sequence heterogeneity in the second internal transcribed spacer (rDNA) within some members of theHypodontus macropi (Nematoda: Strongyloidea) complex

Single-strand conformation polymorphism analysis was employed to investigate sequence variation in the second internal transcribed spacer (ITS-2) of nuclear ribosomal DNA within and among individuals representing three operational taxonomic units (OTUs) of Hypodontus macropi from different species of Australian macropodid marsupials. Of the 96 nematodes analysed, totals of 3 (OTU1 from Petrogale persephone), 10 (OTU2 from Macropus robustus) and 7 (OTU9 from Macropus rufus) representative individuals were selected for DNA sequencing to characterise and estimate the magnitude of nucleotide variation in the ITS-2. While no unequivocal nucleotide difference in the ITS-2 was detectable within OTU1, most sequence variation (3/44.7%) detected within OTU2 and OTU9 was related chiefly to dinucleotide (CA, TA, or a combination of both) differences. This microsatellite variability in some H. macropi OTUs suggests that the ITS-2 rDNA may be subjected to slippage events during DNA replication, resulting in short dinucleotide repeat tracts being dispersed throughout the ITS-2 lineages, or possibly transposition and/or crossing-over events. Nucleotide variation in the ITS-2 of individual OTUs was related to the proposed secondary structure for the precursor ribosomal RNAs. Most of the sequence heterogeneity or polymorphism within OTU2 and OTU9 occurred in loops or bulges of the predicted secondary structure, which appear not to be under functional constraint. The findings of this study have implications for investigating speciation events and population differentiation in nematodes at the molecular level.

Genomic and genetic research on bursate nematodes: significance, implications and prospects

Molecular genetic research on parasitic nematodes (order Strongylida) is of major significance for many fundamental and applied areas of medical and veterinary parasitology. The advent of gene technology has led to some progress for this group of nematodes, particularly in studying parasite systematics, drug resistance and population genetics, and in the development of diagnostic assays and the characterisation of potential vaccine and drug targets. This paper gives an account of the molecular biology and genetics of strongylid nematodes, mainly of veterinary socio-economic importance, indicates the implications of such research and gives a perspective on genome research for this important parasite group, in light of recent technological advances and knowledge of the genomes of other metazoan organisms.