Genetic control of horizontal virus transmission in the chestnut blight fungus, Cryphonectria parasitica

Vegetative incompatibility in fungi has long been known to reduce the transmission of viruses between individuals, but the barrier to transmission is incomplete. In replicated laboratory assays, we showed conclusively that the transmission of viruses between individuals of the chestnut blight fungus Cryphonectria parasitica is controlled primarily by vegetative incompatibility (vic) genes. By replicating vic genotypes in independent fungal isolates, we quantified the effect of heteroallelism at each of six vic loci on virus transmission. Transmission occurs with 100% frequency when donor and recipient isolates have the same vic genotypes, but heteroallelism at one or more vic loci generally reduces virus transmission. Transmission was variable among single heteroallelic loci. At the extremes, heteroallelism at vic4 had no effect on virus transmission, but transmission occurred in only 21% of pairings that were heteroallelic at vic2. Intermediate frequencies of transmission were observed when vic3 and vic6 were heteroallelic (76 and 32%, respectively). When vic1, vic2, and vic7 were heteroallelic, the frequency of transmission depended on which alleles were present in the donor and the recipient. The effect of heteroallelism at two vic loci was mostly additive, although small but statistically significant interactions (epistasis) were observed in four pairs of vic loci. A logistic regression model was developed to predict the probability of virus transmission between vic genotypes. Heteroallelism at vic loci, asymmetry, and epistasis were the dominant factors controlling transmission, but host genetic background also was statistically significant, indicating that vic genes alone cannot explain all the variation in virus transmission. Predictions from the logistic regression model were highly correlated to independent transmission tests with field isolates. Our model can be used to estimate horizontal transmission rates as a function of host genetics in natural populations of C. parasitica.     

Mitotic stability and nuclear inheritance of integrated viral cDNA in engineered hypovirulent strains of the chestnut blight fungus.

Transmissible hypovirulence is a novel form of biological control in which virulence of a fungal pathogen is attenuated by an endogenous RNA virus. The feasibility of engineering hypovirulence was recently demonstrated by transformation of the chestnut blight fungus, Cryphonectria parasitica, with a full-length cDNA copy of a hypovirulence-associated viral RNA. Engineered hypovirulent transformants were found to contain both a chromsomally integrated cDNA copy of the viral genome and a resurrected cytoplasmically replicating double-stranded RNA form. We now report stable maintenance of integrated viral cDNA through repeated rounds of asexual sporulation and passages on host plant tissue. We also demonstrate stable nuclear inheritance of the integrated viral cDNA and resurrection of the cytoplasmic viral double-stranded RNA form in progeny resulting from the mating of an engineered hypovirulent C. parasitica strain and a vegetatively incompatible virulent strain. Mitotic stability of the viral cDNA ensures highly efficient transmission of the hypovirulence phenotype through conidia. Meiotic transmission, a mode not observed for natural hypovirulent strains, introduces virus into ascospore progeny representing a spectrum of vegetative compatibility groups, thereby circumventing barriers to anastomosis-mediated transmission imposed by the fungal vegetative incompatibility system. These transmission properties significantly enhance the potential of engineered hypovirulent C. parasitica strains as effective biocontrol agents.     

Approximating the coalescent with recombination

The coalescent with recombination describes the distribution of genealogical histories and resulting patterns of genetic variation in samples of DNA sequences from natural populations. However, using the model as the basis for inference is currently severely restricted by the computational challenge of estimating the likelihood. We discuss why the coalescent with recombination is so challenging to work with and explore whether simpler models, under which inference is more tractable, may prove useful for genealogy-based inference. We introduce a simplification of the coalescent process in which coalescence between lineages with no overlapping ancestral material is banned. The resulting process has a simple Markovian structure when generating genealogies sequentially along a sequence, yet has very similar properties to the full model, both in terms of describing patterns of genetic variation and as the basis for statistical inference.     

Spatial but not temporal co-divergence of a virus and its mammalian host

