Synthesizing environmental, epidemiological and vector and parasite genetic data to assist decision making for disease elimination.

We present a framework for identifying when conditions are favourable for transmission of vector-borne diseases between communities by incorporating predicted disease prevalence mapping with landscape analysis of sociological, environmental and host/parasite genetic data. We explored the relationship between environmental features and gene flow of a filarial parasite of humans, Onchocerca volvulus, and its vector, blackflies in the genus Simulium. We generated a baseline microfilarial prevalence map from point estimates from 47 locations in the ecological transition separating the savannah and forest in Ghana, where transmission of O. volvulus persists despite onchocerciasis control efforts. We generated movement suitability maps based on environmental correlates with mitochondrial population structure of 164 parasites from 15 communities and 93 vectors from only four sampling sites, and compared these to the baseline prevalence map. Parasite genetic distance between sampling locations was significantly associated with elevation (r = .793, p = .005) and soil moisture (r = .507, p = .002), while vector genetic distance was associated with soil moisture (r = .788, p = .0417) and precipitation (r = .835, p = .0417). The correlation between baseline prevalence and parasite resistance surface maps was stronger than that between prevalence and vector resistance surface maps. The centre of the study area had high prevalence and suitability for parasite and vector gene flow, potentially contributing to persistent transmission and suggesting the importance of re-evaluating transmission zone boundaries. With suitably dense sampling, this framework can help delineate transmission zones for onchocerciasis and would be translatable to other vector-borne diseases.

Characterisation of seryl tRNA synthetase (srs-2) in Haemonchus contortus and Teladorsagia circumcincta.

The aim of the study was to purify and characterise recombinant proteins with the potential as an anti-parasite vaccine. Full-length cDNAs encoding seryl-tRNA synthetase (srs-2) were cloned from Haemonchus contortus (HcSRS-2) and Teladorsagia circumcincta (TcSRS-2). TcSRS-2 and HcSRS-2 cDNA (1458bp) encoded proteins of 486 amino acids, each of which was present as a single band of about 55 kDa on SDS-PAGE. Multiple alignments of the protein sequences showed homology of 94% between TcSRS-2 and HcSRS-2, 76-93% with SRS-2s of eight nematodes and 68% with Mus musculus SRS-2. The predicted three-dimensional structures revealed an overall structural homology of TcSRS-2 and HcSRS-2, highly conserved binding and catalytic sites, and minor differences in the tautomerase binding site residues in other nematode SRS-2 homologues. A phylogenetic tree was constructed using helminth and mammalian SRS-2 sequences. Soluble C-terminal SRS-2 proteins were expressed in Escherichia coli strain AY2.4 and purified. Recombinant HcSRS-2 assay shows that the recombinant enzyme was active and stable. The Km and Vmax for ATP were 3.9 ± 1.0 µM and 2.7 ± 0.1 µmol min-1 mg-1 protein, respectively. Antibodies in serum and saliva from field-immune, but not nematode-naïve, sheep recognised recombinant HcSRS-2 and TcSRS-2 in enzyme-linked immunosorbent assays. Recognition of the recombinant proteins by antibodies generated by exposure of sheep to the native enzyme indicates similar antigenicity of the two proteins.

Chemosensory mechanisms of host seeking and infectivity in skin-penetrating nematodes.

Approximately 800 million people worldwide are infected with one or more species of skin-penetrating nematodes. These parasites persist in the environment as developmentally arrested third-stage infective larvae (iL3s) that navigate toward host-emitted cues, contact host skin, and penetrate the skin. iL3s then reinitiate development inside the host in response to sensory cues, a process called activation. Here, we investigate how chemosensation drives host seeking and activation in skin-penetrating nematodes. We show that the olfactory preferences of iL3s are categorically different from those of free-living adults, which may restrict host seeking to iL3s. The human-parasitic threadworm Strongyloides stercoralis and hookworm Ancylostoma ceylanicum have highly dissimilar olfactory preferences, suggesting that these two species may use distinct strategies to target humans. CRISPR/Cas9-mediated mutagenesis of the S. stercoralis tax-4 gene abolishes iL3 attraction to a host-emitted odorant and prevents activation. Our results suggest an important role for chemosensation in iL3 host seeking and infectivity and provide insight into the molecular mechanisms that underlie these processes.

