Clinal resistance structure and pathogen local adaptation in a serpentine flax-flax rust interaction

Because disease resistance is a hallmark signature of pathogen-mediated selection pressure on hosts, studies of resistance structure (the spatial distribution of disease resistance genes among conspecific host populations) can provide valuable insights into the influence of pathogens on host evolution and spatial variation in the magnitude of their effects. To date few studies of wild plant-pathogen interactions have characterized resistance structure by sampling across the host's biogeographic range, and only a handful have paired such investigations with studies of disease levels under natural conditions. I used a greenhouse cross-inoculation experiment to characterize genetic resistance of 16 populations of California dwarf flax (Hesperolinon californicum) to attack by multiple samples of the rust fungus Melampsora lini. I documented a latitudinal cline in resistance structure, manifest across the host's biogeographic range, which mirrored almost identically a cline in infection prevalence documented through field surveys of disease in study populations. These results provide empirical evidence for clinal patterns of antagonistic selection pressure, demonstrate that such patterns can be manifest across broad biogeographic scales, and suggest that rates of disease prevalence in wild plant populations may be tightly linked to the distribution of host resistance genes. Tests for local adaptation of the fungus revealed evidence of the phenomenon (significantly greater infection in sympatric plant-fungal pairings) as well as the potential for substantial bias to be introduced into statistical analyses by spatial patterns of host resistance structure.     

Population genetic structure of a California endemic branchiopod,Branchinecta sandiegonensis

Branchinecta sandiegonensis (Crustacea: Anostraca) is a narrow range endemic fairy shrimp discontinuously distributed in ephemeral pools on coastal mesas in San Diego County, USA. Ten populations across the range of the species were subjected to allozyme analysis for eleven loci. The species exhibits low variability (P₉₅ =9.1–45.5) and one third of the loci tested did not conform to Hardy-Weinberg equilibrium expectations. The species also exhibited a high degree of genetic differentiation between populations. F _ST values (fixation index) for most pairs of populations were above 0.25 (0.036–0.889).Low genetic variability and high genetic structure may result from low gene flow and founder effects due to habitat fragmentation and the lack of potential vectors for cyst dispersal. The unpredictable rainfall of the region also creates potential for variable population sizes which could affect structure and variability. This revised version was published online in August 2006 with corrections to the Cover Date.     

Population structure of the invasive forest pathogen Hymenoscyphus pseudoalbidus

Understanding the genetic diversity and structure of invasive pathogens in source and in introduced areas is crucial to the revelation of hidden biological features of an organism, to the reconstruction of the course of invasions and to the establishment of effective control measures. Hymenoscyphus pseudoalbidus (anamorph: Chalara fraxinea) is an invasive and highly destructive fungal pathogen found on common ash Fraxinus excelsior in Europe and is native to East Asia. To gain insights into its dispersal mechanisms and history of invasion, we used microsatellite markers and characterized the genetic structure and diversity of H. pseudoalbidus populations at three spatial levels: (i) between Europe and Japan, (ii) in Europe and (iii) at the epidemic's front in Switzerland. Phylogenetic and network analysis demonstrated that individuals from both regions are conspecific. However, populations from Japan harboured a higher genetic diversity and were genetically differentiated from European ones. No evident population structure was found among the 1208 European strains using Bayesian and multivariate clustering analysis. Only the distribution of genetic diversity in space, pairwise population differentiation (G_ST) and the spatial analysis of principal components revealed a faint geographical pattern around Europe. A significant allele deficiency in most European populations pointed to a recent genetic bottleneck, whereas no pattern of isolation by distance was found. Our data suggest that H. pseudoalbidus was introduced just once by at least two individuals. The potential source region of H. pseudoalbidus is vast, and further investigations are required for a more accurate localization of the source population.     

Population structure of loggerhead shrikes in the California Channel Islands

