Host-pathogen coevolution promotes the evolution of general, broad-spectrum resistance and reduces foreign pathogen spillover risk.

Genetic variation for disease resistance within host populations can strongly impact the spread of endemic pathogens. In plants, recent work has shown that within-population variation in resistance can also affect the transmission of foreign spillover pathogens if that resistance is general. However, most hosts also possess specific resistance mechanisms that provide strong defenses against coevolved endemic pathogens. Here we use a modeling approach to ask how antagonistic coevolution between hosts and their endemic pathogen at the specific resistance locus can affect the frequency of general resistance, and therefore a host's vulnerability to foreign pathogens. We develop a two-locus model with variable recombination that incorporates both general resistance (effective against all pathogens) and specific resistance (effective against endemic pathogens only). With coevolution, when pathogens can evolve to evade specific resistance, we find that the regions where general resistance can evolve are greatly expanded, decreasing the risk of foreign pathogen invasion. Furthermore, coevolution greatly expands the conditions that maintain polymorphisms at both resistance loci, thereby driving greater genetic diversity within host populations. This genetic diversity often leads to positive correlations between host resistance to foreign and endemic pathogens, similar to those observed in natural populations. However, if resistance loci become linked, the resistance correlations can shift to negative. If we include a third linkage-modifying locus in our model, we find that selection often favors complete linkage. Our model demonstrates how coevolutionary dynamics with an endemic pathogen can mold the resistance structure of host populations in ways that affect its susceptibility to foreign pathogen spillovers, and that the nature of these outcomes depends on resistance costs, as well as the degree of linkage between resistance genes.

Quantitative disease resistance in wild Silene vulgaris to its endemic pathogen Microbotryum silenes-inflatae.

The evolution of disease resistances is an expected feature of plant-pathogen systems, but whether the genetics of this trait most often produces qualitative or quantitative phenotypic variation is a significant gap in our understanding of natural populations. These two forms of resistance variation are often associated with differences in number of underlying loci, the specificities of host-pathogen coevolution, as well as contrasting mechanisms of preventing or slowing the infection process. Anther-smut disease is a commonly studied model for disease of wild species, where infection has severe fitness impacts, and prior studies have suggested resistance variation in several host species. However, because the outcome of exposing the individual host to this pathogen is binary (healthy or diseased), resistance has been previously measured at the family level, as the proportion of siblings that become diseased. This leaves uncertain whether among-family variation reflects contrasting ratios of segregating discrete phenotypes or continuous trait variation among individuals. In the host Silene vulgaris, plants were replicated by vegetative propagation in order to quantify the infection rates of the individual genotype with the endemic anther-smut pathogen, Microbotryum silenes-inflatae. The variance among field-collected families for disease resistance was significant, while there was unimodal continuous variation in resistance among genotypes. Using crosses between genotypes within ranked resistance quartiles, the offspring infection rate was predicted by the parental resistance values. While the potential remains in this system for resistance genes having major effects, as there were suggestions of such qualitative resistance in a prior study, here the quantitative disease resistance to the endemic anther-smut pathogen is indicated for S. vulgaris. The variation in natural populations and strong heritability of the trait, combined with severe fitness consequences of anther-smut disease, suggests that resistance in these host populations is highly capable of responding to disease-induced selection.

Dynamics of transposable element accumulation in the non-recombining regions of mating-type chromosomes in anther-smut fungi.

In the absence of recombination, the number of transposable elements (TEs) increases due to less efficient selection, but the dynamics of such TE accumulations are not well characterized. Leveraging a dataset of 21 independent events of recombination cessation of different ages in mating-type chromosomes of Microbotryum fungi, we show that TEs rapidly accumulated in regions lacking recombination, but that TE content reached a plateau at ca. 50% of occupied base pairs by 1.5 million years following recombination suppression. The same TE superfamilies have expanded in independently evolved non-recombining regions, in particular rolling-circle replication elements (Helitrons). Long-terminal repeat (LTR) retrotransposons of the Copia and Ty3 superfamilies also expanded, through transposition bursts (distinguished from gene conversion based on LTR divergence), with both non-recombining regions and autosomes affected, suggesting that non-recombining regions constitute TE reservoirs. This study improves our knowledge of genome evolution by showing that TEs can accumulate through bursts, following non-linear decelerating dynamics.