Genomic introgression mapping of field-derived multiple-anthelmintic resistance in Teladorsagia circumcincta.

Preventive chemotherapy has long been practiced against nematode parasites of livestock, leading to widespread drug resistance, and is increasingly being adopted for eradication of human parasitic nematodes even though it is similarly likely to lead to drug resistance. Given that the genetic architecture of resistance is poorly understood for any nematode, we have analyzed multidrug resistant Teladorsagia circumcincta, a major parasite of sheep, as a model for analysis of resistance selection. We introgressed a field-derived multiresistant genotype into a partially inbred susceptible genetic background (through repeated backcrossing and drug selection) and performed genome-wide scans in the backcross progeny and drug-selected F2 populations to identify the major genes responsible for the multidrug resistance. We identified variation linking candidate resistance genes to each drug class. Putative mechanisms included target site polymorphism, changes in likely regulatory regions and copy number variation in efflux transporters. This work elucidates the genetic architecture of multiple anthelmintic resistance in a parasitic nematode for the first time and establishes a framework for future studies of anthelmintic resistance in nematode parasites of humans.

Dominant negative mutations of Caenorhabditis elegans daf-7 confer a novel developmental phenotype.

BACKGROUND: TGF-ß signaling pathways are involved in the control of development in every member of the animal kingdom. As such, TGF-ß ligands are widely divergent yet retain a set of core conserved features, specifically, a pre-protein cleavage site and several conserved ligand domain residues, the disruption of which produces a dominant negative phenotype. RESULTS: We have extended these observations into an invertebrate system by creating a series of loss-of-function Caenorhabditis elegans daf-7 transgenes. When we tested these mutant transgenes in a daf-7/+ background, we saw a molting and excretory canal phenotype. Members of both pathways downstream of daf-4 were required for this phenotype. CONCLUSIONS: Our results show that the basic mechanisms of TGF-ß function are conserved across the animal kingdom. A subset of our daf-7 mutations also produced an unexpected and novel phenotype. Epistasis experiments demonstrated that both daf-3/-5 and sma-4/-9 were downstream of our mutant daf-7 transgenes, which suggests not only a role for DAF-7 in the control of molting and the development of the excretory system but also that daf-7 and dbl-1 signaling may converge downstream of their shared Type II receptor, daf-4. Our approach may unveil new roles in development for other invertebrate TGF-ß ligands.

Levels of Amyloid Beta (Aß) Expression in the Caenorhabditis elegans Neurons Influence the Onset and Severity of Neuronally Mediated Phenotypes.

A characteristic feature of Alzheimer's disease (AD) is the formation of neuronal extracellular senile plaques composed of aggregates of fibrillar amyloid ß (Aß) peptides, with the Aß1-42 peptide being the most abundant species. These Aß peptides have been proposed to contribute to the pathophysiology of the disease; however, there are few tools available to test this hypothesis directly. In particular, there are no data that establish a dose-response relationship between Aß peptide expression level and disease. We have generated a panel of transgenic Caenorhabditis elegans strains expressing the human Aß1-42 peptide under the control of promoter regions of two pan-neuronal expressed genes, snb-1 and rgef-1. Phenotypic data show strong age-related defects in motility, subtle changes in chemotaxis, reduced median and maximum lifespan, changes in health span indicators, and impaired learning. The Aß1-42 expression level of these strains differed as a function of promoter identity and transgene copy number, and the timing and severity of phenotypes mediated by Aß1-42 were strongly positively correlated with expression level. The pan-neuronal expression of varying levels of human Aß1-42 in a nematode model provides a new tool to investigate the in vivo toxicity of neuronal Aß expression and the molecular and cellular mechanisms underlying AD progression in the absence of endogenous Aß peptides. More importantly, it allows direct quantitative testing of the dose-response relationship between neuronal Aß peptide expression and disease for the first time. These strains may also be used to develop screens for novel therapeutics to treat Alzheimer's disease.