The loggerhead shrike (Lanius ludovicianus), a songbird that hunts like a small raptor, maintains breeding populations on seven of the eight California Channel Islands. One of the two subspecies, L. l. anthonyi, was described as having breeding populations on six of the islands while a second subspecies, L. l. mearnsi, was described as being endemic to San Clemente Island. Previous genetic studies have demonstrated that the San Clemente Island loggerhead shrike is well differentiated genetically from both L. l. anthonyi and mainland populations, despite the fact that birds from outside the population are regular visitors to the island. Those studies, however, did not include a comparison between San Clemente Island shrikes and the breeding population on Santa Catalina Island, the closest island to San Clemente. Here we use mitochondrial control region sequences and nuclear microsatellites to investigate the population structure of loggerhead shrikes in the Channel Islands. We confirm the genetic distinctiveness of the San Clemente Island loggerhead shrike and, using Bayesian clustering analysis, demonstrate the presence and infer the source of the nonbreeding visitors. Our results indicate that Channel Island loggerhead shrikes comprise three distinct genetic clusters that inhabit: (i) San Clemente Island, (ii) Santa Catalina Island and (iii) the Northern Channel Islands and nearby mainland; they do not support a recent suggestion that all Channel Island loggerhead shrikes should be managed as a single entity.     

Host dispersal shapes the population structure of a tick‐borne bacterial pathogen

Birds are hosts for several zoonotic pathogens. Because of their high mobility, especially of longdistance migrants, birds can disperse these pathogens, affecting their distribution and phylogeography. We focused on Borrelia burgdorferi sensu lato, which includes the causative agents of Lyme borreliosis, as an example for tick‐borne pathogens, to address the role of birds as propagation hosts of zoonotic agents at a large geographical scale. We collected ticks from passerine birds in 11 European countries. B. burgdorferi s.l. prevalence in Ixodes spp. was 37% and increased with latitude. The fieldfare Turdus pilaris and the blackbird T. merula carried ticks with the highest Borrelia prevalence (92 and 58%, respectively), whereas robin Erithacus rubecula ticks were the least infected (3.8%). Borrelia garinii was the most prevalent genospecies (61%), followed by B. valaisiana (24%), B. afzelii (9%), B. turdi (5%) and B. lusitaniae (0.5%). A novel Borrelia genospecies “Candidatus Borrelia aligera” was also detected. Multilocus sequence typing (MLST) analysis of B. garinii isolates together with the global collection of B. garinii genotypes obtained from the Borrelia MLST public database revealed that: (a) there was little overlap among genotypes from different continents, (b) there was no geographical structuring within Europe, and (c) there was no evident association pattern detectable among B. garinii genotypes from ticks feeding on birds, questing ticks or human isolates. These findings strengthen the hypothesis that the population structure and evolutionary biology of tick‐borne pathogens are shaped by their host associations and the movement patterns of these hosts.     

Population structure and diversity of an invasive pine needle pathogen reflects anthropogenic activity

Dothistroma septosporum is a haploid fungal pathogen that causes a serious needle blight disease of pines, particularly as an invasive alien species on Pinus radiata in the Southern Hemisphere. During the course of the last two decades, the pathogen has also incited unexpected epidemics on native and non‐native pine hosts in the Northern Hemisphere. Although the biology and ecology of the pathogen has been well documented, there is a distinct lack of knowledge regarding its movement or genetic diversity in many of the countries where it is found. In this study we determined the global population diversity and structure of 458 isolates of D. septosporum from 14 countries on six continents using microsatellite markers. Populations of the pathogen in the Northern Hemisphere, where pines are native, displayed high genetic diversities and included both mating types. Most of the populations from Europe showed evidence for random mating, little population differentiation and gene flow between countries. Populations in North America (USA) and Asia (Bhutan) were genetically distinct but migration between these continents and Europe was evident. In the Southern Hemisphere, the population structure and diversity of D. septosporum reflected the anthropogenic history of the introduction and establishment of plantation forestry, particularly with Pinus radiata. Three introductory lineages in the Southern Hemisphere were observed. Countries in Africa, that have had the longest history of pine introductions, displayed the greatest diversity in the pathogen population, indicating multiple introductions. More recent introductions have occurred separately in South America and Australasia where the pathogen population is currently reproducing clonally due to the presence of only one mating type.     

