# Sentinel Trees as a Tool to Forecast Invasions of Alien Plant Pathogens

Anna-Maria Vettraino, Alain Roques, Annie Yart, Jian-Ting Fan, Jiang-Hua Sun, Andrea Vannini

► To cite this version:

Anna-Maria Vettraino, Alain Roques, Annie Yart, Jian-Ting Fan, Jiang-Hua Sun, et al. Sentinel Trees as a Tool to Forecast Invasions of Alien Plant Pathogens. PLoS ONE, 2015, 10 (3), 15 p. 10.1371/journal.pone.0120571. [DOI]

HAL Id: hal-02640514

https://hal.inrae.fr/hal-02640514

Submitted on 28 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Sentinel Trees as a Tool to Forecast Invasions of Alien Plant Pathogens

AnnaMaria Vettraino¹, Alain Roques², Annie Yart², Jian-ting Fan³, Jiang-hua Sun⁴, Andrea Vannini¹*

1 DIBAF, University of Tuscia, Viterbo, Italy, 2 INRA-UR633, Zoologie Forestière, Centre de recherche d'Orléans, Orléans, France, 3 School of Forestry and Bio-technology, Zhejiang A & F University, Lin'an, China, 4 State key laboratory of Integrated Management of pest Insects and Rodents, Institute of Zoology, Chinese Academy of Sciences, Beijing, China

* vannini@unitus.it

<figure><img src="image_2.png" /></figure>

CrossMark
click for updates

## Abstract

Recent disease outbreaks caused by alien invasive pathogens into European forests posed a serious threat to forest sustainability with relevant environmental and economic effects. Many of the alien tree pathogens recently introduced into Europe were not previously included on any quarantine lists, thus they were not subject to phytosanitary inspections. The identification and description of alien fungi potentially pathogenic to native European flora before their introduction in Europe, is a paramount need in order to limit the risk of invasion and the impact to forest ecosystems. To determine the potential invasive fungi, a sentinel trees plot was established in Fuyang, China, using healthy seedlings of European tree species including *Quercus petreae*, *Q. suber*, and *Q. ilex*. The fungal assemblage associated with symptomatic specimens was studied using the tag-encoded 454 pyrosequencing of the nuclear ribosomal internal transcribed spacer-1 (ITS 1). Taxa with probable Asiatic origin were identified and included plant pathogenic genera. These results indicate that sentinel plants may be a strategic tool to improve the prevention of bioinvasions.

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Vettraino A, Roques A, Yart A, Fan J-t, Sun J-h, Vannini A (2015) Sentinel Trees as a Tool to Forecast Invasions of Alien Plant Pathogens. PLoS ONE 10(3): e0120571. doi:10.1371/journal.pone.0120571

**Academic Editor:** Boris Alexander Vinatzer, Virginia Tech, UNITED STATES

**Received:** November 28, 2014

**Accepted:** January 24, 2015

**Published:** March 31, 2015

**Copyright:** © 2015 Vettraino et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper.

**Funding:** This work was supported by ISEFOR (Increasing Sustainability of European Forests: Modelling for Security Against Invasive Pests and Pathogens under Climate Change - collaborative project 245268) and COST PERMIT (Pathway Evaluation and pest Risk Management InTransport- FP1002-181110-06882). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

## Introduction

In the past 200yrs, the number of fungal and fungal-like Invasive Alien Species (IAS) challenging plants in Europe has increased exponentially, especially in regards to ascomycetes [1]. Most of them (57%) have been introduced to Europe through living plants and have caused serious epidemics. The recent outbreak of ash dieback in Central Europe and the UK [2], the spread of pitch canker of pine in plantations in Spain [3], the sudden appearance of *Phytophthora* spp. outbreaks in the UK [4]-[5] challenging both plantations and natural forests are only recent examples of new threats to be added to 'old' invasions still threatening European forests.
**Competing Interests:** The authors have declared that no competing interests exist.

The temporal trend of invasion in Europe reflects the history of trading living plants and associated commodities between Europe and other countries. After World War I, the commercial exchanges with North America increased as well as the number of invasive pathogens from USA [1]. The proliferation of pathogens from Asia reflects the trend in trading of living plants, especially between Europe and China [6]. A common feature of invasive pathogens responsible for epidemics in forests is that most of them were first described (in Europe or elsewhere) after they had escaped from their geographic centre of origin and begun causing damage in a susceptible plant community. In other words, they were unknowns to science before invasion.

Most of the National Plant Protection Organizations (NPPOs) regulate the inspections of living plants and commodities based on lists of known organisms described as invasive elsewhere.

Furthermore, most of the diagnostic protocols (e.g. EPPO standards) are based on visual assessment of symptoms or signs on selected susceptible hosts, while alien pathogens often enter new environments through non-hosts, infected asymptomatic native hosts, or associated materials of commerce (e.g. soil, packaging). In fact, as Darwinian evolution predicts, being adapted to and co-evolved with their original native hosts, many of these 'unknowns to science' are unlikely to do noticeable damage on original hosts in their native ecosystems, and are less likely to be detected even when introduced in a new environment with their original native hosts.

