The complete mitochondrial genome of the New Zealand parasitic roundworm Teladorsagia circumcincta (Trichostrongyloidea: Haemonchidae) field strain NZ_Teci_NP

Multi-Omic Profiling, Structural Characterization, and Potent Inhibitor Screening of Evasion-Related Proteins of a Parasitic Nematode, Haemonchus contortus, Surviving Vaccine Treatment

The emergence of drug-resistant parasitic nematodes in both humans and livestock calls for development of alternative and cost-effective control strategies. Barbervax^® is the only registered vaccine for the economically important ruminant strongylid Haemonchus contortus. In this study, we compared the microbiome, genome-wide diversity, and transcriptome of H. contortus adult male populations that survived vaccination with an experimental vaccine after inoculation in sheep. Our genome-wide SNP analysis revealed 16 putative candidate vaccine evasion genes. However, we did not identify any evidence for changes in microbial community profiling based on the 16S rRNA gene sequencing results of the vaccine-surviving parasite populations. A total of fifty-eight genes were identified as significantly differentially expressed, with six genes being long non-coding (lnc) RNAs and none being putative candidate SNP-associated genes. The genes that highly upregulated in surviving parasites from vaccinated animals were associated with GO terms belonging to predominantly molecular functions and a few biological processes that may have facilitated evasion or potentially lessened the effect of the vaccine. These included five targets: astacin (ASTL), carbonate dehydratase (CA2), phospholipase A2 (PLA2), glutamine synthetase (GLUL), and fatty acid-binding protein (FABP3). Our tertiary structure predictions and modelling analyses were used to perform in silico searches of all published and commercially available inhibitor molecules or substrate analogs with potential broad-spectrum efficacy against nematodes of human and veterinary importance.

Sequencing and Reconstructing Helminth Mitochondrial Genomes Directly from Genomic Next-Generation Sequencing Data

We present a detailed method for extraction of high-molecular weight genomic DNA suitable for numerous DNA sequencing applications, and a straightforward in silico approach for reconstructing novel mitochondrial (mt) genomes directly from total genomic DNA extracts derived from next-generation sequencing (NGS) data sets. The in silico post-sequencing pipeline described is fast, accurate, and highly efficient, with modest memory requirements that can be performed using a standard desktop computer. The approach is particularly effective for obtaining mitochondrial genomes for species with little or no mitochondrial sequence information currently available and overcomes many of the limitations of traditional strategies. The described methodologies are also applicable for metagenomics sequencing from mixed or pooled samples containing multiple species and subsequent specific assembly of specific mitochondrial genomes.

Comparative mitogenomics of Spirocerca lupi from South Africa and China: Variation and possible heteroplasmy

The complete mitochondrial genome of Spirocerca lupi isolated from a dog in South Africa was sequenced using next generation sequencing (NGS) technology and the 12 protein coding genes along with the two rRNA genes were compared to 18 other nematode species as well as S. lupi from China. The mitochondrial genome of S. lupi South Africa had a mean genetic diversity of 6.1 % compared to S. lupi China with some variation in nucleotide composition, gene positioning and size. Pairwise distance results indicated slightly higher variation when compared to the pairwise distances of other closely related species, however, this variation was not high enough for it to be considered a cryptic species. Phylogenetic analysis indicated that S. lupi from the two continents are very similar. In addition, single nucleotide polymorphisms were detected in the nad2 gene with ten sequence variants identified from 10 clones from a single nematode, suggesting possible heteroplasmy. The origin of the heteroplasmy is currently unknown but it is speculated to have arisen from accumulated mutations in the mitochondria during somatic replication.

Comparative evaluation of different molecular methods for DNA extraction from individual Teladorsagia circumcincta nematodes

Background

The purpose of this study was to develop a reliable DNA extraction protocol to use on individual Teladorsagia circumcincta nematode specimens to produce high quality DNA for genome sequencing and phylogenetic analysis. Pooled samples have been critical in providing the groundwork for T. circumcincta genome construction, but there is currently no standard method for extracting high-quality DNA from individual nematodes. 11 extraction kits were compared based on DNA quality, yield, and processing time.

Results

11 extraction protocols were compared, and the concentration and purity of the extracted DNA was quantified. Median DNA concentration among all methods measured on NanoDrop 2000™ ranged between 0.45–11.5 ng/μL, and on Qubit™ ranged between undetectable – 0.962 ng/μL. Median A260/280 ranged between 0.505–3.925, and median A260/230 ranged − 0.005 – 1.545. Larval exsheathment to remove the nematode cuticle negatively impacted DNA concentration and purity.