Host-pathogen coevolution promotes the evolution of general, broad-spectrum resistance and reduces foreign pathogen spillover risk.

Genetic variation for disease resistance within host populations can strongly impact the spread of endemic pathogens. In plants, recent work has shown that within-population variation in resistance can also affect the transmission of foreign spillover pathogens if that resistance is general. However, most hosts also possess specific resistance mechanisms that provide strong defenses against coevolved endemic pathogens. Here we use a modeling approach to ask how antagonistic coevolution between hosts and their endemic pathogen at the specific resistance locus can affect the frequency of general resistance, and therefore a host's vulnerability to foreign pathogens. We develop a two-locus model with variable recombination that incorporates both general (resistance to all pathogens) and specific (resistance to endemic pathogens only). We find that introducing coevolution into our model greatly expands the regions where general resistance can evolve, decreasing the risk of foreign pathogen invasion. Furthermore, coevolution greatly expands which conditions maintain polymorphisms at both resistance loci, thereby driving greater genetic diversity within host populations. This genetic diversity often leads to positive correlations between host resistance to foreign and endemic pathogens, similar to those observed in natural populations. However, if resistance loci become linked, the resistance correlations can shift to negative. If we include a third, linkage modifying locus into our model, we find that selection often favors complete linkage. Our model demonstrates how coevolutionary dynamics with an endemic pathogen can mold the resistance structure of host populations in ways that affect its susceptibility to foreign pathogen spillovers, and that the nature of these outcomes depends on resistance costs, as well as the degree of linkage between resistance genes.

Specific resistance prevents the evolution of general resistance and facilitates disease emergence.

Host-shifts, where pathogens jump from an ancestral host to a novel host, can be facilitated or impeded by standing variation in disease resistance, but only if resistance provides broad-spectrum general resistance against multiple pathogen species. Host resistance comes in many forms and includes both general resistance, as well as specific resistance, which may only be effective against a single pathogen species or even genotype. However, most evolutionary models consider only one of these forms of resistance, and we have less understanding of how these two forms of resistance evolve in tandem. Here, we develop a model that allows for the joint evolution of specific and general resistance and asks if the evolution of specific resistance drives a decrease in the evolution of general resistance. We also explore how these evolutionary outcomes affect the risk of foreign pathogen invasion and persistence. We show that in the presence of a single endemic pathogen, the two forms of resistance are strongly exclusionary. Critically, we find that specific resistance polymorphisms can prevent the evolution of general resistance, facilitating the invasion of foreign pathogens. We also show that specific resistance polymorphisms are a necessary condition for the successful establishment of foreign pathogens following invasion, as they prevent the exclusion of the foreign pathogen by the more transmissible endemic pathogen. Our results demonstrate the importance of considering the joint evolution of multiple forms of resistance when evaluating a population's susceptibility to foreign pathogens.

Multimodal pathogen transmission as a limiting factor in host distribution.

Theoretical models suggest that infectious diseases could play a substantial role in determining the spatial extent of host species, but few studies have collected the empirical data required to test this hypothesis. Pathogens that sterilize their hosts or spread through frequency-dependent transmission could have especially strong effects on the limits of species' distributions because diseased hosts that are sterilized but not killed may continue to produce infectious stages and frequency-dependent transmission mechanisms are effective even at very low population densities. We collected spatial pathogen prevalence data and population abundance data for alpine carnations infected by the sterilizing pathogen Microbotryum dianthorum, a parasite that is spread through both frequency-dependent (vector-borne) and density-dependent (aerial spore transmission) mechanisms. Our 13-year study reveals rapid declines in population abundance without a compensatory decrease in pathogen prevalence. We apply a stochastic, spatial model of parasite spread that accommodates spatial habitat heterogeneity to investigate how the population dynamics depend on multimodal (frequency-dependent and density-dependent) transmission. We found that the observed rate of population decline could plausibly be explained by multimodal transmission, but is unlikely to be explained by either frequency-dependent or density-dependent mechanisms alone. Multimodal pathogen transmission rates high enough to explain the observed decline predicted that eventual local extinction of the host species is highly likely. Our results add to a growing body of literature showing how multimodal transmission can constrain species distributions in nature.

