Genetic assignment methods for gaining insight into the management of infectious disease by understanding pathogen, vector, and host movement

Genetic analysis of the tomato russet mite provides evidence of oligophagy and a widespread pestiferous haplotype

Worldwide, the tomato russet mite (TRM), Aculops lycopersici (Eriophyidae), is a key pest on cultivated tomato in addition to infesting other cultivated and wild Solanaceae; however, basic information on TRM supporting effective control strategies is still lacking, mainly regarding its taxonomic status and genetic diversity and structure. As A. lycopersici is reported on different species and genera of host plants, populations associated with different host plants may constitute specialized cryptic species, as shown for other eriophyids previously considered generalists. The main aims of this study were to (i) confirm the TRM taxonomic unity of populations from different host plants and localities as well as the species' oligophagy, and (ii) to advance the understanding of TRM host relationship and invasion history. For this purpose, we evaluated the genetic variability and structure of populations from different host plants along crucial areas of occurrence, including the area of potential origin, based on DNA sequences of mitochondrial (cytochrome c oxidase subunit I) and nuclear (internal transcribed spacer, D2 28S) genomic regions. Specimens from South America (Brazil) and Europe (France, Italy, Poland, The Netherlands) were collected from tomato and other solanaceous species from the genera Solanum and Physalis. Final TRM datasets were composed of 101, 82 and 50 sequences from the COI (672 bp), ITS (553 bp) and D2 (605 bp) regions, respectively. Distributions and frequencies of haplotypes (COI) and genotypes (D2 and ITS1) were inferred; pairwise genetic distance comparisons, and phylogenetic analysis were performed, including Bayesian Inference (BI) combined analysis. Our results showed that genetic divergences for mitochondrial and nuclear genomic regions from TRM associated with different host plants were lower than those observed in other eriophyid taxa, confirming conspecificity of TRM populations and oligophagy of this eriophyid mite. Four haplotypes (cH) were identified from the COI sequences with cH1 being the most frequent, representing 90% of all sequences occurring in all host plants studied (Brazil, France, The Netherlands); the other haplotypes were present exclusively in Brazilian populations. Six variants (I) were identified from the ITS sequences: I-1 was the most frequent (76.5% of all sequences), spread in all countries and associated with all host plants, except S. nigrum. Just one D2 sequence variant was found in all studied countries. The genetic homogeneity among populations highlights the occurrence of a highly invasive and oligophagous haplotype. These results failed to corroborate the hypothesis that differential symptomatology or damage intensity among tomato varieties and solanaceous host plants could be due to the genetic diversity of the associated mite populations. The genetic evidence, along with the history of spread of cultivated tomato, corroborates the hypothesis of a South American origin of TRM.

Can genetic assignment tests provide insight on the influence of captive egression on the epizootiology of chronic wasting disease?

Identifying the sources of ongoing and novel disease outbreaks is critical for understanding the diffusion of epizootic diseases. Identifying infection sources is difficult when few physical differences separate individuals with different origins. Genetic assignment procedures show great promise for assessing transmission dynamics in such situations. Here, we use genetic assignment tests to determine the source of chronic wasting disease infections in free-ranging white-tailed deer (Odocoileus virginianus) populations. Natural dispersal is thought to facilitate the geographic diffusion of chronic wasting disease, but egression from captive cervid populations represents an alternative source of infection that is difficult to detect due to physical similarities with wild deer. Simulated reference populations were created based on allele frequencies from 1,912 empirical microsatellite genotypes collected in four sampling subregions and five captive facilities. These reference populations were used to assess the likelihood of ancestry and assignment of 1,861 free-ranging deer (1,834 noninfected and 27 infected) and 51 captive individuals to captive or wild populations. The ancestry (Q) and assignment scores (A) for free-ranging deer to wild populations were high (average Q wild = 0.913 and average A wild = 0.951, respectively), but varied among subregions (Q wild = 0.800-0.947, A wild = 0.857-0.976). These findings suggest that captive egression and admixture are rare, but risk may not be spatially uniform. Ancestry and assignment scores for two free-ranging deer with chronic wasting disease sampled in an area where chronic wasting disease was previously unobserved in free-ranging herds indicated a higher likelihood of assignment and proportion of ancestry attributable to captive populations. While we cannot directly assign these individuals to infected facilities, these findings suggest that rare egression events may influence the epizootiology of chronic wasting disease in free-ranging populations. Continued disease surveillance and genetic analyses may further elucidate the relative disease risk attributable to captive and wild sources.