In this context it is necessary to develop strategies to describe and study new potential plant pathogens and evaluate the risk they might pose to European flora before their introduction to the old continent. The Sentinel plant concept seems to properly respond to this need. According to this hypothesis, living plant collections cultivated in an exotic environment are exposed to the inoculum pressure of the native pathogens. If a compatible pathogenic interaction occurs, due to the lack of coevolution with the native pathogens, the plants likely will show evident symptoms. The identification of the causal agents of the different symptomatology provides a list of new pathogens potentially harmful to those plants. An International Sentinel Plant Network (ISPN) has been recently developed where living plant collections at botanic gardens and arboreta around the world are networked and serve as early warning systems to help predict and prevent the invasion of new pests (insects, pathogens, or plants) [7]. However, ISPN is dependent on already established exotic plant collections while the efficiency of the strategy relies on the accurate selection of plants to cultivate according to climatic requirements, and the presence of the same genus or family in the exotic area where plantations are established. In order to explore the efficiency of sentinel plantations as an early warning system to prevent invasion of insects and pathogens potentially harmful to EU flora, two experimental plots with European trees were established in China in 2008 [8]. Results regarding colonization of European trees by Chinese insects are presented in a separate paper [9]. This paper reports on the results of the pathogen survey of one of the two plantations by mean of symptoms records, detection of organisms by molecular (mass sequencing) techniques, their taxonomic collocation, and the potential of pathogenicity to EU hosts.

# Materials and Methods

## Sample area

The study was carried out in a sentinel trees plot established near Fuyang (30,003333 N; 119,799722 E; 110m elevation), ca. 40 km southwest of Hangzhou (Zhejiang province, South eastern China), on a small parcel of agricultural land completely surrounded by a mixed
conifer- broadleaved forest. The forest was mostly composed of *Pinus massonnianna*, *Cunning-hammia lanceolata*, various species of *Fagaceae* and bamboos. The surrounding agricultural land included field crops with an emphasis on rice. The European trees planted in China considered in the study included *Quercus petraea*, *Q. Suber* and *Q. ilex*. The Fuyang site was planted in May 2008 from one-year old bare-rooted seedlings shipped from Europe to China and kept 3 weeks in quarantine at the port of entry before being planted. Before export to China, the bare-rooted seedlings were individually submitted to thorough phytosanitary inspections by quarantine services to ensure the absence of any insect or pathogen damage, and an additional insecticide and fungicide treatment was applied. The seedlings were planted immediately after being allowed to enter China. The initial height of the seedlings was ca.1-1,5 m-tall at the plantation. Twenty-five seedlings per species were planted in randomized blocks in order to allow statistical analyses. Twice per year from 2008 to 2010, the plot was surveyed for the presence of symptoms and signs of diseases on the base of a visual and, where necessary, microscopic observation of specimens.

At the end of the experiment in autumn 2010, all the seedlings were removed, including roots, and destroyed by burning.

## Characterization of the fungal community of leaves

The detection of the fungal community was performed only on *Q. petraea*, *Q. suber* and *Q. ilex* symptomatic tissues of mature leaves collected in 2010, during the growing season. Samples from each tree were bulked and processed, 1sample/host species.

Genomic DNA was extracted from symptomatic tissue samples in China using the Qiagen DNeasy Plant Mini Kit (Qiagen, Milan, Italy) according to the manufacturer's protocols. DNA extracts were quantified using the Qubit Quantitation Kit (Life Technologies, Monza, Italy). DNA extraction was carried out in China and sent to DIBAF-University of Tuscia in Italy for library preparation. The non-coding nuclear rDNA ITS region was amplified using the fungal primer pair ITS1F (5'-AxxxCTTGGTCATTTAGAGGAAGTAA- 3') and ITS2 (5'-BGCT GCGTTCTTCATCGATGC- 3'), where A and B represent the two pyrosequencing adaptors (CCATCTCATCCCTGCGTGTCTCCGACTCAG and CCTATCCCCTGTGTGCCTT GGCAGTCTCAG) and xxx the sample identification barcoding key. The following PCR conditions were used: 94°C for 4 min, 30 cycles of 30 s at 94°C (denaturation), 50°C for 1 min (annealing) and 72°C for 90 s (extension), followed by 10 min at 72°C.

Positive and negative controls were included in the PCR. Amplicons were purified using the Agencourt AMPure XP system (Beckman Coulter Inc., USA) and quantified with the Qubit Quantitation Kit (Invitrogen, USA). The library was sequenced in 1/16 of a standard PicoTitre Plate with the 454 Life Sciences GS-FLX System (Macrogen).

## Bioinformatics and statistical analyses

Raw data were trimmed and analysed using the software CLOTU [10]. Reads with undetermined nucleotide, with any mismatch against tags and primers, and shorter than 150 bp were removed from the data set. Sequence reads from each sample were clustered into similarity-based OTUs (operational taxonomic units) using BLASTClust (single linkage clustering) at 97% similarity and 75% overlap in sequence length between reads in the pairwise alignments [11–14].
One representative sequence from each OTU was submitted to BLASTn [15] for compari-
son against the GenBank non-redundant (NCBI-nr) database.

