Sampling and isolation of Fusarium
Fusarium isolates were obtained from sunflower seeds of eight cultivars (Aguará 04, Aguará 06, BRS 321, BRS 322, BRS G26, GF 101, Hélio 250, and Olisun 3) exploited in Alagoinha, Paraíba State, Brazil (06º 57' 00'' S, 35º 32' 42'' W), during the 2017/18 crop season. Seeds were provided by ‘Empresa Brasileira de Pesquisa Agropecuária – EMBRAPA’, Campina Grande, Paraíba State, Brazil. After discarding the damaged ones, a set of 200 seeds were surface disinfected by submersing in a 1.0% aqueous solution of NaClO for 3 min, rinsed three times with sterile water, and left to dry on sterilized filter paper prior subjection to the blotter test (BRASIL, 2009). Seeds were individually transferred to Petri dishes (15.0 cm diameter) double lined with sterilized filter paper moistened with 10 mL of sterile water and incubated for seven days into Biological Oxigen Demand (B.O.D.) chamber set to 25 ± 2°C, 12:12 L:D daylength. By using a sterile histological pin, typical structures of Fusarium that grew over the seeds were transferred to Petri dishes (9.0 cm diameter) previously filled with potato dextrose agar (PDA) media, and incubated for seven days into B.O.D. under similar conditions. Based on morphological features of Fusarium spp. (Leslie & Summerell, 2006), monosporic cultures were established, and the isolates were deposited on a collection at the Phytopathology Laboratory of the ‘Universidade Federal da Paraíba (UFPB)’, Areia, Paraíba State, Brazil.
Morphological characterization
Representative isolates of Fusarium spp. were morphologically characterized following Leslie and Summerell (2006). The pigmentation of 10–14 days old colonies grown at 25 ± 2 ºC on PDA was determined using the mycological color charts of Rayner (1970). In contrast, 10–14 days old colonies grown at 25 ± 2 ºC on synthetic low-nutrient agar (SNA) media added with carnation leaves were used to observe length, width, shape, and septation of micro- and macroconidia (n = 15 of each); disposition of the microconidia (i.e., false head or chains); nature of conidiogenous cells (i.e., mono- and/or polyphialides); and the presence or absence of chlamydospores. Both micro- and macroconidia were measured using the program Motic Images Plus 2.0 in a computer coupled with an Olympus BX53 microscope. Photographs were taken using an Olympus DP73 camera and the program CellSens Dimension at 200, 400 and 1000x magnification. A composite photo plate was assembled from individual photographs to enable the comparison of relevant micromorphological features.
To evaluate the mycelial growth of the representative isolates, a disc (5.0 mm diameter) from each colony was transferred to four Petri dishes (9.0 cm diameter) filled with PDA incubated in the dark at 25 ± 2 ºC. For seven consecutive days, the diameter of the colonies was measured with a digital caliper rule, and the Mycelial Growth Index (MGI) was estimated based in Oliveira (1991) as follows: MGI = (D - Dp)/N, where D is the colony diameter in the current day; Dp is the colony diameter from the previous day; and N is the number of days after inoculation, respectively. The final colony size (FCS) was measured on the 7th day post inoculation.
To estimate the conidia production, a spore suspension was prepared by adding 20 mL of sterile water to 7 days old colonies grown on PDA at 25 ± 2 ºC, which were scraped with a soft paint brush to dislodge the spores. This suspension was filtered through double folded sterile gauze and the spores were counted using a hemocytometer (Alfenas & Mafia, 2016).
All experiments followed a randomized design.
DNA extraction, PCR and DNA sequencing
Genomic DNA extraction from 5–10 days old colonies of each Fusarium isolate were performed using the cetyltrimethylammonium ammonium bromide (CTAB) method (Doyle & Doyle, 1987), with some modifications as described previously (Materatski et al., 2019). The quality and concentration of DNA was determined by using a Quawell Q9000 micro spectrophotometer (Quawell Technology, Beijing, China).
The elongation factor gene (EF-1α) was amplified with Fa + 7 and Ra + 6 primers (Karlsson et al., 2016) to estimate the initial genetic diversity and select representative isolates for amplification of the β-tubulin 2 gene (TUB2). Different haplotypes were identified with DnaSP 4.0 (Rozas et al., 2003), and a haplotype network was generated through PopART v. 1.7 (Leigh & Bryant, 2015) by using the algorithm TCS (Clement et al., 2000; Clement et al., 2002). The representative isolates were chosen considering the sunflower cultivars, and their β-tubulin 2 (TUB2) partial gene region was amplified with T1 and T22 primers (O'Donnel & Cizelnik, 1997).
Polymerase chain reactions (PCR) were performed in a final volume of 25 µL containing a mixture of 7–20 ng of genomic DNA, Tris-HCl 10 mM (pH 8.6), KCl 50 mM, MgCl2 1.5 mM, and dNTPs 0.2 mM, 0.2 µM of each primer, and 2.5 U of DreamTaq DNA polymerase. PCR cycles were performed in a thermocycler (BioRad) as follow: initial denaturation at 95 ºC for 2 min, followed by 40 cycles of denaturation at 95 ºC for 30 s, annealing at 67 ºC (EF-1α) or 58 ºC (TUB2) for 50 s, extension at 72 ºC for 60 s and a final extension at 72 ºC for 10 min. The PCR products were separated by electrophoresis in a 1.5% agarose gel with Tris-acetate-EDTA (TAE) 1.0X buffer and were revealed under UV light. These PCR products were purified with DNA Clean & Concentrator (Zymo Research) according to the manufacturer guidelines. DNA sequencing was performed by Macrogen (Madrid, Spain).
Sequence alignment and phylogenetic analysis
Consensus were obtained using the Staden Package (Staden et al., 1998). All obtained consensus sequences were used as queries to the Fusarium-ID database (Torres-Cruz et al., 2022) and NCBI nucleotide database using the blast algorithm (Johnson et al., 2008) to confirm the taxonomic assignments of the isolates. Sequences representing ex-types and related sequences were retrieved from GenBank. Sequences from the present study were deposited in GenBank (Supplementary table S1).
Multiple sequence alignments for each individual locus were estimated online using the G-INS-i strategy in MAFFT version 7 (Katok & Standley, 2013; Katoh et al., 2019), with default parameters for gap opening and extension and 200PAM / κ = 2 nucleotide scoring matrix. Wherever needed, they were manually adjusted in MEGA6 (Tamura et al., 2013).
Phylogenetic analyses were performed using the Maximum Likelihood (ML) and Bayesian Inference (BI) methods for both individual and concatenated genes. ML and BI analyses were performed using RAXML-HCP2 v.7.0.4 (Stamatakis, 2014) and MrBayes v 3.2.1 (Ronquist et al., 2012), respectively, implemented in the CIPRES cluster (https://www.phylo.org/portal2/home.action). ML analyses were carried out with 1000 pseudo replicates (-m GTRGAMMA -p 12345 -k -f a -N 1000 -x 12345) under the GTR-GAMMA model. With respect to BI, evolution models were estimated in MrModeltest 2.3 (Nylander, 2004) using the Akaike information criterion (AIC) for each gene. The combined data set was partitioned to reflect the most appropriate nucleotide substitution model for each of the single gene data sets for the BI of the combined data set. Four Markov Chain Monte Carlo (MCMC) chains were conducted for 107 generations, with samplings every 1000 generations. All parameters convergence were checked through Tracer v 1.5 (Rambaut & Drummond, 2010), and the first 25% generations were discarded as burn-in.
Clades were considered well supported when ML bootstrap support was ≥ 70, and BI posterior probability was ≥ 0.95. Phylogenetic trees were visualized in FigTree v1.4.3 (Rambaut, 2012) and edited in Adobe Illustrator CS6 software (Adobe Systems, USA).
Prevalence of Fusarium species and Seed-plant transmission disease bioassay
The prevalence of Fusarium species obtained was determined by calculating the Isolation Rate (IR) as proposed by Veloso et al. (2018): IR (%) = (Cx / Ct) x 100, where Cx is the number of isolates belonging to one species and Ct is the number of isolates per cultivar.
A randomized design bioassay was carried out to access the transmission of Fusarium-associated diseases from sunflower seeds to growing plants. Sunflower seeds, cv BRS 324, from the 2018/19 crop season, were previously tested to measure seed germination and sanity status according to Brasil (2009). They exhibited 93% germination with no incidence of Fusarium species, thus being eligible for the transmission bioassay performed following the methodology described by Sousa et al. (2008). Briefly, 100 surface disinfected (as above) sunflower seeds were distributed over 7 days old colonies of Fusarium representative isolates grown at 25 ± 2°C and 12:12 L:D daylength on PDA amended with mannitol to a -1.0 MPa osmotic potential. Another set of 100 disinfected seeds distributed over PDA + mannitol without fungi served as the control treatment. After exposure for 48 h, the seeds were sown 2.0 cm deep into plastic vases (5.0 L of volume) filled with the sterile commercial substrate Brasplant® and vermiculite (2:1 v/v). They were watered twice a day for 35 days inside a greenhouse with night temperature about 16–18 ºC and day temperature about 28–32 ºC.
Symptomatic growing plants were tallied throughout the experiment, whereas ungerminated seeds were surface disinfected (as described above) and kept into a humid chamber for seven days to observe whether Fusarium species were present. By the 35th day after sowing, all plants were harvested and roots were cleaned with tap water. From each plant, there were taken three fragments of roots, stems, and leaves, which were surface disinfected by submersing in 70% ethanol for 30 s, 1.5% NaClO aqueous solution for 1 min, rinsed three times with sterile water, and dried on sterilized filter paper. These fragments were distributed on PDA incubated at 25 ± 2°C and 12:12 L:D daylength for seven days, and the transmission of Fusarium-associated seed-borne diseases was confirmed whenever typical structures of Fusarium grew over them.
The incidence of Fusarium-associated seed-borne diseases on sunflower was calculated as follows: Disease Incidence (%) = (n / N) × 100, where n is the number of symptomatic plants + the number of ungerminated seeds contaminated with Fusarium spp., and N is the total number of examined plants and seeds. The transmission rate (TR) of Fusarium spp. from sunflower seeds to growing plants was calculated after adapting an equation proposed by Teixeira and Machado (2003): TR (%) = (IR / IS) x 100, where IR is the infection rate of Fusarium spp. on sunflower seeds, roots, stems, and/or leaves fragments and IS is the incidence of each isolate on artificially inoculated seeds.
Statistical analysis
Data from mycelial growth, conidia production, and transmission bioassays were checked for normality and homoscedasticity of variance using the Shapiro-Wilk’s and Levene’s tests, respectively. Dependent variables that passed these tests were analyzed by one-way ANOVA and means were grouped by Scott-Knott test (p < 0.05) using the R software statistical package (R Core Team, 2020).