Target screening using RNA interference in the sheep abomasal nematode parasite Haemonchus contortus.

RNA interference (RNAi) on parasitic nematodes has been described as a valuable tool for screening putative targets that could be used as novel drug and/or vaccine candidates. This study aimed to set up a pipeline to identify potential targets using RNAi for vaccine/anti-parasite therapy development against Haemonchus contortus, a blood-feeding abomasal nematode parasite. The available H. contortus sequence data was mined for targets, which were tested for essentiality using RNAi electroporation assays. A total of 56 genes were identified and tested for knockdown using electroporation of first-stage larvae (L1) H. contortus with the target double-stranded RNA. Electroporation of L1 proved to be effective overall; 17 targets had a strong phenotype and significant reduction in alive H. contortus, and another 24 had a moderate phenotype with a significant reduction in larvae development. A total of 28 targets showed a significant reduction in the development of H. contortus larvae to the infective stage (L3) following the RNAi assay. Down-regulation of target transcript levels was evaluated in some targets by semi-quantitative PCR. Four out of five genes tested showed complete knockdown of mRNA levels via semi-quantitative PCR, whereas the knockdown was partial for one. In conclusion, the results indicate that the RNAi pathway is confirmed in H. contortus and that several target genes have the potential to be investigated further as possible vaccine candidates.

Genome-based tools for onchocerciasis elimination: utility of the mitochondrial genome for delineating Onchocerca volvulus transmission zones.

National programs in Africa have expanded their objectives from control of onchocerciasis (river blindness) as a public health problem to elimination of parasite transmission, motivated by the reduction of Onchocerca volvulus infection prevalence in many African meso- and hyperendemic areas due to mass drug administration of ivermectin (MDAi). Given the large, contiguous hypo-, meso-, and hyperendemic areas, sustainable elimination of onchocerciasis in sub-Saharan Africa requires delineation of geographic boundaries for parasite transmission zones, so that programs can consider the risk of parasite re-introduction through vector or human migration from areas with ongoing transmission when making decisions to stop MDAi. We propose that transmission zone boundaries can be delineated by characterising the parasite genetic population structure within and between potential zones. We analysed whole mitochondrial genome sequences of 189 O. volvulus adults to determine the pattern of genetic similarity across three West African countries: Ghana, Mali, and Côte d'Ivoire. Population genetic structure indicates that parasites from villages near the Pru, Daka, and Black Volta rivers in central Ghana belong to one parasite population, indicating that the assumption that river basins constitute individual transmission zones is not supported by the data. Parasites from Mali and Côte d'Ivoire are genetically distinct from those from Ghana. This research provides the basis for developing tools for elimination programs to delineate transmission zones, to estimate the risk of parasite re-introduction via vector or human movement when intervention is stopped in one area while transmission is ongoing in others, to identify the origin of infections detected post-treatment cessation, and to investigate whether persisting prevalence despite ongoing interventions in one area is due to parasites imported from others.

Cytotaxonomic characterization and estimation of migration patterns of onchocerciasis vectors (Simulium damnosum sensu lato) in northwestern Ethiopia based on RADSeq data.