Population genetic diversity and structure of two rare vernal pool grasses in central California

Vernal pool ecosystems are declining throughout California, with only 10% of historic habitat remaining. This has endangered many specialist endemic plant species, leaving extant populations fragmented, isolated, and threatened or endangered. Recovery plans for the increasing number of endangered vernal pool species require information on their genetic and ecological status to guide conservation and restoration efforts. Federally threatened Neostapfia colusana (Colusa grass) and federally endangered Tuctoria greenei (Greene’s tuctoria) are two endemic vernal pool grasses of high conservation concern in central California. Remaining populations are highly fragmented due to range-wide habitat destruction. Using five polymorphic microsatellite markers for each species, we performed genetic surveys of 240 individuals from eight vernal pools for N. colusana, and 317 individuals from 13 vernal pools for T. greenei. We detected high within-population genetic diversity for both species, with average allelic diversities of 24 alleles/locus (mean Hobs = 0.68, mean Hexp = 0.71) for N. colusana, and 19 alleles/locus (mean Hobs = 0.77, and mean Hexp = 0.79) for T. greenei. Bayesian clustering and AMOVA indicated two genetically distinct population groups for N. colusana (Fst = 0.268, P < 0.0001), and three for T. greenei (Fst = 0.11, P < 0.0001). We found very slight temporal genetic structure at one N. colusana (Fst = 0.013, P < 0.05) and two T. greenei (Fst = 0.015, Fst = 0.018, P < 0.05) pools. These estimates of population genetic diversity and structure are critical measures for both species that will help inform recovery management actions.     

Perfect and near-perfect adaptation in a model of bacterial chemotaxis

The signaling apparatus mediating bacterial chemotaxis can adapt to a wide range of persistent external stimuli. In many cases, the bacterial activity returns to its prestimulus level exactly, and this perfect adaptability is robust against variations in various chemotaxis protein concentrations. We model the bacterial chemotaxis signaling pathway, from ligand binding to CheY phosphorylation. By solving the steady-state equations of the model analytically, we derive a full set of conditions for the system to achieve perfect adaptation. The conditions related to the phosphorylation part of the pathway are discovered for the first time, while other conditions are generalizations of the ones found in previous works. Sensitivity of the perfect adaptation is evaluated by perturbing these conditions. We find that, even in the absence of some of the perfect adaptation conditions, adaptation can be achieved with near-perfect precision as a result of the separation of scales in both chemotaxis protein concentrations and reaction rates, or specific properties of the receptor distribution in different methylation states. Since near-perfect adaptation can be found in much larger regions of the parameter space than that defined by the perfect adaptation conditions, their existence is essential to understand robustness in bacterial chemotaxis.     

Population structure of California coyotes corresponds to habitat-specific breaks and illuminates species history

Little is known about the relationship between animal movements and the emergent structure of populations, especially for species occupying large continuous distributions. Some such mammals disperse disproportionately into habitat similar to their natal habitat, a behavioural bias that might be expected to lead to habitat-conforming genetic structure. We hypothesized that coyotes (Canis latrans) would exhibit such natal-biased dispersal, and used 13 microsatellite loci to test, correspondingly, whether genetic structure conformed to major habitat breaks. First, we used a model-based approach to assign coyote genotypes to distinct genetic clusters irrespective of geographical location. Visualization on a geographical information system revealed a strong concordance between the locations of cluster assignments and habitat bioregions, not explainable in terms of physical dispersal barriers or intervening low-quality habitat. Next, we used a multiple Mantel test, which controlled for effects of geographical distance (itself, marginally significant; P = 0.06), to statistically determine that genetic distance was indeed higher between than within bioregions (P < 0.001). Whereas previously published examples of landscape effects on gene flow have typically been explainable in terms of species-wide habitat affinities or dispersal barriers, our finding that genetic subdivisions were associated with unobstructed boundaries between contiguous habitats suggests a role for intraspecific variability in habitat affinities as a factor underlying genetic structure. In addition, our data combined with previously published data suggest a pattern of genetic isolation-by-distance throughout western North America, consistent with independent evidence that the western half of the coyote range predates European settlement.     

