INTRODUCTION
Barley (Hordeum vulgare) is the fourth most cultivated cereal worldwide, after corn (Zea mays ssp. mays), common wheat (Triticum aestivum), and rice (Oryza sativa) (FAO, 2020; Producción Agrícola Mundial, 2020). Climate change has led to the implementation of new cultivation practices, which has resulted in changes in the profiles and ranges of pathogens that affect these crops (Pautasso et al., 2012). Barley plants are affected by the presence of several types of foliar pathogenic fungi, such as Ramularia collo-cygni, Pyrenophora teres, Cochliobolus sativus, Puccinia hordei and Rhynchosporium secalis, which are responsible for causing diseases (Carretero et al., 2015; Romero et al., 2020). In fact, net blotch (NB) caused by P. teres (anamorph Drechslera teres) can induce yield losses greater than 50% in barley (Statkeviciute et al., 2010; McLean and Hollaway, 2018) due to a reduction in the number of kernels per spike, grain size, and total photosynthetic leaf area (Dumalasová et al., 2012). In template regions, high rainfall (between 153 and 154 mm) aggravates the damage caused by NB (Steffenson and Webster, 1992; Wu et al., 2020). However, epidemics have also occurred in areas with lower rainfall (Steffenson and Webster, 1992; Ellwood et al., 2018), where barley seeds with a high fungal load constitute an inoculum for P. teres survival from one season to another, becoming a very efficient way of transmitting the pathogen to the seedling, coleoptile and roots (Couretot et al., 2017). NB has two forms differentiated by the symptoms: net form (NFNB) and spot form (SFNB). NFNB, caused by P. teres f. teres, induces horizontal and vertical crisscrossed dark brown venation on barley leaves that can become chlorotic (Dumalasová et al., 2012), whereas SFNB, caused by P. teres f. maculata, produces circular or elliptical dark brown spots surrounded by chlorotic leaf tissue (Carlsen et al., 2017).
The pathogen agent of NFNB can be easily identified in barley plants based on symptomatology. However, SFNB presents non-cross-linked symptoms without a net-spot (Marshall et al., 2015), which closely resemble those caused by C. sativus (Rehman et al., 2020). When spikes are attacked by any of these fungi, the infected seeds become of a dark color (Carmona and Sautua, 2015). Therefore, macroscopic characteristics observed in barley seeds are not useful for the identification of the fungus responsible for the infection (Carmona and Sautua, 2015) because there are several fungi that can cause discoloration and seed darkening, such as Alternaria spp., Drechslera verticillata, D. campanulata, and Bipolaris zeicola (Medd et al., 2003; Manamgoda et al., 2014; Romero-Cortes et al., 2019). Consequently, microscopic observations of the conidia are required to distinguish the pathogen agent of SFNB (Poudel et al., 2019), and also to differentiate between C. sativus and P. teres net-spot blotch-causing fungi since they present no macroscopic morphological differences (Statkeviciute et al., 2010).
The genus Cochliobolus is another group of fungi that affects barley, including C. sativus (Sivanesan, 1990), (anamorph Bipolaris sorokiniana, syn. Helminthosporium sativum), which presents a wide range of hosts in the Poaceae family (Musaed et al., 2013). C. sativus is widely distributed (Bashyal et al., 2011) and constitutes a continuous genetic group of isolates that vary in virulence and aggressiveness against several cereals and grasses (Duveiller and Altamirano, 2000). Affected crops include barley, spring wheat, rye, weeds, and other grass species (Jones et al., 1983; Wu et al., 2020). Of these cereal species, wheat and barley are the most highly produced crops (Murray et al., 2009; Mata-Santoyo et al., 2018). In barley, C. sativus causes different diseases such as black point, common root rot, and spot blotch, which can affect both plants and/or seeds (Kumar et al., 2002; Manamgoda et al., 2014; Romero et al., 2020; Al-Sadi, 2021). Severity of attack and damage caused to the crop is greater in soils with poor fertility, and high temperature and humidity conditions (Duveiller and Altamirano, 2000; Savary et al., 2011; Kumar et al., 2020).
Two different pathological stages can be distinguished: the interference of the photosynthesis process, followed by the interference in the absorption of nutrients and water by the roots (Forcelini, 1991; Carmona and Sautua, 2015; Carmona and Sautua, 2017). Affected plants result in fewer tillers and grains per spike, as well as reduced seed germination, plant growth, and yield. In particular, yield reduction by C. sativus can reach up to 30% (Kumar et al., 2002; Ghazvini and Tekauz, 2007; Al-Sadi and Deadman, 2010; Sultana et al., 2018).
The damage caused by the members of the Pleosporaceae family in barley results in severe economic losses worldwide. In this sense, some farmers are unable to use technology and rely on family labor. In the highlands of Mexico, a recent study focused on the morphological identification of fungi belonging to the leaf spot complex in barley (Romero et al., 2020). However, there is a need for the molecular identification of the pathogenic fungi associated with this crop in the area. Therefore, the aim of this research was to determine, morphologically and molecularly, strains of pathogenic fungi isolated from barley infected seeds in the highlands of Mexico, in the states of Hidalgo, Puebla and Guanajuato.
MATERIALS AND METHODS
Fungal strain isolation
Fungal strains were isolated from five varieties of barley seeds with dark color symptoms collected during the years 2008 to 2014 in the highlands of Mexico. Fungi were isolated in Petri dishes, containing 20 seeds previously disinfected by immersion in a 2% solution of aqueous sodium hypochlorite for 120 s, followed by three washes with sterile water. Afterwards, the seeds were placed in moist chambers and incubated at 15°C under visible light to induce sporulation. Fungi were isolated from barley seeds using potato dextrose agar (PDA) (Sigma, St. Louis, MO, USA) at 25°C for 7 days. Subsequently, growing colonies were selected based on their macroscopic characteristics, cut using a scalpel, and transferred to a new PDA plate for incubation under the same aforementioned conditions. Conidia and hyphae were maintained at -84°C until analysis.
Morphological characterization
The morphological description of C. sativus (Kumar et al., 2002; Manamgoda et al., 2014) and P. teres (Louw et al., 1995; Crous et al., 1995) was conducted as described. Growth rate and morphological descriptions (conidia and conidiophores) were performed using a PDA culture medium. A 0.5 cm-diameter mycelium disk containing isolated conidia was placed in the center of each Petri dish. Fungi were grown at 25°C for one week or until they filled the entire Petri dish. Micrographs were taken with a scanning electron microscope (SEM) (JEOL, Model IT300, Boston, MA, USA).
Molecular characterization
The DNA of the strains identified as C. sativus and P. teres was obtained from mycelium grown in 250 ml Erlenmeyer flasks containing 50 ml of potato dextrose broth (PDB) (20% potato, 15% dextrose (w/v)), using a modified method described by Aamir et al. (2015). Strains were grown with shaking (220 rpm) at 25 °C for a week. Biomass was recovered and washed twice with TE buffer (10 mM Tris-HCl, 1 mM EDTA). DNA extraction of each isolate was performed according to Cuervo-Parra et al. (2011). The ITS regions (ITS1 and ITS2) of the ribosomal DNA (rDNA) gene were amplified by PCR using the primers ITS1 and ITS4 (Ramírez-Cariño et al., 2020).
The PCR system was the same as described by Cuervo-Parra et al. (2014). PCR reactions were placed in an automatic thermal cycler (Techne PRIME3, Series No. 31309, Cole Parmer, Staffordshire, UK) under the following conditions: 5 min initial denaturation at 95°C, followed by 30 cycles of 1 min denaturation at 95°C; 1 min primer annealing at 57°C; 1 min extension at 72°; and a final extension period of 12 min at 72°C. The resulting products were purified with the GeneClean® II kit (MP Biomedicals, Carlsbad, CA, USA) according to the manufacturer's protocol. Amplicons were observed on a 2% agarose gel stained with ethidium bromide. Sequence pair alignment was carried out using BLAST software (Zhang et al., 2000), and the DNA sequences of C. sativus and P. teres were aligned with other related sequences from GenBank (Andrie et al., 2008; Aggarwal et al., 2010a; Aggarwal et al., 2010b; Kleczewski, 2010; Liu and Friesen, 2010; Lartey et al., 2013; Moya et al., 2013; Manamgoda et al., 2014; Shim, 2014; Phillips et al., 2016) through ClustalW (Larkin et al., 2007) to identify variable regions or base sequences. Distance matrices were calculated between all pairs of sequences considered in the multiple alignment, and two phylogenetic trees were generated based on Neighbor-Joining (NJ) and Maximum Likelihood (ML) statistical methods using MEGA6 software (Tamura et al., 2013). Same as for the statistical analysis, the phylogenetic analysis was performed based on: bootstrap test of phylogeny using 1000 random bootstrap replications, nucleotide substitutions, Kimura's two-parameter model, uniform rates among sites and homogeneous pattern among lineages. Additionally, the following statistical parameters were used for the ML statistical method: ML Heuristic method, Nearest-Neighbor-Interchange (NNI) and branch swap filter, very strong. All parameters were calculated using MEGA6 software (Tamura et al., 2013).
RESULTS AND DISCUSSION
Isolation of strains
Sixty-nine fungal strains were isolated and morphologically characterized from five barley varieties (Esmeralda, Forrajera, Alina, Adabella, and Gaviota) (data not shown). Based on the morphological characterization and the importance of some of these fungi as pathogens, 16 strains were identified as C. sativus and P. teres (Table 1). It is observed that these fungi have been present in the field over several years. C. sativus strains have been present in Puebla and in all the municipalities of Hidalgo, showing a constant presence in barley plantations.
Morphological characterization
Colonies grown on PDA medium (5-7 d, 25 °C) showed a diameter of 36 to 42 mm, with an initial whitish growth that later acquired a gray-brown
to black color, with irregularly shaped mycelial structures of white color, velvety, with an irregular or wavy margin. Historically, the genus Cochliobolus has been characterized by having brown conidiophores and conidia that can be fusoid, straight, or curved, and germinating from a germinal tube at each end (Manamgoda et al., 2014). Abundant sporulation was observed after one week. The reverse was dark or black with concentric rings. The hyphae were gray to brown, smooth or verrucous. Asexual morphological characteristics on PDA medium showed abundant conidiophores, straight to flexuous, with septa, arising singly or grouped, simple or branched (Fig. 1A-C), with size of 130.59-154.53 * 6.18-7.76 |jm (av. = 142.56 |jm, SD = 11.97, n = 32, long; av. = 6.94 |jm, SD = 0.77, n = 32, wide). Conidia were dark olivaceous brown to golden brown or slightly pale at ends, curved or straight, oval, cylindrical, fusiform or with broadly ellipsoidal shape, with rounded ends and a basal scar (Fig. 1D-G). A slight curvature was observed in some conidia. The size of the conidia was 77.78-146.4 * 15.66-19.64 |jm, (av. = 112.09 |jm, SD = 34.31, n = 32, long; av. = 17.65 |jm, SD = 1.99, n = 32, wide) with 3 to 10 pseudosepta. The conidial wall was smooth and distinctly thick in the pseudosepta. The morphological values obtained in the present study are within the range observed by other authors for other strains of C. sativus (Müller et al., 2005; Manamgoda et al., 2014). Therefore, based on the morphological characterization, the strains were identified as C. sativus (Kumar et al., 2002; Manamgoda et al., 2014; Romero et al., 2020).
Eight strains (JCP2, JCP7, JCP10, JCP16, JCP19, JCP50, JCP51 and JCP57) grown on PDA showed an initial hyaline growth that over time acquired an olivaceous brown to white-cream coloration, with cottony colonies with a 5.2 mm day-1 growth rate. The reverse showed a crenulated, dentate, serrate or an irregular growth with a coloration between yellow-creamy to light brown. The conidiophores observed on PDA medium showed different shapes, septate, being in the range of 3.64-79.48 * 10.56-16.82 |jm (av. = 41.56 |jm, SD = 37.92, n = 32; av. = 13.69 |jm, SD = 3.13, n = 32).
Conidia of isolates were 36.37-102.55 * 13.95 25.61 |jm (av. = 69.46 |jm, SD = 33.09, n = 32; av. = 19.78 |jm, SD = 5.83, n = 32), olivaceous brown, cylindrical, tapering to subtruncate bases, with transverse septa (dictyoseptate), constricted at the level of the septum (Fig. 2). The morphological results obtained for the strains of P. teres were in agreement with those previously reported in the literature (Crous et al., 1995; Louw et al., 1995; Romero et al., 2020). Although the teleomorph obtained in this study had minimum values for the conidiophores and conidia, the results agree with values reported by Crous et al. (1995) and Romero et al. (2020). Conidia are more important for the propagation of NB disease because they are generated throughout the growing season, meanwhile sexual reproduction can lead to the formation of new pathotypes (Statkeviciute et al., 2010).
Molecular characterization
The analysis of the ITS regions of rDNA is a simple and reproducible molecular tool
to identify and characterize fungal species (Mohammadi and Amini, 2015). A total of 16 new sequences of phytopathogenic fungi were generated in this study with additional sequences downloaded from GenBank. Phytopathogenic fungi of this study were identified using the ITS1 and ITS4 primers, by analyzing the DNA region containing the ITS1/5.8s/ITS2 sequence. The amplified regions of the isolated fungal strains had a size between 507 and 617 bp and were deposited into the NCBI GenBank databases with the following accession numbers KX951412, KX066030, KX066031, KX066032, KX066033, KX066034, KX066035, KX066036 for C. sativus, and KX115407, KX115408, KX115409, KX115410, KX115411, KX115412, KX115413, KX115414 for P. teres.
From ITS regions, the data obtained for all the sequences of Cochliobolus and Pyrenophora were consistent with the results reported in previous studies for Pyrenophora graminea, P. teres, and P. teres f. maculate (Carlsen et al., 2017; Poudel et al., 2019) and for several species of Bipolaris (Manamgoda et al., 2014). The genetic differences observed between isolates of C. sativus and P. teres, collected from different cereal-growing regions in Mexico, could explain the difference in virulence (Christensen, 1926; Sultana et al., 2018). In this sense, a study that evaluated a barley core collection for SFNB reaction revealed distinct genotype-specific pathogen virulence and host susceptibility (Neupane et al., 2015). Wu et al. (2020) performed a PCR assay with the ToxA-specific marker, a significant virulence factor related to wheat spot blotch disease, for 196 Mexican strains of C. sativus isolated from wheat and barley, and reported the presence of ToxA in 20 of the isolates, but none of the eight B. sorokiniana isolated from barley was ToxA-positive. Compared with reports for other geographical regions (McDonald et al., 2018; Friesen et al., 2018; Navathe et al., 2020), these results suggest that the Mexican populations of C. sativus are not carriers of ToxA.
Phylogenetic analysis of the Mexican barley sequences and other 64 related sequences from GenBank generated a tree by the NJ method (Fig. 3), and another tree by the ML method (Fig. 4). The trees were rooted with the GenBank sequences AB474001 and KM659218 of Bacillus subtilis as outgroup to infer phylogenetic relationships (El-Sayed et al., 2008; You and Zhang, 2014). In the present study, two distinctive fungal taxa isolated from barley seeds were > 96% similar to Cochliobolus and Pyrenophora. Of the ITS regions, the ITS2 region was the most conserved for all strains of Cochliobolus and Pyrenophora. The optimal NJ tree with the sum of branch length = 1.62618503 and the ML tree with the Log Likelihood = -1809.97 are shown. The NJ and ML analyses showed that the DNA sequences of P. teres and C. sativus, lined up at separate branches, formed two major clusters (I and II). Cluster I contained the JCP50, JCP19, JCP16, JCP7, JCP57, JCP2, JCP51 and JCP10 strains, and related Pyrenophora/Drechslera isolates from GenBank. Cluster II contained the JCP12, JCP14, JCP65, JCP30, JCP3, JCP43, JCP6 and JCP44 strains, and related strains of Cochliobolus also from GenBank. Each node of the phylogenetic trees was well supported by high bootstrap values (Figs. 3 and 4).
At the top of the phylogenetic trees, it can be seen that the Mexican strains isolated from barley as P. teres were phylogenetically related to the NCBI GenBank sequences of P. teres from South Korea (Shim, 2014), Argentina (Moya et al., 2013), and United States of America (Liu and Friesen, 2010; Moya et al., 2013), and of P. teres f. teres (Liu and Friesen, 2010) and P. teres f. maculata, both from the United States of America (Lartey et al., 2013). There were some variations between the sequences of the strains of P. teres isolated from Mexican barley in terms of deletions, insertions, substitutions and base transversions during multiple sequence alignment using ClustalW software (Larkin et al., 2007). However, the differences do not represent enough variation within the NCBI GenBank Pyrenophora sequences as to indicate that the Mexican strains of P. teres are a different species.
The bootstrap value obtained for the strains of P. teres by the NJ method was 96 during the phylogenetic analysis, while values generated by the ML method were higher. Since values were high, these strains were easily distinguishable from other sequences of Pyrenophora species from GenBank grouped in other branches of the same Cluster I of the phylogenetic trees (e.g., Pyrenophora tetrarrhenae, Drechslera poae, and Pyrenophora tricini-repentis). Nevertheless, more gene region information is needed to establish genetic divergence from other species or continuity within the P. teres complex. In 1967, Pyrenophora japonica was considered as a mutant strain of P. teres, consisting of numerous clones that are tailored to environments (McDonald, 1967). Later studies have confirmed that P. japonica was not a separate species (Crous et al., 1995; Campbell et al., 1999). Other studies have described that P. graminea morphologically resembles P. teres but causes a leaf stripe disease in the barley plant, while it cannot infect leaves directly through conidia (Smedegard-Petersen, 1977; Liu et al., 2011). Furthermore, phylogenetic analyzes show that groups of P. graminea
sequences are the closest neighbors to groups of P. teres sequences, being more closely related to P. teres f. maculata than to P. teres f. teres (Bakonyi and Justesen, 2007; Rau et al., 2007).
In Australia and South Africa, the anamorph state of Pyrenophora hordei shows morphological characteristics similar to those of P. teres, but its teleomorph state is different in both morphology and disease symptoms in barley (Scott, 1994). Another study suggests that the sexual recombination between the two forms of NB of barley occurs naturally in the field, leading to the appearance of a new SFNB differential set (McLean et al., 2014). Furthermore, to improve the resolution of the phylogenetic tree, a second gene such as gpd could be incorporated to make the phylogenetic analysis more robust (Zhang and Berbee, 2001; Lepoint et al., 2010; Vasighzadeh et al., 2019). In addition, it would be useful to detect intraspecific variability in P. teres by using specific primers to detect and discriminate between net and spot forms (Williams et al., 2001; Svobodova et al., 2005; Leisova et al., 2006; Poudel et al., 2019).
The fungal sequences of the genus Cochliobolus isolated from barley grains were aligned in the upper part of Cluster II in the two phylogenetic trees, together with other C. sativus sequences from GenBank from several geographical regions (Figs. 3 and 4). For the NJ method, the JCP12, JCP14, and JCP44 strains were located near C. sativus sequences from India (Aggarwal et al., 2010a) and New Zealand (Phillips et al., 2016). For the ML method, these sequences were aligned in other branches of the same clade, along with the other C. sativus sequences isolated from Mexican barley. The JCP30, JCP3, and JCP43 strains lined up next to other strains of C. sativus isolated in the United States of America (Kleczewski, 2010), supported by a high bootstrap value for the NJ and ML methods (Figs. 3 and 4). Additionally, the JCP65 strain was aligned near two sequences of C. sativus isolated in Italy and the United States of America (Manamgoda et al., 2014), respectively, with a bootstrap support value of 57 for the NJ method. For the ML method, this sequence was aligned near the JCP30, JCP3, JCP43 strains isolated from barley samples in the present study and other strain of C. sativus isolated in the United States of America (Kleczewski, 2010). Finally, the JCP6 strain was aligned between two others C. sativus sequences isolated in India (Aggarwal et al., 2010b) and Japan (Manamgoda et al., 2014) by the NJ method. Furthermore, the phylogenetic analysis conducted by the ML method allowed differentiating the C. sativus sequences isolated from barley in Mexico from the GenBank database from other geographical regions (Fig. 4).
Although the bootstrap values for the C. sativus sequences in the two phylogenetic trees were not very high (67 and 68 for the NJ and ML methods, respectively), it was possible to clearly differentiate the isolates from the other sequences belonging to other species within the genus Cochliobolus. In addition, the genetic variation of C. sativus strains isolated from barley seeds from different cereal-growing regions in Mexico is consistent with what reported in the literature (Duveiller and Altamirano, 2000; Kumar et al., 2002; Bashyal et al., 2011; Manamgoda et al., 2012; Manamgoda et al., 2014).
The obtained molecular results are in agreement with previous studies (Manamgoda et al., 2012; Manamgoda et al., 2014). The large number of fungal strains considered in various studies conducted worldwide shows that C. sativus derives from a continuous isolation with variations in virulence and aggressiveness with specific and non-specific interactions (Maraite et al., 1998; Ghazvini and Tekauz, 2007). On the other hand, the distinction of the genus Cochliobolus from its sister group observed in the phylogenetic analysis generated for Cochliobolus and Pyrenophora sequences by the NJ and ML methods (based on the ITS alignment), agrees with molecular data available in the literature (Berbee et al., 1999), then redefined by using phylogenetic analyses of ITS, GPDH, TEF, and LSU sequences (Manamgoda et al., 2012).
CONCLUSIONS
The fungi P. teres and C. sativus that cause spot blotch disease of barley and the pathogenesis involving these fungi were subjected to detailed studies at the molecular level. Morphological and ITS sequence-based identification showed that the fungal strains isolated from barley grown in the highlands of Mexico were between 96 and 99% similar to C. sativus and P. teres. Strains of C. sativus can be identified with high precision using BLAST software for ITS sequences against the sequences deposited in the NCBI database. However, this is not the case for P. teres as the genus can be determined but not the intragenic species. The data obtained in this study will allow for the subsequent development of biological control strategies for these pathogens with other fungal agents, such as Trichoderma spp.