Co-divergence between host and parasites suggests that evolutionary processes act across similar spatial and temporal scales. Although there has been considerable work on the extent and correlates of co-divergence of RNA viruses and their mammalian hosts, relatively little is known about the extent to which virus evolution is determined by the phylogeographic history of host species. To test hypotheses related to co-divergence across a variety of spatial and temporal scales, we explored phylogenetic signatures in Andes virus (ANDV) sampled from Chile and its host rodent, Oligoryzomys longicaudatus. ANDV showed strong spatial subdivision, a phylogeographic pattern also recovered in the host using both spatial and genealogical approaches, and despite incomplete lineage sorting. Lineage structure in the virus seemed to be a response to current population dynamics in the host at the spatial scale of ecoregions. However, finer scale analyses revealed contrasting patterns of genetic structure across a latitudinal gradient. As predicted by their higher substitution rates, ANDV showed greater genealogical resolution than the rodent, with topological congruence influenced by the degree of lineage sorting within the host. However, despite these major differences in evolutionary dynamics, the geographic structure of host and virus converged across large spatial scales.     

Evidence for interspecies transmission of viruses in natural populations of filamentous fungi in the genus Cryphonectria

Interspecies transmission is a significant evolutionary event that has allowed a variety of pathogens to invade new host species. We investigated interspecies transmission of viruses between the chestnut blight fungus, Cryphonectria parasitica, and a sympatric unidentified Cryphonectria species in Japan. Two isolates of Cryphonectria sp. were found to contain Cryphonectria hypovirus 1 (CHV-1), which has been typically found in C. parasitica. Three lines of evidence support the hypothesis of interspecies transmission of CHV-1. First, host species occur sympatrically and therefore have the opportunity to come into physical contact. Second, we transmitted CHV-1 between species experimentally in the laboratory. Third, phylogenetic analysis of 476 bp of the ORF B region of CHV-1 showed that sequences from Cryphonectria sp. were more closely related to those from C. parasitica than to each other. Local geographical subdivision of virus sequences from both host species argues against the alternative hypothesis of independent evolution of CHV-1 since speciation of their hosts. Based on these findings, we rule out the hypotheses that CHV-1 diverged from viruses in a common ancestor of the hosts, or that ancestral polymorphisms in CHV-1 persisted in the two host taxa. Estimating the direction and frequency of interspecies transmission in nature will require more extensive samples of CHV-1 from both host species.     

Spatial-temporal stratifications in natural populations and how they affect understanding and estimation of effective population size

The concept of effective population size (N(e) ) is based on an elegantly simple idea which, however, rapidly becomes very complex when applied to most real-world situations. In natural populations, spatial and temporal stratifications create different classes of individuals with different vital rates, and this in turn affects (generally reduces) N(e) in complex ways. I consider how these natural stratifications influence our understanding of effective size and how to estimate it, and what the consequences are for conservation and management of natural populations. Important points that emerge include the following: 1 The relative influences of local vs metapopulation N(e) depend on a variety of factors, including the time frame of interest. 2 Levels of diversity in local populations are strongly influenced by even low levels of migration, so these measures are not reliable indicators of local N(e) . 3 For long-term effective size, obtaining a reliable estimate of mutation rate is the most important consideration; unless this is accomplished, estimates can be biased by orders of magnitude. 4 At least some estimators of contemporary N(e) appear to be robust to relatively high (approximately 10%) equilibrium levels of migration, so under many realistic scenarios they might yield reliable estimates of local N(e) . 5 Age structure probably has little effect on long-term estimators of N(e) but can strongly influence contemporary estimates. 6 More research is needed in several key areas: (i) to disentangle effects of selection and drift in metapopulations connected by intermediate levels of migration; (ii) to elucidate the relationship between N(b) (effective number of breeders per year) and N(e) per generation in age-structured populations; (iii) to perform rigorous sensitivity analyses of new likelihood and coalescent-based methods for estimating demographic and evolutionary histories.     

Flies on the move: an inherited virus mirrors Drosophila melanogaster's elusive ecology and demography