Population genetic correlates of declining transmission in a human pathogen

Pathogen control programs provide a valuable, but rarely exploited, opportunity to directly examine the relationship between population decline and population genetics. We investigated the impact of an ~12‐fold decline in transmission on the population genetics of Plasmodium falciparum infections (n = 1731) sampled from four clinics on the Thai–Burma border over 10 years and genotyped using 96 genome‐wide SNPs. The most striking associated genetic change was a reduction in the frequency of infections containing multiple parasite genotypes from 63% in 2001 to 14% in 2010 (P = 3 × 10^?15). Two measures of the clonal composition of populations (genotypic richness and the β‐parameter of the Pareto distribution) declined over time as more people were infected by parasites with identical multilocus genotypes, consistent with increased selfing and a reduction in the rate at which multilocus genotypes are broken apart by recombination. We predicted that the reduction in transmission, multiple clone carriage and outbreeding would be mirrored by an increased influence of genetic drift. However, geographical differentiation and expected heterozygosity remained stable across the sampling period. Furthermore, N_e estimates derived from allele frequencies fluctuation between years remained high (582 to ∞) and showed no downward trend. These results demonstrate how genetic data can compliment epidemiological assessments of infectious disease control programs. The temporal changes in a single declining population parallel to those seen in comparisons of parasite genetics in regions of differing endemicity, strongly supporting the notion that reduced opportunity for outbreeding is the key driver of these patterns.     

Population genetic structure, local adaptation, and conservation genetics of Kandelia obovata

Topographic changes during the Pleistocene glacial/interglacial cycles affected the distribution of coastline mangroves and influenced their population genetic structure. The submergence of the continental shelf off southeast China during the postglacial age caused coastline expansions and resulted in the colonization of mangroves. Here, we performed multilocus genome scans using amplified fragment-length polymorphisms to explore the effects of topography and natural selection in structuring Kandelia obovata populations. Long-term isolation by the Taiwan Strait since the end of the last glacial maximum, which obstructed gene flow, differentiated the Taiwanese and Chinese populations. Founders that colonized from both outlets of the Taiwan Strait were sourced from the northern South China Sea and the Ryukyus, thereby creating a melting pot in the Taiwan Strait. Inner-strait currents played roles as vectors for propagule dispersal among populations. Upon examination of the allele-frequency distributions of outlier loci, most negative outliers reflected the widespread polymorphisms shared by common ancestors. Furthermore, significant differentiation in the genetic components of positive outliers between this and other populations and the negative correlation with geographic distance suggested the presence of geography- or latitude-independent population divergence. Restored populations were compared with their sources and revealed biased sampling of nursery seedlings, which caused within-population substructures and reduced effective population sizes. This study indicated that multiple factors affect the population structure of the mangroves off southeast China.     

Genetic and geographical structure in grapevines from northwestern Spain

