Variation in phenotype and virulence among enzootic and panzootic amphibian chytrid lineages

The Global Panzootic Lineage of the fungus Batrachochytrium dendrobatidis (Bd-GPL) is threatening amphibians worldwide. In contrast, four lineages (Bd-Brazil, Bd-CH, Bd-Cape, and Bd-Korea) that diverged early in the history of Bd have not yet been directly linked to amphibian declines. Bd likely evolves in response to strong selective pressure imposed by hosts and the environment, leading to differences among pathogen phenotypes and genotypes that may directly affect virulence. Here, we report on variation in phenotype, genotype, and virulence of Bd-Brazil and Bd-GPL. Specifically, we (i) used a controlled infection experiment to compare virulence between one Bd-Brazil and three Bd-GPL isolates on a North American amphibian host (Lithobates sylvaticus), (ii) tested for relative phenotypic and genotypic differentiation among Bd isolates from Brazil, and (iii) tested for possible correlations between environmental variables and Bd phenotypes. We found substantial variation in virulence among Bd-GPL isolates and found that our Bd-Brazil isolate showed virulence comparable to an average North American Bd-GPL. North American hosts infected with a Bd-GPL isolate from Panama did not show significant mortality. Bd phenotypes varied significantly across sampling locations; these phenotypes were neither spatially clustered nor correlated with any environmental variables. Additionally, we found a surprising lack of correlation between genotypic divergence and zoospore and zoosporangium sizes in our sample. Although Bd-Brazil was less virulent infecting L. sylvaticus than one Bd-GPL isolate, this endemic lineage still caused ∼50% mortality in our experimental North American hosts. This indicates that Bd-Brazil has the potential to kill amphibians if introduced to naïve wild populations. Our findings underscore that characterizing virulence of multiple Bd isolates and lineages is important for understanding the evolutionary history and diversity of Bd.     

Population structure and diversity in sexual and asexual populations of the pathogenic fungus Melampsora lini

Many pathogens undergo both sexual and asexual reproduction to varying degrees, yet the ecological, genetic and evolutionary consequences of different reproductive strategies remain poorly understood. Here we investigate the population genetic structure of wild populations of the plant pathogen Melampsora lini on its host Linum marginale , using amplified fragment length polymorphism (AFLP) markers, two genes underlying pathogen virulence, and phenotypic variation in virulence. In Australia, M. lini occurs as two genetically and geographically divergent lineages (AA and AB), one of which is completely asexual (AB), and the other able to reproduce both clonally and sexually (AA). To quantify the genetic and evolutionary consequences of these different life histories, we sampled five populations in each of two biogeographical regions. Analysis of AFLP data obtained for 275 isolates revealed largely disjunct geographical distributions for the two different lineages, low genetic diversity within lineages, and strong genetic structure among populations within each region. We also detected significant divergence among populations for both Avr genes and virulence phenotypes, although generally these values were lower than those obtained with AFLP markers. Furthermore, isolates belonging to lineage AA collectively harboured significantly higher genotypic and phenotypic diversity than lineage AB isolates. Together these results illustrate the important roles of reproductive modes and geographical structure in the generation and maintenance of virulence diversity in populations of M. lini .     

四川省马铃薯晚疫病菌群体表型和遗传变异的分析

致病疫霉Phytophthora infestans引起的晚疫病是马铃薯的一种毁灭性病害。为了对四川省近年马铃薯晚疫病菌进行系统多样性分析,本研究从表现型和基因型两个方面鉴定了四川省马铃薯晚疫病菌的群体多态性。交配型、甲霜灵敏感性和生理小种的鉴定结果发现A1交配型菌株只有2份,A2型有29份,自育型12份;甲霜灵敏感菌株3株,中抗菌株22株,高抗菌株18株;共测定出11个生理小种,其中全毒力小种发生频率最高,占供试菌株的25.58%。基因型鉴定结果表明,43份材料共发现了10种SSR基因型,5个SSR标记在该种群上共产生了16个SSR位点,每个标记平均有3.2个位点,多态性信息含量平均值为0.46,SSR4是多态性最丰富的标记。本试验中A2菌株和自育型的基因型SSR相同,其中全毒力小种是当地的优势致病菌,该研究为今后四川马铃薯晚疫病的有效防控提供理论依据。     