Vertically transmitted parasites rely on their host's reproduction for their transmission, leading to the evolutionary histories of both parties being intimately entwined. Parasites can thus serve as a population genetic magnifying glass for their host's demographic history. Here, we study the fruitfly Drosophila melanogaster's vertically transmitted sigma virus DMelSV. The virus has a high mutation rate and low effective population size, allowing us to reconstruct at a fine scale how the combined forces of the movement of flies and selection on the virus have shaped its migration patterns. We found that the virus is likely to have spread to Europe from Africa, mirroring the colonization route of Drosophila. The North American DMelSV population appears to be the result of a recent single immigration from Europe, invading together with its host in the late 19th century. Across Europe, DMelSV migration rates are low and populations are highly genetically structured, likely reflecting limited fly movement. Despite being intolerant of extreme cold, viral diversity suggests that fly populations can persist in harsh continental climates and that recolonization from the warmer south plays a minor role. In conclusion, studying DMelSV can provide insights into the poorly understood ecology of D. melanogaster, one of the best‐studied organisms in biology.     

Utility of DNA viruses for studying human host history: case study of JC virus

Microbial pathogens, and viruses in particular, can serve as important complements to traditional genetic markers when investigating the population histories of their human host. The range of mutation rates for DNA viruses suggests that DNA viruses can be useful markers of both recent and ancient events in their host histories. Here, we assess the utility of a well known DNA virus, JC virus (JCV), for investigating human history and demography. Using complete coding viral genomes, we confirm the phylogeographic structure of JCV in populations worldwide and provide coalescent estimates of its evolutionary rate under two alternative models of its history. Using these rate estimates, we compare Bayesian skyline plots of population size changes for JCV to those of its human host as estimated with coding mitochondrial genomes of the latter. These comparisons, when combined with other evidence including a log Bayes Factor model test, show that JCV is evolving rapidly and is therefore tracking the recent history of its human host. These results support the hypothesis that post-World War II societal changes are most likely responsible for the recent demographic patterns observed among different regional JCV populations. In sum, fast evolving DNA viruses, such as JCV, can complement RNA viruses to provide novel insights about the recent history and demography of their human host.     

Geographical patterns of adaptation within a species' range: interactions between drift and gene flow

We use individual-based stochastic simulations and analytical deterministic predictions to investigate the interaction between drift, natural selection and gene flow on the patterns of local adaptation across a fragmented species' range under clinally varying selection. Migration between populations follows a stepping-stone pattern and density decreases from the centre to the periphery of the range. Increased migration worsens gene swamping in small marginal populations but mitigates the effect of drift by replenishing genetic variance and helping purge deleterious mutations. Contrary to the deterministic prediction that increased connectivity within the range always inhibits local adaptation, simulations show that low intermediate migration rates improve fitness in marginal populations and attenuate fitness heterogeneity across the range. Such migration rates are optimal in that they maximize the total mean fitness at the scale of the range. Optimal migration rates increase with shallower environmental gradients, smaller marginal populations and higher mutation rates affecting fitness.     

PLEIOTROPY CAUSES LONG-TERM GENETIC CONSTRAINTS ON LIFE-HISTORY EVOLUTION IN BRASSICA RAPA

Fundamental, long-term genetic trade-offs constrain life-history evolution in wild crucifer populations. I studied patterns of genetic constraint in Brassica rapa by estimating genetic correlations among life-history components by quantitative genetic analyses among ten wild populations, and within four populations. Genetic correlations between age and size at first reproduction were always greater than +0.8 within and among all populations studied. Although quantitative genetics might provide insight about genetic constraints if genetic parameters remain approximately constant, little evidence has been available to determine the constancy of genetic correlations. I found strong and consistent estimates of genetic correlations between life-history components, which were very similar within four natural populations. Population differentiation also showed these same trade-offs, resulting from long-term genetic constraint. For some traits, evolutionary changes among populations were incompatible with a model of genetic drift. Historical patterns of natural selection were inferred from population differentiation, suggesting that correlated response to selection has caused some traits to evolve opposite to the direct forces of natural selection. Comparison with Arabidopsis suggests that these life-history trade-offs are caused by genes that regulate patterns of resource allocation to different components of fitness. Ecological and energetic models may correctly predict these trade-offs because there is little additive genetic variation for rates of resource acquisition, but resource allocation is genetically variable.     

Differential propagation of the metazoan parasite Myxobolus cerebralis by Limnodrilus hoffmeisteri, Ilyodrilus templetoni, and genetically distinct strains of Tubifex tubifex