OTUs were only assigned to species level if i) the query sequence matched database sequences from fungal isolates (including at least one vouchered specimen) with E-values ≤e-100 and percentage sequence identity ≥97% and ii) at the lowest E-values, when there were no contradictions among different species within the same genus [16]. At lower percentage sequence identity, the query read was identified to genus (identity 95-96%) family (identity 90-94%) order (identity 80-88%) and class level (identity <80%).

The quality of the BLAST-based taxonomic assignment of the OTUs was assessed using two approaches. First, GenBank was queried using the BLASTn algorithm for each sequence, and the top hits recorded. Blast results were inspected manually to remove inconsistencies: sequences were not considered when there was a contradiction in genus assignment in the reference database, since we suspected misidentification of the database sequence. Second, an alignment was generated for each genus using our query sequences, selected top BLASTn hits and the sequence of the first taxonomically named BLASTn hit, if the top hits were unidentified. These alignments were assembled automatically in CLUSTALW [17] with manual editing. Neighbour joining (NJ) analysis, based on K2P distances and tested by 1000 bootstrap replicates (BS), were conducted for each alignment to confirm the identification of query sequence to the genus or species level, using PHYLIP 3.6 [18]. Bayesian analysis was conducted on the same individual dataset used in the NJ analysis. The best nucleotide substitution model was evaluated by MrModeltest v 2.5 [19]. Phylogenetic analysis was performed with Mrbayes v3.1 [20] applying a general time reversible (GTR) substitution model with gamma (G) and proportion of invariable (I) site parameters to accommodate variable rates across sites. Two independent runs of Markov chain Monte Carlo (MCMC) were performed using four chains for 2x10⁶generation and trees were saved each 1000 generations. The initial 250 trees in each run were discarded as "burn-in". The remaining trees were used to construct majority-rule consensus trees. Bayesian post probabilities (PPs) for each clade were derived from trees remaining after the discarding the burn-in sample. Convergence of runs and parameters was controlled by the diagnostics implemented in MrBayes and in Tracer 1.5 [21]. Trees were drawn using TREEVIEW [22].

To assess sampling efficiency, a rarefaction curve was generated by means of the Analytic Rarefaction 1.3 software [23], using 10 specimens as a step parameter for calculation. Statistical analysis of data was performed with Prism version 5.00 software (GraphPad Software, California, USA) and Systat 11 (Systat Software, Inc.).

Ethics Statement

Specific permissions were not required for these locations and activities and the field studies
did not involve endangered or protected species.

Results

Characterization of the fungal community of leaves

No symptoms or signs were observed on sentinel oaks in 2008–2009 surveys. However during the 2010 growing season a powdery mildew infection developed with signs consisting of white, talcum powder-like coating on the leaves of *Q. petraea*. Sometimes the infected leaves became distorted (Fig. 1A). Leaf spots were also detected on the three *Quercus* species (Fig. 1B). Signs of *Leptosphaeria* cankers were observed on *Q. petraea* branches on which perithecia were
<figure><img src="image_4.png" /><figcaption>Fig 1. Signs of powdery mildew on leaves of *Q. petraea* (A); necrotic spots on leaves of *Q. petraea* (B) *Leptosphaeria* sp. perytheciaon branches of *Q. petraea* (C); ascopores of *Leptosphaeria* sp. (D); conidia of *Alternaria* sp. observed on specimens of all oak hosts (E); conidia of *Pestalotiopsis* sp. from conidiomata on necrotic spots of *Q. petraea* leaves (F). Bars 20 µm.</figcaption></figure>

doi:10.1371/journal.pone.0120571.g001

present (Fig. 1C) containing bitunicate asci with typically3-septate ascospores (Fig. 1D). *Alternaria*-like conidia were also observed on all *Quercus* spp. specimens (Fig. 1E). Asexual fruiting bodies containing *Pestalotiopsis*-like conidia were found in leaf spots of *Q. petraea* (Fig. 1F).

After the quality filtering and removal of singletons, considered to be mostly artefacts [14] a total of 14.069 reads were retained for supplementary analyses. Numbers of reads were 746, 2.959, and 10.364 for *Q. ilex*, *Q. suber* and *Q. petraea*, respectively. A total of 106 OTUs distributed among *Q. ilex* (30), *Q. suber* (59) and *Q. Petraea* (93) were considered for further analyses. Four OTUs accounted for the largest fraction of reads (59%) (Fig. 2). Among these, OTU 1 to 3 were ubiquitous on the 3 host species and assigned to *Alternaria* sp., unidentified Basidiomycota and *Epicoccum nigrum* respectively, while OTU4 was specific to *Q. petraea* and assigned to *Erysiphe quercicola* (see below).

The plot of OTUs vs the number of samples resulted in a rarefaction curve that approached a plateau (Fig. 3) indicating that the majority of taxa present in the site were covered by the analysis.
<figure><img src="image_5.png" /><figcaption>**Fig 2.** Distribution of number of reads within the 106 identified operational taxonomic units (OTUs), singletons excluded.</figcaption></figure>

doi:10.1371/journal.pone.0120571.g002

