Evaluation of genome skimming to detect and characterise human and livestock helminths.

The identification of gastrointestinal helminth infections of humans and livestock almost exclusively relies on the detection of eggs or larvae in faeces, followed by manual counting and morphological characterisation to differentiate species using microscopy-based techniques. However, molecular approaches based on the detection and quantification of parasite DNA are becoming more prevalent, increasing the sensitivity, specificity and throughput of diagnostic assays. High-throughput sequencing, from single PCR targets through to the analysis of whole genomes, offers significant promise towards providing information-rich data that may add value beyond traditional and conventional molecular approaches; however, thus far, its utility has not been fully explored to detect helminths in faecal samples. In this study, low-depth whole genome sequencing, i.e. genome skimming, has been applied to detect and characterise helminth diversity in a set of helminth-infected human and livestock faecal material. The strengths and limitations of this approach are evaluated using three methods to characterise and differentiate metagenomic sequencing data based on (i) mapping to whole mitochondrial genomes, (ii) whole genome assemblies, and (iii) a comprehensive internal transcribed spacer 2 (ITS2) database, together with validation using quantitative PCR (qPCR). Our analyses suggest that genome skimming can successfully identify most single and multi-species infections reported by qPCR and can provide sufficient coverage within some samples to resolve consensus mitochondrial genomes, thus facilitating phylogenetic analyses of selected genera, e.g. Ascaris spp. Key to this approach is both the availability and integrity of helminth reference genomes, some of which are currently contaminated with bacterial and host sequences. The success of genome skimming of faecal DNA is dependent on the availability of vouchered sequences of helminths spanning both taxonomic and geographic diversity, together with methods to detect or amplify minute quantities of parasite nucleic acids in mixed samples.

Can avian flyways reflect dispersal barriers of clinostomid parasites? First evidence from the mitogenome of Clinostomum complanatum.

Clinostomum complanatum (Rudolphi, 1814) is an economically important parasitic flatworm (Trematoda, Digenea), yet little is known on the population structure of these animals. We characterise a new mitochondrial genome for C. complanatum, derived from an Iranian specimen. The newly obtained sequence is used to position the species in the digenean tree of life. The first-ever intraspecific comparison at mitogenome scale within C. complanatum revealed a high degree of similarity to the previously sequenced mitogenome of a distant (Italian) population. Avian migratory routes mirror phylogenetic clustering, and hence we suggest that infection of a flying host enables genetic exchange between parasites across large geographic distances. Comparative mitogenomic work in Clinostomum spp. at both the intra- (C. complanatum) and interspecific (C. complanatum-C. sinensis) level further shows that usage of new and/or additional mitochondrial markers is preferred over single-gene methods for high-resolution diagnostics and population biology.

Combining host and vector data informs emergence and potential impact of an Usutu virus outbreak in UK wild birds.

Following the first detection in the United Kingdom of Usutu virus (USUV) in wild birds in 2020, we undertook a multidisciplinary investigation that combined screening host and vector populations with interrogation of national citizen science monitoring datasets to assess the potential for population impacts on avian hosts. Pathological findings from six USUV-positive wild passerines were non-specific, highlighting the need for molecular and immunohistochemical examinations to confirm infection. Mosquito surveillance at the index site identified USUV RNA in Culex pipiens s.l. following the outbreak. Although the Eurasian blackbird (Turdus merula) is most frequently impacted by USUV in Europe, national syndromic surveillance failed to detect any increase in occurrence of clinical signs consistent with USUV infection in this species. Furthermore, there was no increase in recoveries of dead blackbirds marked by the national ringing scheme. However, there was regional clustering of blackbird disease incident reports centred near the index site in 2020 and a contemporaneous marked reduction in the frequency with which blackbirds were recorded in gardens in this area, consistent with a hypothesis of disease-mediated population decline. Combining results from multidisciplinary schemes, as we have done, in real-time offers a model for the detection and impact assessment of future disease emergence events.