Infections on the move: how transient phases of host movement influence disease spread

Animal movement impacts the spread of human and wildlife diseases, and there is significant interest in understanding the role of migrations, biological invasions and other wildlife movements in spatial infection dynamics. However, the influence of processes acting on infections during transient phases of host movement is poorly understood. We propose a conceptual framework that explicitly considers infection dynamics during transient phases of host movement to better predict infection spread through spatial host networks. Accounting for host transient movement captures key processes that occur while hosts move between locations, which together determine the rate at which hosts spread infections through networks. We review theoretical and empirical studies of host movement and infection spread, highlighting the multiple factors that impact the infection status of hosts. We then outline characteristics of hosts, parasites and the environment that influence these dynamics. Recent technological advances provide disease ecologists unprecedented ability to track the fine-scale movement of organisms. These, in conjunction with experimental testing of the factors driving infection dynamics during host movement, can inform models of infection spread based on constituent biological processes.

Global mtDNA genetic structure and hypothesized invasion history of a major pest of citrus, Diaphorina citri (Hemiptera: Liviidae)

The Asian citrus psyllid Diaphorina citri Kuwayama is a key pest of citrus as the vector of the bacterium causing the “huanglongbing” disease (HLB). To assess the global mtDNA population genetic structure, and possible dispersal history of the pest, we investigated genetic variation at the COI gene collating newly collected samples with all previously published data. Our dataset consists of 356 colonies from 106 geographic sites worldwide. High haplotype diversity (H‐mean = 0.702 ± 0.017), low nucleotide diversity (π‐mean = 0.003), and significant positive selection (Ka/Ks = 32.92) were observed. Forty‐four haplotypes (Hap) were identified, clustered into two matrilines: Both occur in southeastern and southern Asia, North and South America, and Africa; lineages A and B also occur in eastern and western Asia, respectively. The most abundant haplotypes were Hap4 in lineage A (35.67%), and Hap9 in lineage B (41.29%). The haplotype network identified them as the ancestral haplotypes within their respective lineages. Analysis of molecular variance showed significant genetic structure (F_ST = 0.62, p < .0001) between the lineages, and population genetic analysis suggests geographic structuring. We hypothesize a southern and/or southeastern Asia origin, three dispersal routes, and parallel expansions of two lineages. The hypothesized first route involved the expansion of lineage B from southern Asia into North America via West Asia. The second, the expansion of some lineage A individuals from Southeast Asia into East Asia, and the third involved both lineages from Southeast Asia spreading westward into Africa and subsequently into South America. To test these hypotheses and gain a deeper understanding of the global history of D. citri, more data‐rich approaches will be necessary from the ample toolkit of next‐generation sequencing (NGS). However, this study may serve to guide such sampling and in the development of biological control programs against the global pest D. citri.

Genetic Evidence of Contemporary Dispersal of the Intermediate Snail Host of Schistosoma japonicum: Movement of an NTD Host Is Facilitated by Land Use and Landscape Connectivity

BACKGROUND