BACKGROUND: While much progress has been made in the control and elimination of onchocerciasis across Africa, the extent to which vector migration might confound progress towards elimination or result in re-establishment of endemism in areas where transmission has been eliminated remains unclear. In Northern Ethiopia, Metema and Metekel-two foci located near the Sudan border-exhibit continuing transmission. While progress towards elimination has been faster in Metema, there remains a problematic hotspot of transmission. Whether migration from Metekel contributes to this is currently unknown. METHODOLOGY/PRINCIPLE FINDINGS: To assess the role of vector migration from Metekel into Metema, we present a population genomics study of 151 adult female vectors using 47,638 RADseq markers and mtDNA CoI sequencing. From additional cytotaxonomy data we identified a new cytoform in Metema, closely related to S. damnosum s.str, here called the Gondar form. RADseq data strongly indicate the existence of two distinctly differentiated clusters within S. damnosum s.l.: one genotypic cluster found only in Metema, and the second found predominantly in Metekel. Because blackflies from both clusters were found in sympatry (in all four collection sites in Metema), but hybrid genotypes were not detected, there may be reproductive barriers preventing interbreeding. The dominant genotype in Metema was not found in Metekel while the dominant genotype in Metekel was found in Metema, indicating that (at the time of sampling) migration is primarily unidirectional, with flies moving from Metekel to Metema. There was strong differentiation between clusters but little genetic differentiation within clusters, suggesting migration and gene flow of flies within the same genetic cluster are sufficient to prevent genetic divergence between sites. CONCLUSIONS/SIGNIFICANCE: Our results confirm that Metekel and Metema represent different transmission foci, but also indicate a northward movement of vectors between foci that may have epidemiological importance, although its significance requires further study.

Strongyloides questions-a research agenda for the future.

The Strongyloides genus of parasitic nematodes have a fascinating life cycle and biology, but are also important pathogens of people and a World Health Organization-defined neglected tropical disease. Here, a community of Strongyloides researchers have posed thirteen major questions about Strongyloides biology and infection that sets a Strongyloides research agenda for the future. This article is part of the Theo Murphy meeting issue 'Strongyloides: omics to worm-free populations'.

Parastrongyloides trichosuri suggests that XX/XO sex determination is ancestral in Strongyloididae (Nematoda).

The parasitic roundworms Strongyloides stercoralis (in man) and Strongyloides ratti (in rats) employ environmentally controlled XX/XO sex determination with a pair of X chromosomes and two pairs of autosomes. Strongyloides papillosus (in sheep) has only two pairs of chromosomes, one of which combines the genetic material homologous to the S. ratti chromosomes X and I. This species creates males through the elimination of one copy of the portion related to the X chromosome (chromatin diminution). It is not clear which one of these two sex-determining mechanisms is ancestral. We demonstrate that Strongyloides vituli (in cattle) has two pairs of chromosomes like its very close relative S. papillosus whereas Parastrongyloides trichosuri, a closely related out-group to Strongyloides spp. in Australian brushtail possums, has three chromosome pairs and employs XX/XO sex determination. The X chromosome of P. trichosuri is homologous to the X chromosome of S. ratti. Our data strongly suggest that the last common ancestor of Strongyloides spp. and Parastrongyloides spp. had two pairs of autosomes along with two or one X chromosome in females and males, respectively. The situation with two pairs of chromosomes is likely derived and occurred through the fusion of the X chromosome with an autosome.

Procedure for Handling and Storage of Onchocerca volvulus Microfilariae Obtained from Skin Snips for Downstream Genetic Work.

WHO and endemic countries target elimination of transmission of Onchocerca volvulus, the parasite causing onchocerciasis. Population genetic analysis of O. volvulus may provide data to improve the evidence base for decisions on when, where, and for how long to deploy which interventions and post-intervention surveillance to achieve elimination. Development of necessary methods and tools requires parasites suitable for genetic analysis. Based on our experience with microfilariae obtained from different collaborators, we developed a microfilariae transfer procedure for large-scale studies in the Democratic Republic of Congo (DRC) comparing safety and efficacy of ivermectin, the mainstay of current onchocerciasis elimination strategies, and moxidectin, a new drug. This procedure is designed to increase the percentage of microfilariae in skin snips suitable for genetic analysis, improve assignment to metadata, and minimize time and materials needed by the researchers collecting the microfilariae. Among 664 microfilariae from South Sudan, 35.7% and 39.5% failed the mitochondrial and nuclear qPCR assay. Among the 576 microfilariae from DRC, 16.0% and 16.7% failed these assays, respectively. This difference may not only be related to the microfilariae transfer procedure but also to other factors, notably the ethanol concentration in the tubes in which microfilariae were stored (64% vs. ≥75%).