High-throughput sequencing of faeces provides evidence for dispersal of parasites and pathogens by migratory waterbirds.

Examination of faecal material has demonstrated how a broad range of organisms are distributed by bird movements. Such research has largely focused on dispersal of plant seeds by frugivores and of freshwater organisms by waterbirds. However, with few exceptions (e.g. avian influenza, Ebola virus), there is a dearth of evidence for transport of parasites and pathogens. High-throughput sequencing methods now provide a powerful means of addressing this knowledge gap by elucidating faecal contents in unprecedented detail. We collected faeces excreted by a range of migratory waterbirds in south-west Spain and pooled faecal DNA to create libraries reflective of feeding behavior. We created sets of libraries using high-throughput metagenomic and amplicon sequencing. For the latter we employed two sets of primers to broadly target the V4 region of the 18S rRNA gene (one set amplifying the region across all eukaryotes, the other excluding amplification of metazoans). Libraries revealed a wide diversity of eukaryotes, including parasites of the faecal producers themselves, parasites of food items, or those incidentally ingested. We also detected novel microbial eukaryotic taxa and found that parasite assemblage profiles were relatively distinct. Comparing the performance of the methods used supports their joint use for future studies of diversity and abundance. Because viable stages of many parasites are likely to be present in faeces, our results suggest significant levels of bird-mediated dispersal of parasites (both from avian and other hosts). Our methods revealed much hidden biodiversity, and allowed identification of the individuals who produced the faecal samples to species level, facilitating the study of interaction networks.

Evolutionary transitions in broad tapeworms (Cestoda: Diphyllobothriidea) revealed by mitogenome and nuclear ribosomal operon phylogenetics.

Broad tapeworms (Diphyllobothriidea) are parasites whose adults are capable of infecting a wide range of freshwater, marine and terrestrial tetrapods including humans. Previous works examining the evolution of habitat and host use in this group have been hampered by the lack of a well-resolved phylogeny. In order to produce a robust phylogenetic framework for diphyllobothriideans, we sequenced the complete mitochondrial genome of 13 representatives, carefully chosen to cover the major clades, and two outgroup species representing the Spathebothriidea and Haplobothriidea. In addition, complementary data from the nuclear ribosomal operon was sequenced for 10 representative taxa. Mitogenomes and ssrDNA and lsrDNA were used towards elucidating the phylogenetic framework for the Diphyllobothriidea. The Cephalochlamydidae is confirmed as the earliest diverging diphyllobothriidean lineage, and Solenophoridae and Diphyllobothriidae are sister groups. We infer a probable freshwater origin of the diphyllobothriideans. The ancestral condition for life cycle complexity could not be unambiguously resolved. However, we infer exclusive use of a three-host life cycle following the origin of the Solenophoridae + Diphyllobothriidae. Regarding definitive host use, although we infer reptiles as the most likely ancestral condition, this result should be revisited with a more densely sampled phylogeny in future studies. Freshwater habitat is used by the early diverging lineages within the Solenophoridae + Diphyllobothriidae clade. For the latter, habitat use shifts between freshwater and marine environments, and definitive host use includes marine and terrestrial mammals and birds. We use mitochondrial genomes to distinguish Schistocephalus species occurring in different species of sticklebacks and demonstrate conspecificity of Ligula cf. intestinalis specimens collected from two Fennoscandian ringed seal subspecies.

The first mitochondrial genomes of endosymbiotic rhabdocoels illustrate evolutionary relaxation of atp8 and genome plasticity in flatworms.

The first three mitochondrial (mt) genomes of endosymbiotic turbellarian flatworms are characterised for the rhabdocoels Graffilla buccinicola, Syndesmis echinorum and S. kurakaikina. Interspecific comparison of the three newly obtained sequences and the only previously characterised rhabdocoel, the free-living species Bothromesostoma personatum, reveals high mt genomic variability, including numerous rearrangements. The first intrageneric comparison within rhabdocoels shows that gene order is not fully conserved even between congeneric species. Atp8, until recently assumed absent in flatworms, was putatively annotated in two sequences. Selection pressure was tested in a phylogenetic framework and is shown to be significantly relaxed in this and another protein-coding gene: cox1. If present, atp8 appears highly derived in platyhelminths and its functionality needs to be addressed in future research. Our findings for the first time allude to a large degree of undiscovered (mt) genomic plasticity in rhabdocoels. It merits further attention whether this variation is correlated with a symbiotic lifestyle. Our results illustrate that this phenomenon is widespread in flatworms as a whole and not exclusive to the better-studied neodermatans.