Do pathogens become more virulent as they spread? Evidence from the amphibian declines in Central America

The virulence of a pathogen can vary strongly through time. While cyclical variation in virulence is regularly observed, directional shifts in virulence are less commonly observed and are typically associated with decreasing virulence of biological control agents through coevolution. It is increasingly appreciated, however, that spatial effects can lead to evolutionary trajectories that differ from standard expectations. One such possibility is that, as a pathogen spreads through a naive host population, its virulence increases on the invasion front. In Central America, there is compelling evidence for the recent spread of pathogenic Batrachochytrium dendrobatidis (Bd) and for its strong impact on amphibian populations. Here, we re-examine data on Bd prevalence and amphibian population decline across 13 sites from southern Mexico through Central America, and show that, in the initial phases of the Bd invasion, amphibian population decline lagged approximately 9 years behind the arrival of the pathogen, but that this lag diminished markedly over time. In total, our analysis suggests an increase in Bd virulence as it spread southwards, a pattern consistent with rapid evolution of increased virulence on Bd''s invading front. The impact of Bd on amphibians might therefore be driven by rapid evolution in addition to more proximate environmental drivers.     

Batrachochytrium dendrobatidis: requirement for further isolate collection and archiving

The fungal pathogen Batrachochytrium dendrobatidis (Bd) causes the disease chytridiomycosis, which is lethal to many species of amphibians worldwide. Many studies have investigated the epidemiology of chytridiomycosis in amphibian populations, but few have considered possible host-pathogen coevolution. More specifically, investigations focused on the evolution of Bd, and the link with Bd virulence, are needed. Such studies, which may be important for conservation management of amphibians, depend on access to Bd isolates. Here we provide a summary of known Bd isolates that have been collected and archived in various locations around the world. Of 257 Bd isolates, we found that 53% originate from ranids in the United States. In many cases, detailed information on isolate origin is unavailable, and it is unknown how many isolates are cryo-archived. We suggest the creation of a centralized database of isolate information, and we urge researchers and managers to isolate and archive Bd to facilitate future research on chytridiomycosis.     

Variation in tolerance and virulence in the chestnut blight fungus-hypovirus interaction

Chestnut blight, caused by the fungus Cryphonectria parasitica, has been effectively controlled with double-stranded RNA hypoviruses in Europe for over 40 years. The marked reduction in the virulence of C. parasitica by hypoviruses is a phenomenon known as hypovirulence. This virus-fungus pathosystem has become a model system for the study of biological control of fungi with viruses. We studied variation in tolerance to hypoviruses in fungal hosts and variation in virulence among virus isolates from a local population in Italy. Tolerance is defined as the relative fitness of a fungal individual when infected with hypoviruses (compared to being uninfected); virulence is defined for each hypovirus as the reduction in fitness of fungal hosts relative to virus-free hosts. Six hypovirus-infected isolates of C. parasitica were sampled from the population, and each hypovirus was transferred into six hypovirus-free recipient isolates. The resulting 36 hypovirus-fungus combinations were used to estimate genetic variation in tolerance to hypoviruses, in hypovirus virulence, and in virus-fungus interactions. Four phenotypes were evaluated for each virus-fungus combination to estimate relative fitness: (i) sporulation, i.e., the number of asexual spores (conidia) produced; (ii) canker area on field-inoculated chestnut trees, (iii) vertical transmission of hypoviruses into conidia, and (iv) conidial germination. Two-way analysis of variance (ANOVA) revealed significant interactions (P < 0.001) between viruses and fungal isolates for sporulation and canker area but not for conidial germination or transmission. One-way ANOVA among hypoviruses (within each fungal isolate) and among fungal isolates (within each hypovirus) revealed significant genetic variation (P < 0.01) in hypovirus virulence and fungal tolerance within several fungal isolates, and hypoviruses, respectively. These interactions and the significant genetic variation in several fitness characters indicate the potential for future evolution of these characters. However, biological control is unlikely to break down due to evolution of tolerance to hypoviruses in the fungus because the magnitudes of tolerance and interactions were relatively small.     

