Introduction

Hazelnut (Corylus avellana L.) is a perennial bushy plant belonging to Betulaceae native to Europe and western Asia. It represents one of the most economically important nut crops worldwide. The Mediterranean and Black sea areas are historical sites of hazelnut cultivation with Turkey covering more than 60% of the world’s production, followed by Italy, USA, Azerbaijan, Georgia, Chile, China, Iran, France and Spain (FAOSTAT 2023). In Italy, 84,526 ha are cultivated with hazelnut and 98,666 tonnes of fruit were harvested in 2022 (ISTAT 2023). Italian hazelnut production has always been present in few geographical areas that are specialized for their environmental and climatic conditions, for the developed technical knowledge and socio-economic context and for the adoption of high quality cultivars (Botta and Valentini 2018). The regions with the highest production are Piedmont (30.5%), Latium (28.7%), Campania (22.5%) and Sicily (12.8%). In Piedmont, hazelnut cultivation covers more than 27,000 ha and the majority (16,558 ha) is located in the Cuneo province (ISTAT 2023). In this region, hazelnut production is mainly based on the cultivar Tonda Gentile (synonyms: ‘Tonda Gentile delle Langhe’, ‘Tonda trilobata’). This cultivar is recognized as one of the best hazelnut cultivar worldwide for its taste and aroma appreciated for fresh consumption and, in particular, for the industry, due to the presence of major international companies specialized in confectionery production (Silvestri et al. 2021). During the last decade, a general expansion in terms of acreage (from 14,375 in 2013 to 24,701 ha in 2022) and yield (from 23,797 in 2013 to 30,180 tons in 2022) in hazelnut cultivation was observed in this region (ISTAT 2023). However, the trend in nut production is highly influenced by climatic conditions (An et al. 2020). In the last 2 years, a slight decrease was observed probably due to water deficiency, late-spring frosts and hot summer temperatures. Moreover, a considerable spread of cankers and twig dieback symptoms was observed in hazelnut orchards in Piedmont increasing the concern of the producers.

Cankers, blight, dieback and decay affecting trunk, branches, twigs and, in severe cases, the entire plant, represent a serious threat for hazelnut worldwide (Teviotdale et al. 2002). In the present study, we will include under the general name of fungal trunk diseases (FTD) the diseases caused by fungal pathogens characterized by the mentioned symptoms and that were reported in the main hazelnut growing countries across the world. In the USA, the most important hazelnut growing area is the Willamette valley in Oregon, where the eastern filbert blight (EFB) caused by Anisogramma anomala represented a severe problem (Pinkerton 1992). The adoption of resistant cultivars provided an effective strategy to control EFB disease (Johnson et al. 1996). However, the emergence of previously unknown hazelnut disease is observed and previously unassociated fungal pathogens, such as Diplodia mutila, Diaporthe eres, Dothiorella omnivora and Valsa cf. eucalypti were isolated from cankers on trunks and branches (Wiman et al. 2019). In Chile, Diplodia coryli, Dip. mutila and Diaporthe australafricana were reported as causal agent of hazelnut dieback and cankers (Guerrero and Pérez 2013a, b; Guerrero et al. 2014). Dothiorella parva was isolated from Corylus avellana in Iran and described as a new species (Abdollahzadeh et al. 2014). Studies conducted in the Guilan province, the main cultivation area of Iran, reported Botryosphaeria dothidea and Diplodia theobromae as pathogen on hazelnut plants and Dia. amygdali as causal agent of hazelnut trees decline (Mir Hosseini Moghaddam and Taherzadeh 2007; Mohammadi and Jabbari Firoozjah 2019; Ghasemi-Doodaran and Davari 2020). Recently, Cytospora sp., Phomopsis sp., Lasiodiplodia sp. and Pestalotiopsis sp. were found in association with dieback of hazelnut in the same region (Houshyarfard 2020). In Turkey, Botryosphaeria dothidea was reported as pathogen on hazelnut plants (Polat et al. 2022). Moreover, recent studies conducted on Corylus heterophylla in China found Dia. eres, Dia. corylicola, Dia. donglingensis and Dia. huairouensis as FTD causal agents (Gao et al. 2021; Bai et al. 2022). In Italy, different studies were carried out in the main hazelnut growing areas across the peninsula. The presence of Cytospora canker, caused by Cytospora corylicola and known as ‘Mal dello stacco’, was reported throughout the most relevant production regions with a significant yield decrease in Campania (Botta and Valentini 2018). Although the spread of C. corylicola, there is a lack of information about the biology and phylogeny of this species and its pathogenicity was long debated (Scortichini 2006; Lamichhane et al. 2014). In Piedmont, sporadic branch cankers were observed in hazelnut orchards and the causal agent were identified as belonging to the genera Phomopsis and Sphaeropsis (Botta and Valentini 2018). A recent study was conducted in Sardinia on twigs and branches of hazelnut trees with exudates and cankers and different fungal species were isolated form symptomatic plant materials and identified: Dip. sapinea, Dip. seriata, Dothiorella iberica, Do. omnivora, Do. parva, Do. symphoricarposicola, Diaporthella cryptica, Gnomoniopsis smithogilvyi and Anthostoma decipiens (Linaldeddu et al. 2016). The same study provided a clarification on the etiology of Cytospora canker on hazelnut caused by A. decipiens (Linaldeddu et al. 2016).

