Insights into the genetic history of Green‐legged Partridgelike fowl: mtDNA and genome‐wide SNP analysis

The Green‐legged Partridgelike (GP) fowl, an old native Polish breed, is characterised by reseda green‐coloured shanks rather than yellow, white, slate or black commonly observed across most domestic breeds of chicken. Here, we investigate the origin, genetic relationships and structure of the GP fowl using mtDNA D‐loop sequencing and genome‐wide SNP analysis. Genome‐wide association analysis between breeds enables us to verify the genetic control of the reseda green shank phenotype, a defining trait for the breed. Two mtDNA D‐loop haplogroups and three autosomal genetic backgrounds are revealed. Significant associations of SNPs on chromosomes GGA24 and GGAZ indicate that the reseda green leg phenotype is associated with recessive alleles linked to the W and Id loci. Our results provide new insights into the genetic history of European chicken, indicating an admixd origin of East European traditional breeds of chicken on the continent, as supported by the presence of the reseda green phenotype and the knowledge that the GP fowl as a breed was developed before the advent of commercial stocks.     

High MHC gene copy number maintains diversity despite homozygosity in a Critically Endangered single‐island endemic bird,but no evidence of MHC‐based mate choice

Small population sizes can, over time, put species at risk due to the loss of genetic variation and the deleterious effects of inbreeding. Losing diversity in the major histocompatibility complex (MHC) could be particularly harmful, given its key role in the immune system. Here, we assess MHC class I (MHC‐I) diversity and its effects on mate choice and survival in the Critically Endangered Raso lark Alauda razae, a species restricted to the 7 km² islet of Raso, Cape Verde, since ~1460, whose population size has dropped as low as 20 pairs. Exhaustively genotyping 122 individuals, we find no effect of MHC‐I genotype/diversity on mate choice or survival. However, we demonstrate that MHC‐I diversity has been maintained through extreme bottlenecks by retention of a high number of gene copies (at least 14), aided by cosegregation of multiple haplotypes comprising 2–8 linked MHC‐I loci. Within‐locus homozygosity is high, contributing to low population‐wide diversity. Conversely, each individual had comparably many alleles, 6–16 (average 11), and the large and divergent haplotypes occur at high frequency in the population, resulting in high within‐individual MHC‐I diversity. This functional immune gene diversity will be of critical importance for this highly threatened species’ adaptive potential.     

Genomic analysis of a cardinalfish with larval homing potential reveals genetic admixture in the Okinawa Islands

Discrepancies between potential and observed dispersal distances of reef fish indicate the need for a better understanding of the influence of larval behaviour on recruitment and dispersal. Population genetic studies can provide insight on the degree to which populations are connected, and the development of restriction site‐associated sequencing (RAD‐Seq) methods has made such studies of nonmodel organisms more accessible. We applied double‐digest RAD‐Seq methods to test for population differentiation in the coral reef‐dwelling cardinalfish, Siphamia tubifer, which based on behavioural studies, have the potential to use navigational cues to return to natal reefs. Analysis of 11,836 SNPs from fish collected at coral reefs in Okinawa, Japan, from eleven locations over 3 years reveals little genetic differentiation between groups of S. tubifer at spatial scales from 2 to 140 km and between years at one location: pairwise F_ST values were between 0.0116 and 0.0214. These results suggest that the Kuroshio Current largely influences larval dispersal in the region, and in contrast to expectations based on studies of other cardinalfishes, there is no evidence of population structure for S. tubifer at the spatial scales examined. However, analyses of outlier loci putatively under selection reveal patterns of temporal differentiation that indicate high population turnover and variable larval supply from divergent source populations between years. These findings highlight the need for more studies of fishes across various geographic regions that also examine temporal patterns of genetic differentiation to better understand the potential connections between early life‐history traits and connectivity of reef fish populations.     

Phylogeography of a pearl oyster (Pinctada maxima) across the Indo‐Australian Archipelago: evidence of strong regional structure and population expansions but no phylogenetic breaks