The grapevine germplasm bank from Estación de Viticultura e Enoloxía de Galicia contains 110 accessions collected from more than 100‐year‐old grapevines from northwestern Spain. We identified 60 different genotypes (among which ‘Pinot Noir’ and ‘Shyrah’ were included as references) with 32 polymorphic microsatellites (SSRs), 21 of which were unlinked SSRs and without null alleles, which were used to detect genetic structure. We found four reconstructed populations that were highly differentiated by analysis of molecular variance (18.3%): (i) western Galician group (‘Caíño Bravo’); (ii) eastern Galician group (‘Merenzao’); (iii) Castilian group (‘Garnacha’) and (iv) southern Spanish group (‘Moscatel’). The greatest genetic differentiation was found between the western Galician/eastern Galician groups and the southern Spanish/Castilian groups, with Fst > 0.13. Western and eastern Galician groups contained most of the ancient Galician cultivars (86%) and included ancestors from central Europe, such as ‘Pinot Noir’ and ‘Shyrah’ in the western Galician group and the introgressant ‘Merenzao’ (French ‘Trousseau’) in the eastern Galician group; the Castilian and southern Spain groups included traditional cultivars from central and southern Spanish and suggest a further introduction of those cultivars in northwestern plantations. Genetic structure corresponded with East‐to‐West geographical differentiation in the northern Iberian Peninsula, following the old path to Santiago de Compostela from France, with the southern Spanish and Castilian groups in the East (Castilian plateau and southern Spain), the eastern Galician group and the western Galician group in the West. Five new relationships were discovered: ‘Caíño Bravo’/‘Caíño Longo2’, ‘Albarín Tinto’/‘Torrontés1’, ‘Treixadura’/‘Treixadura Francesa’, ‘Moscatel’/‘Moscatel de Bago Miúdo’ and ‘Moscatel’/‘Moscatel de Hamburgo’, indicating that hybridisation and further selection by growers have been the origin of most of the genetic diversity found in northwestern Spain.     

Population genetic structure and ancestry of Oncorhynchus mykiss populations above and below dams in south-central California

Genetic analyses of coastal Oncorhynchus mykiss, commonly known as steelhead/rainbow trout, at the southern extreme of their geographic range in California are used to evaluate ancestry and genetic relationships of populations both above and below large dams. Juvenile fish from 20 locations and strains of rainbow trout commonly planted in reservoirs in the five study basins were evaluated at 24 microsatellite loci. Phylogeographic trees and analysis of molecular variance demonstrated that populations within a basin, both above and below dams, were generally each other’s closest relatives. Absence of hatchery fish or their progeny in the tributaries above dams indicates that they are not commonly spawning and that above-barrier fish are descended from coastal steelhead trapped at dam construction. Finally, no genetic basis was found for the division of populations from this region into two distinct biological groups, contrary to current classification under the US and California Endangered Species Acts.     

Population genetic structure and adaptation of malaria parasites on the edge of endemic distribution

To determine whether the major human malaria parasite Plasmodium falciparum exhibits fragmented population structure or local adaptation at the northern limit of its African distribution where the dry Sahel zone meets the Sahara, samples were collected from diverse locations within Mauritania over a range of ~1000 km. Microsatellite genotypes were obtained for 203 clinical infection samples from eight locations, and Illumina paired‐end sequences were obtained to yield high coverage genomewide single nucleotide polymorphism (SNP) data for 65 clinical infection samples from four locations. Most infections contained single parasite genotypes, reflecting low rates of transmission and superinfection locally, in contrast to the situation seen in population samples from countries further south. A minority of infections shared related or identical genotypes locally, indicating some repeated transmission of parasite clones without recombination. This caused some multilocus linkage disequilibrium and local divergence, but aside from the effect of repeated genotypes there was minimal differentiation between locations. Several chromosomal regions had elevated integrated haplotype scores (|iHS|) indicating recent selection, including those containing drug resistance genes. A genomewide F_ST scan comparison with previous sequence data from an area in West Africa with higher infection endemicity indicates that regional gene flow prevents genetic isolation, but revealed allele frequency differentiation at three drug resistance loci and an erythrocyte invasion ligand gene. Contrast of extended haplotype signatures revealed none to be unique to Mauritania. Discrete foci of infection on the edge of the Sahara are genetically highly connected to the wider continental parasite population, and local elimination would be difficult to achieve without very substantial reduction in malaria throughout the region.     

Secrets of actin-based motility revealed by a bacterial pathogen

Actin-based cell motility is a complex process involving a dynamic, self-organizing cellular system. Experimental problems initially limited our understanding of this type of motility, but the use of a model system derived from a bacterial pathogen has led to a breakthrough. Now, all the molecular components necessary for dynamic actin self-organization and motility have been identified, setting the stage for future mechanistic studies.