A total of 45 out of 106 OTUs (42%) belonged to Basidiomycota, 59 OTUs (56%) to Ascomycota and only 2 OTUs (2%) to Glomeromycota. Dothideomycetes (19%) and Tremellomyces (26%) were the most frequent classes of Ascomycota and Basidiomycota, respectively.

A total of 9 OTUs have been identified of possible Asiatic origin, based on clustering in NJ trees within groups of sequences of taxa mostly reported from Asia.

Among these species only 2 belong to Basidiomycota:

OTU91, based on the phylogenetic analysis (BS = 88 and PPs = 95) was ascribed to Fungal sp. NW-AN003, belonging to the genus *Cryptococcus* and previously recorded as endophyte of larch in China (NCBI accession number JN809912);

OTU 48 was assigned to *Cryptococcus* sp.; the best blast hits of this read consisted in 8 iso-
lates of *Cryptococcus* sp., from plants leaves in China (NCBI accessions numbers: HQ890368,

<figure><img src="image_6.png" /><figcaption>**Fig 3.** Rarefaction curves describing the observed number of taxa as a function of the number of samples processed.</figcaption></figure>

doi:10.1371/journal.pone.0120571.g003
HQ890367, HQ890366, HQ890365, HQ890364, HQ890363, HQ890362, HQ890369) (BS = 99 and PPs = 99);

Six OTUs were assigned to Ascomycota

OTU88 and OTU38 were ascribed to the genus *Mycosphaerella* (Fig. 4):

OTU88, was detected on *Q. petraea* and *Q. suber*. Itforms a monophyletic clade with Cf. *Ramichloridium* (NCBI accession number GU214688), *Mycosphaerella* sp. M16 (NCBI accession number HM595520) and *Mycosphaerella* sp.C 2–2 (NCBI accession number JX290142) (BS = 100 and PPs = 100), isolated from *Phellodendron amurense* in South Korea *Abies beshan-zuensis* in China and *Toona sinensis* in Korea, respectively [24], [25]. OTU38 was detected on *Q. ilex*, representing 3.2% of the total reads, and, at lesser extent, on *Q. petraea* and *Q. suber*. It has been identified as *Mycosphaerella* sp. Ponipodef 06 (NCBI accession number HQ731642) isolated from *Populus nigra* x *P. deltoides* in China (BS = 88 and PPs = 60).

OTU25 was detected on *Q. petraea* and *Q. suber* (and other hosts, data not shown) and as-
signed to *Strelitziana albiziae* (NCBI accession number HQ599584) isolated from leaves of
*Albizia julibrissin* infected with *Camptomeris albiziae* in South Korea [26].

<figure><img src="image_7.png" /><figcaption>**Fig 4. Phylogenetic relationship between OTU88, OTU38 and top BLAST hits species (with GenBank accession numbers and geographic distribution).** The tree was constructed using the neighbor-joining cluster algorithm. Bootstrap support (%) and posterior probabilities are provided on the branches (BS/PP). Bootstrap values below 50 are not shown. Scale bar, 0.002 substitutions/site.</figcaption></figure>

doi:10.1371/journal.pone.0120571.g004
<figure><img src="image_8.png" /><figcaption>**Fig 5. Phylogenetic relationship between OTU102 and top BLAST hits species (with GenBank accession numbers and geographic distribution).**
The tree was constructed using the neighbor-joining cluster algorithm. Bootstrap support (%) and posterior probabilities are given on the branches (BS/PP). Bootstrap values below 50 are not shown. Scale bar, 0.002 substitutions/site.</figcaption></figure>

doi:10.1371/journal.pone.0120571.g005

OTU102 was detected on *Q. petraea* and identified as *Leptosphaeria* sp. It formed a well-supported clade (BS = 93 and PPs = 97) including only isolates from China (NCBI accession numbers:HM537061,KF143798,FI537121,KF143787) (Fig. 5).

The blast and the phylogenic analysis assigned OTU4 to *Erysiphe quercicola*. No additional OTUs related to Erysiphales have been identified. From the phylogenetic analysis it is clear that OTU4 is closely related to two isolates from Asia (Thailand and Japan). The division between OTU4 in one clade and the isolate from France on another clade is well supported (BB = 66 and PPs = 81) (Fig. 6).

OTU77 showed 97% (E-value = 3e-95) pairwise similarity to *Lalaria* sp. JS-42 (NCBI accession number JF706660) isolated from flowers in Korea (BB = 100 and PPs = 97).

In addition, OTU16 was assigned to *Pestalotiopsis* sp. (97% similarity; E-value = 3e-100). It accounts for about 2.5% of the total reads on *Q. petraea*. Although the Asiatic origin of this OTU cannot be demonstrated, it probably refers to a new species.
<figure><img src="image_9.png" /><figcaption>**Fig 6.** Phylogenetic relationship between OTU4 and top BLAST hits species (with GenBank accession numbers and geographic distribution). The tree was constructed using the neighbor-joining cluster algorithm. Bootstrap support (%) and posterior probabilities are given on the branches (BS/PP). Bootstrap values below 50 are not shown. Scale bar, 0.002 substitutions/site.</figcaption></figure>