Genetic diversity in Hemileia vastatrix based on RAPD markers

Random amplified polymorphic DNA (RAPD) was used to assess the genetic structure of Hemileia vastatrix populations. Forty-five rust isolates with different virulence spectra and from different hosts and geographical regions were analyzed. Out of 45 bands, generated with three RAPD primers, 35 (78%) were polymorphic and scored as molecular markers. Cluster analysis exhibits unstructured variability of this pathogen with regard to physiological race, geographical origin or host. The genotypic diversity (H') inferred from Shannon's index was higher than gene diversity (Ht), suggesting that diversity is distributed among clonal lineages. Estimates of gene diversity in Africa and Asia populations were higher in total (Ht) as compared to within population diversity (Hs). Genetic differentiation was considerable among coffee rust isolates from Africa (Gst = 0.865) and Asia (Gst = 0.768) but not among isolates from South America (Gst = 0.266). We concluded that genetic diversity in H. vastatrix was moderately low and that the genetic differentiation among populations shows that asexual reproduction is likely to play an important role in the population biology of this fungus. This should be taken into account for the development of breeding programs.     

Variation in Tolerance and Virulence in the Chestnut Blight Fungus-Hypovirus Interaction

Chestnut blight, caused by the fungus Cryphonectria parasitica, has been effectively controlled with double-stranded RNA hypoviruses in Europe for over 40 years. The marked reduction in the virulence of C. parasitica by hypoviruses is a phenomenon known as hypovirulence. This virus-fungus pathosystem has become a model system for the study of biological control of fungi with viruses. We studied variation in tolerance to hypoviruses in fungal hosts and variation in virulence among virus isolates from a local population in Italy. Tolerance is defined as the relative fitness of a fungal individual when infected with hypoviruses (compared to being uninfected); virulence is defined for each hypovirus as the reduction in fitness of fungal hosts relative to virus-free hosts. Six hypovirus-infected isolates of C. parasitica were sampled from the population, and each hypovirus was transferred into six hypovirus-free recipient isolates. The resulting 36 hypovirus-fungus combinations were used to estimate genetic variation in tolerance to hypoviruses, in hypovirus virulence, and in virus-fungus interactions. Four phenotypes were evaluated for each virus-fungus combination to estimate relative fitness: (i) sporulation, i.e., the number of asexual spores (conidia) produced; (ii) canker area on field-inoculated chestnut trees, (iii) vertical transmission of hypoviruses into conidia, and (iv) conidial germination. Two-way analysis of variance (ANOVA) revealed significant interactions (P < 0.001) between viruses and fungal isolates for sporulation and canker area but not for conidial germination or transmission. One-way ANOVA among hypoviruses (within each fungal isolate) and among fungal isolates (within each hypovirus) revealed significant genetic variation (P < 0.01) in hypovirus virulence and fungal tolerance within several fungal isolates, and hypoviruses, respectively. These interactions and the significant genetic variation in several fitness characters indicate the potential for future evolution of these characters. However, biological control is unlikely to break down due to evolution of tolerance to hypoviruses in the fungus because the magnitudes of tolerance and interactions were relatively small.     

Allelic exchange and site-directed mutagenesis probe the contribution of ActA amino-acid variability to phosphorylation and virulence-associated phenotypes among Listeria monocytogenes strains

To test the hypothesis that actA allelic variation contributes to virulence differences among Listeria monocytogenes strains, cell-to-cell spread and intracellular ActA phosphorylation patterns were characterized for 14 wild-type isolates and selected isogenic mutants. Our data show that (i) while actA allelic variation is not responsible for enhanced cell-to-cell spread observed in epidemic clone I strains, actA allelic variation may contribute to reduced plaque size observed in some isolates, (ii) actA sequence alone determines phosphorylation-dependent ActA banding patterns, and (iii) sequence variation at the positively selected ActA residue 498 does not contribute to ActA phosphorylation patterns or to differences in cell-to-cell spread.     

Phenotypic variation and genomic variation in insect virulence traits reveal patterns of intraspecific diversity in a locust-specific fungal pathogen

Intraspecific pathogen diversity is crucial for understanding the evolution and maintenance of adaptation in host–pathogen interactions. Traits associated with virulence are often a significant source of variation directly impacted by local selection pressures. The specialist fungal entomopathogen, Metarhizium acridum, has been widely implemented as a biological control agent of locust pests in tropical regions of the world. However, few studies have accounted for natural intraspecific phenotypic and genetic variation. Here, we examine the diversity of nine isolates of M. acridum spanning the known geographic distribution, in terms of (1) virulence towards two locust species, (2) growth rates on three diverse nutrient sources, and (3) comparative genomics to uncover genomic variability. Significant variability in patterns of virulence and growth was shown among the isolates, suggesting intraspecific ecological specialization. Different patterns of virulence were shown between the two locust species, indicative of potential host preference. Additionally, a high level of diversity among M. acridum isolates was observed, revealing increased variation in subtilisin-like proteases from the Pr1 family. These results culminate in the first in-depth analysis regarding multiple facets of natural variation in M. acridum, offering opportunities to understand critical evolutionary drivers of intraspecific diversity in pathogens.     

Genotypic and phenotypic diversity of a baculovirus population within an individual insect host

It is becoming increasingly apparent that many pathogen populations, including those of insects, show high levels of genotypic variation. Baculoviruses are known to be highly variable, with isolates collected from the same species in different geographical locations frequently showing genetic variation and differences in their biology. More recent studies at smaller scales have also shown that virus DNA profiles from individual larvae can show polymorphisms within and between populations of the same species. Here, we investigate the genotypic and phenotypic variation of an insect baculovirus infection within a single insect host. Twenty four genotypically distinct nucleopolyhedrovirus (NPV) variants were isolated from an individual pine beauty moth, Panolis flammea, caterpillar by in vivo cloning techniques. No variant appeared to be dominant in the population. The PaflNPV variants have been mapped using three restriction endonucleases and shown to contain three hypervariable regions containing insertions of 70-750 bp. Comparison of seven of these variants in an alternative host, Mamestra brassicae, demonstrated that the variants differed significantly in both pathogenicity and speed of kill. The generation and maintenance of pathogen heterogeneity are discussed.     

Be flexible and adapt easily—The great role of plasticity relative to genetic variation for aggressiveness of Fusarium culmorum isolates

The phenotypic variation in an array of pathogen isolates in natural environments can be partitioned into genotypic variation and environmental plasticity. The present study uses a mixed-model approach to partition the relative contribution of both factors among isolates of Fusarium culmorum from natural field populations in various environments. Twenty-eight and 38 isolates from an international collection were phenotyped for aggressiveness and deoxynivalenol (DON) accumulation across two locations during the years 2015 and 2016, respectively, on four winter type cereals as hosts: bread wheat, durum wheat, triticale and rye, thus providing 16 environments. Aggressiveness, measured as Fusarium head blight (FHB) severity, was assessed by visually rating the symptoms of all isolates on infected hosts, and for 10 isolates, additionally the mycotoxin deoxynivalenol (DON) was measured in the grain after harvest. Despite significant genotypic variation among the isolates, the interactions with years and locations explained the largest proportion of variance which disentangled the overwhelming role of plasticity. Host-by-isolate interaction was not significant and no significant (p < .001) change in the ranking of isolates from one host to another was detected. As the main factor of plasticity was isolate-by-year interaction, this implies that seasonal changes might be an important evolutionary driver in F. culmorum populations.     

Genotypic diversity and migration patterns of Phytophthora infestans in the Nordic countries

In this study we investigated the genotypic diversity and the migration patterns of Phytophthora infestans in the Nordic countries. Isolates of P. infestans from outbreaks in 43 fields sampled in 2008 were collected using stratified sampling with country, field, and disease foci as the different strata. Microsatellites were used as markers to determine the genotypic variation in the sampled material. The results show a high genotypic variation of P. infestans in the Nordic countries with most of the genotypes found only once among the collected isolates. The major part of the genotypic variation was observed within the fields, with low differentiation between the fields. The observed low association of alleles among loci is consistent with frequent sexual reproduction of P. infestans in the Nordic countries. Coalescence analyses did not support a single common population for the four countries, thus indicating some degree of geographic differentiation. The analyses of migration patterns showed differing levels of gene flow among the Nordic countries. No correlation between migration rates and geographical distance could be seen. This could be explained by different degrees of genetic similarity between the pathogen populations in the different countries.     

Comparison of genotypic and phenotypic techniques for assessing the variability of the fungus Epicoccum nigrum

AIMS: The diversity within a collection of worldwide isolates of Epicoccum nigrum has been studied using several phenotypic approaches. In addition, the abilities of phenotypic and genotypic techniques for the differentiation of a set of isolates are compared. METHODS AND RESULTS: The methodology used include the study of isozymes (acetyl esterase and alkaline phosphatase), HPLC profile of metabolites and antibiotic activities against a panel of bacteria, yeasts and filamentous fungi, and cytotoxicity against three mammalian cell lines. Two procedures for assessing the relationships within a collection of isolates, using a combination of the techniques, were evaluated, comparing the advantages and disadvantages of each method. CONCLUSIONS: The results showed that each individual technique allows differentiation of the isolates studied to some degree and that the information provided by each technique could be considered as complementary. Genotypic techniques were more powerful than the phenotypic ones to discriminate among the strains. SIGNIFICANCE AND IMPACT OF THE STUDY: This work evaluates the predictive value of several phenotypic techniques on a collection of fungal isolates, and compares the results obtained with genotypic techniques performed on the same strains.     

Chromosomal Copy Number Variation,Selection and Uneven Rates of Recombination Reveal Cryptic Genome Diversity Linked to Pathogenicity

Pathogenic fungi constitute a growing threat to both plant and animal species on a global scale. Despite a clonal mode of reproduction dominating the population genetic structure of many fungi, putatively asexual species are known to adapt rapidly when confronted by efforts to control their growth and transmission. However, the mechanisms by which adaptive diversity is generated across a clonal background are often poorly understood. We sequenced a global panel of the emergent amphibian pathogen, Batrachochytrium dendrobatidis (Bd), to high depth and characterized rapidly changing features of its genome that we believe hold the key to the worldwide success of this organism. Our analyses show three processes that contribute to the generation of de novo diversity. Firstly, we show that the majority of wild isolates manifest chromosomal copy number variation that changes over short timescales. Secondly, we show that cryptic recombination occurs within all lineages of Bd, leading to large regions of the genome being in linkage equilibrium, and is preferentially associated with classes of genes of known importance for virulence in other pathosystems. Finally, we show that these classes of genes are under directional selection, and that this has predominantly targeted the Global Panzootic Lineage (BdGPL). Our analyses show that Bd manifests an unusually dynamic genome that may have been shaped by its association with the amphibian host. The rates of variation that we document likely explain the high levels of phenotypic variability that have been reported for Bd, and suggests that the dynamic genome of this pathogen has contributed to its success across multiple biomes and host-species.     

Collagenase production and hemolytic activity related to 16S rRNA variability among Parvimonas micra oral isolates

Parvimonas micra are gram positive anaerobic cocci isolated from the oral cavity and frequently related to polymicrobial infections in humans. Despite reports about phenotypic differences, the genotypic variation of P. micra and its role in virulence are still not elucidated. The aim of this study was to determine the genotypic diversity of P. micra isolates obtained from the subgingival biofilm of subjects with different periodontal conditions and to correlate these findings with phenotypic traits. Three reference strains and 35 isolates of P. micra were genotyped by 16S rRNA PCR-RFLP and phenotypic traits such as collagenase production, elastolytic and hemolytic activities were evaluated. 16S rRNA PCR-RFLP showed that P. micra could be grouped into two main clusters: C1 and C2; cluster C1 harbored three genotypes (HG1259-like, HG1467-like and ICBMO583-like) while cluster C2 harbored two genotypes (ATCC33270-like and ICBMO36). A wide variability in collagenolytic activity intensities was observed among all isolates, while elastolytic activity was detected in only two isolates. There was an association between hemolytic activity in rabbit erythrocytes and cluster C2. There was an association between hemolytic activity in rabbit erythrocytes and cluster C1. Although these data suggest a possible association between P. micra genetic diversity and their pathogenic potential, further investigations are needed to confirm this hypothesis.     

Molecular analyses suggest monospecificity of the genus Sarcoptes (Acari: Sarcoptidae).

To clarify the taxonomic status of mites of the genus Sarcoptes, the second internal transcribed spacer (ITS-2) of the rRNA gene, as well as phenotypic characters, were investigated in 23 isolates from nine host species in four continents. Phenotypic differences among isolates were observed, but the range of variation within each isolate precluded the differentiation of individual mites. Genotypically, there was no delimitation between distinct genotypic groups and no correlation with host species or geographic origin was evident. These results support the conspecificity of the mites investigated and confirm the view that the genus Sarcoptes consists of a single, heterogenous species.     

THE IMPACT OF SPATIAL SCALE AND HABITAT CONFIGURATION ON PATTERNS OF TRAIT VARIATION AND LOCAL ADAPTATION IN A WILD PLANT PARASITE

Theory indicates that spatial scale and habitat configuration are fundamental for coevolutionary dynamics and how diversity is maintained in host–pathogen interactions. Yet, we lack empirical data to translate the theory to natural host–parasite systems. In this study, we conduct a multiscale cross‐inoculation study using the specialist wild plant pathogen Podosphaera plantaginis on its host plant Plantago lanceolata. We apply the same sampling scheme to a region with highly fragmented (Åland) and continuous (Saaremaa) host populations. Although theory predicts higher parasite virulence in continuous regions, we did not detect differences in traits conferring virulence among the regions. Patterns of adaptation were highly scale dependent. We detected parasite maladaptation among regions, and among populations separated by intermediate distances (6.0–40.0 km) within the fragmented region. In contrast, parasite performance did not vary significantly according to host origin in the continuous landscape. For both regions, differentiation among populations was much larger for genetic variation than for phenotypic variation, indicating balancing selection maintaining phenotypic variation within populations. Our findings illustrate the critical role of spatial scale and habitat configuration in driving host–parasite coevolution. The absence of more aggressive strains in the continuous landscape, in contrast to theoretical predictions, has major implications for long‐term decision making in conservation, agriculture, and public health.     

Quantifying the genetic component of phenotypic variation in unpedigreed wild plants: tailoring genomic scan for within-population use

The study of adaptive genetic variation in natural populations is central to evolutionary biology. Quantitative genetics methods, however, are hardly applicable to long-lived organisms, and current knowledge on adaptive genetic variation in wild plants mostly refers to annuals and short-lived perennials. Studies on long-lived species are essential to explore possible life-history correlates of genetic variation, selection, and trait heritability. In this paper, we propose a method based on molecular markers to quantify the genetic basis of individual phenotypic differences in wild plants under natural conditions. Rather than focusing on inferring individual relatedness to estimate the heritability of phenotypic traits, we directly estimate the proportion of observed phenotypic variance that is statistically accounted for by genotypic differences between individuals. This is achieved by (i) identifying loci that are correlated across individuals with the phenotypic trait of interest by means of an amplified fragment length polymorphism (AFLP)-based explorative genomic scan, and (ii) fitting multiple regression and linear random effect models to estimate the effects of genotype, environment and genotype × environment on phenotypes. We apply this method to estimate genotypic and environmental effects on cumulative maternal fecundity in a wild population of the long-lived Viola cazorlensis monitored for 20 years. Results show that between 56–63% (depending on estimation method) of phenotypic variance in fecundity is accounted for by genotypic differences in 11 AFLP loci that are significantly related to fecundity. Genotype × environment effects accounted for 38% of fecundity variance, which may help to explain the unexpectedly high levels of genetic variance for fecundity found.