Assessing Onchocerca volvulus Intensity of Infection and Genetic Diversity Using Mitochondrial Genome Sequencing of Single Microfilariae Obtained before and after Ivermectin Treatment.

Onchocerciasis is a neglected tropical disease targeted for elimination using ivermectin mass administration. Ivermectin kills the microfilariae and temporarily arrests microfilariae production by the macrofilariae. We genotyped 436 microfilariae from 10 people each in Ituri, Democratic Republic of the Congo (DRC), and Maridi County, South Sudan, collected before and 4-5 months after ivermectin treatment. Population genetic analyses identified 52 and 103 mitochondrial DNA haplotypes among the microfilariae from DRC and South Sudan, respectively, with few haplotypes shared between people. The percentage of genotype-based correct assignment to person within DRC was ~88% and within South Sudan ~64%. Rarefaction and extrapolation analysis showed that the genetic diversity in DRC, and even more so in South Sudan, was captured incompletely. The results indicate that the per-person adult worm burden is likely higher in South Sudan than DRC. Analyses of haplotype data from a subsample (n = 4) did not discriminate genetically between pre- and post-treatment microfilariae, confirming that post-treatment microfilariae are not the result of new infections. With appropriate sampling, mitochondrial haplotype analysis could help monitor changes in the number of macrofilariae in a population as a result of treatment, identify cases of potential treatment failure, and detect new infections as an indicator of continuing transmission.

Geospatial modeling of pre-intervention nodule prevalence of Onchocerca volvulus in Ethiopia as an aid to onchocerciasis elimination.

BACKGROUND: Onchocerciasis is a neglected tropical filarial disease transmitted by the bites of blackflies, causing blindness and severe skin lesions. The change in focus for onchocerciasis management from control to elimination requires thorough mapping of pre-control endemicity to identify areas requiring interventions and to monitor progress. Onchocerca volvulus nodule prevalence in sub-Saharan Africa is spatially continuous and heterogeneous, and highly endemic areas may contribute to transmission in areas of low endemicity or vice-versa. Ethiopia is one such onchocerciasis-endemic country with heterogeneous O. volvulus nodule prevalence, and many districts are still unmapped despite their potential for onchocerciasis transmission. METHODOLOGY/PRINCIPLE FINDINGS: A Bayesian geostatistical model was fitted for retrospective pre-intervention nodule prevalence data collected from 916 unique sites and 35,077 people across Ethiopia. We used multiple environmental, socio-demographic, and climate variables to estimate the pre-intervention prevalence of O. volvulus nodules across Ethiopia and to explore their relationship with prevalence. Prevalence was high in southern and northwestern Ethiopia and low in Ethiopia's central and eastern parts. Distance to the nearest river (RR: 0.9850, 95% BCI: 0.9751-0.995), precipitation seasonality (RR: 0.9837, 95% BCI: 0.9681-0.9995), and flow accumulation (RR: 0.9586, 95% BCI: 0.9321-0.9816) were negatively associated with O. volvulus nodule prevalence, while soil moisture (RR: 1.0218, 95% BCI: 1.0135-1.0302) was positively associated. The model estimated the number of pre-intervention cases of O. volvulus nodules in Ethiopia to be around 6.48 million (95% BCI: 3.53-13.04 million). CONCLUSIONS/SIGNIFICANCE: Nodule prevalence distribution was correlated with habitat suitability for vector breeding and associated biting behavior. The modeled pre-intervention prevalence can be used as a guide for determining priorities for elimination mapping in regions of Ethiopia that are currently unmapped, most of which have comparatively low infection prevalence.

Enhanced annealing of mismatched oligonucleotides using a novel melting curve assay allows efficient in vitro discrimination and restriction of a single nucleotide polymorphism.