doi:10.1371/journal.pone.0120571.g006
<figure><img src="image_10.png" /><figcaption>**Fig 7.** Hierarchically clustered heat map of the taxa distribution with Asian origin. The relationship among samples was determined by Bray-Curtis distance. The heat map plot depicts the relative percentage of each fungal taxa (variables clustering on the vertical-axis) within each tree host (horizon-axis clustering). The relative abundance for fungal taxa are indicated by color intensity with the legend indicated at the bottom of the figure.</figcaption></figure>

doi:10.1371/journal.pone.0120571.g007

A graphical representation of OTUs abundance and distribution among hosts is showed in Fig. 7.

In Table 1, the information about identification of taxa and association with signs and symptoms are summarized.

<table><caption>**Table 1.** Summary table on identification and association to signs and symptoms of the OTUs of possible Asiatic origin.</caption><thead><tr><th rowspan="2">N°</th><th rowspan="2">OTUs</th><th rowspan="2">Taxa</th><th colspan="2">Identification</th><th rowspan="2">Association with signs and symptoms</th></tr><tr><th>Molecular</th><th>Morphological</th></tr></thead><tbody><tr><td>1</td><td>4</td><td>Erysiphe quercicola</td><td>YES</td><td>YES</td><td>YES</td></tr><tr><td>2</td><td>16</td><td>Pestalotiopsis sp.</td><td>YES</td><td>YES</td><td>YES</td></tr><tr><td>3</td><td>25</td><td>Strelitziana albiziae</td><td>YES</td><td>NO</td><td>NO</td></tr><tr><td>4</td><td>38</td><td>Mycosphaerella sp. 1</td><td>YES</td><td>NO</td><td>NO</td></tr><tr><td>5</td><td>48</td><td>Cryptococcus sp.1</td><td>YES</td><td>NO</td><td>NO</td></tr><tr><td>6</td><td>77</td><td>Lalaria sp.</td><td>YES</td><td>NO</td><td>NO</td></tr><tr><td>7</td><td>88</td><td>Mycosphaerella sp.</td><td>YES</td><td>NO</td><td>NO</td></tr><tr><td>8</td><td>91</td><td>Cryptococcus sp.2</td><td>YES</td><td>NO</td><td>NO</td></tr><tr><td>9</td><td>102</td><td>Leptosphaeria sp.</td><td>YES</td><td>YES</td><td>YES</td></tr></tbody></table>

doi:10.1371/journal.pone.0120571.t001
Discussion

Our study reports for the first time the colonization of a selected collection of EU trees by na-
tive fungi while growing in China, and the potential of sentinel plantations to detect taxa path-
ogenic to EU trees that might pose a biosecurity risk for Europe. The sentinel plantation of
Fuyang represents an improvement of the ISPN hypothesis based on the use of existing botani-
cal gardens worldwide. Infact, in this study, the EU trees were selected according to ecological
criteria. They were planted in China in sites comparable in climate and plant community com-
position to those of their area of origin in Europe. Presence of native species belonging to the
same family or genus (e.g. the *Fagaceae* family in the Fuyang plot) and having similar ecologi-
cal requirements enhances the probability of exposure of sentinel trees to natural inoculum
pressure. The placing of the European trees in these specific sites was thought to favours "host
jump" by native pathogens facilitating the evaluation of their aggressiveness.

Four out of the nine taxa of Asiatic origin, are related to plant pathogenic genera and were
associated with symptoms and signs typical of infection on the sentinel trees.

The OTU 38 and 88 refer to the *Mycosphaerella* complex (Fig. 4) and specifically to fungi isolated from woody hosts other than *Fagaceae* in Asia. Specifically OTU88 refers to taxa isolated from *Abies beshanzuensis* and *Toona sinensis* host species present in the Fuyang forest areas, and plantations (*Phellodendron amurense*). *Phellodendron amurense* is native of North-eastern provinces of China but it is widely cultivated in the Zhejiang provinces. Furthermore, this host is reported to be impacted by at least 3 anamorphic *Mycosphaerellaceae* in China [27] although no sequence data are available for comparison. OTU38 refers to taxa isolated from leaves of *Populus x Euroamericana* also cultivated in the Zhejiang province. These 2 OTUs were identified from leaf spots of sentinel *Quercus petraea*, *Q. suber* and *Q. ilex* plants. Additional OTUs (37 and 44) referring to the *Mycosphaerella* anamorphs *Cercospora* and *Pseudocercosporella* were identified and also associated to leaf spots of sentinel oak trees. However their Asiatic origin could not be proved. *Mycosphaerella* spp. are among the most common and destructive plant pathogens known, causing considerable economic losses on a wide variety of host plants worldwide, including economically important horticultural and ornamental crops [28], [29]. The term *Mycosphaerella* is linked to more than 30 anamorphic form genera and account for several thousand species (ca 10.000) most of which are plant pathogens [24]. Interestingly, Crous [24] in his review on *Mycosphaerella* complex, argues that most of the investigated host plants are associated to several undescribed taxa in the *Mycosphaerella* complex. Many of these species can colonise the same lesion, have the ability to appear in non-host tissue, and to actively undergo host shifts while infecting plant tissue. In such context *Mycosphaerella* spp. are perfect candidates to become invasive in new environments on new hosts. In addition *Quercus* is already reported as a common host of several *Mycosphaerella* species; Gilman and Wadley [30] listed at least 8 species pathogenic on oaks in the USA.

