This paper presents the European distribution of the understudied, host-specific rust fungus, Uromyces erythronii (Pucciniomycetes, Pucciniales, Pucciniaceae). Distribution data were derived from the survey of herbarium materials of its European host plant, Erythronium dens-canis. We demonstrate the presence of this rust fungus in 14 countries within the distribution area of its host. The temporal trend of emergence of the two rust fungus generations (aecia and telia) is presented. Based on the study of 1700 E. dens-canis individuals, we conclude that the overall infection rate has not changed significantly over the last 200 years. During field surveys, U. erythronii infection was detectable in most of the studied Erythronium populations (88.5%). A high similarity in the nrITS region was detected among samples from Europe (Croatia, Romania) and Asia (Japan).
Citation: Nagy T., Pfliegler W. P., Takács A., Tökölyi J. & Molnár V. A. 2019: Distribution, infection rates and DNA barcoding of Uromyces erythronii (Pucciniaceae), a parasite of Erythronium (Liliaceae) in Europe. – Willdenowia 49: 13–20. doi: https://doi.org/10.3372/wi.49.49103
Version of record first published online on 14 March 2019 ahead of inclusion in April 2019 issue.
Despite (or just due to) their antagonistic nature, parasitic relationships play an important role in shaping ecosystems and biodiversity (Gómez & Nichols 2013). Diseases constrain hosts to evolve resistance or tolerance strategies (Roy & Kirchner 2000), thus the hosts and parasites are always in antagonistic co-evolutionary race (Stahl & Bishop 2000), unless the interaction eventually evolves towards commensalism (Miller & al. 2006).
The autoecious rust fungus Uromyces erythronii (DC.) Pass. parasitizes some spring geophyte taxa of the vascular plant family Liliaceae Juss., e.g. Erythronium japonicum Decne., Amana edulis (Mig.) Honda and A. latifolia (Makino) Honda (Fukuda & Nakamura 1987). In Europe, the only known host plant of this fungus is E. dens-canis L., a Red List species in several countries (Witkowski & al. 2003). According to the criteria of the IUCN Red List, E. dens-canis is categorized as Vulnerable (“VU”) in the entire Carpathians and in Hungary, Critically Endangered (“CR”) in Slovakia, and Endangered (“EN”) in Ukraine. The species is present in Romania but not threatened (“+”) (Witkowski & al. 2003). Even though threat assessment and conservation of parasites may be just as important as that of their hosts (Gómez & Nichols 2013), co-endangerment remains poorly investigated, currently being restricted to a few published works (e.g. Dunn & al. 2009; Mihalca & al. 2011 on parasitic insects).
The distribution and abundance of parasites are dependent on host populations, as shown in the case of Uromyces erythronii and Erythronium japonicum (Fukuda & Nakamura 1987). In Europe, the presence of U. erythronii in several countries was only recorded recently, e.g. in 2007 from Croatia (Miličević & al. 2008) and in 2016 in Ukraine (Tykhonenko & al. 2017). The earliest deposited specimens of U. erythronii in the Fungi Collection of the Royal Botanical Gardens, Kew, were collected in 1946 from France, in 1949 from Italy and in 1962 from Romania (with further items without collection date from Switzerland, Germany and Romania; Kew s.d.+). From Hungary, six samples of U. erythronii (collected in 1965, 1987, 1988 and 2007 from three localities) were deposited in the Fungi Collection of the Hungarian Natural History Museum (cf. Jandrasits & Fischl 2014). Outside the range of the indigenous distribution of its host, the presence of U. erythronii has been confirmed in England (Henderson 2000) and in Germany (Kruse & al. 2014), where it infects cultivated populations of Erythronium. The low number of known occurrence data, and the fact that all of these were recorded in the last c. 60 years, indicates a poor documentation of the distribution of this rust in Europe, or may refer to recent colonization and spreading.
Our aims in this study were (1) to map the distribution range and the frequency of Uromyces erythronii in Europe, predominantly the Carpathian Basin, and (2) to document possible temporal trends in the presence of the aecial and telial generations of the rust. We also aimed (3) to study whether the infection rate of the rust fungus has changed over time and (4) to generate ITS DNA barcodes for U. erythronii samples from the Carpathian Basin to facilitate future genetic comparisons of the populations of the species in different geographical areas and on different host plants.
Material and methods
We reviewed 525 herbarium sheets (altogether more than 2000 individuals) of Erythronium dens-canis in eight herbaria (BP, CL, DE, GK, LJU, RRM, ZA, ZAHO – herbarium codes according to Thiers 2018+, with the exception of the unregistered RRM = Rippl-Rónai Museum, Kaposvár, Hungary) and 24 items in the fungi collection of BP (Table 1). We screened leaves and petioles of herbarium specimens for the presence of Uromyces erythronii with a 15× hand lens. We recorded the total number of individuals and the number of infected individuals on each sheet (with the exception of the fungi collection items, because the paper bags usually contain fragmented leaf materials instead of entire individuals, therefore we treated every item as an individual sample). Information from the herbarium labels (site, date [year, month, day] and name of collector for each collection) were registered. Number of aecial- and telial-infected individuals were quantified separately only in the materials of BP and DE. Digital photographs of all herbarium sheets were taken. We used QGIS 2.18 (Quantum GIS Development Team 2017) software for generating a distribution map.
Collection period, sample sizes and infection rates of the screened Erythronium dens-canis herbarium specimens, summarized for each studied herbarium collection.
We used a generalized linear model (GLM) with binomial error distribution to evaluate the change in infection rate of Uromyces erythronii over the last 200 years. Analysis was done based on 1700 Erythronium dens-canis individuals from 427 herbarium sheets, which were labelled with at least the year of collection. All sheets with uncertain year of collection were disregarded in the analysis. The number of infected and uninfected individuals in each year entered the model as a dependent variable, while year of collection was included as an explanatory variable. Analysis was done in the R Statistical Environment (R Core Team 2017).
We checked the presence of Uromyces erythronii aecia in 26 (16 Romanian, seven Hungarian, three Croatian) populations of Erythronium dens-canis during flowering and fruiting phenological states. Fieldwork was carried out between 16 April and 2 May 2015, between 6 March and 9 April 2016, and between 9 and 10 March 2017. We checked 16–135 individuals (mean ± SD = 73 ± 30) in each population, depending on the local population size of the host plant (Table 3).
Infection rates derived from the studied Erythronium dens-canis herbarium sheets. Number of screened sheets, the total number of individuals on these, as well as the year of collection of the oldest infected specimen are given for each country separately.
Main data of field survey: geographic location, elevation (alt.), date of observation, number of screened host individuals (n) and infection rate observed. Localities are listed alphabetically, first by country then by settlement. Abbreviations of countries studied: Cro – Croatia, Hu – Hungary, Ro – Romania.
DNA was isolated from aeciospores scraped off of minute pieces of aecia removed from the leaf of a Romanian sample (collection year: 2016) and the leaf of a Croatian sample (collection year: 2017) of Erythronium dens-canis. Samples were deposited in the herbarium of Debrecen University (DE-soo-43284, DE-soo-43285). Following the microwave extraction protocol described in Haelewaters & al. (2015), the spores were placed in 0.5 mL PCR tubes and microwave-treated (750 W for 5 min). Subsequently 50 Î¼L ddH20 was added to the tube, and the spores were manually crushed with a sterile pipette tip while viewed under a dissecting microscope. The PCR tubes were then incubated at –20°C for 20 min. For each PCR reaction, 5 Î¼L of the extracted material was used. PCR amplification of the internal transcribed spacer (ITS) region was carried out using GoTaq polymerase (Promega, Madison, WI, U.S.A.) and the primer pair ITS1f and ITS4 (White & al. 1990; Gardes & Bruns 1993). The PCR protocol was as follows: 95°C for 5 min, 30× (94°C 50 s, 55°C 50 s, 72°C 50 s), 72°C for 5 min. For amplification, an Applied Biosystems (Foster City, CA, U.S.A.) 2720 thermal cycler and a final volume of 50 Î¼L were used. PCR products were loaded onto 1.2% agarose gels for electrophoresis at 100 V for 15 min and UV transillumination was used to check the product size. PCR products were cleaned with a Geneaid (New Taipei City, Taiwan) DF100 PCR cleaning kit and sequenced in both directions using the same primers (Microsynth AG, Switzerland). Sequences were trimmed of ambiguous bases at both ends in Chromas 2.6.5. (Technelysium Pty. Ltd.). Sequences compiled from the reads were deposited in ENA (European Nucleotide Archive; accession numbers: MH205916 and MH205917) and blasted in NCBI GenBank (NCBI s.d.+). Boundaries of ITS1, the 5.8S rDNA, and ITS2 were identified using ITSx (Bengtsson-Palme & al. 2013).
Light microscopy images were taken of material mounted in HPVA medium, with an Olympus BD40 microscope (Debrecen, Hungary) equipped with an Olympus 40× lens and digital microscope camera, using the Olympus DP Controller software. Macrophotographs were taken with a Pentax k7 DSLR camera with macrophotography setup and flash.
Uromyces erythronii infection was detected on 199 herbarium sheets (335 individuals) of Erythronium dens-canis (Table 1). Altogether 81.2 % of herbarium specimens (446 items) had exact locality data, and 62.1% of herbarium specimens (341 items) had an exact date (year, month and day). Herbarium sheets (including items from the BP fungi collection) were collected between 1811–2017 (Table 1). Among herbarium specimens, the earliest collection date across all years was 27 February, and the latest was 22 June. Additionally, an exceptionally late date, 26 October, was registered for a specimen with a vague locality description, but based on the date, a high mountain locality in the Carpathians is to be suspected.
Most of the studied herbarium sheets of Erythronium (198), as well as most rust-infected sheets (70) and individuals (115), originated from Romania (Table 2, Fig. 2). Numerous further infected individuals were collected in Hungary (62), Croatia (55) and Ukraine (29) (Table 2, Fig. 2). Infection rate of the specimens from the three most represented countries were 56.0% (Hungary), 50.0% (Croatia) and 35.4% (Romania). The infection rate of individuals was 28.5%, 18.8% and 14.7%, respectively (Table 2).
Aecial and telial infection were separately quantified on a total of 153 specimens (251 individuals). Aecia were found on 141 specimens (225 individuals) collected between 27 February and 28 May, and telia were found on 22 specimens (32 individuals) collected between 7 March and 22 June (Fig. 3). Simultaneous occurrence of both developmental stages of the fungus were detected only in the case of six individuals.
According to the GLM, the year of collection had no significant effect on the infection rate of Erythronium dens-canis individuals by Uromyces erythronii (DF = 135, E = 0.002, SE = 0.002, t-value = 0.945, p-value = 0.344).
Uromyces erythronii infection was detectable for 88.5% of the studied Erythronium populations. Average infection rate of populations was mean ± SD = 25 ± 21% (range: 0–69%). The highest infection rates were documented in Tömörd, Hungary (69%) and Feleacu, Romania (65%) (Table 3). The fungus remained undetected in only three screened host populations (Hungary: Miskolc; Romania: Haţeg and Juliţa).
We generated ITS (ITS1-5.8S-ITS2) DNA barcodes for two samples collected during field surveys in Romania (MH205916) and Croatia (MH205917). The quick heat extraction protocol (Haelewaters & al. 2015) was successfully applied, with minute pieces of aecia being sufficient for subsequent PCR amplification. The obtained sequences were identical in their 589 bp overlapping regions (of which 517 bases corresponded to the ITS15.8S-ITS2). Comparison with the single available DNA barcode sequence of the species in GenBank (accession number LC203755), generated for a sample collected in Japan from the host Erythronium japonicum Decne., revealed a difference of only two nucleotides in the ITS15.8S-ITS2 region (= 99.61% similarity). Both of these, a substitution and an additional nucleotide in the European samples, were located in the ITS1.
Our new records of Uromyces erythronii revealed a distribution overlapping with much of its European host, Erythronium dens-canis (Fig. 1). Formal reports did not exist on the presence of the rust in Bosnia and Herzegovina, Bulgaria, Montenegro, Slovakia and Slovenia, although possible presence was suggested (e.g. Kruse & al. 2014; Tykhonenko & al. 2017). Our study confirms the presence of the rust in these countries (Table 2) based on frequently unintentional collections of the fungus in the form of Erythronium herbarium specimens. These specimens occasionally predate the formal description of the fungus (Table 2, Fig. 2).
The often remarkably high infection rate in herbarium material was reflected in the results of the field surveys conducted within the framework of this study. An overwhelming majority of studied Croatian, Hungarian and Romanian Erythronium populations were found to be infected. Based on our experiences, examining relatively few Erythronium individuals during the flowering and fruiting seasons may be sufficient to detect Uromyces erythronii in the field. On the other hand, infection rates vary widely across populations (0–69%; Table 3). Variable infection rate is a known phenomenon in case of rust infections (Ericson & al. 1999), and it has to be noted that the highest infection rate was detected in a population (Tömörd, Hungary) that most probably originated from a deliberate planting (Molnár V., ined.). Such inbred populations are generally more susceptible to pathogens than outcrossed populations (e.g. Burdon & al. 1999).
The connected life cycles of Uromyces erythronii and Erythronium japonicum were illustrated and detailed by Fukuda & Nakamura (1987). They found that the growing period of E. japonicum is restricted to one or two months (from late February to early May). However, we documented a somewhat longer period for E. dens-canis: 10% of the infected individuals were collected in the second and third 10-day periods of May. Fukuda & Nakamura highlighted that the early-emerging individuals of E. japonicum are more rusty than later-growing ones. It seems that this phenomenon is also true in the case of E. dens-canis, because individuals collected in late February bear aecia in higher frequency (29%) than individuals collected from March to May (8–14%) (Fig. 3). Fukuda & Nakamura also stated that aecia may develop only on the larger leaves in case of two-leaved individuals of Erythronium. Contrary to this, we observed aecial infection on both leaves in some cases, which may be caused by independent infection events by basidiospores. Otherwise, the less tightly rolled leaves of E. dens-canis in the emerging stage may also explain the higher infection rate of these. Secondary infection by aeciospores generating the telial generation affected only a small proportion of individuals (Fig. 3) in line with the observations of Fukuda & Nakamura.
Invasions and emerging infectious diseases caused by representatives of various fungus families have been extensively discussed (e.g. Anderson & al. 2004; Parker & Gilbert 2004; Desprez-Loustau & al. 2007). Fungal invasions often originate from host shift, even across continents (Palm & Rossman 2003). According to our results, the infection rate of Uromyces erythronii on Erythronium dens-canis did not change during the documented last two centuries (GLM, DF= 135, E = 0.002, SE = 0.002, t-value = 0.945, p-value = 0.344). Consequently, it is apparently not a novel host-pathogen relationship and the current wide distribution of the host emerged much earlier.
DNA barcode sequencing revealed a high similarity in the ITS region among samples from Croatia, Romania and Japan. The 100% identity of the two European sequences and their 99.61% similarity to a Japanese sample from a different host species, Erythronium japonicum, in the ITS region suggests a low intraspecific geographic or host-related diversity and no indication of cryptic speciation in Uromyces erythronii.
In conclusion, Uromyces erythronii is more widespread and abundant in Europe than previously reported. The species potentially occurs in the whole distribution area of its host. The apparent expansion of U. erythronii is not a consequence of its recent spread, because its earlier distribution is well documented but previously unstudied by herbarium specimens collected decades or even centuries earlier. Our results also signify the importance of incorporating vascular plant collections in the study of microscopic fungi (see also Denchev & Denchev 2016).
We thank the contribution of herbarium material for study in BP, CL, DE, GK, LJU, RRM, ZA and ZAHO. The authors are grateful to Vedran Šegota for his kind help in localization of Croatian herbarium data in the ZA and ZAHO collections. We are much obliged to Danny Haelewaters and an anonymous reviewer for the professional comments on the manuscript of this paper. Linguistic corrections by Orsolya Vincze (Debrecen–Kolozsvár) are also acknowledged. The research was supported by the OTKA K108992 grant and EFOP-3.6.3-VEKOP-16-2017-00008 project, which is co-financed by the European Union and the European Social Fund.
Anderson P. K., Cunningham A. A., Patel N. G., Morales F. J., Epstein P. R. & Daszak P. 2004: Emerging infectious diseases of plants: pathogen pollution, climate change and agrotechnology drivers. – Trends Ecol. Evol.19: 535–544.
Bengtsson-Palme J., Veldre V., Ryberg M., Hartmann M., Branco S., Wang Z., Godhe A., Bertrand Y, De Wit P., Sanchez M., Ebersberger I., Sanli K, de Souza F., Kristiansson E., Abarenkov K., Eriksson K. M. & Nilsson R. H. 2013: ITSx: improved software detection and extraction of ITS1 and ITS2 from ribosomal ITS sequences of fungi and other eukaryotes for use in environmental sequencing. – Meth. Ecol. Evol.4: 914–919.
Denchev T. T. & Denchev C. M. 2016: Contribution to the smut fungi of Greece. – Willdenowia46: 233–244.
Desprez-Loustau M. L., Robin C., Buee M., Courtecuisse R., Garbaye J., Suffert F., Sache I. & Rizzo D. M. 2007: The fungal dimension of biological invasions. – Trends Ecol. Evol.22: 472–480.
Dunn R. R., Harris N. C., Colwell R. K., Koh L. P. & Sodhi N. S. 2009: The sixth mass coextinction: are most endangered species parasites and mutualists? – Proc. Roy. Soc. London, Ser. B, Biol. Sci.276: 3037–3045.
Ericson L., Burdon J. J. & Müller W. J. 1999: Spatial and temporal dynamics of epidemics of the rust fungus Uromyces valerianae on populations of its host Valeriana salina. – J. Ecol.87: 649–658.
Fukuda T. & Nakamura S. 1987: Biotic interaction between a rust fungus, Uromyces erythronii Pass., and its host plant, Erythronium japonicum Decne. (Liliaceae). – Pl. Spec. Biol.2: 75–83.
Gardes M. & Bruns T. D. 1993: ITS primers with enhanced specificity for basidiomycetes-application to the identification of mycorrhizae and rusts. – Molec. Ecol.2: 113–118.
Gómez A. & Nichols E. 2013: Neglected wild life: parasitic biodiversity as a conservation target. – Int. J. Parasitol. Parasites & Wildlife2: 222–227.
Haelewaters D., Gorczak M., Pfliegler W. P., Tartally A., Tischer M., Wrzosek M. & Pfister D. H. 2015: Bringing Laboulbeniales into the 21st century: enhanced techniques for extraction and PCR amplification of DNA from minute ectoparasitic fungi. – IMA Fungus6: 363–372.
Kews.d.+ [continuously updated]: HerbIMI fungarium database.Royal Botanic Gardens, Kew. – Published at http://www.herbimi.info [accessed 26 Oct 2018].
MihalcaA. D., GhermanC. M. & CozmaV.2011: Co-endangered hard-ticks: threatened or threatening? – Parasites & Vectors4(1): 71.
MiličevićT., IvićD., KaliternaJ. & CvjetkovićB.2008: Occurrence of phytopathogenic fungi on some endemic and protected plants in natural ecosystems of Croatia. – P. 87in: Abstract book of XII International Congress of Mycology, August 2008, Istanbul. – Istanbul: International Union of Microbiological Societies.
Miller M. R., White A. & Boots M. 2006: The evolution of parasites in response to tolerance in their hosts: the good, the bad, and apparent commensalism. – Evolution60: 945–956.
NCBIs.d.+ [continuously updated]: BLAST Basic Local Alignment Search Tool. – Published at http://ncbi.nlm.nih.gov/blast/Blast.cgi [accessed 26 Oct 2018].
Parker I. M. & Gilbert G. S. 2004: The evolutionary ecology of novel plant-pathogen interactions. – Annual Rev. Ecol. Evol. Syst.35: 675–700.
Quantum GIS Development Team2017: Quantum GIS Development Team Quantum GIS Geographic Information System. 2017. – Published at https://qgis.org/hu/site/index.html
R Core Team2017: R: a language and environment for statistical computing. – Vienna: R Foundation for Statistical Computing. – Published at https://www.r-project.org/
Roy B. A. & Kirchner J. W. 2000: Evolutionary dynamics of pathogen resistance and tolerance. – Evolution54: 51–63.
Stahl E. A. & Bishop J. G. 2000: Plant-pathogen arms races at the molecular level. – Curr. Opin. Pl. Biol.3: 299–304.
ThiersB.2018+ [continuously updated]: Index Herbariorum: a global directory of public herbaria and associated staff.New York Botanical Garden's Virtual Herbarium. – Published at http://sweetgum.nybg.org/science/ih/
Tykhonenko Y., Sytschak N. N., Kagalo A. A. & Orlov O. O. 2017: New records of Uromyces erythronii (Pucciniales) from Ukraine. – Ukrayins'k. Bot. Zhurn.74: 184–188.
White T. J., Bruns T., Lee S. J. W. T. & Taylor J. L. 1990: Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. – P. C. R. Meth. Applic.18: 315–322.