Characterization of 21 microsatellite loci for the precocious, grass-shrimp trematode Alloglossidium renale.

We developed microsatellite markers to use in studying the population genetics of the trematode Alloglossidium renale, a fluke with a precocious life cycle where sexual maturation occurs in a grass shrimp. Among 21 tested loci in a Mississippi population sample, 14 were polymorphic, 12 of which significantly deviated from Hardy-Weinberg Equilibrium (HWE). We estimated identity disequilibrium (ID) to confirm whether the deviations from HWE were due to significant amounts of selfing or due to technical factors. The selfing rate derived from FIS was 86.6%, whereas the selfing rate obtained by ID was 83.9%, indicating that the deviation in HWE was due to a high amount of selfing within the population. These markers will be useful for ecological and evolutionary studies of A. renale especially in relation to the interplay of hermaphroditic mating systems, inbreeding depression, and transmission dynamics.

TREMATODE CLONE ABUNDANCE DISTRIBUTIONS: AN ECO-EVOLUTIONARY LINK BETWEEN PARASITE TRANSMISSION AND PARASITE MATING SYSTEMS.

Most trematodes and some cestodes have obligate life history features that include an asexual developmental stage that can produce genetically-identical individuals (clonemates) followed by an adult stage with sexual reproduction. These life history features can influence the evolutionary mechanism of inbreeding in parasites, especially among self-compatible hermaphroditic endoparasites whose mating opportunities are restricted to within hosts. As clonemate mating in hermaphroditic species produces a genetic inbreeding signature identical to that of self-mating, it is important to understand how clonemates are transmitted through their life stages. A handful of prior studies compared clone richness (number of clones) across life cycle stages to infer transmission processes and to characterize clone abundance distributions (CADs) among hosts. Here we illustrate the use of the proportion of clonemate dyads (PC) within hosts to describe the CADs. PC has several advantages as an ecological metric in that it is unbiased by sample size, takes into account relative parasite burdens, and has a direct transmission interpretation, i.e., the probability of cotransmitting clonemates. Moreover, PC is also an evolutionary metric as it can be used to estimate a potential clonemate mating rate. We demonstrate the use of PC in comparing CADs within and across 2 trematode developmental stages in the lancet fluke Dicrocoelium dendriticum. Also, we show how genetic estimates of apparent selfing (true selfing plus clonemate mating) at larval and juvenile stages can be compared to PC estimated at the adult stage to assess the contribution of clonemate mating to apparent selfing. The eco-evolutionary links presented are generalizable to assess sibling cotransmission as well. Thus, the framework presented herein will facilitate future field-based studies on the transmission and mating systems of parasitic flatworms.

QUANTIFYING BILATERAL INFECTION PATTERNS IN THE TREMATODE ALLOGLOSSIDIUM RENALE.

Within-host distributions of parasites can have relevance to parasite competition, parasite mating, transmission, and host health. We examined the within-host distribution of the adult trematode Alloglossidium renale infecting the paired antennal glands of grass shrimp. There are 4 possible parasite distributions for infections of paired organs: random, uniform, biased aggregation to 1 particular organ (e.g., left vs. right), or inconsistently biased (aggregated, but does not favor 1 side). Previous work has shown that morphological asymmetries in hosts can lead to biased infections of paired organs. Apparent symmetry between the antennal glands of grass shrimp leads to the prediction that there would be no bias for 1 particular organ. However, an alternative prediction stems from the fact that A. renale is hermaphroditic: aggregation between glands would increase outcrossing opportunities and thus, avoid inbreeding via self-mating. Existing methods to test for an overall pattern did not apply to the A. renale system because of low-intensity infections as well as many 0 values for abundance per unit of the antennal gland. Hence, we used Monte Carlo simulations to determine if the observed overall patterns differed from those expected by randomly allocating parasites into groups of 2. We found that in 3 of 4 data sets, A. renale infections did not deviate from random distributions. The fourth data set had a more uniform pattern than expected by chance. As there was no aggregation between glands and the proportion of worms in single gland infections did not differ from that expected by chance alone, we found no evidence of inbreeding avoidance as might be manifested via a within-host distribution. Given the large proportion of worms in single infections, we predict as a major evolutionary outcome that populations of A. renale will be largely inbred.