OTU 102 refers most probably to a Chinese *Leptosphaeria* sp. as supported by phylogenetic
analysis (Fig. 5). The closest sequences were assigned to an undescribed *Leptosphaeria* sp. iso-
lated from plant specimens in China [31] The genus *Leptosphaeria* includes a large number of
species mainly considered saprobic or necrotrophic on stems and leaves. The genus also com-
prises highly pathogenic species such as *L. maculans* and *L. biglobosa*, the causal agents of
‘blackleg’ of *Brassicae* crops [32], and *L. coniothyrium* causal agent of stem canker of raspberry
and roses [33]. The association of typical symptoms and signs with molecular detection by
mass sequencing suggests the interaction of *Leptosphaeria* sp. with *Q. petraea*, although the
pathogenic behaviour must be demonstrated.

The finding of the biotrophic pathogen *E. quercicola* on *Q. petraea* is of interest because this
species has been described from *Q. phillyraeoides* by Takamatsu et al. [34]. It was described as a
new species genetically and morphologically distinct from *E. alphitoides* and having an Asiatic and Australian distribution. Most recently *E. quercicola* has been detected in France and associated with flag-shoot symptoms of *Q. robur* and *Q. petraea* [35]. Mougou et al. [36] suggested that *E. quercicola*, similar to other oak powdery mildews in Europe, was introduced on European oaks through host shift from traded exotic host species such as *Q. phillyraeoides* (present in many nursery catalogues in France). Furthermore Kirshner and Liu [37] reported anamorphic *E. quercicola* on 2 tropical hosts commonly traded to Europe as ornamentals, *Cinnamomum camphora* and *Murraya paniculata* highlighting the risk of host shift from taxonomic distant species. These native hosts, including *Q. phillyraeoides*, are common in the Fuyang area supporting the presence of source of natural inoculum of *E. quercicola*. In the sentinel plantation studied, signs and symptoms of powdery mildew were observed in 2010 (third year from establishment of the plantation) on *Q. petraea* but not on *Q. suber* and *Q. ilex*. The molecular detection confirmed the presence of *E. quercicola* limited to *Q. petraea* (1214 reads). Moreover, the phylogenetic analysis confirms that OTU4 is clearly distinguished from the sequence detected in France (Fig. 5). The recovery of *E. quercicola* on *Q. petraea* in the Fuyang sentinel plantation represents an "a posteriori" demonstration of the efficacy of such an approach to detect plant pathogens cable of invading and causing negative impacts to European forests.

The detection of *Pestalotiopsis* sp. is particularly interesting. While no specific signs of *Mycosphaerella* (reproductive structures, spores) have been reported from specimen observations, conidiomata of *Pestalotiopsis* sp. and associated conidia were observed on *Q. petraea* leaf spots. Species of *Pestalotiopsis* commonly cause diseases on a variety of plants. They are also frequently isolated as endophytes or occur as saprobes [38]. In such context the presence of conidiomata of *Pestalotiopsis* on necrotic tissues does not necessarily demonstrate the pathogenic behaviour of this taxa to *Q. petraea*. However, the low value of identity of OTU16 (97%) with the sequences of *Pestalotiopsis* spp. in the NCBI database and its considerable presence on *Q. petraea* (2.5% of the total reads) support the hypothesis of a new species not previously recorded.

Although the potential of sentinel plantations to intercept local plant pathogens seems high, several bottlenecks need to be resolved before this approach can be routinely utilized. The Fuyang plantation was established with 1–2 year old seedlings shipped directly from Europe. This sentinel plantation was primarily designed to facilitate the detection of damaging pests (insects). In this context, the choice of propagation material for the plantation dealt only with the need of having medium sized trees in place for insect pests attack. However, the type of propagation material gains importance when fungal (or other plant pathogens) detection is considered. The use of seedlings imported from Europe to establish the sentinel plantation does not eliminate the risk of their colonization by fungal taxa at the area of production. To circumvent this, the use of certified seeds to establish the plantation is recommended in order to minimize the risk of endophytic infection by European taxa. In any case, a characterization of possible microbial community residents in the propagation material is mandatory before the plantation is established to provide a baseline for colonization by native organisms in the sentinel plantation. In this study the detection of putative Chinese fungi infecting the European sentinel trees was performed by symptom/sign assessment and molecular methods. Both approaches are appropriate and necessary, but certainly not sufficient to prove pathogenicity, especially for taxa behaving both as saprophyte and pathogens. While the assessment of specific signs and symptoms coupled with the molecular detection are exhaustive for biotrophic pathogens, like *E. quercicola*, the presence of different *Mycosphaerella* spp. associated with leaf spots could be explained by both a pathogenic and saprophytic interaction. Biological characterization of the organisms associated with the disease and the reproduction of symptoms through inoculation trials on healthy hosts is required to prove pathogenicity and should
always follow the molecular characterization. Obviously, sentinel plant material suspected to be infected by native plant pathogens cannot be shipped abroad and should be processed in specialised laboratories "in loco". In this study symptomatic material was processed 'in loco' for DNA extraction and purification; unfortunately biological detection and isolation of candidate pathogenic fungi failed and could not be completed due to bottlenecks in collaborative agreements. For the complete success of a sentinel tree strategy, the quality of reciprocal agreement between the collaborating countries including shared protocols for plantation establishment and samples processing is of importance.

