Migration patterns among global populations of the pathogenic fungus Mycosphaerella graminicola

DNA sequences from five nuclear loci and data from three microsatellites were collected from 360 isolates representing 14 globally distributed populations of the plant pathogenic fungus Mycosphaerella graminicola. Haplotype networks were constructed for the five sequence loci and population subdivision was assessed using Hudson's permutation test. Migration estimates were calculated using six regional populations for both the sequence and microsatellite loci. While subdivision was detected among the six regional populations, significant gene flow was indicated among some of the populations. The European and Israeli populations contributed the majority of historical immigrants to the New World. Migration estimates for microsatellite loci were used to infer more recent migration events among specific New World populations. We conclude that gene flow was an important factor in determining the demographic history of Mycosphaerella graminicola.     

Origin and domestication of the fungal wheat pathogen Mycosphaerella graminicola via sympatric speciation

The Fertile Crescent represents the center of origin and earliest known place of domestication for many cereal crops. During the transition from wild grasses to domesticated cereals, many host-specialized pathogen species are thought to have emerged. A sister population of the wheat-adapted pathogen Mycosphaerella graminicola was identified on wild grasses collected in northwest Iran. Isolates of this wild grass pathogen from 5 locations in Iran were compared with 123 M. graminicola isolates from the Middle East, Europe, and North America. DNA sequencing revealed a close phylogenetic relationship between the pathogen populations. To reconstruct the evolutionary history of M. graminicola, we sequenced 6 nuclear loci encompassing 464 polymorphic sites. Coalescence analyses indicated a relatively recent origin of M. graminicola, coinciding with the known domestication of wheat in the Fertile Crescent around 8,000-9,000 BC. The sympatric divergence of populations was accompanied by strong genetic differentiation. At the present time, no genetic exchange occurs between pathogen populations on wheat and wild grasses although we found evidence that gene flow may have occurred since genetic differentiation of the populations.     

Variation for neutral markers is correlated with variation for quantitative traits in the plant pathogenic fungus Mycosphaerella graminicola

We compared genetic variation and population differentiation at RFLP marker loci with seven quantitative characters including fungicide resistance, temperature sensitivity, pycnidial size, pycnidial density, colony size, percentage of leaves covered by pycnidia (PLACP) and percentage of leaves covered by lesions (PLACL) in Mycosphaerella graminicola populations sampled from four regions. Wide variation in population differentiation was found across the quantitative traits assayed. Fungicide resistance, temperature sensitivity, and PLACP displayed a significantly higher Q(ST) than G(ST), consistent with selection for local adaptation, while pycnidial size, pycnidial density and colony size displayed a lower or significantly lower Q(ST) than G(ST), consistent with constraining selection. There was not a statistical difference between Q(ST) and G(ST) in PLACL. We also found a positive and significant correlation between genetic variation in molecular marker loci and quantitative traits at the multitrait scale, suggesting that estimates of overall genetic variation for quantitative traits in M. graminicola could be derived from analysis of the molecular genetic markers.     

Genetic basis of adaptation in Arabidopsis thaliana: local adaptation at the seed dormancy QTL DOG1

Local adaptation provides an opportunity to study the genetic basis of adaptation and investigate the allelic architecture of adaptive genes. We study delay of germination 1 (DOG1), a gene controlling natural variation in seed dormancy in Arabidopsis thaliana and investigate evolution of dormancy in 41 populations distributed in four regions separated by natural barriers. Using F(ST) and Q(ST) comparisons, we compare variation at DOG1 with neutral markers and quantitative variation in seed dormancy. Patterns of genetic differentiation among populations suggest that the gene DOG1 contributes to local adaptation. Although Q(ST) for seed dormancy is not different from F(ST) for neutral markers, a correlation with variation in summer precipitation supports that seed dormancy is adaptive. We characterize dormancy variation in several F(2) -populations and show that a series of functionally distinct alleles segregate at the DOG1 locus. Theoretical models have shown that the number and effect of alleles segregatin at quantitative trait loci (QTL) have important consequences for adaptation. Our results provide support to models postulating a large number of alleles at quantitative trait loci involved in adaptation.     

Genetic stability in a population of a plant pathogenic fungus over time