Whirling disease, caused by the parasite Myxobolus cerebralis, has infected rainbow trout (Oncorhynchus mykiss) and other salmonid fish in the western United States, often with devastating results to native populations but without a discernible spatial pattern. The parasite develops in a complex 2-host system in which the aquatic oligochaete Tubifex tubifex is an obligate host. Because substantial differences in whirling disease severity in different areas of North America did not seem explainable by environmental factors or features of the parasite or its fish host, we sought to determine whether ecological or genetic variation within oligochaete host populations may be responsible. We found large differences in compatibility between the parasite and various laboratory strains of T. tubifex that were established from geographic regions with different whirling disease histories. Moreover, 2 closely related species of tubificids, Limnodrilus hoffmeisteri and Ilyodrilus templetoni, which occur naturally in mixed species assemblages with T. tubifex, were incompatible with M. cerebralis. Virulence of the parasite was directly correlated with the numbers of triactinomyxon spores that developed within each strain of T. tubifex. Thus, the level of virulence was directly related to the compatibility between the host strain and the parasite. Genetic analyses revealed relationships that were in agreement with the level of parasite production. Differences in compatibilities between oligochaetes and M. cerebralis may contribute to the spatial variance in the severity of the disease among salmonid populations.     

Genotype-environment interactions of spontaneous mutations for vegetative fitness in the human pathogenic fungus Cryptococcus neoformans

Spontaneous mutation is the ultimate source of all genetic variation. By interacting with environmental factors, genetic variation determines the phenotype and fitness of individuals in natural populations. However, except in a few model organisms, relatively little is known about the patterns of genotype-environment interactions of spontaneous mutations. Here I examine the rates of spontaneous mutation and the patterns of genotype-environment interaction of mutations affecting vegetative growth in the human fungal pathogen Cryptococcus neoformans. Eight mutation accumulation (MA) lines were established from a single clone on the nutrient-rich medium YEPD for each of two temperatures, 25 degrees and 37 degrees. Cells from generations 100, 200, 400, and 600 for each of the 16 MA lines were stored and assayed for vegetative growth rates under each of four conditions: (i) 25 degrees on SD (a synthetic dextrose minimal medium); (ii) 25 degrees on YEPD; (iii) 37 degrees on SD; and (iv) 37 degrees on YEPD. Both MA conditions and assay environments for vegetative growth showed significant influence on the estimates of genomic mutation rates, average effect per mutation, and mutational heritability. Significant genotype-environment interactions were detected among the newly accumulated spontaneous mutations. Overall, clones from MA lines maintained at 37 degrees showed less decline in vegetative fitness than those maintained at 25 degrees. The result suggests that a high-temperature environment might be very important for the maintenance of the ability to grow at a high temperature. Results from comparisons between clinical and environmental samples of C. neoformans were consistent with laboratory experimental population analyses. This study calls into question our long-standing view that warm-blooded mammals were only occasional and accidental hosts of this human fungal pathogen.     

Host heterogeneity in susceptibility and disease dynamics: tests of a mathematical model

Most mathematical models of disease assume that transmission is linearly dependent on the densities of host and pathogen. Recent data for animal diseases, however, have cast doubt on this assumption, without assessing the usefulness of alternative models. In this article, we use a combination of laboratory dose-response experiments, field transmission experiments, and observations of naturally occurring populations to show that virus transmission in gypsy moths is a nonlinear function of virus density, apparently because of heterogeneity among individual gypsy moth larvae in their susceptibility to the virus. Dose-response experiments showed that larvae from a laboratory colony of gypsy moths are substantially less heterogeneous in their susceptibility to the virus than are larvae from feral populations, and field experiments showed that there is a more strongly nonlinear relationship between transmission and virus density for feral larvae than for lab larvae. This nonlinearity in transmission changes the dynamics of the virus in natural populations so that a model incorporating host heterogeneity in susceptibility to the virus gives a much better fit to data on virus dynamics from large-scale field plots than does a classical model that ignores host heterogeneity. Our results suggest that heterogeneity among individuals has important effects on the dynamics of disease in insects at several spatial and temporal scales and that heterogeneity in susceptibility may be of general importance in the ecology of disease.     

Invasion history and demographic pattern of Cryphonectria hypovirus 1 across European populations of the chestnut blight fungus