Can disease resistance evolve independently at different ages? Genetic variation in age-dependent resistance to disease in three wild plant species.

Juveniles are typically less resistant (more susceptible) to infectious disease than adults, and this difference in susceptibility can help fuel the spread of pathogens in age-structured populations. However, evolutionary explanations for this variation in resistance across age remain to be tested.One hypothesis is that natural selection has optimized resistance to peak at ages where disease exposure is greatest. A central assumption of this hypothesis is that hosts have the capacity to evolve resistance independently at different ages. This would mean that host populations have (a) standing genetic variation in resistance at both juvenile and adult stages, and (b) that this variation is not strongly correlated between age classes so that selection acting at one age does not produce a correlated response at the other age.Here we evaluated the capacity of three wild plant species (Silene latifolia, S. vulgaris and Dianthus pavonius) to evolve resistance to their anther-smut pathogens (Microbotryum fungi), independently at different ages. The pathogen is pollinator transmitted, and thus exposure risk is considered to be highest at the adult flowering stage.Within each species we grew families to different ages, inoculated individuals with anther smut, and evaluated the effects of age, family and their interaction on infection.In two of the plant species, S. latifolia and D. pavonius, resistance to smut at the juvenile stage was not correlated with resistance to smut at the adult stage. In all three species, we show there are significant age × family interaction effects, indicating that age specificity of resistance varies among the plant families. Synthesis. These results indicate that different mechanisms likely underlie resistance at juvenile and adult stages and support the hypothesis that resistance can evolve independently in response to differing selection pressures as hosts age. Taken together our results provide new insight into the structure of genetic variation in age-dependent resistance in three well-studied wild host-pathogen systems.

Onset and stepwise extensions of recombination suppression are common in mating-type chromosomes of Microbotryum anther-smut fungi.

Sex chromosomes and mating-type chromosomes can display large genomic regions without recombination. Recombination suppression often extended stepwise with time away from the sex- or mating-type-determining genes, generating evolutionary strata of differentiation between alternative sex or mating-type chromosomes. In anther-smut fungi of the Microbotryum genus, recombination suppression evolved repeatedly, linking the two mating-type loci and extended multiple times in regions distal to the mating-type genes. Here, we obtained high-quality genome assemblies of alternative mating types for four Microbotryum fungi. We found an additional event of independent chromosomal rearrangements bringing the two mating-type loci on the same chromosome followed by recombination suppression linking them. We also found, in a new clade analysed here, that recombination suppression between the two mating-type loci occurred in several steps, with first an ancestral recombination suppression between one of the mating-type locus and its centromere; later, completion of recombination suppression up to the second mating-type locus occurred independently in three species. The estimated dates of recombination suppression between the mating-type loci ranged from 0.15 to 3.58 million years ago. In total, this makes at least nine independent events of linkage between the mating-type loci across the Microbotryum genus. Several mating-type locus linkage events occurred through the same types of chromosomal rearrangements, where similar chromosome fissions at centromeres represent convergence in the genomic changes leading to the phenotypic convergence. These findings further highlight Microbotryum fungi as excellent models to study the evolution of recombination suppression.

Tempo of Degeneration Across Independently Evolved Nonrecombining Regions.