Collections of the plant pathogenic fungus Mycosphaerella graminicola were made from the same field of wheat over a 3-year period. The field was planted with small plots containing four varieties of wheat grown in pure stand and in all possible two-, three- and four-way mixtures. In each year, the wheat field was recolonized by a local source of inoculum of unknown origin. Allele frequencies at 10 RFLP loci were compared at two different times within a growing season and over the 3-year period. No significant differences in allele frequencies were found for any of the RFLP loci over any of the time periods. DNA fingerprints were used to identify clones produced by asexual reproduction. Genotypic diversity based on the frequency of each clone was compared for each collection. No significant changes in genotypic diversity were found within a year or between years. Identical genotypes were found in the field at different times within a season, but no clones were conserved between years. No clone existed in a high frequency in any year, suggesting that selection for particular asexual lineages was weak. The founding population each year probably originated from wind-borne ascospores of the teleomorph, which may exist as an indigenous population on alternative hosts, such as Poa annua (annual bluegrass).     

Agronomic and meteorological factors affecting the severity of leaf blotch caused by Mycosphaerella graminicola in commercial wheat crops in England

Factors affecting the severity of leaf blotch on the two upper leaves of wheat plants in crops at the milky ripe growth stage (GS 73–75) were investigated using survey data from 3513 randomly selected wheat crops sampled during 1985–1996. Year‐to‐year variation in disease severity was greater than spatial variability at county level, although both showed significant differences. The presence of disease above a 5% severity threshold was modelled using random effects logistic regression (Generalised Linear Mixed Model), which enabled risk variables measured at the field level to be combined with meteorological variables estimated at county level. The final model included terms for the fixed effects of disease resistance rating, date of sowing, high risk septoria periods in May and June, number of fungicide sprays and number of days with frost (≤‐2°C) in November. The percentage of crops above the threshold decreased with later sowing, increased number of November frost days and increased number of fungicide sprays. In contrast, high risk septoria periods (rain splash events) in May and June showed a positive correlation with the percentage of crops above the threshold. There were benefits from using resistant cultivars. The model showed that a range of risk variables were of broadly equivalent importance in determining the development of leaf blotch. These risk variables should be integrated in any scheme designed to support fungicide use decisions.     

MgHog1 regulates dimorphism and pathogenicity in the fungal wheat pathogen Mycosphaerella graminicola

The dimorphic ascomycete pathogen Mycosphaerella graminicola switches from a yeastlike form to an infectious filamentous form that penetrates the host foliage through stomata. We examined the biological function of the mitogen-activated protein kinase-encoding gene MgHog1 in M. graminicola. Interestingly, MgHog1 mutants were unable to switch to filamentous growth on water agar that mimics the nutritionally poor conditions on the foliar surface and, hence, exclusively developed by a yeastlike budding process. Consequently, due to impaired initiation of infectious germ tubes, as revealed by detailed in planta cytological analyses, the MgHog1 mutants failed to infect wheat leaves. We, therefore, conclude that MgHog1 is a new pathogenicity factor involved in the regulation of dimorphism in M. graminicola. Furthermore, MgHog1 mutants are osmosensitive, resistant to phenylpyrrole and dicarboximide fungicides, and do not melanize.     

Divergent selection and phenotypic plasticity during incipient speciation in Lake Victoria cichlid fish

Divergent selection acting on several different traits that cause multidimensional shifts are supposed to promote speciation, but the outcome of this process is highly dependent on the balance between the strength of selection vs. gene flow. Here, we studied a pair of sister species of Lake Victoria cichlids at a location where they hybridize and tested the hypothesis that divergent selection acting on several traits can maintain phenotypic differentiation despite gene flow. To explore the possible role of selection we tested for correlations between phenotypes and environment and compared phenotypic divergence (P_ST) with that based on neutral markers (F_ST). We found indications for disruptive selection acting on male breeding colour and divergent selection acting on several morphological traits. By performing common garden experiments we also separated the environmental and heritable components of divergence and found evidence for phenotypic plasticity in some morphological traits contributing to species differences.     

Contrasting introduction scenarios among continents in the worldwide invasion of the banana fungal pathogen Mycosphaerella fijiensis