We reconstructed the invasion history of the fungal virus Cryphonectria hypovirus 1 (CHV‐1) in Europe, which infects the chestnut blight fungus Cryphonectria parasitica. The pattern of virus evolution was inferred based on nucleotide sequence variation from isolates sampled across a wide area in Europe at different points in time. Phylogeny and time estimates suggested that CHV‐1 was introduced together with its fungal host to Europe and that it rapidly colonized the central range along the south facing slopes of the Alps and the north‐east facing slopes of the Dinaric Alps. These central populations were the source for two waves of simultaneous invasions toward the southern Balkans and Turkey, as indicated by migration rates. Our results showed that the evolutionary scenarios for CHV‐1 and C. parasitica were spatially congruent. As infection with CHV‐1 reduces the pathogenicity of C. parasitica toward the chestnut tree, CHV‐1 invasions of the newly established C. parasitica populations probably prevented the development of devastating chestnut blight epidemics in Europe. We propose that in this, and supposedly in other pathosystems, geographic, vegetation‐related, demographic, economic, and political factors may help explain the correlated invasion pattern of a parasite and its host.     

Virulence not only costs but also benefits the transmission of a fungal virus

Current theory suggests that cost-benefit relationships govern the evolution of parasite virulence. The cost of virulence is expected to be high for fungal viruses, which are obligate parasites and completely dependent on their hosts. The majority of fungal viruses infect their hosts without any apparent symptoms. Cryphonectria hypovirus 1 (CHV-1), in contrast, is virulent and debilitates its host, Cryphonectria parasitica. However, the virulence of CHV-1 is associated with high costs for virus transmission, such as an attenuated fungal growth and reduced production of the fungal spores spreading the virus. In this study, we tested the hypothesis that virulence may not only have costs but also benefits for transmitting CHV-1 across vegetative incompatibility barriers between fungi. We investigated viruses with low, medium, and high virulence, and determined their transmission rate per host-to-host contact (transmissibility). The average transmission rate across all combinations tested was 53% for the most virulent virus, 37% for the virus with intermediate virulence, and 20% for the virus with lowest virulence. These results showed that increased virulence was strongly correlated with increased transmissibility, potentially counterbalancing virulence costs. This association of virulence and transmissibility may explain why CHV-1 spread widely and evolved higher virulence than most other fungal viruses.     

Intra-population comparison of vegetative and floral trait heritabilities estimated from molecular markers in wild Aquilegia populations

Measuring heritable genetic variation is important for understanding patterns of trait evolution in wild populations, and yet studies of quantitative genetic parameters estimated directly in the field are limited by logistic constraints, such as the difficulties of inferring relatedness among individuals in the wild. Marker-based approaches have received attention because they can potentially be applied directly to wild populations. For long-lived, self-compatible plant species where pedigrees are inadequate, the regression-based method proposed by Ritland has the appeal of estimating heritabilities from marker-based estimates of relatedness. The method has been difficult to implement in some plant populations, however, because it requires significant variance in relatedness across the population. Here, we show that the method can be readily applied to compare the ability of different traits to respond to selection, within populations. For several taxa of the perennial herb genus Aquilegia, we estimated heritabilities of floral and vegetative traits and, combined with estimates of natural selection, compared the ability to respond to selection of both types of traits under current conditions. The intra-population comparisons showed that vegetative traits have a higher potential for evolution, because although they are as heritable as floral traits, selection on them is stronger. These patterns of potential evolution are consistent with macroevolutionary trends in the European lineage of the genus.     

Likelihood analysis of ongoing gene flow and historical association

Abstract.— We develop a Monte Carlo-based likelihood method for estimating migration rates and population divergence times from data at unlinked loci at which mutation rates are sufficiently low that, in the recent past, the effects of mutation can be ignored. The method is applicable to restriction fragment length polymorphisms (RFLPs) and single nucleotide polymorphisms (SNPs) sampled from a subdivided population. The method produces joint maximum-likelihood estimates of the migration rate and the time of population divergence, both scaled by population size, and provides a framework in which to test either for no ongoing gene flow or for population divergence in the distant past. We show the method performs well and provides reasonably accurate estimates of parameters even when the assumptions under which those estimates are obtained are not completely satisfied. Furthermore, we show that, provided that the number of polymorphic loci is sufficiently large, there is some power to distinguish between ongoing gene flow and historical association as causes of genetic similarity between pairs of populations.     