Clonemate cotransmission supports a role for kin selection in a puppeteer parasite.

Host manipulation by parasites is a fascinating evolutionary outcome, but adaptive scenarios that often accompany even iconic examples in this popular field of study are speculative. Kin selection has been invoked as a means of explaining the evolution of an altruistic-based, host-manipulating behavior caused by larvae of the lancet fluke Dicrocoelium dendriticum in ants. Specifically, cotransmission of larval clonemates from a snail first host to an ant second host is presumed to lead to a puppeteer parasite in the ant's brain that has clonemates in the ant abdomen. Clonal relatedness between the actor (brain fluke) and recipients (abdomen flukes) enables kin selection of the parasite's host-manipulating trait, which facilitates transmission of the recipients to the final host. However, the hypothesis that asexual reproduction in the snail leads to a high abundance of clonemates in the same ant is untested. Clonal relationships between the manipulator in the brain and the nonmanipulators in the abdomen are also untested. We provide empirical data on the lancet fluke's clonal diversity within its ant host. In stark contrast to other trematodes, which do not exhibit the same host-manipulating behavioral trait, the lancet fluke has a high abundance of clonemates. Moreover, our data support existing theory that indicates that the altruistic behavior can evolve even in the presence of multiple clones within the same ant host. Importantly, our analyses conclusively show clonemate cotransmission into ants, and, as such, we find support for kin selection to drive the evolution and maintenance of this iconic host manipulation.

Recognizing the Causes of Parasite Morphological Variation to Resolve the Status of a Cryptogenic Pentastome.

Exotic species can threaten biodiversity by introducing parasites to native hosts. Thus, it is critical to identify if the same parasite species infects both native and exotic hosts. However, developmental- or environmental-induced morphological variation may render species identification ambiguous. Our study reports a range expansion in the southern United States of the pentastome Raillietiella indica from the Mediterranean gecko, Hemidactylus turcicus, as well as a host expansion into the green anole, Anolis carolinensis, in the anole's native range. Species identification was based on sequence data and male spicule shape. In agreement with a study from Australia, we found that much of the morphological variation in hook measurements, the primary diagnostic traits of raillietiellid pentastomes, was due to development. Here, we explicitly link this developmental variation to instar stage by incorporating experimental infection data obtained from the literature. We also show that the various hook traits are themselves highly correlated and, thus, likely not independent. Taking instar stage and correlated hook variables into account, we directly controlled for development on a composite hook size measurement. Using a large sample size from H. turcicus, we did not find any consistent effects of potential factors (host sex, host snout-vent-length, or parasite intensity) that may result in environmental-induced variation in relative hook size (corrected for body length). However, relative male spicule size tended to be negatively correlated with parasite intensity. In contrast, both pentastome body length and relative hook size significantly varied among host species whereas relative male spicule size was not significantly different among host species. Our study independently supports the conclusions that developmental- and host-induced morphological variations need to be accounted for to accurately identify pentastome species.

Little to no inbreeding depression in a tapeworm with mixed mating.

Meta-studies on hermaphrodites have found a negative relationship between primary selfing rates and levels of inbreeding depression (ID) and, thus, generally support purging in inbred systems. However, in plants, high among-taxa variance in ID results in no difference in the mean ID between outcrossing and mixed-mating taxa. Selective interference likely explains high ID among mixed-mating taxa, whereas low levels of ID among mixed-mating taxa are not as stressed. Among animal hermaphrodites, primarily molluscs, there are little data on mixed-mating systems. To fill a taxonomic and mating system gap, we tested for ID in a mixed-mating tapeworm, Oochoristica javaensis. We provide a direct estimate of ID across infection of an intermediate host by comparing selfing rates at two life history stages. We found little to no evidence for ID, and the level of ID falls in line with what is reported for highly selfing species even though O. javaensis has mixed mating. We discuss this result within the context of kin mating in O. javaensis. Our results emphasize that primary selfing rates alone may be insufficient to classify the inbreeding history in all species when testing for a relationship to ID. Mixed-mating taxa, and possibly some outcrossing taxa, may exhibit low levels of ID if biparental inbreeding is also driving purging. We advocate that ID studies report estimates of inbreeding history (e.g. FIS or identity disequilibrium) from nature-derived adult samples to provide context rather than relying on primary selfing rates alone.