Reconstructing and characterizing introduction routes is a key step towards understanding the ecological and evolutionary factors underlying successful invasions and disease emergence. Here, we aimed to decipher scenarios of introduction and stochastic demographic events associated with the global spread of an emerging disease of bananas caused by the destructive fungal pathogen Mycosphaerella fijiensis. We analysed the worldwide population structure of this fungus using 21 microsatellites and 8 sequence‐based markers on 735 individuals from 37 countries. Our analyses designated South‐East Asia as the source of the global invasion and supported the location of the centre of origin of M. fijiensis within this area. We confirmed the occurrence of bottlenecks upon introduction into other continents followed by widespread founder events within continents. Furthermore, this study suggested contrasting introduction scenarios of the pathogen between the African and American continents. While potential signatures of admixture resulting from multiple introductions were detected in America, all the African samples examined seem to descend from a single successful founder event. In combination with historical information, our study reveals an original and unprecedented global scenario of invasion for this recently emerging disease caused by a wind‐dispersed pathogen.     

Evolutionary divergence in life history traits among populations of the Lake Malawi cichlid fish Astatotilapia calliptera

During the early stages of adaptive radiation, populations diverge in life history traits such as egg size and growth rates, in addition to eco‐morphological and behavioral characteristics. However, there are few studies of life history divergence within ongoing adaptive radiations. Here, we studied Astatotilapia calliptera, a maternal mouthbrooding cichlid fish within the Lake Malawi haplochromine radiation. This species occupies a rich diversity of habitats, including the main body of Lake Malawi, as well as peripheral rivers and shallow lakes. We used common garden experiments to test for life history divergence among populations, focussing on clutch size, duration of incubation, egg mass, offspring size, and growth rates. In a first experiment, we found significant differences among populations in average clutch size and egg mass, and larger clutches were associated with smaller eggs. In a second experiment, we found significant differences among populations in brood size, duration of incubation, juvenile length when released, and growth rates. Larger broods were associated with smaller juveniles when released and shorter incubation times. Although juvenile growth rates differed between populations, these were not strongly related to initial size on release. Overall, differences in life history characters among populations were not predicted by major habitat classifications (Lake Malawi or peripheral habitats) or population genetic divergence (microsatellite‐based F_ST). We suggest that the observed patterns are consistent with local selective forces driving the observed patterns of trait divergence. The results provide strong evidence of evolutionary divergence and covariance of life history traits among populations within a radiating cichlid species, highlighting opportunities for further work to identify the processes driving the observed divergence.     

The landscape of transposable elements in the finished genome of the fungal wheat pathogen Mycosphaerella graminicola

Background

In addition to gene identification and annotation, repetitive sequence analysis has become an integral part of genome sequencing projects. Identification of repeats is important not only because it improves gene prediction, but also because of the role that repetitive sequences play in determining the structure and evolution of genes and genomes. Several methods using different repeat-finding strategies are available for whole-genome repeat sequence analysis. Four independent approaches were used to identify and characterize the repetitive fraction of the Mycosphaerella graminicola (synonym Zymoseptoria tritici) genome. This ascomycete fungus is a wheat pathogen and its finished genome comprises 21 chromosomes, eight of which can be lost with no obvious effects on fitness so are dispensable.

Results

Using a combination of four repeat-finding methods, at least 17% of the M. graminicola genome was estimated to be repetitive. Class I transposable elements, that amplify via an RNA intermediate, account for about 70% of the total repetitive content in the M. graminicola genome. The dispensable chromosomes had a higher percentage of repetitive elements as compared to the core chromosomes. Distribution of repeats across the chromosomes also varied, with at least six chromosomes showing a non-random distribution of repetitive elements. Repeat families showed transition mutations and a CpA → TpA dinucleotide bias, indicating the presence of a repeat-induced point mutation (RIP)-like mechanism in M. graminicola. One gene family and two repeat families specific to subtelomeres also were identified in the M. graminicola genome. A total of 78 putative clusters of nested elements was found in the M. graminicola genome. Several genes with putative roles in pathogenicity were found associated with these nested repeat clusters. This analysis of the transposable element content in the finished M. graminicola genome resulted in a thorough and highly curated database of repetitive sequences.

Conclusions

This comprehensive analysis will serve as a scaffold to address additional biological questions regarding the origin and fate of transposable elements in fungi. Future analyses of the distribution of repetitive sequences in M. graminicola also will be able to provide insights into the association of repeats with genes and their potential role in gene and genome evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1132) contains supplementary material, which is available to authorized users.     

Local adaptation with high gene flow: temperature parameters drive adaptation to altitude in the common frog (Rana temporaria)