Controlling population evolution in the laboratory to evaluate methods of historical inference

Natural populations of known detailed past demographic history are extremely valuable to evaluate methods of historical inference, yet are extremely rare. As an alternative approach, we have generated multiple replicate microsatellite data sets from laboratory-cultured populations of a gonochoric free-living nematode, Caenorhabditis remanei, that were constrained to pre-defined demographic histories featuring different levels of migration among populations or bottleneck events of different magnitudes. These data sets were then used to evaluate the performances of two recently developed population genetics methods, BayesAss+, that estimates recent migration rates among populations, and Bottleneck, that detects the occurrence of recent bottlenecks. Migration rates inferred by BayesAss+ were generally over-estimates, although these were often included within the confidence interval. Analyses of data sets simulated in-silico, using a model mimicking the laboratory experiments, produced less biased estimates of the migration rates, and showed increased efficiency of the program when the number of loci and sampled genotypes per population was higher. In the replicates for which the pre-bottleneck laboratory-cultured populations did not significantly depart from a mutation/drift equilibrium, an important assumption of the program Bottleneck, only a portion of the bottleneck events were detected. This result was confirmed by in-silico simulations mirroring the laboratory bottleneck experiments. More generally, our study demonstrates the feasibility, and highlights some of the limits, of the approach that consists in generating molecular genetic data sets by controlling the evolution of laboratory-reared nematode populations, for the purpose of validating methods inferring population history.     

Clonal population structure of the chestnut blight fungus in expanding ranges in southeastern Europe

Expanding populations are often less genetically diverse at their margins than at the centre of a species' range. Established, older populations of the chestnut blight fungus, Cryphonectria parasitica, are more variable for vegetative compatibility (vc) types than in expanding populations in southeastern Europe where C. parasitica has colonized relatively recently. To test whether vc types represent clones, we genotyped 373 isolates of C. parasitica from southern Italy, Romania, Bulgaria, Macedonia, Greece and Turkey using 11 sequence-characterized amplified region (SCAR) markers. Ten SCAR loci and six vegetative incompatibility (vic) loci were polymorphic in these samples. These populations are clonal by all criteria tested: (i) among 373 isolates, we found only eight multilocus haplotypes, and the same haplotypes were found in multiple countries, sometimes separated in time by as much as 12 years; (ii) the number of haplotypes observed was significantly less than expected under random mating; (iii) populations are in linkage disequilibrium; (iv) the two sets of independent markers, SCARs and vc types, are highly correlated; and (v) sexual structures of C. parasitica were found only in Bulgaria and Romania. One mating type (MAT-1) was found in 98% of the isolates sampled. In contrast, a population in northern Italy, in the central part of the range in Europe, had 12 multilocus haplotypes among 19 isolates. The spread of a few clones could be the result either of founder effect and restricted migration, or these clones have greater fitness than others and spread because they are better adapted to conditions in southeastern Europe.     

The impact of parasitism on resource allocation in a fungal host: the case of Cryphonectria parasitica and its mycovirus,Cryphonectria Hypovirus 1

Parasites are known to profoundly affect resource allocation in their host. In order to investigate the effects of Cryphonectria Hypovirus 1 (CHV1) on the life‐history traits of its fungal host Cryphonectria parasitica, an infection matrix was completed with the cross‐infection of six fungal isolates by six different viruses. Mycelial growth, asexual sporulation, and spore size were measured in the 36 combinations, for which horizontal and vertical transmission of the viruses was also assessed. As expected by life‐history theory, a significant negative correlation was found between host somatic growth and asexual reproduction in virus‐free isolates. Interestingly this trade‐off was found to be positive in infected isolates, illustrating the profound changes in host resource allocation induced by CHV1 infection. A significant and positive relationship was also found in infected isolates between vertical transmission and somatic growth. This last relationship suggests that in this system, high levels of virulence could be detrimental to the vertical transmission of the parasite. Those results underscore the interest of studying host–parasite interaction within the life‐history theory framework, which might permit a more accurate understanding of the nature of the modifications triggered by parasite infection on host biology.