Recombination is beneficial over the long term, allowing more effective selection. Despite long-term advantages of recombination, local recombination suppression can evolve and lead to genomic degeneration, in particular on sex chromosomes. Here, we investigated the tempo of degeneration in nonrecombining regions, that is, the function curve for the accumulation of deleterious mutations over time, leveraging on 22 independent events of recombination suppression identified on mating-type chromosomes of anther-smut fungi, including newly identified ones. Using previously available and newly generated high-quality genome assemblies of alternative mating types of 13 Microbotryum species, we estimated degeneration levels in terms of accumulation of nonoptimal codons and nonsynonymous substitutions in nonrecombining regions. We found a reduced frequency of optimal codons in the nonrecombining regions compared with autosomes, that was not due to less frequent GC-biased gene conversion or lower ancestral expression levels compared with recombining regions. The frequency of optimal codons rapidly decreased following recombination suppression and reached an asymptote after ca. 3 Ma. The strength of purifying selection remained virtually constant at dN/dS = 0.55, that is, at an intermediate level between purifying selection and neutral evolution. Accordingly, nonsynonymous differences between mating-type chromosomes increased linearly with stratum age, at a rate of 0.015 per My. We thus develop a method for disentangling effects of reduced selection efficacy from GC-biased gene conversion in the evolution of codon usage and we quantify the tempo of degeneration in nonrecombining regions, which is important for our knowledge on genomic evolution and on the maintenance of regions without recombination.

Resistance Correlations Influence Infection by Foreign Pathogens.

AbstractReciprocal selection promotes the specificity of host-pathogen associations and resistance polymorphisms in response to disease. However, plants and animals also vary in response to pathogen species not previously encountered in nature, with potential effects on new disease emergence. Using anther smut disease, we show that resistance (measured as infection rates) to foreign pathogens can be correlated with standing variation in resistance to an endemic pathogen. In Silene vulgaris, genetic variation in resistance to its endemic anther smut pathogen correlated positively with resistance variation to an anther smut pathogen from another host, but the relationship was negative between anther smut and a necrotrophic pathogen. We present models describing the genetic basis for assessing resistance relationships between endemic and foreign pathogens and for quantifying infection probabilities on foreign pathogen introduction. We show that even when the foreign pathogen has a lower average infection ability than the endemic pathogen, infection outcomes are determined by the sign and strength of the regression of the host's genetic variation in infection rates by a foreign pathogen on variation in infection rates by an endemic pathogen as well as by resistance allele frequencies. Given that preinvasion equilibria of resistance are determined by factors including resistance costs, we show that protection against foreign pathogens afforded by positively correlated resistances can be lessened or even result in elevated infection risk at the population level, depending on local dynamics. Therefore, a pathogen's emergence potential could be influenced not only by its average infection rate but also by resistance variation resulting from prior selection imposed by endemic diseases.

From generalist to specialists: Variation in the host range and performance of anther-smut pathogens on Dianthus.

Determining the processes that drive the evolution of pathogen host range can inform our understanding of disease dynamics and the potential for host shifts. In natural populations, patterns of host range could be driven by genetically based differences in pathogen infectivity or ecological differences in host availability. In northwestern Italy, four reproductively isolated lineages of the fungal plant-pathogen Microbotryum have been shown to co-occur on several species in the genus Dianthus. We carried out cross-inoculation experiments to determine whether patterns of realized host range in these four lineages were driven by differences in infectivity and to test whether there was evidence of a trade-off between host range and within-host reproduction. We found strong concordance between field patterns of host range and pathogen infectivity on different Dianthus species using experimental inoculation, indicating that infection ability is a major driving force of host range. However, we found no evidence of a trade-off between the ability to infect a wider range of host species and spore production on a shared host.

Recombination suppression and evolutionary strata around mating-type loci in fungi: documenting patterns and understanding evolutionary and mechanistic causes.

Genomic regions determining sexual compatibility often display recombination suppression, as occurs in sex chromosomes, plant self-incompatibility loci and fungal mating-type loci. Regions lacking recombination can extend beyond the genes determining sexes or mating types, by several successive steps of recombination suppression. Here we review the evidence for recombination suppression around mating-type loci in fungi, sometimes encompassing vast regions of the mating-type chromosomes. The suppression of recombination at mating-type loci in fungi has long been recognized and maintains the multiallelic combinations required for correct compatibility determination. We review more recent evidence for expansions of recombination suppression beyond mating-type genes in fungi ('evolutionary strata'), which have been little studied and may be more pervasive than commonly thought. We discuss testable hypotheses for the ultimate (evolutionary) and proximate (mechanistic) causes for such expansions of recombination suppression, including (1) antagonistic selection, (2) association of additional functions to mating-type, such as uniparental mitochondria inheritance, (3) accumulation in the margin of nonrecombining regions of various factors, including deleterious mutations or transposable elements resulting from relaxed selection, or neutral rearrangements resulting from genetic drift. The study of recombination suppression in fungi could thus contribute to our understanding of recombination suppression expansion across a broader range of organisms.