Both environmental and genetic influences can result in phenotypic variation. Quantifying the relative contributions of local adaptation and phenotypic plasticity to phenotypes is key to understanding the effect of environmental variation on populations. Identifying the selective pressures that drive divergence is an important, but often lacking, next step. High gene flow between high‐ and low‐altitude common frog (Rana temporaria) breeding sites has previously been demonstrated in Scotland. The aim of this study was to assess whether local adaptation occurs in the face of high gene flow and to identify potential environmental selection pressures that drive adaptation. Phenotypic variation in larval traits was quantified in R. temporaria from paired high‐ and low‐altitude sites using three common temperature treatments. Local adaptation was assessed using Q_ST–F_ST analyses, and quantitative phenotypic divergence was related to environmental parameters using Mantel tests. Although evidence of local adaptation was found for all traits measured, only variation in larval period and growth rate was consistent with adaptation to altitude. Moreover, this was only evident in the three mountains with the highest high‐altitude sites. This variation was correlated with mean summer and winter temperatures, suggesting that temperature parameters are potentially strong selective pressures maintaining local adaptation, despite high gene flow.     

Global migration patterns in the fungal wheat pathogen Phaeosphaeria nodorum

The global migration patterns of the fungal wheat pathogen Phaeosphaeria nodorum were analysed using 12 microsatellite loci. Analysis of 693 isolates from nine populations indicated that the population structure of P. nodorum is characterized by high levels of genetic diversity and a low degree of subdivision between continents. To determine whether genetic similarity of populations was a result of recent divergence or extensive gene flow, the microsatellite data were analysed using an isolation-with-migration model. We found that the continental P. nodorum populations diverged recently, but that enough migration occurred to reduce population differentiation. The migration patterns of the pathogen indicate that immigrants originated mainly from populations in Europe, China and North America.     

The effect of Trichoderma harzianum and T. koningii on the control of tan spot (Pyrenophora tritici-repentis) and leaf blotch (Mycosphaerella graminicola) of wheat under field conditions in Argentina

The effect of six isolates of Trichoderma harzianum and one isolate of T. koningii on the incidence and severity of tan spot (Pyrenophora tritici-repentis) and leaf blotch of wheat (Mycosphaerella graminicola) was evaluated under field conditions. Significant differences between wheat cultivars, inoculum types and growth stages were found. Three of the isolates tested (T2 for M. graminicola, T7 for P. tritici-repentis and T5 for both of them) showed the best performance in controlling leaf blotch and tan spot when coated onto seed or sprayed onto wheat leaves at different growth stages, with significant severity reduction up to 56%. At tillering, six of the isolates reduced the severity of P. tritici-repentis and M. graminicola compared to the control by up to 39% and 12-53%, respectively. In some experiments, the biocontrol preparation (T2 and T5) gave a level of disease control similar to that obtained with Tebuconazole (70 and 48%, respectively). The effect of Trichoderma against P. tritici-repentis was also observed at the heading stage, when six of the treatments reduced disease severity by 16-35%. This is the first report on the efficacy of Trichoderma spp. against wheat necrotrophic pathogens under field conditions in Argentina.     

Polymerase chain reaction–restriction fragment length polymorphism markers for the fungal banana pathogen Mycosphaerella fijiensis

Polymerase chain reaction–restriction fragment length polymorphism markers were developed to study populations of the fungal banana pathogen Mycosphaerella fijiensis. Twelve markers were defined, 11 in anonymous and single‐copy nuclear DNA sequences and one in the internal transcribed spacer and 5.8S rDNA sequence. The polymerase chain reaction products obtained with locus‐specific primer pairs were digested with restriction enzymes to reveal polymorphism. Between five and 12 markers were polymorph in M. fijiensis populations from different geographical origins (Papua New Guinea, the Philippines, Cameroon and Latin America). The mean of allele number and gene diversity (expected heterozygosity) per locus in the different geographical populations ranged between 1.4 and 2.7 and 0.17 and 0.45, respectively.     

The impact of bottlenecks on microbial survival,adaptation, and phenotypic switching in host–pathogen interactions