Hurdles in the evolutionary epidemiology of Angiostrongylus cantonensis: Pseudogenes, incongruence between taxonomy and DNA sequence variants, and cryptic lineages.

Angiostrongylus cantonensis, the rat lungworm, is a zoonotic pathogen that is one of the leading causes of eosinophilic meningitis worldwide. This parasite is regarded as an emerging pathogen with a global range expansion out of southeastern Asia post-WWII. To date, molecular systematic/phylogeographic studies on A. cantonensis have mainly used two mitochondrial (mtDNA) markers, cytochrome c oxidase 1 (CO1) and cytochrome b (CYTB), where the focus has largely been descriptive in terms of reporting local patterns of haplotype variants. In order to look for more global evolutionary patterns, we herein provide a collective phylogenetic assessment using the six available whole mtDNA genome samples that have been tagged as A. cantonensis, A. malaysiensis, or A. mackerrasae along with all other GenBank CO1 and CYTB partial sequences that carry these species identifiers. The results reveal three important complications that researchers will need to be aware of, or will need to resolve, prior to conducting future molecular evolutionary studies on A. cantonensis. These three problems are (i) incongruence between taxonomic identifications and mtDNA variants (haplotypes or whole mtDNA genome samples), (ii) the presence of a CYTB mtDNA pseudogene, and (iii) the need to verify A. mackerrasae as a species along with other possible cryptic lineages, of which there is suggestive evidence (i.e., A. cantonensis could be a species complex). We provided a discussion of how these complications are hurdles to our understanding of the global epidemiology of angiostrongyliasis. We call for future studies to be more explicit in morphological traits used for identifications (e.g., provide measurements). Moreover, it will be necessary to repeat prior morphological and life-history studies while simultaneously using sequence data in order to assess possible associations between critical epidemiological data (e.g., biogeography, virulence/pathology, host species use) and specific lineages.

Resolving evolutionary changes in parasite life cycle complexity: Molecular phylogeny of the trematode genus Alloglossidium indicates more than one origin of precociousness.

The evolutionary causes and consequences of changes in complex life cycles are of central importance in parasitology. However, data remain limited because in part, knowledge on phylogenetic relationships among species that differ in life cycle patterns remains scarce. We present a molecular phylogeny of the trematode genus Alloglossidium, which contains several species that display precocious (a.k.a., progenetic) life cycles (i.e., maturation in what is typically regarded as an intermediate host). The molecular phylogeny contrasts with previous morphological and life-history based phylogenetic hypotheses. In particular, a precocious life cycle wherein leeches are used as final hosts evolved early in the history of the genus. Among the remaining species, which are a separate clade, a three-host life cycle using ictalurid catfishes is ancestral. Furthermore, there are at least two additional independent evolutionary events that lead to a precocious life cycle where a catfish host is lost and a crustacean is used as a final host. We conclude with a discussion on how existing hypotheses on the evolution of precociousness, and parasite life cycle complexity in general, may or may not relate to the patterns observed in genus Alloglossidium.

Role of parasite transmission in promoting inbreeding: I. Infection intensities drive individual parasite selfing rates.