Integrated Approach Reveals Role of Mitochondrial Germ-Line Mutation F18L in Respiratory Chain, Oxidative Alterations, Drug Sensitivity, and Patient Prognosis in Glioblastoma.

Glioblastoma is the most common and malignant primary brain tumour in adults, with a dismal prognosis. This is partly due to considerable inter- and intra-tumour heterogeneity. Changes in the cellular energy-producing mitochondrial respiratory chain complex (MRC) activities are a hallmark of glioblastoma relative to the normal brain, and associate with differential survival outcomes. Targeting MRC complexes with drugs can also facilitate anti-glioblastoma activity. Whether mutations in the mitochondrial DNA (mtDNA) that encode several components of the MRC contribute to these phenomena remains underexplored. We identified a germ-line mtDNA mutation (m. 14798T > C), enriched in glioblastoma relative to healthy controls, that causes an amino acid substitution F18L within the core mtDNA-encoded cytochrome b subunit of MRC complex III. F18L is predicted to alter corresponding complex III activity, and sensitivity to complex III-targeting drugs. This could in turn alter reactive oxygen species (ROS) production, cell behaviour and, consequently, patient outcomes. Here we show that, despite a heterogeneous mitochondrial background in adult glioblastoma patient biopsy-derived cell cultures, the F18L substitution associates with alterations in individual MRC complex activities, in particular a 75% increase in MRC complex II_III activity, and a 34% reduction in CoQ10, the natural substrate for MRC complex III, levels. Downstream characterisation of an F18L-carrier revealed an 87% increase in intra-cellular ROS, an altered cellular distribution of mitochondrial-specific ROS, and a 64% increased sensitivity to clomipramine, a repurposed MRC complex III-targeting drug. In patients, F18L-carriers that received the current standard of care treatment had a poorer prognosis than non-carriers (373 days vs. 415 days, respectively). Single germ-line mitochondrial mutations could predispose individuals to differential prognoses, and sensitivity to mitochondrial targeted drugs. Thus, F18L, which is present in blood could serve as a useful non-invasive biomarker for the stratification of patients into prognostically relevant groups, one of which requires a lower dose of clomipramine to achieve clinical effect, thus minimising side-effects.

The first next-generation sequencing approach to the mitochondrial phylogeny of African monogenean parasites (Platyhelminthes: Gyrodactylidae and Dactylogyridae).