This study investigates the genetic structure and phylogeography of a broadcast spawning bivalve mollusc, Pinctada maxima, throughout the Indo‐West Pacific and northern Australia. DNA sequence variation of the mitochondrial cytochrome oxidase subunit I (COI) gene was analysed in 367 individuals sampled from nine populations across the Indo‐West Pacific. Hierarchical AMOVA indicated strong genetic structuring amongst populations (Φ_ST = 0.372, P < 0.001); however, sequence divergence between the 47 haplotypes detected was low (maximum 1.8% difference) and no deep phylogenetic divergence was observed. Results suggest the presence of genetic barriers isolating populations of the South China Sea and central Indonesian regions, which, in turn, show patterns of historical separation from northern Australian regions. In P. maxima, historical vicariance during Pleistocene low sea levels is likely to have restricted planktonic larval transport, causing genetic differentiation amongst populations. However, low genetic differentiation is observed where strong ocean currents are present and is most likely due to contemporary larval transport along these pathways. Geographical association with haplotype distributions may indicate signs of early lineage sorting arising from historical population separations, yet an absence of divergent phylogenetic clades related to geography could be the consequence of periodic pulses of high genetic exchange. We compare our results with previous microsatellite DNA analysis of these P. maxima populations, and discuss implications for future conservation management of this species. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 107 , 632–646.     

Morphometric and genetic differentiation among populations of flat‐headed cusimanse (Crossarchus platycephalus) in Nigeria

Geographic barriers can partition genetic diversity among populations and drive evolutionary divergence between populations, promoting the speciation process and affecting conservation goals. We integrated morphological and genomic data to assess the distribution of variation in the flat‐headed cusimanse (Crossarchus platycephalus), a species of least conservation concern, on either side of the River Niger in Nigeria. Ecological disturbances affect the conservation status of many other animals in this region. The two populations were differentiated in the snout and fore limbs, with greater morphological diversity in the western population. We used Restriction site Associated DNA sequencing (RAD‐seq) and identified two genotypic clusters in a STRUCTURE analysis. Individuals from the eastern population are almost entirely assigned to one cluster, whereas genotypes from the western population are a mixture of the two clusters. The population from west of the River Niger also had higher heterozygosity. The morphological and population genetic data are therefore in agreement that the population from west of the River Niger is more diverse than the eastern population, and the eastern population contains a subset of the genetic variation found in the western population. Our results demonstrate that combining morphological and genotypic measures of diversity can provide a congruent picture of the distribution of intraspecific variation. The results also suggest that future work should explore the role of the River Niger as a natural barrier to migration in Nigeria.     

Habitat continuity and stepping‐stone oceanographic distances explain population genetic connectivity of the brown alga Cystoseira amentacea

Effective predictive and management approaches for species occurring in a metapopulation structure require good understanding of interpopulation connectivity. In this study, we ask whether population genetic structure of marine species with fragmented distributions can be predicted by stepping‐stone oceanographic transport and habitat continuity, using as model an ecosystem‐structuring brown alga, Cystoseira amentacea var. stricta. To answer this question, we analysed the genetic structure and estimated the connectivity of populations along discontinuous rocky habitat patches in southern Italy, using microsatellite markers at multiple scales. In addition, we modelled the effect of rocky habitat continuity and ocean circulation on gene flow by simulating Lagrangian particle dispersal based on ocean surface currents allowing multigenerational stepping‐stone dynamics. Populations were highly differentiated, at scales from few metres up to thousands of kilometres. The best possible model fit to explain the genetic results combined current direction, rocky habitat extension and distance along the coast among rocky sites. We conclude that a combination of variable suitable habitat and oceanographic transport is a useful predictor of genetic structure. This relationship provides insight into the mechanisms of dispersal and the role of life‐history traits. Our results highlight the importance of spatially explicit modelling of stepping‐stone dynamics and oceanographic directional transport coupled with habitat suitability, to better describe and predict marine population structure and differentiation. This study also suggests the appropriate spatial scales for the conservation, restoration and management of species that are increasingly affected by habitat modifications.     

Genetic assignment of recruits reveals short‐ and long‐distance larval dispersal in Pocillopora damicornis on the Great Barrier Reef