In Piedmont, the knowledge on FTD of hazelnut is still limited and no studies were recently conducted to investigate the etiology and epidemiology of this disease. Thus, considering the diversity of the symptoms observed in the field and the relevant economic value of this crop, this study was conducted with the aim to investigate the etiology of hazelnut FTD in this region, in detail: (i) to identify the fungal species in association with FTD of hazelnut trees using molecular tools and phylogeny, (ii) to assess the morphological features of the identified species and (iii) to test the pathogenicity of the species found and to fulfil Koch’s postulates.

Materials and methods

Field survey and fungal isolation

Field surveys were conducted from March 2021 to September 2022 in seven hazelnut orchards in Piedmont (Table 1). Samples of trunks, branches and twigs were collected from symptomatic plants of ‘Tonda Gentile’ showing Cytospora cankers (‘Mal dello stacco’) and dieback symptoms. Wood samples (5–10 mm) were surface sterilised in 1% sodium hypochlorite for 1 min, rinsed in sterile distilled water (SDW) for 1 min and dried on sterile absorbent paper. Small fragments (2–3 mm) were cut from the edge between healthy and necrotic tissues and plated on potato dextrose agar (PDA, VWR Chemicals, Leuven, Belgium) amended with 25 mg l−1 of streptomycin sulphate (PDA-S, AppliChem GmbH, Darmstadt, Germany). The plates were incubated at 25 ± 1 °C under a 12 h photoperiod. Following 48 to 72 h of incubation, mycelial plugs were taken from the margin of developing colony and placed on new PDA-S and water agar (WA, Microbiol Diagnostici, Cagliari, Italy) plates. After 4 to 5 days, pure cultures were established by transferring single hyphal tips on new PDA-S plates. A total of 35 isolates were obtained and used for characterization. Stock cultures of these isolates are kept at −80 °C in the culture collection of the University of Torino, Italy.

Table 1 Coordinates and locations of the orchards surveyed in this study

DNA extraction, polymerase chain reaction (PCR) amplification, and sequencing

Total genomic DNA was extracted from 0.1 g of mycelium grown on PDA-S, using the E.Z.N.A. Fungal DNA Mini Kit (Omega Bio-Tek), following the manufacturer’s instructions. Four different genomic loci were amplified and sequenced for species identification. The nuclear ribosomal internal transcribed spacer (ITS) region of each isolate was amplified using the universal primers ITS1 and ITS4 (White et al. 1990). The primers EF1-728F and EF1-986R (Carbone and Kohn 1999) were used to amplify partial region of translation elongation factor-1α (tef1) in isolates of Botryosphaeriaceae and Diaporthe species (Guarnaccia et al. 2020, 2022a; Aiello et al. 2022). The partial β-tubulin (tub2) gene was amplified with primers T1–Bt2b (Glass and Donaldson 1995; O’Donnell and Cigelnik 1997) for isolates belonging to Botryosphaeriaceae family and Diaporthe genus (Guarnaccia et al. 2020, 2022a), whilst primers Tub2fd–Tub4fd (Woudenberg et al. 2009) were used to amplify the same region in isolates of the genus Nothophoma (Chen et al. 2015). The partial RNA polymerase second largest subunit (rpb2) gene was amplified with the primers: Rpb2-5f2–Rpb2-7cr (Liu et al. 1999; Reeb et al. 2004) for isolates identified as Nothophoma sp. (Chen et al. 2015). PCR mixtures and conditions were followed as described in the above-cited references. PCR amplification was examined by electrophoresis on 1% agarose (VWR Life Science AMRESCO® biochemicals) gels stained with GelRedTM. Eurofins Genomics Service (Cologne, Germany) sequenced PCR products in both directions. The DNA sequences generated were analysed and consensus sequences were computed using the program Geneious v. 11.1.5 (Auckland, New Zealand).