BACKGROUND: Many SNP discrimination strategies employ natural restriction endonucleases to discriminate between allelic states. However, SNPs are often not associated with a restriction site and therefore, a number of attempts have been made to generate sequence-adaptable restriction endonucleases. In this study, a simple, sequence-adaptable SNP discrimination mechanism between a 'wild-type' and 'mutant' template is demonstrated. This model differs from other artificial restriction endonuclease models as cis- rather than trans-orientated regions of single stranded DNA were generated and cleaved, and therefore, overcomes potential issues of either inefficient or non-specific binding when only a single variant is targeted. RESULTS: A series of mismatch 'bubbles' that spanned 0-5-bp surrounding a point mutation was generated and analysed for sensitivity to S1 nuclease. In this model, generation of oligonucleotide-mediated ssDNA mismatch 'bubbles' in the presence of S1 nuclease resulted in the selective degradation of the mutant template while maintaining wild-type template integrity. Increasing the size of the mismatch increased the rate of mutant sequence degradation, until a threshold above which discrimination was lost and the wild-type sequence was degraded. This level of fine discrimination was possible due to the development of a novel high-resolution melting curve assay to empirically determine changes in Tm (~5.0°C per base-pair mismatch) and to optimise annealing conditions (~18.38°C below Tm) of the mismatched oligonucleotide sets. CONCLUSIONS: The in vitro 'cleavage bubble' model presented is sequence-adaptable as determined by the binding oligonucleotide, and hence, has the potential to be tailored to discriminate between any two or more SNPs. Furthermore, the demonstrated fluorometric assay has broad application potential, offering a rapid, sensitive and high-throughput means to determine Tm and annealing rates as an alternative to conventional hybridisation detection strategies.

Intraspecific epitopic variation in a carbohydrate antigen exposed on the surface of Trichostrongylus colubriformis infective L3 larvae.

The carbohydrate larval antigen, CarLA, is present on the exposed surface of all strongylid nematode infective L3 larvae tested, and antibodies against CarLA can promote rapid immune rejection of incoming Trichostrongylus colubriformis larvae in sheep. A library of ovine recombinant single chain Fv (scFv) antibody fragments, displayed on phage, was prepared from B cell mRNA of field-immune sheep. Phage displaying scFvs that bind to the surface of living exsheathed T. colubriformis L3 larvae were identified, and the majority of worm-binding scFvs recognized CarLA. Characterization of greater than 500 worm surface binding phage resulted in the identification of nine different anti-CarLA scFvs that recognized three distinct T. colubriformis CarLA epitopes based on blocking and additive ELISA. All anti-CarLA scFvs were specific to the T. colubriformis species of nematode. Each of the three scFv epitope classes displayed identical Western blot recognition patterns and recognized the exposed surface of living T. colubriformis exsheathed L3 larvae. Surprisingly, each of the anti-CarLA scFvs was able to bind to only a subset of worms. Double-labelling indirect immunofluorescence revealed that the three classes of anti-CarLA scFvs recognize distinct, non-overlapping, T. colubriformis sub-populations. These results demonstrate that individual T. colubriformis L3 larvae display only one of at least three distinct antigenic forms of CarLA on their surface at any given time, and suggest that antigenic variation within CarLA is likely a mechanism of immune evasion in strongylid nematodes.

The genome sequence of the Australian filarial nematode, Cercopithifilaria johnstoni.

We present a genome assembly and annotation of an individual female Cercopithifilaria johnstoni, a parasitic filarial nematode that is transmitted by hard ticks (Ixodidae) to infect a broad range of native Australian murid and marsupial hosts. The genome sequence is 76.9 Mbp in length, and although in draft form (N50 = 99 kbp, N50[n] = 232), is largely complete based on universally conserved orthologs (BUSCOs; genome = 94.9%, protein = 96.5%) and relative to other related filarial species. These data represent the first genomic resources for the genus Cercopithifilaria, a group of parasites with a broad host range, and form the basis for comparative analysis with the human-infective parasite, Onchocerca volvulus, both of which are responsible for similar eye and skin pathologies in their respective hosts.

Spread of anthelmintic resistance in intestinal helminths of dogs and cats is currently less pronounced than in ruminants and horses - Yet it is of major concern.