Understanding connectivity of coral populations among and within reefs over ecologically significant timescales is essential for developing evidence‐based management strategies, including the design of marineprotected areas. Here, we present the first assessment of contemporary connectivity among populations of two Molecular Operational Taxonomic Units (MOTUs) of the brooding coral Pocillopora damicornis. We used individual‐based genetic assignment methods to identify the proportions of philopatric and migrant larval recruits, settling over 12 months at sites around Lizard Island (northern Great Barrier Reef [GBR]) and over 24 months at sites around the Palms Islands (central GBR). Overall, we found spatially and temporally variable rates of self‐recruitment and dispersal, demonstrating the importance of variation in local physical characteristics in driving dispersal processes. Recruitment patterns and inferred dispersal distances differed between the two P. damicornis MOTUs, with type α recruits exhibiting predominantly philopatric recruitment, while the majority of type β recruits were either migrants from identified putative source populations or assumed migrants based on genetic exclusion from all known populations. While P. damicornis invests much energy into brooding clonal larvae, we found that only 15% and 7% of type α and type β recruits, respectively, were clones of sampled adult colonies or other recruits, challenging the hypothesis that reproduction is predominantly asexual in this species on the GBR. We explain high rates of self‐recruitment and low rates of clonality in these MOTUs by suggesting that locally retained larvae originate predominantly from spawned gametes, while brooded larvae are mainly vagabonds.     

Pedigree analysis for the genetic management of group‐living species

Captive breeding programs are an important tool for the conservation of endangered species. These programs are commonly managed using pedigrees containing information about the history of each individual's family, such as breeding pairs and parentage. However, there are some species that are kept in groups where it is hard to distinguish between particular individuals within the group, making it very difficult to record any information at an individual level. Currently, software and methods commonly used for registering and analyzing pedigrees to help manage populations at an individual level are not adequate for managing these group‐living species. Therefore, there is a need to further develop these tools and methodologies for pedigree analysis to better manage group‐living species. PMx is a program used for the management of ex situ populations in zoos and aquariums. We adapted the pedigree analysis method implemented in PMx to analyze pedigrees (records of descendant lineages) of group‐living species. In addition, we developed a group pedigree data entry sheet and group2PMx, a converter program that enables group datasets to be imported into PMx. We show how pedigree analysis of a group‐living species can be used for population management using the studbook of the endangered Texas blind cave salamander Eurycea rathbuni. Such analyses of the pedigree of groups can improve the management of group‐living species in ex situ breeding programs. Firstly, it enables better management decisions based on more accurate genetic measures between groups, allowing for greater control of inbreeding. Secondly, it can improve the conditions in which group‐living species are held by adapting husbandry practices to better reflect conditions of these species living in the wild. The use of the spreadsheet and group2PMx extends the application of PMx_, allowing conservation managers and other institutions outside the zoo and aquarium community to easily import and analyze their pedigree data.     

New view of population genetics of zooplankton: RAD‐seq analysis reveals population structure of the North Atlantic planktonic copepod Centropages typicus

Detection of population genetic structure of zooplankton at medium‐to‐small spatial scales in the absence of physical barriers has remained challenging and controversial. The large population sizes and high rates of gene flow characteristic of zooplankton have made resolution of geographical differentiation very difficult, especially when using few genetic markers and assuming equilibrium conditions. Next‐generation sequencing now allows simultaneous sampling of hundreds to thousands of genetic markers; new analytical approaches allow studies under nonequilibrium conditions and directional migration. Samples of the North Atlantic Ocean planktonic copepod, Centropages typicus, were analysed using restriction site‐associated DNA (RAD) sequencing on a PROTON platform. Although prior studies revealed no genetic differentiation of populations across the geographical range of the species, analysis of RAD tags showed significant structure across the North Atlantic Ocean. We also compared the likelihood for models of connectivity among NW Atlantic populations under various directional flow scenarios that replicate oceanographic conditions of the sampled domain. High‐density marker sampling with RAD sequencing markedly outperformed other technical and analytical approaches in detection of population genetic structure and characterization of connectivity of this high geneflow zooplankton species.     

Range‐wide analysis of genetic structure in a widespread,highly mobile species (Odocoileus hemionus) reveals the importance of historical biogeography