Phylogenetic analyses

The sequences obtained were blasted against the NCBI’s GenBank nucleotide database to determine the closest relatives for a taxonomic framework of the studied isolates. Alignments of different genomic regions, including sequences obtained from this study and sequences downloaded from GenBank, were initially performed with the MAFFT v. 7 online server (http://mafft.cbrc.jp/alignment/server/index.html) (Katoh and Standley 2013), and then manually adjusted in MEGA v. 7 (Kumar et al. 2016). Analyses were conducted individually for each locus (data not shown) and as multilocus sequence analyses using the following loci combinations: ITS, tub2 and tef1 for members of Botryosphaeriaceae and Diaporthe spp. (Guarnaccia et al. 2020; Zhang et al. 2021) and ITS, tub2 and rpb2 for Nothophoma spp. (Chen et al. 2015). For isolates belonging to Anthostoma a single locus analysis was performed on ITS region (Linaldeddu et al. 2016). Lasiodiplodia theobromae (CBS 164.96) (Zhang et al. 2021) was used as outgroup for species belonging to Botryosphaeriaceae. Diaporthella corylina (CBS 121124) (Guarnaccia et al. 2020) was selected as outgroup for Diaporthe spp. Allophoma minor (CBS 325.82) (Chen et al. 2015) was used as outgroup for Nothophoma spp. and Cryptovalsa ampelina (STEU 8113) (Moyo et al. 2018) was selected as outgroup for Anthostoma sp. Phylogenies were based on Bayesian Inference (BI) and Maximum Parsimony (MP) analyses. Regarding BI, the best evolutionary model for each partition was determined with MrModeltest v. 2.3 (Nylander 2004) and incorporated in the analyses. MrBayes v. 3.2.5 (Ronquist et al. 2012) was used to generate phylogenetic trees under optimal criteria per partition. The Markov Chain Monte Carlo (MCMC) analysis used four chains and started from a random tree topology. The heating parameter was set with the value of 0.2 and trees were sampled every 1000 generations. Analyses stopped at the moment which the average standard deviation of split frequencies was below 0.01. The MP analyses were performed using PAUP (Phylogenetic Analysis Using Parsimony, v. 4.0b10) (Swofford 2003). Phylogenetic relationships were estimated by heuristic searches with 100 random addition sequences. Tree bisection-reconnection was used, with the branch swapping option set on ‘best trees’ only with all characters weighted equally and alignment gaps treated as fifth state. Tree length (TL), consistency index (CI), retention index (RI) and rescaled consistence index (RC) were calculated for parsimony and the bootstrap analyses (Hillis and Bull 1993) were based on 1000 replications.

Morphology

Based on molecular characterization, representative isolates were selected to assess their morphological features. Agar plugs (5-mm-diam) were taken from the edge of actively growing cultures and transferred to the center of 9-cm-diam Petri dishes containing PDA-S. Isolates belonging to Botryosphaeriaceae were transferred onto the centre of 9 cm diam Petri dishes containing 2% water agar supplemented with sterile pine needles (PNA) (Smith et al. 1996) to induce sporulation. Isolates belonging to Diaporthe were placed both on PNA and on malt extract agar (MEA; Merck, Darmstadt, Germany) to induce sporulation. Plates were then incubated at 25 ± 1 °C under a 12 h photoperiod. Colony characters were observed after 7 days and culture colours were determined according to Rayner (1970). Cultures were examined periodically for the development of conidiomata. Conidia characteristics were examined by mounting fungal structures in SDW and the length and width of 30 conidia were measured for each isolate using an optic microscope at ×40 magnification. The average and standard deviation were calculated.

Pathogenicity

The pathogenicity of representative isolates of the identified species was determined to fulfil Koch’s postulates. Two representative isolates of each species were used to inoculate 1-year-old healthy hazelnut plants of ‘Tonda Gentile’. Seven plants per fungal isolate were inoculated. For each plant, one inoculation point was considered. The inoculation area was surface disinfected with 70% ethanol solution. A sterile scalpel was used to remove the outer bark to expose the vascular tissues. Mycelium plugs (4 mm diam.) were taken from 7-days-old cultures on PDA-S and placed with the mycelium in contact with the internal plant tissues. The same number of plants were treated with sterile PDA-S plugs as controls. The inoculation points were sealed with Parafilm®. All the plants were placed in a completely randomized design under a shade canopy for 3 months. After this period, the bark was removed, and the internal lesion length was measured. Small wood pieces (2–3 mm) of symptomatic tissue from the margin of the lesions were surface disinfected and placed on PDA-S to re-isolate the inoculated fungal species to fulfill Koch’s postulate. The identification of the pathogens was confirmed by colony characteristics. The data obtained were tested for normality, homogeneity of variances, and residual patterns. Logarithmic transformation of the data was performed. ANOVA was conducted with lesion length as dependent variable and fungal isolates as independent variable. Treatment means of lesion length were compared according to Tukey’s HSD test at α = 0.05. The data analysis was conducted using SPSS software 26 (IBM Corporate).

Results

Field survey and fungal isolation

In all the surveyed orchards (Table 1), with the exception of orchard 3, dieback symptoms, observed with a variable disease incidence from 30 to 40%, seriously compromised entire portions of the tree canopy, mainly in the higher part (Fig. 1a-b). Longitudinal streaks were observed on the bark (Fig. 1c-d) and internal dark brown vascular necrosis were found looking at transversal section (Fig. 1e) and, after bark removal, at longitudinal section (Fig. 1f-g). In four out of seven surveyed orchards (orchard 1, 2, 3, and 6), symptoms of Cytospora canker, known as ‘Mal dello stacco’, were observed with a disease incidence from 10 to 20%. Typical reddish-brown irregular spots and elongated cankers on the bark were found leading, in severe cases, to the break of branches (Fig. 1h-i). Oozing conidia in a reddish-brown matrix were observed on the outer surface of the trunk in late summer (Fig. 1j). Tissues under the bark showed necrotic dark brown vascular discoloration (Fig. 1k). Both the described symptoms were observed more frequently on mature plants (>10 year-old). Fungal isolates obtained from symptomatic plant samples were grouped as Botryosphaeriaceae, Diaporthe, Cytospora-like and Phoma-like species according to their culture characteristics. The species recovered from each orchard are reported in Table 1.

Fig. 1
figure 1

Symptoms of FTD of hazelnut trees observed in the field. a, b Dieback symptoms affecting the higher part of the tree canopy. c, d Longitudinal streaks on the bark. eg Internal dark brown vascular necrosis. h, i Cytospora canker (‘Mal dello stacco’) symptoms causing the break of branches. j Oozing conidia in a reddish-brown matrix. k Necrotic dark brown internal vascular discoloration

Phylogenetic analyses

Sequences generated in this study are deposited in GenBank (Table 2). The combined-locus phylogeny of Botryosphaeriaceae consisted of 79 sequences, including the outgroup, and a total of 1444 characters (ITS: 1–549, tef1: 556–958 and tub2: 965–1444) was included. A total of 31 sequences, including the outgroup, was included in the Diaporthe phylogenetic analyses that consisted of 1690 characters (ITS: 1–612, tef1: 619–1036 and tub2: 1043–1690). The analyses conducted on the Nothophoma group consisted of 25 sequences, including the outgroup, with a total of 1437 characters (ITS: 1–493, tub2: 500–835 and rpb2: 842–1437). The single-locus phylogeny of Anthostoma based on ITS consisted of 28 taxa, including the outgroup, and it consisted of 573 characters. A maximum of 1000 equally most parsimonious trees were saved and characteristics of the combined gene partitions used for each phylogenetic analysis are reported in Table 3. Bootstrap support values from the Maximum Parsimony (MP) analysis were plotted on the Bayesian Inference (BI) phylogenies presented in Figs. 2 to 4. For the BI analyses, the models recommended by MrModeltest are reported in Table 4. Unique site patterns for each partition and all the parameters of the Bayesian analyses are reported in Table 3. In the Botryosphaeriaceae analyses, two isolates clustered with the epitype and one reference isolate of Botryosphaeria. dothidea, four isolates clustered with the epitype and four reference isolates of Diplodia seriata, two isolates clustered with the epitype and one reference isolate of Dip. subglobosa and two isolates clustered with the ex-type and two reference isolates of Dothiorella parva (Fig. 2). Regarding phylogenies of Diaporthe, two isolates clustered with five reference isolates of Dia. rudis, whilst four isolates clustered with four reference isolates of Dia. eres (Fig. 3). The final tree generated for Nothophoma showed that six isolated clustered with the ex-type and two reference isolates of N. brennandiae (Fig. 4). Thirteen isolates clustered with three representative sequences of A. decipiens (Fig. 5).

Table 2 Information on fungal isolates used in the phylogenetic analyses and their corresponding GenBank accession numbers
Table 3 Parsimony and Bayesian analyses characteristics in this study
Fig. 2
figure 2

Consensus phylogram of 3462 trees resulting from a Bayesian analysis of the combined ITS, tef1 and tub2 sequences of Botryosphaeriaceae isolates. Bayesian posterior probability values and bootstrap support values are indicated at the nodes. The tree was rooted to Lasiodiplodia theobromae (CBS 164.96). Isolates from the current study are indicated in red

Fig. 4
figure 3

Consensus phylogram of 1922 trees resulting from a Bayesian analysis of the combined ITS, tub2 and rpb2 sequences of Nothophoma isolates. Bayesian posterior probability values and bootstrap support values are indicated at the nodes. The tree was rooted to Allophoma minor (CBS 325.82). Isolates from the current study are indicated in red

Table 4 Evolutionary models as determined by MrModeltest (Nylander 2004)
Fig. 3
figure 4

Consensus phylogram of 1002 trees resulting from a Bayesian analysis of the combined ITS, tef1 and tub2 sequences of Diaporthe isolates. Bayesian posterior probability values and bootstrap support values are indicated at the nodes. The tree was rooted to Diaporthella corylina (CBS 121124). Isolates from the current study are indicated in red

Fig. 5
figure 5

Consensus phylogram of 1732 trees resulting from a Bayesian analysis of the ITS sequences of Anthostoma isolates. Bayesian posterior probability values and bootstrap support values are indicated at the nodes. The tree was rooted to Cryptovalsa ampelina (STEU 8113). Isolates from the current study are indicated in red

Morphology

Morphological features, supported by phylogenetic inference, were observed and used to describe the eight known species (Fig. 6). Colony characters and colours were observed on plates of PDA-S. Botryosphaeriaceae isolates were characterized by the presence of abundant fast-growing aerial mycelium that covered the entire PDA-S petri dishes after 7 days. Botryosphaeria dothidea colonies were white to pale grey. Colony reverse color was pale grey. Conidia were hyaline, aseptate, thin-walled, fusiform to subclavate, with dimensions of 20.7–28.9 × 4.6–6.7 μm, mean ± SD = 24.8 ± 2 × 5.6 ± 0.6 μm. Colonies of Dip. seriata were light grey on the front side and light grey turning dark in the centre on the reverse side. Conidia were light brown, ovoid with truncated or rounded base and obtuse apex, aseptate, with dimensions of 20.2–26.8 × 8.2–11.2 μm, mean ± SD = 23 ± 1.6 × 9.5 ± 0.8 μm. Colonies of Dip. subglobosa were light grey turning dark grey in the center both on the obverse and reverse sides. Conidia were hyaline, aseptate, smooth, thick-walled, oblong to ovoid, straight, ends broadly rounded and measured 26.9–33.2 × 13.3–17.1 μm, mean ± SD = 30.2 ± 1.3 × 15.4 ± 0.9 μm. Dothiorella parva colonies were pale olivaceous grey on the obverse side and dull green- to dark olivaceous- grey on the reverse side. Conidia were ellipsoid to ovoid, brown, 1-septate, moderately thick-walled, ends rounded, often with a truncate base, with dimensions of 16.3–21.5 × 7.5–10.3 μm, mean ± SD = 20.1 ± 1 × 8.9 ± 0.7 μm. Isolates of Diaporthe showed fluffy aerial mycelium. Diaporthe eres colonies were white to pale grey on the front side and pale grey turning dark in the center on the reverse side. Alpha conidia were aseptate, hyaline, smooth, ovoid to ellipsoid, guttulate and measured 5.4–7.5 × 1.5–2.9 μm, mean ± SD = 6.3 ± 0.6 × 2.2 ± 0.3 μm. Beta conidia were hyaline, aseptate, smooth, spindle shaped, slightly curved, measuring 20.4–30.3 × 1.0–1.4 μm, mean ± SD = 25.1 ± 2.5 × 1.2 ± 0.1. Diaporthe rudis colonies were white to beige with a brownish halo around the margin on the front side and pale beige in the margin to buff honey in the center on the reverse side. Alpha conidia were hyaline, aseptate, smooth, biguttulate and ellipsoid with subtruncate bases and measured 4.2–6.9 × 1.4–2.6 μm, mean ± SD = 5.7 ± 0.7 × 1.9 ± 0.3 μm. Beta conidia were not observed. Colonies of N. brennandiae were characterized by cottony moderate aerial mycelium dark brick to sepia on the obverse side and dark brick to cinnamon on the reverse. Conidia were hyaline becoming brown, ellipsoidal to oblong, straight and measured 3.9–6.2 × 2–4.3 μm, mean ± SD = 5.1 ± 0.6 × 3.2 ± 0.5 μm. Anthostoma decipiens colonies showed cottony scarce aerial mycelium white to pale grey on the front side and honey in the margin to buff in the center on the reverse side. Conidia were hyaline, lunate and unicellular measuring 5.9–9.7 × 1.2–2.4 μm, mean ± SD = 7.4 ± 0.8 × 1.7 ± 0.3 μm.

Fig. 6
figure 6figure 6

Morphological characteristics of the different fungal species grown 7 days on PDA-S. a, b, c Botryosphaeria dothidea. d, e, f Diplodia seriata. g, h, i Diplodia subglobosa. j, k, l Dothiorella parva. m, n, o Diaporthe eres. p, q, r Diaporthe rudis. s, t, u Nothophoma brennandiae. v, w, x Anthostoma decipiens. Scale bar = (c, f, i, l) 20 µm. o, r, u, x 10 µm

Pathogenicity

All the tested isolates caused symptoms similar to those found in the field (Fig. 7). The plants inoculated with the fungal isolates showed external longitudinal streaks or dark brown lesions on the bark. Plants inoculated with representative isolates of A. decipiens were characterized by elongate canker development at the inoculation point. Internal necrosis and vascular discolouration were observed. Isolate CVG1374 of A. decipiens showed the highest lesion length (76.4 ± 22.6 mm), followed by A. decipiens isolate CVG1380 (lesion length = 76.4 ± 22.6 mm) and Dia. eres isolate CVG1334 (lesion length = 30.3 ± 7.3 mm). The values of mean lesion length for nine strains ranged from 13.5 ± 4.0 to 4.8 ± 1.8 mm for B. dothidea strain CVG2219 and Dia. rudis strain CVG1333, respectively, however their aggressiveness in hazelnut plants was not significantly different among them. The four remaining strains caused lesions which showed no significant differences among them, with lesion length ranging from 4.1 ± 0.7 to 2.5 ± 0.6 mm for Dip. subglobosa strain CVG1367 and Do. parva strain CVG1415 (Fig. 8). All the inoculated species were re-isolated and the morphological characteristics (color, shape, mycelium texture and conidia) were assessed to confirm their identity. Koch’s postulates were fulfilled. Weak necrosis observed on control plants (mean lesion length = 2.5 ± 0.6 mm) were restricted to the inoculation point and considered as a reaction to wounding. No fungal colonies were retrieved after re-isolation from control plants.

Fig. 7
figure 7

Internal lesion of hazelnut branches of cv. ‘Tonda Gentile’ at 3 months after inoculation with mycelial plugs of the species: a Anthostoma decipiens. b Botryosphaeria dothidea. c Diaporthe eres. d Dothiorella parva. e Diplodia seriata. f Nothophoma brennandiae. g Diaporthe rudis. h Diplodia subglobosa. i control

Fig. 8
figure 8

Disease severity (lesion length, mm) on hazelnut plants of ‘Tonda Gentile’ at 3 months after inoculation of representative fungal isolates of A. decipiens, B. dothidea, Dia. eres, Dia. rudis, Dip. seriata, Dip. subglobosa, Do. parva and N. brennandiae. Columns represent the mean data of seven replicate plants per strain. Columns with common letters do not differ significantly according to Tukey’s HSD test (P = 0.05) for lesion length. Mean comparison test was applied to log-transformed lesion length data. Horizontal lines on the columns are the standard error of the mean

Discussion

The present study represents the first investigation on the fungal species diversity in association with FTD of hazelnut in Piedmont. Currently, the spread of different FTD symptoms across this area is representing an increasing concern to hazelnut producers. The surveys, conducted in seven orchards, focused on sample collection from plants showing different symptoms, thus, to clarify the identification of the fungal species associated with the main two wood diseases reported. Isolates recovered from plants showing twig and branch dieback were grouped in Botryosphaeriaceae, Diaporthe spp. and Phoma-like spp. according to their colony morphology (Phillips et al. 2013; Udayanga et al. 2014; Chen et al. 2015). Isolates from plants showing symptoms of Cytospora canker, known also as ‘Mal dello stacco’, were grouped in Botryosphaeriaceae and Cytospora-like spp., based on their culture characteristic (Phillips et al. 2013; Lawrence et al. 2018). Sequencing of molecular loci and phylogenetic analyses allowed the identification of eight different species associated with FTD on hazelnut in Piedmont: Anthostoma decipiens, Botryosphaeria dothidea, Diaporthe eres, Dia. rudis, Diplodia seriata, Dip. subglobosa, Dothiorella parva and Nothophoma brennandiae. Particularly, Diaporthe spp., Diplodia spp., Do. parva and N. brennandiae were isolated in association with twig and branch dieback, whilst A. decipiens and B. dothidea were found in association with Cytospora canker (‘Mal dello stacco’). This study is also the first report of B. dothidea and Dia. eres as causal agents of FTD on hazelnut in Italy. Moreover, it is the first report of Dia. rudis, Dip. subglobosa and N. brennandiae causing twig and branch dieback on hazelnut trees worldwide. Species belonging to Botryosphaeriaceae family are well known as polyphagous pathogens for their wide distribution and virulence on multiple plant hosts (Batista et al. 2021). Botryosphaeria dothidea was previously reported on hazelnut trees in Iran and Turkey (Mohammadi and Jabbari Firoozjah 2019; Polat et al. 2022) and Dip. seriata was reported as wood pathogen on hazelnut in Sardinia (Linaldeddu et al. 2016). Hereby, B. dothidea was isolated in association with A. decipiens in the surveyed orchard n°6 (Monteu Roero-CN), whilst Dip. seriata was found in two different sites, orchard n°4 (Feisoglio-CN) and orchard n°7 (Monteu Roero-CN). Diplodia subglobosa was described as a new species in association with Fraxinus spp. and Lonicera nigra and it was recently reported on Fraxinus excelsior in Italy and Slovenia (Alves et al. 2014; Linaldeddu et al. 2020, 2022). In this study, this species was found in co-occurrence with A. decipiens in orchard n°2 (Diano d’Alba-CN). Dothiorella parva was isolated from Corylus avellana in Spain and reported as Dothiorella sp. (Phillips et al. 2008). This species was described as sp. nov. considering isolates from hazelnut in Iran and it was later reported as pathogen on Ostrya carpinifolia, a forest tree within the Betulaceae, in Slovenia and Italy (Abdollahzadeh et al. 2014; Pavlic-Zupanc et al. 2015). Recently, its name was proposed as a synonym of Do. sarmentorum (Zhang et al. 2021). In the present survey, it was found in co-occurrence with N. brennandiae in the orchard n°5 (Feisoglio-CN). All the Botryosphaeriaceae species found in the present study, B. dothidea, Dip. seriata, Dip. subglobosa and Do. parva, were pathogenic when inoculated on healthy hazelnut plants with a different virulence depending on both the tested species and isolates. Particularly, Dip. subglobosa showed minor symptoms with respect to the other tested isolates within Botryosphaeriaceae. The different virulence of the two isolates of Do. parva suggests the presence of intraspecific variability. Diaporthe spp. are included among the most globally relevant causal agents of FTD on different fruit crops (Lawrence et al. 2015; Guarnaccia et al. 2022b). Diaporthe eres was reported as wood pathogen on Corylus avellana in Oregon and Chile (Guerrero and Pérez 2013a; Wiman et al. 2019) and on Corylus heterophylla in China (Gao et al. 2021; Bai et al. 2022). Diaporthe rudis was reported on hazelnut only as causal agent of kernel defects, as well as Dia. eres (Pscheidt et al. 2019; Arciuolo et al. 2020, 2022). In this study, both species were found in co-occurrence with Botryosphaeriaceae in two out of seven orchards, particularly in orchard n°2 (Diano d’Alba-CN) and n°7 (Monteu Roero-CN). The pathogenicity test on these species showed a high virulence of one isolate of Dia. eres, whilst the other isolates caused minor lesions as both the isolates of Dia. rudis. The difference found suggests an intraspecific variability for pathogenicity in Dia. eres that was already reported for this species (Lawrence et al. 2015). Nothophoma brennandiae was originally isolated from Ulmus × hollandica in Italy and reported as N. quercina and then described as a new species isolated from garden soil in The Netherlands (Hou et al. 2020). It was found in co-occurrence with A. decipiens and Do. parva in orchard n°1 (Albaretto della Torre-CN) and orchard n°5 (Feisoglio-CN), respectively. This species caused minor symptoms when inoculated on hazelnut plants. The identification of A. decipiens represents the clarification of the causal agent of Cytospora canker (‘Mal dello stacco’) in the surveyed area and it is in line with the previous report in Sardinia (Linaldeddu et al. 2016). The doubts about the fungal species causing the Cytospora canker on hazelnut was long debated in Piedmont as well as in other Italian regions (Scortichini 2006; Botta and Valentini 2018). The cultural characteristics of A. decipiens, which have a high similarity with Cytospora spp. in terms of size and shape of conidia, have so far led the attribution of Cytospora canker to Cytospora corylicola or Cytospora spp. (Linaldeddu et al. 2016). Anthostoma decipiens was reported as pathogen on Carpinus betulus, another plant of the Betulaceae, in Italy (Saracchi et al. 2008; Rocchi et al. 2010) and Iran (Mirabolfathy et al. 2018). Moreover, it was found in association with grapevine trunk disease in Spain (Luque et al. 2012) and reported as pathogen of different plant hosts such as Alnus glutinosa, Betula pendula, Castanea sativa, Corylus avellana, Vagus sylvatica and Ostrya carpinifolia (Saracchi et al. 2015). In the present study, both isolates used in the pathogenicity test showed high virulence and caused the most severe lesions with respect to the other tested species. The ability of this pathogen to cause severe internal lesion with consequent cracks in the bark represents a serious threat to hazelnut, especially in late summer when conidia evasion occurs promoting the inoculum spread. Conidia of A. decipiens are thought to be dispersed through wind and rains as those of Botryosphaeriaceae and Diaporthe spp., representing a primary source of inoculum (Phillips et al. 2013; Lawrence et al. 2015). Pruning wounds, along with natural openings, are the main pathway for the colonization (Moral et al. 2019). The use of non-contaminated tools, the protection of wounds and the removal of pruning debris must be adopted. Proper agronomic practices contribute to avoid the inoculum spread, especially considering that several pathogens hereby found were reported in the same area on other fruit crops such as apple, grapevine and blueberry (Martino et al. 2022, 2023; Dardani et al. 2023). Additional surveys in more orchards as well as pathogenicity tests are planned, to investigate the incidence of the disease in North-western Italy as well as the infection rate caused by spore suspensions of each fungal pathogen. Moreover, the influence of inoculum density on the disease development will be assessed along with the seasonal dynamic of spore production, then to provide a correct pruning timing. Further studies on co-infection of different fungal species will be addressed to provide useful information on the complexity of FTD on hazelnut. Several fungal species found in the present study caused minor symptoms on healthy hazelnut plants, suggesting that they could act as secondary pathogens. Insights on the biology and ecology of these species are needed to establish effective disease management strategies. Moreover, understanding the interaction with abiotic stress and wood microbial community could contribute to the knowledge on FTD, promoting the high value of the cultivar Tonda Gentile produced in Piedmont to which the Protected Geographical Indication (PGI) was recognized by the European Union with the name ‘Nocciola Piemonte’ (‘Nocciola Piemonte delle Langhe’, from 2019).