While the dispersal of hosts and vectors-through active or passive movement-is known to facilitate the spread and re-emergence of certain infectious diseases, little is known about the movement ecology of Oncomelania spp., intermediate snail host of the parasite Schistosoma japonicum, and its consequences for the spread of schistosomiasis in East and Southeast Asia. In China, despite intense control programs aimed at preventing schistosomiasis transmission, there is evidence in recent years of re-emergence and persistence of infection in some areas, as well as an increase in the spatial extent of the snail host. A quantitative understanding of the dispersal characteristics of the intermediate host can provide new insights into the spatial dynamics of transmission, and can assist public health officials in limiting the geographic spread of infection.

METHODOLOGY/PRINCIPAL FINDINGS

Oncomelania hupensis robertsoni snails (n = 833) were sampled from 29 sites in Sichuan, China, genotyped, and analyzed using Bayesian assignment to estimate the rate of recent snail migration across sites. Landscape connectivity between each site pair was estimated using the geographic distance distributions derived from nine environmental models: Euclidean, topography, incline, wetness, land use, watershed, stream use, streams and channels, and stream velocity. Among sites, 14.4% to 32.8% of sampled snails were identified as recent migrants, with 20 sites comprising >20% migrants. Migration rates were generally low between sites, but at 8 sites, over 10% of the overall host population originated from one proximal site. Greater landscape connectivity was significantly associated with increased odds of migration, with the minimum path distance (as opposed to median or first quartile) emerging as the strongest predictor across all environmental models. Models accounting for land use explained the largest proportion of the variance in migration rates between sites. A greater number of irrigation channels leading into a site was associated with an increase in the site's propensity to both attract and retain snails.

CONCLUSIONS/SIGNIFICANCE

Our findings have important implications for controlling the geographic spread of schistosomiasis in China, through improved understanding of the dispersal capacity of the parasite's intermediate host.

The genetic structure of an invasive pest, the Asian citrus psyllid Diaphorina citri (Hemiptera: Liviidae)

The Asian citrus psyllid Diaphorina citri is currently the major threat to the citrus industry as it is the vector of Candidatus Liberibacter, the causal agent of huanglongbing disease (HLB). D. citri is native to Asia and now colonizes the Americas. Although it has been known in some countries for a long time, invasion routes remain undetermined. There are no efficient control methods for the HLB despite the intensive management tools currently in use. We investigated the genetic variability and structure of populations of D. citri to aid in the decision making processes toward sustainable management of this species/disease. We employed different methods to quantify and compare the genetic diversity and structure of D. citri populations among 36 localities in Brazil, using an almost complete sequence of the cytochrome oxidase I (COI) gene. Our analyses led to the identification of two geographically and genetically structured groups. The indices of molecular diversity pointed to a recent population expansion, and we discuss the role of multiple invasion events in this scenario. We also argue that such genetic diversity and population structure may have implications for the best management strategies to be adopted for controlling this psyllid and/or the disease it vectors in Brazil.

Simulating population genetics of pathogen vectors in changing landscapes: guidelines and application with Triatoma brasiliensis

Background

Understanding the mechanisms that influence the population dynamics and spatial genetic structure of the vectors of pathogens infecting humans is a central issue in tropical epidemiology. In view of the rapid changes in the features of landscape pathogen vectors live in, this issue requires new methods that consider both natural and human systems and their interactions. In this context, individual-based model (IBM) simulations represent powerful yet poorly developed approaches to explore the response of pathogen vectors in heterogeneous social-ecological systems, especially when field experiments cannot be performed.

Methodology/Principal Findings

We first present guidelines for the use of a spatially explicit IBM, to simulate population genetics of pathogen vectors in changing landscapes. We then applied our model with Triatoma brasiliensis, originally restricted to sylvatic habitats and now found in peridomestic and domestic habitats, posing as the most important Trypanosoma cruzi vector in Northeastern Brazil. We focused on the effects of vector migration rate, maximum dispersal distance and attraction by domestic habitat on T. brasiliensis population dynamics and spatial genetic structure. Optimized for T. brasiliensis using field data pairwise fixation index (FST) from microsatellite loci, our simulations confirmed the importance of these three variables to understand vector genetic structure at the landscape level. We then ran prospective scenarios accounting for land-use change (deforestation and urbanization), which revealed that human-induced land-use change favored higher genetic diversity among sampling points.