Among parasitic organisms, inbreeding has been implicated as a potential driver of host-parasite co-evolution, drug-resistance evolution and parasite diversification. Yet, fundamental topics about how parasite life histories impact inbreeding remain to be addressed. In particular, there are no direct selfing-rate estimates for hermaphroditic parasites in nature. Our objectives were to elucidate the mating system of a parasitic flatworm in nature and to understand how aspects of parasite transmission could influence the selfing rates of individual parasites. If there is random mating within hosts, the selfing rates of individual parasites would be an inverse power function of their infection intensities. We tested whether selfing rates deviated from within-host random mating expectations with the tapeworm Oochoristica javaensis. In doing so, we generated, for the first time in nature, individual selfing-rate estimates of a hermaphroditic flatworm parasite. There was a mixed-mating system where tapeworms self-mated more than expected with random mating. Nevertheless, individual selfing rates still had a significant inverse power relationship to infection intensities. The significance of this finding is that the distribution of parasite infection intensities among hosts, an emergent property of the transmission process, can be a key driver in shaping the primary mating system, and hence the level of inbreeding in the parasite population. Moreover, we demonstrated how potential population selfing rates can be estimated using the predicted relationship of individual selfing rates to intensities and showed how the distribution of parasites among hosts can indirectly influence the primary mating system when there is density-dependent fecundity.

Role of parasite transmission in promoting inbreeding: II. Pedigree reconstruction reveals sib-transmission and consequent kin-mating.

Even though parasitic flatworms are one of the most species-rich groups of hermaphroditic organisms, we know virtually nothing of their mating systems (selfing or kin-mating rates) in nature. Hence, we lack an understanding of the role of inbreeding in parasite evolution. The natural mating systems of parasitic flatworms have remained elusive due to the inherent difficulty in generating progeny-array data in many parasite systems. New developments in pedigree reconstruction allow direct inference of realized selfing rates in nature by simply using a sample of genotyped individuals. We built upon this advancement by utilizing the closed mating systems, that is, individual hosts, of endoparasites. In particular, we created a novel means to use pedigree reconstruction data to estimate potential kin-mating rates. With data from natural populations of a tapeworm, we demonstrated how our newly developed methods can be used to test for cosibling transmission and inbreeding depression. We then showed how independent estimates of the two mating system components, selfing and kin-mating rates, account for the observed levels of inbreeding in the populations. Thus, our results suggest that these natural parasite populations are in inbreeding equilibrium. Pedigree reconstruction analyses along with the new companion methods we developed will be broadly applicable across a myriad of parasite species. As such, we foresee that a new frontier will emerge wherein the diverse life histories of flatworm parasites could be utilized in comparative evolutionary studies to broadly address ecological factors or life history traits that drive mating systems and hence inbreeding in natural populations.

Cryptic lineage diversity in the zoonotic pathogen Angiostrongylus cantonensis.

Delimitation of species is still a necessity among parasitic pathogens especially where morphological characters provide limited discernibility. Identification of cryptic lineages (independently evolving lineages that are morphologically similar) is critical as there could be lineage-specific traits that are of epidemiological importance. Angiostrongylus cantonensis is a zoonotic pathogen that can cause eosinophilic meningitis in humans. Recent reports of single marker sequence divergence hint at the potential for cryptic diversity in this lungworm. However, to definitively address if single marker divergence corresponds to independent evolving lineages, a multilocus approach is necessary. Using multilocus data, our goal was to determine if there were cryptic lineages within Thailand, a country plagued by several outbreaks and isolated cases of A. cantonensis infection. We analyzed the genetic structure of A. cantonensis samples collected from snails, Achatina fulica, across provinces in Thailand. Multilocus data (mitochondrial sequence data and 12 nuclear microsatellites) and individual based analyses were used to test for cryptic lineages. We found strong linkage disequilibrium patterns between mitochondrial haplotypes and nuclear-identified genetic clusters. There were clearly two divergent and independent clades. Moreover, within each clade, the data suggested additional substructure where individual provinces were likely to harbor unique genetic clusters. The results indicate there are at minimum two and possibly up to eight cryptic lineages within the assumed single species of A. cantonensis. Importantly, the two main clades do not show geographic affiliation and can be found in sympatry. With recent studies highlighting A. cantonensis strain diversity in pathogenicity and infectivity, it will be important to determine if these critical epidemiological traits are associated with specific lineages.