BACKGROUND: Monogenean flatworms are the main ectoparasites of fishes. Representatives of the species-rich families Gyrodactylidae and Dactylogyridae, especially those infecting cichlid fishes and clariid catfishes, are important parasites in African aquaculture, even more so due to the massive anthropogenic translocation of their hosts worldwide. Several questions on their evolution, such as the phylogenetic position of Macrogyrodactylus and the highly speciose Gyrodactylus, remain unresolved with available molecular markers. Also, diagnostics and population-level research would benefit from the development of higher-resolution genetic markers. We aim to offer genetic resources for work on African monogeneans by providing mitogenomic data of four species (two belonging to Gyrodactylidae, two to Dactylogyridae), and analysing their gene sequences and gene order from a phylogenetic perspective. RESULTS: Using Illumina technology, the first four mitochondrial genomes of African monogeneans were assembled and annotated for the cichlid parasites Gyrodactylus nyanzae, Cichlidogyrus halli, Cichlidogyrus mbirizei (near-complete mitogenome) and the catfish parasite Macrogyrodactylus karibae (near-complete mitogenome). Complete nuclear ribosomal operons were also retrieved, as molecular vouchers. The start codon TTG is new for Gyrodactylus and for Dactylogyridae, as is the incomplete stop codon TA for Dactylogyridae. Especially the nad2 gene is promising for primer development. Gene order was identical for protein-coding genes and differed between the African representatives of these families only in a tRNA gene transposition. A mitochondrial phylogeny based on an alignment of nearly 12,500 bp including 12 protein-coding and two ribosomal RNA genes confirms that the Neotropical oviparous Aglaiogyrodactylus forficulatus takes a sister group position with respect to the other gyrodactylids, instead of the supposedly 'primitive' African Macrogyrodactylus. Inclusion of the African Gyrodactylus nyanzae confirms the paraphyly of Gyrodactylus. The position of the African dactylogyrid Cichlidogyrus is unresolved, although gene order suggests it is closely related to marine ancyrocephalines. CONCLUSIONS: The amount of mitogenomic data available for gyrodactylids and dactylogyrids is increased by roughly one-third. Our study underscores the potential of mitochondrial genes and gene order in flatworm phylogenetics, and of next-generation sequencing for marker development for these non-model helminths for which few primers are available.

Aedes nigrinus (Eckstein, 1918) (Diptera, Culicidae), a new country record for England, contrasted with Aedes sticticus (Meigen, 1838).

We report the discovery of Aedes nigrinus (Eckstein, 1918) in the New Forest of southern England, bringing to 36 the number of mosquito species recorded in Britain. Because it seems that this species has been misidentified previously in Britain as the morphologically similar Aedes sticticus (Meigen, 1838), the two species are contrasted and distinguished based on distinctive differences exhibited in the adult and larval stages. The pupa of Ae. nigrinus is unknown, but the pupa of Ae. sticticus is distinguished from the pupae of other species of Aedes by modification of the most recent key to British mosquitoes. The history of the mosquito fauna recorded in the UK is summarized and bionomical information is provided for the two species.

The complete mitochondrial genome of Dixella aestivalis (Diptera: Nematocera: Dixidae).

Dixidae, meniscus midges, belong to the suborder Nematocera of the order Diptera. The family includes 197 known species classified in nine genera. The complete mitochondrial genome of the Dixella aestivalis (Meigen) from the United Kingdom is reported here, along with its annotation and comparison with the genome of an unidentified species of Dixella from China. The circular genome consists of 16 465 bp and has a gene content consisting of 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes and a non-coding, A + T-rich, control region. The mitochondrial genome of D. aestivalis can be used to identify genetic markers for species identification, and will be valuable for resolving phylogenetic relationships within the genus, family Dixidae and suborder Nematocera.

Assessing myxozoan presence and diversity using environmental DNA.

Amplicon sequencing on a High Throughput Sequencing platform (custom barcoding) was used to detect and characterise myxosporean communities in environmental DNA samples from marine and freshwater environments and in faeces of animals that may serve as hosts or whose prey may host myxosporean infections. A diversity of myxozoans in filtered water samples and in faeces of piscivores (otters and great cormorants) was detected, demonstrating the suitability of lineage-specific amplicons for characterising otherwise difficult to sample parasite communities. The importance of using this approach was highlighted by the lack of myxosporean detection using commonly employed, broadly targeted eukaryotic primers. These results suggest that, despite being frequently present in eDNA samples, myxozoans have been generally overlooked in "eukaryote-wide" surveys. Lineage-specific primers, in contrast, detected 107 operational taxonomic units that were assigned to both the "freshwater" and "marine" myxosporean lineages. Only 7% of these OTUs clustered with sequences in GenBank, providing evidence for substantial undescribed myxosporean diversity. Many new operational taxonomic units, including those found in otter faeces, clustered with a clade of myxosporeans previously characterised by sequences from invertebrate hosts and water samples only. Because myxozoan species identification is heavily reliant on molecular signatures, lineage-specific amplicon sequencing offers an effective and non-destructive means of improving our knowledge of myxozoan diversity. In addition, the analysis of myxozoan DNA in faeces of piscivores offers a potentially efficient method of sampling for diversity and revealing life cycles as piscivore activities may integrate myxozoan infections in fish over relatively broad spatial scales.