Conclusions/Significance

Our work shows that mechanistic models may be useful tools to link observed patterns with processes involved in the population genetics of tropical pathogen vectors in heterogeneous social-ecological landscapes. Our hope is that our study may provide a testable and applicable modeling framework to a broad community of epidemiologists for formulating scenarios of landscape change consequences on vector dynamics, with potential implications for their surveillance and control.

Worldwide, humans are modifying landscapes at an unprecedented rate. These modifications have an influence on the ecology of pathogen vectors, yet this issue has received relatively little input from modeling research. The current study presents guidelines for the use of a modeling framework for the representation of the dynamics and spatial genetic structure of pathogen vectors. It allows considering spatiotemporal landscape modifications explicitly, to represent human-altered modifications and consequences. We applied this modeling framework to Triatoma brasiliensis, vector of the pathogen Trypanosoma cruzi responsible for the Chagas disease, in the semi-arid Northeastern Brazil. Using field data of pairwise fixation index (FST) from microsatellite loci, we found that migration rate, maximum dispersal distance and attraction by domestic habitat were all key parameters to understand vector spatial genetic structure at the landscape level. At the interface across disciplines, this study provides to the community of epidemiologists a testable and applicable framework to foresee landscape modification consequences on vector dynamics and genetic structure, with potential implications for their surveillance and control.

Juxtaposition between host population structures: implications for disease transmission in a sympatric cervid community

Sympatric populations of phylogenetically related species are often vulnerable to similar communicable diseases. Although some host populations may exhibit spatial structure, other hosts within the community may have unstructured populations. Thus, individuals from unstructured host populations may act as interspecific vectors among discrete subpopulations of sympatric alternate hosts. We used a cervid-bovine tuberculosis (Mycobacterium bovis) system to investigate the landscape-scale potential for bovine tuberculosis transmission within a nonmigratory white-tailed deer (Odocoileus virginianus) and elk (Cervus canadensis) community. Using landscape population genetics, we tested for genetic and spatial structure in white-tailed deer. We then compared these findings with the sympatric elk population that is structured and which has structure that correlates spatially and genetically to physiognomic landscape features. Despite genetic structure that indicates the white-tailed deer population forms three sympatric clusters, the absence of spatial structure suggested that intraspecific pathogen transmission is not likely to be limited by physiognomic landscape features. The potential for intraspecific transmission among subpopulations of elk is low due to spatial population structure. Given that white-tailed deer are abundant, widely distributed, and exhibit a distinct lack of spatial population structure, white-tailed deer likely pose a greater threat as bovine tuberculosis vectors among elk subpopulations than elk.

Approaches to genotyping individual miracidia of Schistosoma japonicum

Molecular genetic tools are needed to address questions as to the source and dynamics of transmission of the human blood fluke Schistosoma japonicum in regions where human infections have reemerged, and to characterize infrapopulations in individual hosts. The life stage that interests us as a target for collecting genotypic data is the miracidium, a very small larval stage that consequently yields very little DNA for analysis. Here, we report the successful development of a multiplex format permitting genotyping of 17 microsatellite loci in four sequential multiplex reactions using a single miracidium held on a Whatman Classic FTA indicating card. This approach was successful after short storage periods, but after long storage (>4 years), considerable difficulty was encountered in multiplex genotyping, necessitating the use of whole genome amplification (WGA) methods. WGA applied to cards stored for long periods of time resulted in sufficient DNA for accurate and repeatable genotyping. Trials and tests of these methods, as well as application to some field-collected samples, are reported, along with the discussion of the potential insights to be gained from such techniques. These include recognition of sibships among miracidia from a single host, and inference of the minimum number of worm pairs that might be present in a host.