Characterization of nine microsatellite loci for Dicrocoelium dendriticum, an emerging liver fluke of ungulates in North America, and their use to detect clonemates and random mating.

This study characterizes polymorphic microsatellite loci from adults of the liver fluke Dicrocoelium dendriticum sampled from a population of sympatric beef cattle and wapiti in a region of emergence in southern Alberta, Canada. We also scrutinized the markers to validate their use in studying the population genetics of this complex life cycle parasite. Among the nine loci described, four deviated significantly from Hardy Weinberg Equilibrium (HWE) due to technical artefacts. The remaining five loci were in HWE. These five provided sufficient resolution to identify clonemates produced from the obligate asexual reproduction phase of the life cycle in snails and to assess the impact of non-random transmission of clonemates on measures of FIS, FST and genotypic disequilibrium. Excluding clonemates, we show that the sub-population of worms was in HWE, that average FIS within hosts was 0.003 (p=0.4922) and that there was no population genetic structure among hosts FST=0.001 (p=0.3243). These markers will be useful for studies of Dicrocoelium dendriticum ecology, transmission, and evolution.

Physiological Status Drives Metabolic Rate in Mediterranean Geckos Infected with Pentastomes.

Negative effects of parasites on their hosts are well documented, but the proximate mechanisms by which parasites reduce their host's fitness are poorly understood. For example, it has been suggested that parasites might be energetically demanding. However, a recent meta-analysis suggests that they have statistically insignificant effects on host resting metabolic rate (RMR). It is possible, though, that energetic costs associated with parasites are only manifested during and/or following periods of activity. Here, we measured CO2 production (a surrogate for metabolism) in Mediterranean geckos (Hemidactylus turcicus) infected with a lung parasite, the pentastome Raillietiella indica, under two physiological conditions: rested and recently active. In rested geckos, there was a negative, but non-significant association between the number of pentastomes (i.e., infection intensity) and CO2 production. In recently active geckos (chased for 3 minutes), we recorded CO2 production from its maximum value until it declined to a stationary phase. We analyzed this decline as a 3 phase function (initial decline, secondary decline, stationary). Geckos that were recently active showed, in the secondary phase, a significant decrease in CO2 production as pentastome intensity increased. Moreover, duration of the secondary phase showed a significant positive association with the number of pentastomes. These results suggest that the intensity of pentastome load exerts a weak effect on the metabolism of resting geckos, but a strong physiological effect on geckos that have recently been active; we speculate this occurs via mechanical constraints on breathing. Our results provide a potential mechanism by which pentastomes can reduce gecko fitness.

Evolutionary consequence of a change in life cycle complexity: A link between precocious development and evolution toward female-biased sex allocation in a hermaphroditic parasite.

The evolutionary consequences of changes in the complex life cycles of parasites are not limited to the traits that directly affect transmission. For instance, mating systems that are altered due to precocious sexual maturation in what is typically regarded as an intermediate host may impact opportunities for outcrossing. In turn, reproductive traits may evolve to optimize sex allocation. Here, we test the hypothesis that sex allocation evolved toward a more female-biased function in populations of the hermaphroditic digenean trematode Alloglossidium progeneticum that can precociously reproduce in their second hosts. In these precocious populations, parasites are forced to self-fertilize as they remain encysted in their second hosts. In contrast, parasites in obligate three-host populations have more opportunities to outcross in their third host. We found strong support that in populations with precocious development, allocation to male resources was greatly reduced. We also identified a potential phenotypically plastic response in a body size sex allocation relationship that may be driven by the competition for mates. These results emphasize how changes in life cycle patterns that alter mating systems can impact the evolution of reproductive traits in parasites.

Species mtDNA genetic diversity explained by infrapopulation size in a host-symbiont system.