Next-Generation Mitogenomics: A Comparison of Approaches Applied to Caecilian Amphibian Phylogeny.

Mitochondrial genome (mitogenome) sequences are being generated with increasing speed due to the advances of next-generation sequencing (NGS) technology and associated analytical tools. However, detailed comparisons to explore the utility of alternative NGS approaches applied to the same taxa have not been undertaken. We compared a 'traditional' Sanger sequencing method with two NGS approaches (shotgun sequencing and non-indexed, multiplex amplicon sequencing) on four different sequencing platforms (Illumina's HiSeq and MiSeq, Roche's 454 GS FLX, and Life Technologies' Ion Torrent) to produce seven (near-) complete mitogenomes from six species that form a small radiation of caecilian amphibians from the Seychelles. The fastest, most accurate method of obtaining mitogenome sequences that we tested was direct sequencing of genomic DNA (shotgun sequencing) using the MiSeq platform. Bayesian inference and maximum likelihood analyses using seven different partitioning strategies were unable to resolve compellingly all phylogenetic relationships among the Seychelles caecilian species, indicating the need for additional data in this case.

Entomological aspects and the role of human behaviour in malaria transmission in a highland region of the Republic of Yemen.

BACKGROUND: The Republic of Yemen has the highest incidence of malaria in the Arabian Peninsula, yet little is known of its vectors or transmission dynamics. METHODS: A 24-month study of the vectors and related epidemiological aspects of malaria transmission was conducted in two villages in the Taiz region in 2004-2005. RESULTS: Cross-sectional blood film surveys recorded an overall malaria infection rate of 15.3 % (250/1638), with highest rates exceeding 30 % in one village in May and December 2005. With one exception, Plasmodium malariae, all infections were P. falciparum. Seven Anopheles species were identified among 3407 anophelines collected indoors using light traps (LT) and pyrethrum knockdown catches (PKD): Anopheles arabiensis (86.9 %), An. sergentii (9 %), An. azaniae, An. dthali, An. pretoriensis, An. coustani and An. algeriensis. Sequences for the standard barcode region of the mitochondrial COI gene confirmed the presence of two morphological forms of An. azaniae, the typical form and a previously unrecognized form not immediately identifiable as An. azaniae. ELISA detected Plasmodium sporozoites in 0.9 % of 2921 An. arabiensis (23 P. falciparum, two P. vivax) confirming this species as the primary malaria vector in Yemen. Plasmodium falciparum sporozoites were detected in An. sergentii (2/295) and a single female of An. algeriensis, incriminating both species as malaria vectors for the first time in Yemen. A vector in both wet and dry seasons, An. arabiensis was predominantly anthropophilic (human blood index = 0.86) with an entomological inoculation rate of 1.58 infective bites/person/year. Anopheles sergentii fed on cattle (67.3 %) and humans (48.3; 20.7 % mixed both species), but only 14.7 % were found in PKDs, indicating predominantly exophilic behaviour. A GIS analysis of geographic and socio-economic parameters revealed that An. arabiensis were significantly higher (P < 0.001) in houses with televisions, most likely due to the popular evening habit of viewing television collectively in houses with open doors and windows. CONCLUSIONS: The predominantly indoor human biting vectors recorded in this study could be targeted effectively with LLINs, indoor residual spraying and/or insecticide-treated window/door curtains reinforced by education to instil a perception that effective and affordable malaria prevention is achievable.

The mitochondrial genome and ribosomal operon of Brachycladium goliath (Digenea: Brachycladiidae) recovered from a stranded minke whale.

Members of the Brachycladiidae are known to cause pathologies implicated in cetacean strandings and it is important to develop accurate diagnostic markers to differentiate these and other helminths found in cetaceans. Brachycladium goliath (van Beneden, 1858) is a large trematode found, as adults, usually in the hepatic (bile) and pancreatic ducts of various cetaceans. Complete sequences were determined for the entire mitochondrial genome, and phylogenetically informative nuclear genes contained within the ribosomal operon, from a small piece of an individual worm taken from a common minke whale Balaenoptera acutorostrata Lacépède, 1804. Genomic DNA was sequenced using an Illumina MiSeq platform. The mtDNA is 15,229 bp in length consisting of 12 protein-coding genes, 22 tRNA genes, 2 rRNA genes and 2 non-coding regions of which the larger is comprised of 4 tandemly repeated units (260 bp each). The ribosomal RNA operon is 9297 bp long. These data provide a rich resource of molecular markers for diagnostics, phylogenetics and population genetics in order to better understand the role, and associated pathology of helminth infections in cetaceans.

Complete mitochondrial genome of the geniculate calcified red alga, Corallina officinalis (Corallinales, Rhodophyta).

We present the first mitochondrial genome of the calcified, geniculate coralline red alga Corallina officinalis (Corallinales). The circular genome consists of 26,504 bp and has a gene content consisting of 23 protein-coding genes, 26 transfer RNA genes and two ribosomal RNA genes, with an overall GC content of 30.1%.

The mitochondrial genome of Parascaris univalens--implications for a "forgotten" parasite.

BACKGROUND: Parascaris univalens is an ascaridoid nematode of equids. Little is known about its epidemiology and population genetics in domestic and wild horse populations. PCR-based methods are suited to support studies in these areas, provided that reliable genetic markers are used. Recent studies have shown that mitochondrial (mt) genomic markers are applicable in such methods, but no such markers have been defined for P. univalens. METHODS: Mt genome regions were amplified from total genomic DNA isolated from P. univalens eggs by long-PCR and sequenced using Illumina technology. The mt genome was assembled and annotated using an established bioinformatic pipeline. Amino acid sequences inferred from all protein-encoding genes of the mt genomes were compared with those from other ascaridoid nematodes, and concatenated sequences were subjected to phylogenetic analysis by Bayesian inference. RESULTS: The circular mt genome was 13,920 bp in length and contained two ribosomal RNA, 12 protein-coding and 22 transfer RNA genes, consistent with those of other ascaridoids. Phylogenetic analysis of the concatenated amino acid sequence data for the 12 mt proteins showed that P. univalens was most closely related to Ascaris lumbricoides and A. suum, to the exclusion of other ascaridoids. CONCLUSIONS: This mt genome representing P. univalens now provides a rich source of genetic markers for future studies of the genetics and epidemiology of this parasite and its congener, P. equorum. This focus is significant, given that there is no published information on the specific prevalence and distribution of P. univalens infection in domestic and wild horse populations.

Can long-range PCR be used to amplify genetically divergent mitochondrial genomes for comparative phylogenetics? A case study within spiders (Arthropoda: Araneae).

The development of second generation sequencing technology has resulted in the rapid production of large volumes of sequence data for relatively little cost, thereby substantially increasing the quantity of data available for phylogenetic studies. Despite these technological advances, assembling longer sequences, such as that of entire mitochondrial genomes, has not been straightforward. Existing studies have been limited to using only incomplete or nominally intra-specific datasets resulting in a bottleneck between mitogenome amplification and downstream high-throughput sequencing. Here we assess the effectiveness of a wide range of targeted long-range PCR strategies, encapsulating single and dual fragment primer design approaches to provide full mitogenomic coverage within the Araneae (Spiders). Despite extensive rounds of optimisation, full mitochondrial genome PCR amplifications were stochastic in most taxa, although 454 Roche sequencing confirmed the successful amplification of 10 mitochondrial genomes out of the 33 trialled species. The low success rates of amplification using long-Range PCR highlights the difficulties in consistently obtaining genomic amplifications using currently available DNA polymerases optimised for large genomic amplifications and suggests that there may be opportunities for the use of alternative amplification methods.