Anthelmintic resistance (AR) has thus far only rarely been reported for intestinal helminths of dogs and cats, in contrast to parasites of livestock and horses. We highlight possible reasons for this striking and important discrepancy, including ecological, biological and genetic factors and/or intervention regimens of key intestinal helminths concerning both host groups. In view of the current knowledge related to the genetics, mechanisms and principles of AR development, we point at issues which in our view contribute to a comparatively lower risk of AR development in intestinal helminths of dogs and cats. Finally, we specify research needs and provide recommendations by which, based on the available information about AR in ruminant and equine helminths, the development of AR in dog and cat helminths may best be documented, prevented or at least postponed.

Genetic epidemiology of lymphatic filariasis in American Samoa after mass drug administration.

Over 892 million people in 48 countries are at risk of infection by nematodes that cause lymphatic filariasis. As part of the Global Programme to Eliminate Lymphatic Filariasis, mass drug administration is distributed to communities until surveillance indicates infection rates are below target prevalence thresholds. In some countries, including American Samoa, lymphatic filariasis transmission persists despite years of mass drug administration and/or has resurged after cessation. Nothing is known about the population genetics of Wuchereria bancrofti worms in Polynesia, or whether local transmission is persisting and/or increasing due to inadequate mass drug administration coverage, expansion from residual hotspots, reintroduction from elsewhere, or a combination. We extracted DNA from microfilariae on blood slides collected during prevalence surveys in 2014 and 2016, comprising 31 pools of five microfilariae from 22 persons living in eight villages. We sequenced 1104 bp across three mitochondrial markers (ND4, COI, CYTB). We quantified parasite genetic differentiation using variant calls and estimated haplotypes using principal components analysis, F-statistics, and haplotype networks. Of the variants called, all but eight were shared across the main island of Tutuila, and three of those were from a previously described hotspot village, Fagali'i. Genotypic data did not support population genetic structure among regions or villages in 2016, although differences were observed between worms collected in Fagali'i in 2014 and those from 2016. Because estimated haplotype frequency varied between villages, these statistics suggested genetic differentiation, but were not consistent among villages. Finally, haplotype networks demonstrated American Samoan sequence clusters were related to previously published sequences from Papua New Guinea. These are, to our knowledge, the first reports of W. bancrofti genetic variation in Polynesia. The resurgent parasites circulating on the main island of American Samoa represent a single population. This study is the first step towards investigating how parasite population structure might inform strategies to manage resurgence and elimination of lymphatic filariasis.

Cryptic species diversity in ticks that transmit disease in Australia.

Ticks are important vectors of a broad range of pathogens in Australia. Many tick species are morphologically similar and are therefore difficult to identify using morphology alone, particularly when collected in the larval and nymphal life stages. We report here the application of molecular methods to examine the species diversity of ixodid ticks at two sites in southern New South Wales, Australia. Our taxon sampling included six morphologically characterised adult stage voucher specimens of Ixodes trichosuri, Ixodes tasmani, Ixodes fecialis and Ixodes holocyclus (the paralysis tick) and ~250 field collected specimens that were in the larva or nymph stage and thus not morphologically identifiable. One nuclear and two mitochondrial amplicons were sequenced using a combination of Sanger and Illumina MiSeq sequencing. Phylogenetic relationships were estimated using both maximum likelihood and Bayesian methods. Two clades with strong bootstrap and Bayesian support were observed across trees estimated from each of three markers and from an analysis of the concatenated sequences. One voucher specimen of I. trichosuri was located in one of these clades, while the other I. trichosuri voucher specimen was in a second clade with the remaining three identified species, suggesting these morphologically similar ticks may represent different cryptic species. Unidentified specimens were found across both clades, and molecular divergence of many of these is equal to or greater than that observed between identified species, suggesting additional unidentified species may exist. Further studies are required to understand the taxonomic status of ticks in Australia, and how this species diversity impacts disease risk for livestock, domestic animals, wildlife and humans.