In conclusion, the sentinel plantations represent an appropriate answer to the growing need of preventive identification and characterization of plant pathogens at risk of introduction and invasion in Europe (and elsewhere). However the preliminary experience carried out in the Fuyang plantation highlighted the logistic and technical difficulties in utilizing sentinel trees for pathogens detection. Efforts should be made to solve regulatory and logistic bottlenecks associated with biosecurity issues and the possibility of full characterization of biological samples, in order to make possible the routinary and reciprocal utilization of this approach between countries.

## Acknowledgments

This is research was supported by the EU project Increasing Sustainability of European Forests: Modelling for Security Against Invasive Pests and Pathogens under Climate Change (ISEFOR, EU FP7). We also acknowledge the COST Action FP1002 Pathway Evaluation and pest Risk Management In Transport-PERMIT that allowed a critical development of the Sentinel tree concept.

## Author Contributions

Conceived and designed the experiments: AMV AV AR. Performed the experiments: AMV AY JF. Analyzed the data: AMV. Contributed reagents/materials/analysis tools: JF JS AV. Wrote the paper: AMV AV.

## References

1. Santini A, Ghelardini L, De Pace C, Desprez Loustau ML, Capretti P, Chandelier A, et al. Biogeographical patterns and determinants of invasion by forest pathogens in Europe. New Phytologist 2013; 197: 238–250. doi: 10.1111/j.1469-8137.2012.04364.x PMID: 23057437

2. Pautasso M, Gregor A, Queloz V, Holdenrieder O. European ash (Fraxinus excelsior) dieback—A conservation biology challenge. Biological Conservation 2013; 158:37–49.

3. Landeras E, García P, Fernández Y, Braña M, Fernández-Alonso O, Méndez-Lodcs S, et al. Outbreak of pitch canker caused by *Fusarium circinatum* on *Pinus* spp in northern Spain. Plant Disease 2005; 89:1015.

4. Brown AV, Brasier CM. Colonization of tree xylem by *Phytophthora ramorum*, *P. kernoviae* and other *Phytophthora* species. Plant Pathology 2007; 56:227–241.

5. Green S, Brasier CM, Schlenzig A, McCracken A, MacAskill GA, Wilson M, et al. The destructive invasive pathogen *Phytophthora lateralis* found on *Chamaecyparis lawsoniana* across the UK. Forest Pathology 2013; 43:19–28.

6. Brasier CM. The biosecurity threat to the UK and global environment from international trade in plants. Plant Pathology 2008; 57:792–808.

7. Britton KO, White P, Kramer A, Hudler G. A new approach to stopping the spread of invasive insects and pathogens: early detection and rapid response via a global network of sentinel plantings. New Zealand Journal of Forestry Science 2010; 40:109–114.

8. Roques A, Fan JT, Courtial B, Yart A, Sun JH. Sentinel Plantings in Other Continents a Novel Method to Identify Potential Insect Invaders—Final Results of a Chinese Trial Using European Trees. ICBI International Congress on biological invasions October 23–27 2013 Qingdao China. Available: http://www.icbi2013.org/admin/down/201401171009122826.pdf.
9. Roques A, Fan JT, Courtial B, Zhang YZ, Yart A, Auger-Rozenberg MA, et al. Planting sentinel European trees in Asia as a novel method to identify potential insect pest invaders. PLOS ONE. 2015; In press.

10. Kumar S, Carlsen T, Mevik BH, Enger P, Blaalid R, Shalchian-Tabrizi K, et al. CLOTU: an online pipeline for processing and clustering of 454 amplicon reads into OTUs followed by taxonomic annotation. BMC Bioinformatics 2011; 12:182. doi: 10.1186/1471-2105-12-182 PMID: 21599929

11. Bjorbækmo MFM, Carlsen T, Brysting AK, Vrålstad T, Høiland K, Ugland KI, et al. High diversity of root associated fungi in both arctic and alpine *Dryas octopetala* BMC Plant Biology 2010; 10:244. doi: 10.1186/1471-2229-10-244 PMID: 21070665

12. Mohamed DJ, Martiny JBH. Patterns of fungal diversity and composition along a salinity gradient. The ISME Journal 2010; 5:379-388. doi: 10.1038/ismej.2010.137 PMID: 20882058

13. Smith ME, Douhan GW, Rizzo DM. Ectomycorrhizal community structure in a xeric *Quercus* woodland based on rDNA sequence analysis of sporocarps and pooled roots. New Phytologist 2007; 174:847-863. PMID: 17504467