Conclusions

A Schistosoma sp. DNA extraction method was determined as most suitable for individual T. circumcincta nematode specimens due to its resulting DNA concentration, purity, and relatively fast processing time.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12896-021-00695-6.

Characterization of the complete mitochondrial genome of the New Zealand parasitic blowfly Calliphora vicina (Insecta: Diptera: Calliphoridae)

In this study, the complete mitochondrial (mt) genome of the New Zealand parasitic blowfly Calliphora vicina (blue bottle blowfly) field strain NZ_CalVic_NP was generated using next-generation sequencing technology and annotated. The 16,518 bp mt genome consists of 13 protein-coding genes, two ribosomal RNAs, 22 transfer RNAs, and a 1689 bp non-coding region, similar to the two other available C. vicina and most metazoan mt genomes. Phylogenetic analysis showed that C. vicina NZ_CalVic_NP forms a monophyletic cluster with the remaining three Calliphorinae species. The complete mt genome sequence of C. vicina NZ_CalVic_NP is a resource to facilitate future species- and strain-level identification research and investigations into the evolutionary provenance within the Calliphoridae.

The complete mitochondrial genome of the New Zealand parasitic blowfly Lucilia sericata (Insecta: Diptera: Calliphoridae)

In this study, the complete mitochondrial (mt) genome of the New Zealand parasitic blowfly Lucilia sericata (green bottle blowfly) field strain NZ_LucSer_NP was generated using next-generation sequencing technology. The length of complete the mt genome is 15,938 bp, with 39.4% A, 13.0% C, 9.3% G, and 38.2% T nucleotide distribution. The complete mt genome consists of 13 protein-coding genes (PCGs), two ribosomal RNAs, 22 transfer RNAs, and a 1124 bp non-coding region, similar to most metazoan mt genomes. Phylogenetic analysis showed that L. sericata NZ_LucSer_NP forms a monophyletic cluster with the remaining six Lucilia species and the Calliphoridae are polyphyletic. This study provides the first complete mt genome sequence for a L. sericata blowfly species derived from New Zealand to facilitate species identification and phylogenetic analysis.

Characterization of the Complete Mitochondrial Genomes of Two Sibling Species of Parasitic Roundworms, Haemonchus contortus and Teladorsagia circumcincta

Haemonchus contortus and Teladorsagia circumcincta are among the two most pathogenic internal parasitic nematodes infecting small ruminants, such as sheep and goats, and are a global animal health issue. Accurate identification and delineation of Haemonchidae species is essential for development of diagnostic and control strategies with high resolution for Trichostrongyloidea infection in ruminants. Here, we describe in detail and compare the complete mitochondrial (mt) genomes of the New Zealand H. contortus and T. circumcincta field strains to improve our understanding of species- and strain-level evolution in these closely related roundworms. In the present study, we performed extensive comparative bioinformatics analyses on the recently sequenced complete mt genomes of the New Zealand H. contortus NZ_Hco_NP and T. circumcincta NZ_Teci_NP field strains. Amino acid sequences inferred from individual genes of each of the two mt genomes were compared, concatenated and subjected to phylogenetic analysis using Bayesian inference (BI), Maximum Likelihood (ML), and Maximum Parsimony (MP). The AT-rich mt genomes of H. contortus NZ_Hco_NP and T. circumcincta NZ_Teci_NP are 14,001 bp (A+T content of 77.4%) and 14,081 bp (A+T content of 77.3%) in size, respectively. All 36 of the typical nematode mt genes are transcribed in the forward direction in both species and comprise of 12 protein-encoding genes (PCGs), 2 ribosomal RNA (rrn) genes, and 22 transfer RNA (trn) genes. The secondary structures for the 22 trn genes and two rrn genes differ between H. contortus NZ_Hco_NP and T. circumcincta NZ_Teci_NP, however the gene arrangements of both are consistent with other Trichostrongylidea sequenced to date. Comparative analyses of the complete mitochondrial nucleotide sequences, PCGs, A+T rich and non-coding repeat regions of H. contortus NZ_Hco_NP and T. circumcincta NZ_Teci_NP further reinforces the high levels of diversity and gene flow observed among Trichostrongylidea, and supports their potential as ideal markers for strain-level identification from different hosts and geographical regions with high resolution for future studies. The complete mt genomes of H. contortus NZ_Hco_NP and T. circumcincta NZ_Teci_NP presented here provide useful novel markers for further studies of the meta-population connectivity and the genetic mechanisms driving evolution in nematode species.