Highly mobile species that thrive in a wide range of habitats are expected to show little genetic differentiation across their range. A limited but growing number of studies have revealed that patterns of broad‐scale genetic differentiation can and do emerge in vagile, continuously distributed species. However, these patterns are complex and often shaped by both historical and ecological factors. Comprehensive surveys of genetic variation at a broad scale and at high resolution are useful for detecting cryptic spatial genetic structure and for investigating the relative roles of historical and ecological processes in structuring widespread, highly mobile species. In this study, we analysed 10 microsatellite loci from over 1900 samples collected across the full range of mule deer (Odocoileus hemionus), one of the most widely distributed and abundant of all large mammal species in North America. Through both individual‐ and population‐based analyses, we found evidence for three main genetic lineages, one corresponding to the ‘mule deer’ morphological type and two to the ‘black‐tailed deer’ type. Historical biogeographic events likely are the primary drivers of genetic divergence in this species; boundaries of the three lineages correspond well with predictions based on Pleistocene glacial cycles, and substructure within each lineage demonstrates island vicariance. However, across large geographic areas, including the entire mule deer lineage, we found that genetic variation fit an isolation‐by‐distance pattern rather than discrete clusters. A lack of genetic structure across wide geographic areas of the continental west indicates that ecological processes have not resulted in restrictions to gene flow sufficient for spatial genetic structure to emerge. Our results have important implications for our understanding of evolutionary mechanisms of divergence, as well as for taxonomy, conservation and management.     

metapop2: Re‐implementation of software for the analysis and management of subdivided populations using gene and allelic diversity

Management programmes often have to make decisions based on the analysis of the genetic properties and diversity of populations. Expected heterozygosity (or gene diversity) and population structure parameters are often used to make recommendations for conservation, such as avoidance of inbreeding or migration across subpopulations. Allelic diversity, however, can also provide complementary and useful information for conservation programmes, as it is highly sensitive to population bottlenecks, and is more related to long‐term selection response than heterozygosity. Here we present a completely revised and updated re‐implementation of the software metapop for the analysis of diversity in subdivided populations, as well as a tool for the management and dynamic estimation of optimal contributions in conservation programmes. This new update includes computation of allelic diversity for population analysis and management, as well as a simulation mode to forecast the consequences of taking different management strategies over time. Furthermore, the new implementation in C++ includes code optimization and improved memory usage, allowing for fast analysis of large data sets including single nucleotide polymorphism markers, as well as enhanced cross‐software and cross‐platform compatibility.     

A novel approach to parasite population genetics: Experimental infection reveals geographic differentiation,recombination and host‐mediated population structure in Pasteuria ramosa,a bacterial parasite of Daphnia

The population structure of parasites is central to the ecology and evolution of host‐parasite systems. Here, we investigate the population genetics of Pasteuria ramosa, a bacterial parasite of Daphnia. We used natural P. ramosa spore banks from the sediments of two geographically well‐separated ponds to experimentally infect a panel of Daphnia magna host clones whose resistance phenotypes were previously known. In this way, we were able to assess the population structure of P. ramosa based on geography, host resistance phenotype and host genotype. Overall, genetic diversity of P. ramosa was high, and nearly all infected D. magna hosted more than one parasite haplotype. On the basis of the observation of recombinant haplotypes and relatively low levels of linkage disequilibrium, we conclude that P. ramosa engages in substantial recombination. Isolates were strongly differentiated by pond, indicating that gene flow is spatially restricted. Pasteuria ramosa isolates within one pond were segregated completely based on the resistance phenotype of the host—a result that, to our knowledge, has not been previously reported for a nonhuman parasite. To assess the comparability of experimental infections with natural P. ramosa isolates, we examined the population structure of naturally infected D. magna native to one of the two source ponds. We found that experimental and natural infections of the same host resistance phenotype from the same source pond were indistinguishable, indicating that experimental infections provide a means to representatively sample the diversity of P. ramosa while reducing the sampling bias often associated with studies of parasite epidemics. These results expand our knowledge of this model parasite, provide important context for the large existing body of research on this system and will guide the design of future studies of this host‐parasite system.     

Large‐scale introduction of the Indo‐Pacific damselfish Abudefduf vaigiensis into Hawai'i promotes genetic swamping of the endemic congener A. abdominalis