Understanding what shapes variation in genetic diversity among species remains a major challenge in evolutionary ecology, and it has been seldom studied in parasites and other host-symbiont systems. Here, we studied mtDNA variation in a host-symbiont non-model system: 418 individual feather mites from 17 feather mite species living on 17 different passerine bird species. We explored how a surrogate of census size, the median infrapopulation size (i.e., the median number of individual parasites per infected host individual), explains mtDNA genetic diversity. Feather mite species genetic diversity was positively correlated with mean infrapopulation size, explaining 34% of the variation. As expected from the biology of feather mites, we found bottleneck signatures for most of the species studied but, in particular, three species presented extremely low mtDNA diversity values given their infrapopulation size. Their star-like haplotype networks (in contrast with more reticulated networks for the other species) suggested that their low genetic diversity was the consequence of severe bottlenecks or selective sweeps. Our study shows for the first time that mtDNA diversity can be explained by infrapopulation sizes, and suggests that departures from this relationship could be informative of underlying ecological and evolutionary processes.

Alloglossidium floridense n. sp. (Digenea: Macroderoididae) from a spring run in North Central Florida.

A new species of Alloglossidium is described from the intestines of 2 madtom species (Noturus leptacanthus and Noturus gyrinus) that were collected from the run of a small, unnamed spring system that drains into the Santa Fe River, Florida. Alloglossidium floridense n. sp. is morphologically very similar to other nonprecocious Alloglossidium spp. that use ictalurids as definitive hosts, but can be distinguished by a combination of its smaller overall size (length and width), large eggs in relation to its small body size, position of the vitellaria, ovary shape, and position of the ovary in relation to the cirrus sac. A comparison of nuclear rDNA sequences (spanning partial 18s, complete ITS1, 5.8s, ITS2, and partial 28s regions) showed that A. floridense n. sp. diverged by 0.70-3.17% from the other 4, nonprecocious species that infect ictalurids (Alloglossidium corti, Alloglossidium fonti, Alloglossidium geminum, and Alloglossidium kenti). The new species of Alloglossidium, described herein, is the first of the genus to be reported from Florida and the first to be recorded from N. leptacanthus . In light of the subtle morphological differences among the nonprecocious species that infect ictalurids, we discuss how previous descriptions of species traits that are not supported with genetic data are difficult to interpret because of the possible past nonrecognition of distinct species.

Testing local-scale panmixia provides insights into the cryptic ecology, evolution, and epidemiology of metazoan animal parasites.

When every individual has an equal chance of mating with other individuals, the population is classified as panmictic. Amongst metazoan parasites of animals, local-scale panmixia can be disrupted due to not only non-random mating, but also non-random transmission among individual hosts of a single host population or non-random transmission among sympatric host species. Population genetics theory and analyses can be used to test the null hypothesis of panmixia and thus, allow one to draw inferences about parasite population dynamics that are difficult to observe directly. We provide an outline that addresses 3 tiered questions when testing parasite panmixia on local scales: is there greater than 1 parasite population/species, is there genetic subdivision amongst infrapopulations within a host population, and is there asexual reproduction or a non-random mating system? In this review, we highlight the evolutionary significance of non-panmixia on local scales and the genetic patterns that have been used to identify the different factors that may cause or explain deviations from panmixia on a local scale. We also discuss how tests of local-scale panmixia can provide a means to infer parasite population dynamics and epidemiology of medically relevant parasites.

Host genetics and population structure effects on parasitic disease.

Host genetic factors exert significant influences on differential susceptibility to many infectious diseases. In addition, population structure of both host and parasite may influence disease distribution patterns. In this study, we assess the effects of population structure on infectious disease in two populations in which host genetic factors influencing susceptibility to parasitic disease have been extensively studied. The first population is the Jirel population of eastern Nepal that has been the subject of research on the determinants of differential susceptibility to soil-transmitted helminth infections. The second group is a Brazilian population residing in an area endemic for Trypanosoma cruzi infection that has been assessed for genetic influences on differential disease progression in Chagas disease. For measures of Ascaris worm burden, within-population host genetic effects are generally more important than host population structure factors in determining patterns of infectious disease. No significant influences of population structure on measures associated with progression of cardiac disease in individuals who were seropositive for T. cruzi infection were found.