Meiotic recombination in the offspring of Microbotryum hybrids and its impact on pathogenicity.

BACKGROUND: Hybridization is a central mechanism in evolution, producing new species or introducing important genetic variation into existing species. In plant-pathogenic fungi, adaptation and specialization to exploit a host species are key determinants of evolutionary success. Here, we performed experimental crosses between the two pathogenic Microbotryum species, M. lychnidis-dioicae and M. silenes-acaulis that are specialized to different hosts. The resulting offspring were analyzed on phenotypic and genomic levels to describe genomic characteristics of hybrid offspring and genetic factors likely involved in host-specialization. RESULTS: Genomic analyses of interspecific fungal hybrids revealed that individuals were most viable if the majority of loci were inherited from one species. Interestingly, species-specific loci were strictly controlled by the species' origin of the mating type locus. Moreover we detected signs of crossing over and chromosome duplications in the genomes of the analyzed hybrids. In Microbotryum, mitochondrial DNA was found to be uniparentally inherited from the a2 mating type. Genome comparison revealed that most gene families are shared and the majority of genes are conserved between the two species, indicating very similar biological features, including infection and pathogenicity processes. Moreover, we detected 211 candidate genes that were retained under host-driven selection of backcrossed lines. These genes and might therefore either play a crucial role in host specialization or be linked to genes that are essential for specialization. CONCLUSION: The combination of genome analyses with experimental selection and hybridization is a promising way to investigate host-pathogen interactions. This study manifests genetic factors of host specialization that are required for successful biotrophic infection of the post-zygotic stage, but also demonstrates the strong influence of intra-genomic conflicts or instabilities on the viability of hybrids in the haploid host-independent stage.

Differential Gene Expression between Fungal Mating Types Is Associated with Sequence Degeneration.

Degenerative mutations in non-recombining regions, such as in sex chromosomes, may lead to differential expression between alleles if mutations occur stochastically in one or the other allele. Reduced allelic expression due to degeneration has indeed been suggested to occur in various sex-chromosome systems. However, whether an association occurs between specific signatures of degeneration and differential expression between alleles has not been extensively tested, and sexual antagonism can also cause differential expression on sex chromosomes. The anther-smut fungus Microbotryum lychnidis-dioicae is ideal for testing associations between specific degenerative signatures and differential expression because 1) there are multiple evolutionary strata on the mating-type chromosomes, reflecting successive recombination suppression linked to mating-type loci; 2) separate haploid cultures of opposite mating types help identify differential expression between alleles; and 3) there is no sexual antagonism as a confounding factor accounting for differential expression. We found that differentially expressed genes were enriched in the four oldest evolutionary strata compared with other genomic compartments, and that, within compartments, several signatures of sequence degeneration were greater for differentially expressed than non-differentially expressed genes. Two particular degenerative signatures were significantly associated with lower expression levels within differentially expressed allele pairs: upstream insertion of transposable elements and mutations truncating the protein length. Other degenerative mutations associated with differential expression included nonsynonymous substitutions and altered intron or GC content. The association between differential expression and allele degeneration is relevant for a broad range of taxa where mating compatibility or sex is determined by genes located in large regions where recombination is suppressed.

Exploring density- and frequency-dependent interactions experimentally: An r program for generating hexagonal fan designs.