Hybridization in the ocean was once considered rare, a process prohibited by the rapid evolution of intrinsic reproductive barriers in a high‐dispersal medium. However, recent genetic surveys have prompted a reappraisal of marine hybridization as an important demographic and evolutionary process. The Hawaiian Archipelago offers an unusual case history in this arena, due to the recent arrival of the widely distributed Indo‐Pacific sergeant (Abudefduf vaigiensis), which is hybridizing with the endemic congener, A. abdominalis. Surveys of mtDNA and three nuclear loci across Hawai'i (N = 396, Abudefduf abdominalis and N = 314, A. vaigiensis) reveal that hybridization is significantly higher in the human‐perturbed southeast archipelago (19.8%), tapering off to 5.9% in the pristine northwest archipelago. While densities of the two species varied throughout Hawai'i, hybridization was highest in regions with similar species densities, contradicting the generalization that the rarity of one species promotes interspecific mating. Our finding of later generation hybrids throughout the archipelago invokes the possibility of genetic swamping of the endemic species. Exaptation, an adaptation with unintended consequences, may explain these findings: the endemic species has transient yellow coloration during reproduction, whereas the introduced species has yellow coloration continuously as adults, in effect a permanent signal of reproductive receptivity. Haplotype diversity is higher in Hawaiian A. vaigiensis than in our samples from the native range, indicating large‐scale colonization almost certainly facilitated by the historically recent surge of marine debris. In this chain of events, marine debris promotes colonization, exaptation promotes hybridization, and introgression invokes the possible collapse of an endemic species.     

Genome‐wide transcript analysis of early maize leaf development reveals gene cohorts associated with the differentiation of C4 Kranz anatomy

Photosynthesis underpins the viability of most ecosystems, with C₄ plants that exhibit ‘Kranz’ anatomy being the most efficient primary producers. Kranz anatomy is characterized by closely spaced veins that are encircled by two morphologically distinct photosynthetic cell types. Although Kranz anatomy evolved multiple times, the underlying genetic mechanisms remain largely elusive, with only the maize scarecrow gene so far implicated in Kranz patterning. To provide a broader insight into the regulation of Kranz differentiation, we performed a genome‐wide comparative analysis of developmental trajectories in Kranz (foliar leaf blade) and non‐Kranz (husk leaf sheath) leaves of the C₄ plant maize. Using profile classification of gene expression in early leaf primordia, we identified cohorts of genes associated with procambium initiation and vascular patterning. In addition, we used supervised classification criteria inferred from anatomical and developmental analyses of five developmental stages to identify candidate regulators of cell‐type specification. Our analysis supports the suggestion that Kranz anatomy is patterned, at least in part, by a SCARECROW/SHORTROOT regulatory network, and suggests likely components of that network. Furthermore, the data imply a role for additional pathways in the development of Kranz leaves.     

Microsatellite‐based analysis of genetic structure and gene flow of Mythimna separata (Walker) (Lepidoptera: Noctuidae) in China

The oriental armyworm, Mythimna separata, is a serious agricultural pest in China. Seasonal and roundtrip migration has recently led to sudden, localized outbreaks and crop losses. To evaluate genetic differentiation between populations in eastern and western China and elucidate gene flow, the genetic structure of 20 natural populations from nine provinces was examined using seven microsatellite markers. The results indicated high genetic diversity. However, little to moderate (0 < F_ST < 0.15) genetic differentiation was detected, and there was no correlation between genetic distance and geographical distance. Bayesian clustering analysis identified three groups whereas discriminant analysis of principal components identified ten clusters that were considered as two clear‐cut clusters and one admixed group. Gene flow occurred frequently in most population pairs, and an asymmetrical migration rate was detected in several pairwise population comparisons. The bottleneck test showed that few populations had experienced recent bottlenecks. Correspondingly, large‐scale and long‐distance migration of M. separata has caused low genetic differentiation and frequent gene exchange. Our findings are important for studying genetic evolution and help to improve predictions of M. separata outbreaks in China.     

Structure and stability of genetic variance–covariance matrices: A Bayesian sparse factor analysis of transcriptional variation in the three‐spined stickleback

The genetic variance–covariance matrix ( G ) is a quantity of central importance in evolutionary biology due to its influence on the rate and direction of multivariate evolution. However, the predictive power of empirically estimated G ‐matrices is limited for two reasons. First, phenotypes are high‐dimensional, whereas traditional statistical methods are tuned to estimate and analyse low‐dimensional matrices. Second, the stability of G to environmental effects and over time remains poorly understood. Using Bayesian sparse factor analysis (BSFG) designed to estimate high‐dimensional G ‐matrices, we analysed levels variation and covariation in 10,527 expressed genes in a large (n = 563) half‐sib breeding design of three‐spined sticklebacks subject to two temperature treatments. We found significant differences in the structure of G between the treatments: heritabilities and evolvabilities were higher in the warm than in the low‐temperature treatment, suggesting more and faster opportunity to evolve in warm (stressful) conditions. Furthermore, comparison of G and its phenotypic equivalent P revealed the latter is a poor substitute of the former. Most strikingly, the results suggest that the expected impact of G on evolvability—as well as the similarity among G ‐matrices—may depend strongly on the number of traits included into analyses. In our results, the inclusion of only few traits in the analyses leads to underestimation in the differences between the G ‐matrices and their predicted impacts on evolution. While the results highlight the challenges involved in estimating G , they also illustrate that by enabling the estimation of large G ‐matrices, the BSFG method can improve predicted evolutionary responses to selection.     