14. Tedersoo L, Nilsson RH, Abarenkov K, Jairus T, Sadam A, Saa I, et al. 454 Pyrosequencing and Sanger sequencing of tropical mycorrhiza fungi provide similar results but reveal substantial methodological biases. New Phytologist 2010; 188:291-301. doi: 10.1111/j.1469-8137.2010.03373.x PMID: 20636324

15. Altschul SF. Evaluating the statistical significance of multiple distinct local alignments. In: Suhai S, editor. Theoretical and Computational Methods in Genome Research. New York: Plenum; 1997. pp. 1-14.

16. Orgiazzi A, Bianciotto V, Bonfante P, Daghino S, Ghignone S, Lazzari A, et al. 454 Pyrosequencing Analysis of Fungal Assemblages from Geographically Distant Disparate Soils Reveals Spatial Patterning and a Core Mycobiome. Diversity 2013; 5:73-98.

17. Thompson JD, Higgins DG, Gibson TJ. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting position-specific gap penalties and weight matrix choice. Nucleic Acids Res 1994; 22:4673-4680. PMID: 7984417

18. Felsenstein J. PHYLIP—phylogeny inference package (version 32). Cladistics 1989; 5:164-166.

19. Nylander JAA. MrModeltest v2 Program distributed by the author. Evolutionary Biology Centre, Uppsala University; 2004.

20. Ronquist F, Huelsenback JP. MrBayes: Bayesian phylogenetic inference under mixed models. Bioinformatics 2003; 19: 1572-1574. PMID: 12912839

21. Rambaut A, Drummond AJ. Tracer v15 2009. Available: http://treebioedacuk/software/tracer/.

22. Page RDM. TreeView version 166 Glasgow: University of Glasgow; 2001. Available: http://taxonomy zoology.gla.ac.uk/ROD/treeview.html.

23. Holland SM. Analytic rarefaction 2003. Available: http://wwwugaedu/strata/software/.

24. Crous PW. Taxonomy and phylogeny of the genus *Mycosphaerella* and its anamorphs. Fungal Diversity 2009; 38: 1-24.

25. Yuan ZL, Rao LB, Chen YC, Zhang CL, Wu YG. From pattern to process: species and functional diversity in fungal endophytes of *Abies beshanzuensis*. Fungal Biol 1 2011; 15: 197-213.

26. Crous PW, Groenewald JZ, Shin H-D. Fungal Planet 56: *Strelitziana albiziae*. Persoonia 2010; 25:132-133.

27. Zheng H. Invasive Plants of Asian Origin Established In the United States and Their Natural Enemies Vol. 2. Forest Service Forest Health Technology Enterprise Team 1; 2005. Available: http://www.invasive.org/weeds/asianv2.pdf.

28. Farr DF, Bills GF, Chamuris GP, Rossman AY. Fungi on plants and plant products in the United States. St Paul, MN: APS Press; 1995.

29. Braun U, Crous PW, Dugan F, Groenewald JZ, Hoog GS. Phylogeny and taxonomy of cladosporium-like hyphomycetes, including *Davidiella* gen. nov., the teleomorph of *Cladosporium* s.str. Mycological Progress 2003; 2: 3-18.

30. Gilman JC, Wadley BN. The Ascigerous Stage of *Septoria querceti* Thuem. Mycologia 1952; 44: 216-220.

31. Zhao ZJ, Pan YZ, Liu QJ, Li XH. Exposure assessment of lovastatin in Pu-erh tea. International Journal of Food Microbiology 2013; 164: 26-31. doi: 10.1016/j.ijfoodmicro.2013.03.018 PMID: 23587710

32. Stonard JF, Latunde-Dada AO, Huang YJ, West JS, Evans N, Fitt BDL. Geographic variation in severity of *Phoma* stem canker and *Leptosphaeria maculans/L. biglobosa* populations on UK winter oilseed rape (*Brassica napus*). European Journal of Plant Pathology 2010; 126: 97-109.
33. Williamson B. Cane blight. In: Ellis MA, Converse RH, Williams RN, Williamson B, editors. Compendium of Raspberry and Blackberry Diseases and Insects. St. Paul, MN: APS Press; 1991. pp. 5-7.

34. Takamatsu S, Braun U, Limkaisang S, Kom-un S, Sato Y, Cunningham JH. Phylogeny and taxonomy of the oak powdery mildew *Erysiphe alphitoides* sensu lato. Mycol Res 2007; 111: 809-826. PMID: 17681225

35. Feau N, Lauron-Moreau A, Piou D, Marçais B, Dutech C, Desprez-Loustau ML. Niche partitioning of genetic lineages involved in the oak powdery mildew complex. Fungal Ecology 2012; 5: 154-162.

36. Mougou A, Dutech C, Desprez-Loustau ML. New insights into the identity and origin of the causal agent of oak powdery mildew in Europe. Forest Pathology 2008; 38: 275-287.

37. Kirschner R, Liu WA. Two new hosts of anamorphic *Erysiphe quercicola*: *Cinnamomum camphora* and *Murraya paniculata*. *Mycoscience* 2014; 55: 190-195.

38. Jeewon R, Liew ECY, Hyde KD. Phylogenetic evaluation of species nomenclature of *Pestalotiopsis* in relation to host association. Fungal Diversity 2004; 17: 39-55.