Species interactions and diversity are strongly impacted by local processes, with both the density of a focal species and its frequency in the community having an impact on its growth, survival and fecundity. Yet, studies that attempt to control for variation in both frequency and density have traditionally required a large number of replicates.Hexagonal fan designs can include a range of both densities and frequencies in a single plot, providing large economies in space and material for studying local interactions such as competition and disease transmission. However, in practice such experiments can be difficult to plan and implement.This study presents an R program whereby the user can rapidly view a variety of designs and determine the configurations that work best with their space and material constraints. Simple instructions for implementing the fan in any design setting are also provided.We illustrate the implementation of a simple form of the hexagonal fan design in a field experiment to assess the impact of host density on pollinator movement and disease transmission.

Higher Gene Flow in Sex-Related Chromosomes than in Autosomes during Fungal Divergence.

Nonrecombining sex chromosomes are widely found to be more differentiated than autosomes among closely related species, due to smaller effective population size and/or to a disproportionally large-X effect in reproductive isolation. Although fungal mating-type chromosomes can also display large nonrecombining regions, their levels of differentiation compared with autosomes have been little studied. Anther-smut fungi from the Microbotryum genus are castrating pathogens of Caryophyllaceae plants with largely nonrecombining mating-type chromosomes. Using whole genome sequences of 40 fungal strains, we quantified genetic differentiation among strains isolated from the geographically overlapping North American species and subspecies of Silene virginica and S. caroliniana. We inferred that gene flow likely occurred at the early stages of divergence and then completely stopped. We identified large autosomal genomic regions with chromosomal inversions, with higher genetic divergence than the rest of the genomes and highly enriched in selective sweeps, supporting a role of rearrangements in preventing gene flow in genomic regions involved in ecological divergence. Unexpectedly, the nonrecombining mating-type chromosomes showed lower divergence than autosomes due to higher gene flow, which may be promoted by adaptive introgressions of less degenerated mating-type chromosomes. The fact that both mating-type chromosomes are always heterozygous and nonrecombining may explain such patterns that oppose to those found for XY or ZW sex chromosomes. The specific features of mating-type chromosomes may also apply to the UV sex chromosomes determining sexes at the haploid stage in algae and bryophytes and may help test general hypotheses on the evolutionary specificities of sex-related chromosomes.

Understanding Adaptation, Coevolution, Host Specialization, and Mating System in Castrating Anther-Smut Fungi by Combining Population and Comparative Genomics.

Anther-smut fungi provide a powerful system to study host-pathogen specialization and coevolution, with hundreds of Microbotryum species specialized on diverse Caryophyllaceae plants, castrating their hosts through manipulation of the hosts' reproductive organs to facilitate disease transmission. Microbotryum fungi have exceptional genomic characteristics, including dimorphic mating-type chromosomes, that make this genus anexcellent model for studying the evolution of mating systems and their influence on population genetics structure and adaptive potential. Important insights into adaptation, coevolution, host specialization, and mating system evolution have been gained using anther-smut fungi, with new insights made possible by the recent advent of genomic approaches. We illustrate with Microbotryum case studies how using a combination of comparative genomics, population genomics, and transcriptomics approaches enables the integration of different evolutionary perspectives across different timescales. We also highlight current challenges and suggest future studies that will contribute to advancing our understanding of the mechanisms underlying adaptive processes in populations of fungal pathogens.

Sympatry and interference of divergent Microbotryum pathogen species.

The impact of infectious diseases in natural ecosystems is strongly influenced by the degree of pathogen specialization and by the local assemblies of potential host species. This study investigated anther-smut disease, caused by fungi in the genus Microbotryum, among natural populations of plants in the Caryophyllaceae. A broad geographic survey focused on sites of the disease on multiple host species in sympatry. Analysis of molecular identities for the pathogens revealed that sympatric disease was most often due to co-occurrence of distinct, host-specific anther-smut fungi, rather than localized cross-species disease transmission. Flowers from sympatric populations showed that the Microbotryum spores were frequently moved between host species. Experimental inoculations to simulate cross-species exposure to the pathogens in these plant communities showed that the anther-smut pathogen was less able to cause disease on its regular host when following exposure of the plants to incompatible pathogens from another host species. These results indicate that multi-host/multi-pathogen communities are common in this system and they involve a previously hidden mechanism of interference between Microbotryum fungi, which likely affects both pathogen and host distributions.