Regional and local patterns of genetic variation and structure in yellow‐necked mice ‐ the roles of geographic distance,population abundance,and winter severity

The goal of this study, conducted in seven large woodlands and three areas with small woodlots in northeastern Poland in 2004–2008, was to infer genetic structure in yellow‐necked mouse Apodemus flavicollis population and to evaluate the roles of environmental and population ecology variables in shaping the spatial pattern of genetic variation using 768 samples genotyped at 13 microsatellite loci. Genetic variation was very high in all studied regions. The primal genetic subdivision was observed between the northern and the southern parts of the study area, which harbored two major clusters and the intermediate area of highly admixed individuals. The probability of assignment of individual mice to the northern cluster increased significantly with lower temperatures of January and July and declined in regions with higher proportion of deciduous and mixed forests. Despite the detected structure, genetic differentiation among regions was very low. Fine‐scale structure was shaped by the population density, whereas higher level structure was mainly shaped by geographic distance. Genetic similarity indices were highly influenced by mouse abundance (which positively correlated with the share of deciduous forests in the studied regions) and exhibited the greatest change between 0 and 1 km in the forests, 0 and 5 km in small woodlots. Isolation by distance pattern, calculated among regions, was highly significant but such relationship between genetic and geographic distance was much weaker, and held the linearity at very fine scale (~1.5 km), when analyses were conducted at individual level.     

A new method for studying population genetics of cyst nematodes based on Pool‐Seq and genomewide allele frequency analysis

Cyst nematodes are important agricultural pests responsible for billions of dollars of losses each year. Plant resistance is the most effective management tool, but it requires a close monitoring of population genetics. Current technologies for pathotyping and genotyping cyst nematodes are time‐consuming, expensive and imprecise. In this study, we capitalized on the reproduction mode of cyst nematodes to develop a simple population genetic analysis pipeline based on genotyping‐by‐sequencing and Pool‐Seq. This method yielded thousands of SNPs and allowed us to study the relationships between populations of different origins or pathotypes. Validation of the method on well‐characterized populations also demonstrated that it was a powerful and accurate tool for population genetics. The genomewide allele frequencies of 23 populations of golden nematode, from nine countries and representing the five known pathotypes, were compared. A clear separation of the pathotypes and fine genetic relationships between and among global populations were obtained using this method. In addition to being powerful, this tool has proven to be very time‐ and cost‐efficient and could be applied to other cyst nematode species.     

Geographic and host‐mediated population genetic structure in a cestode parasite of the three‐spined stickleback

Comparative studies of genetic diversity and population structure can shed light on the ecological and evolutionary factors that influence host–parasite interactions. Here we examined whether geography, time and genetic variation in Alaskan three‐spined stickleback (Gasterosteus aculeatus Linneaus) hosts shape the population genetic structure of the diphyllobothridean cestode parasite Schistocephalus solidus (Müller, 1776). Host lineages and haplotypes were identified by sequencing the mitochondrial cytochrome b gene, and host population structure was assessed by Bayesian clustering analysis of allelic variation at 11 microsatellite loci. Parasite population structure was characterized according to allelic variation at eight microsatellite loci. Mantel tests and canonical redundancy analysis were conducted to evaluate the proportion of parasite genetic variation attributable to time and geography vs. host lineage, haplotype, and genotypic cluster. Host and parasite population structure were largely discordant across the study area, probably reflecting differences in gene flow, environmental influences external to the host, and genomic admixture among host lineages. We found that geography explained the greatest proportion of parasite genetic variation, but that variation also reflects time, host lineage, and host haplotype. Associations with host haplotypes suggest that one parasite genotypic cluster exhibits a narrower host range, predominantly infecting the most common host haplotypes, whereas the other parasite cluster infects all haplotypes equally, including rare haplotypes. Although experimental infection trials might prove otherwise, distributional differences in hosts preferentially infected by S. solidus could underlie the observed pattern of population structure.