Microbial pathogens and viruses can often maintain sufficient population diversity to evade a wide range of host immune responses. However, when populations experience bottlenecks, as occurs frequently during initiation of new infections, pathogens require specialized mechanisms to regenerate diversity. We address the evolution of such mechanisms, known as stochastic phenotype switches, which are prevalent in pathogenic bacteria. We analyze a model of pathogen diversification in a changing host environment that accounts for selective bottlenecks, wherein different phenotypes have distinct transmission probabilities between hosts. We show that under stringent bottlenecks, such that only one phenotype can initiate new infections, there exists a threshold stochastic switching rate below which all pathogen lineages go extinct, and above which survival is a near certainty. We determine how quickly stochastic switching rates can evolve by computing a fitness landscape for the evolutionary dynamics of switching rates, and analyzing its dependence on both the stringency of bottlenecks and the duration of within‐host growth periods. We show that increasing the stringency of bottlenecks or decreasing the period of growth results in faster adaptation of switching rates. Our model provides strong theoretical evidence that bottlenecks play a critical role in accelerating the evolutionary dynamics of pathogens.     

Chromosomal location of a gene for resistance to septoria tritici blotch (Mycosphaerella graminicola)in the hexaploid wheat ’Synthetic 6x’

Septoria tritici blotch, caused by the fungus Mycosphaerella graminicola,is currently the major foliar disease of wheat world-wide, and new sources of resistance and knowledge about the genetics of resistance are needed to improve breeding for resistance to this disease. Sears’s ’Synthetic 6x’ hexaploid wheat, derived from a hybrid of Triticum dicoccoides and Triticum tauschii, was resistant to 12 of 13 isolates of M. graminicola tested. Chromosome 7D of ’Synthetic 6x’ was identified as carrying resistance to all 12 isolates in tests of seedlings of inter-varietal chromosome substitution lines of ’Synthetic 6x’ into ’Chinese Spring’ and to two isolates in tests of adult plants. A septoria tritici blotch resistance gene, named Stb5, was identified using the M. graminicola isolate IPO94269 and mapped on the short arm of chromosome 7D, near the centromere, in a population of single homozygous chromosome-recombinant lines for the 7D chromosome. Received: 1 February 2001 / Accepted: 17 April 2001     

Fine-scale thermal adaptation in a green turtle nesting population

The effect of climate warming on the reproductive success of ectothermic animals is currently a subject of major conservation concern. However, for many threatened species, we still know surprisingly little about the extent of naturally occurring adaptive variation in heat-tolerance. Here, we show that the thermal tolerances of green turtle (Chelonia mydas) embryos in a single, island-breeding population have diverged in response to the contrasting incubation temperatures of nesting beaches just a few kilometres apart. In natural nests and in a common-garden rearing experiment, the offspring of females nesting on a naturally hot (black sand) beach survived better and grew larger at hot incubation temperatures compared with the offspring of females nesting on a cooler (pale sand) beach nearby. These differences were owing to shallower thermal reaction norms in the hot beach population, rather than shifts in thermal optima, and could not be explained by egg-mediated maternal effects. Our results suggest that marine turtle nesting behaviour can drive adaptive differentiation at remarkably fine spatial scales, and have important implications for how we define conservation units for protection. In particular, previous studies may have underestimated the extent of adaptive structuring in marine turtle populations that may significantly affect their capacity to respond to environmental change.     

The global genetic structure of the wheat pathogen Mycosphaerella graminicola is characterized by high nuclear diversity,low mitochondrial diversity,regular recombination,and gene flow

A total of 1673 Mycosphaerella graminicola strains were assayed for DNA fingerprints and restriction fragment length polymorphism (RFLP) markers in the nuclear and mitochondrial genomes. The isolates were collected from 17 wheat fields located in 11 countries on five continents over a six year period (1989-1995). Our results indicate that genetic diversity in the nuclear genome of this fungus was high for all but three of the field populations surveyed and that populations sampled from different continents had similar frequencies for the most common RFLP alleles. Hierarchical analysis revealed that more than 90% of global gene diversity was distributed within a wheat field, while approximately 5% of gene diversity was distributed among fields within regions and approximately 3% was distributed among regions on different continents. These findings suggest that gene flow has occurred on a global scale. On average, each leaf was colonized by a different nuclear genotype. In contrast, only seven mtDNA haplotypes were detected among the 1673 isolates and the two most common mtDNA haplotypes represented approximately 93% of the world population, consistent with a selective sweep. Analysis of multilocus associations indicated that all field populations were in gametic equilibrium, suggesting that sexual recombination is a regular occurrence globally.