Skip to main content

ORIGINAL RESEARCH article

Front. Plant Sci., 25 November 2022
Sec. Plant Pathogen Interactions
This article is part of the Research Topic Disease and Pest Resistance in Rice View all 6 articles

Genetic diversity and structure of Bipolaris oryzae and Exserohilum rostratum populations causing brown spot of rice in Burkina Faso based on genotyping-by-sequencing

Kouka Hilaire Kabor,*Kouka Hilaire Kaboré1,2*Abalo Itolou KassankognoAbalo Itolou Kassankogno3Henri Adreit,Henri Adreit2,4Joëlle Milazzo,Joëlle Milazzo2,4Sonia Guillou,Sonia Guillou2,4Laurence Blondin,Laurence Blondin2,4Laurie Chopin,Laurie Chopin2,4Sbastien Ravel,Sébastien Ravel2,4Florian Charriat,Florian Charriat2,4Mariam Barro,Mariam Barro3,4Charlottte TollenaereCharlottte Tollenaere4Marc-Henri LebrunMarc-Henri Lebrun1Didier Tharreau,Didier Tharreau2,4
  • 1Université Paris-Saclay, INRAE, UR BIOGER, Palaiseau, France
  • 2CIRAD, UMR PHIM, Montpellier, France
  • 3CNRST/INERA, Laboratoire de Phytopathologie, Bobo-Dioulasso, Burkina Faso
  • 4PHIM, Univ Montpellier, CIRAD, INRAE, IRD, Montpellier SupAgro, Montpellier, France

In recent years, Brown spot disease of rice (BSR) has been observed on leaves and seeds of rice in all rice-growing areas of Burkina Faso. Bipolaris oryzae and Exserohilum rostratum are the main fungal species isolated from BSR infected tissues and they are frequently observed in the same field. However, we are lacking information on the genetic diversity and population structure of these fungi in Burkina Faso. The mode of reproduction is also unknown. The genetic diversity of isolates of B. oryzae (n=61) and E. rostratum (n=151), collected from major rice-growing areas of Burkina Faso, was estimated using genotyping-by-sequencing (GBS). The mean values for nucleotide diversity (π) were 1.9 x10-4 for B. oryzae and 4.8 x10-4 for E. rostratum. There is no genetic differentiation between the geographical populations of each species. The analysis of molecular variance revealed that 89% and 94% of the genetic variances were within the populations of B. oryzae and E. rostratum, respectively. For each species, four genetic clusters were identified by two clustering methods (DAPC and sNMF). The distribution of these genetic groups was independent of the geographical origin of the isolates. Evidence of recombination was detected in the populations of B. oryzae and E. rostratum. For B. oryzae balanced mating type ratios were supporting sexual reproduction. For E. rostratum overrepresentation of MAT1-2 isolates (79%) suggested a predominant asexual reproduction. This study provides important information on the biology and genetics of the two major fungi causing brown spot disease of rice in Burkina Faso.

Introduction

Rice is the second most important source of food after corn in Sub-Saharan Africa (SSA, Nigatu et al., 2017). However, SSA rice cropping is characterized by a low average yield of 1.4 t/ha, compared to over 4 t/ha in Asia (Nwanze et al., 2006). Local rice production cannot meet the needs, leading to import up to half 50% of the rice consumed. The major constraints significantly reducing rice production in SSA are drought, iron toxicity, salinity, weeds, rodents, birds and diseases (Diagne et al., 2013; van Oort, 2018).

Among diseases, Brown spot of rice (BSR) is becoming a threat in SSA. This disease is expanding, and is currently estimated to be responsible for the loss of 3.3% of the regional rice production (Savary et al., 2019). BSR occurs as brown spots on coleoptiles, leaf blades, sheaths, stems, glumes and seeds (Ou, 1985; Sunder et al., 2014). The disease is favored by mineral deficiencies in the soil, especially nitrogen, high temperatures and lack of water (Zadoks, 1974; Ou, 1985; Savary et al., 2000; Savary et al., 2005).

In Burkina Faso, the importance of BSR in rice nurseries and fields was already highlighted in 1962 (Davet, 1962). An outbreak of BSR was recorded in 1994 in the most important rice-growing plain of the country (Sourou) and, since then, the disease has spread to other rice-growing areas (Ouédraogo, 2008). Under natural disease pressure, grain yield losses could reach 16% (Ouédraogo, 2008). A recent study has shown that 80% of rice fields in western Burkina Faso displayed BSR leaf symptoms (Barro et al., 2021). Furthermore, up to 76% of rice grains produced in Burkina Faso are infected (Ouédraogo et al., 2016). This high prevalence leads to reconsider the status of BSR and justifies to improve knowledge of the pathogen biology to ameliorate control methods.

Several species of fungi such as Bipolaris oryzae, B. bicolor, B. indica, B. victoriae, B. zeicola (Motlagh and Kaviani, 2008; Aslam et al., 2021) and Exserohilum rostratum (Majeed et al., 2016; Toher et al., 2016; Kaboré et al., 2022) are responsible for BSR. In Burkina Faso, Bipolaris oryzae (syn. Cochliobolus miyabeanus) and Exserohilum rostratum (syn. Setosphaera rostrata), are the two species most frequently isolated from BSR infected tissues (Kaboré, 2022). These species are frequently found in the same rice fields or in the same seed lots (Ouédraogo et al., 2016; Kaboré, 2022).

Bipolaris oryzae and E. rostratum are heterothallic fungi, that is, their sexual reproduction requires strains of opposite mating types (MAT1-1 and MAT1-2) (Ueyama and Tsuda, 1976). Most studies on B. oryzae populations have reported a high level of genetic and genotypic diversity and low geographical population structure (Castell-Miller and Samac, 2012; Archana et al., 2014a; Ahmadpour et al., 2018). Sexual reproduction of B. oryzae was suspected based on equal frequencies of mating types, but recombination was not confirmed by estimating linkage disequilibrium (Castell-Miller and Samac, 2012; Ahmadpour et al., 2018). Bipolaris oryzae isolates from Burkina Faso showed significant morphological and pathological variability (Ouédraogo et al., 2004; Dembélé, 2014; Boka et al., 2018). Investigations based on mitochondrial DNA RFLP markers revealed the division of B. oryzae isolates into two haplotypes (Ouédraogo et al., 2004).

In the field, the infection cycle of B. oryzae starts with the deposition of asexual spores (conidia), on aerial parts of rice. Conidia produce a germ tube that adhere to the plant surface thanks to a mucilaginous substance. Then, an appressorium is formed that allows direct penetration of cells (Hau and Rush, 1982; Ou, 1985). After tissue colonization by mycelial hyphae, lesions appear 18 to 36 hours post-infection (Tullis, 1935; Johnson and Percich, 1992; Dallagnol et al., 2009). In cultivated rice, a peak of conidia production from new lesions is observed 6 days after infection (Ou, 1985; Barnwal et al., 2013). Less is known about the infection cycle of E. rostratum on rice, but a recent histopathological study revealed that on contact with rice leaves, E. rostratum conidia produce appressoria within 24 h that directly penetrate the epidermal cells with the onset of symptoms within 3 days (Korra et al., 2022).

To date, the genomics studies of B. oryzae and E. rostratum are also limited. The genome sequence of a strain of E. rostratum isolated from maize (BioSample: SAMN19599129) was released recently. The description of the complete circular mitochondrial genome of E. rostratum (ZM170581) isolated from maize followed (Ma et al., 2022). Genomes of E. rostratum strains (BioSamples: SAMN14931704 and SAMN30522240) isolated from rice have also been sequenced recently (https://www.ncbi.nlm.nih.gov/biosample/). For B. oryzae a reference genomes are available from a strain isolated from cultivated rice (Condon et al., 2013) and a strain isolated from wild rice (Zizania palustris) (Castell-Miller et al., 2016). These data were used for genomic structure studies (Condon et al., 2013; Ma et al., 2022).

However, there is a lack of information on the genetic structure and diversity, as well as the mode of reproduction of B. oryzae in rice fields in Burkina Faso. Our knowledge of E. rostratum in Burkina Faso is limited since no population study was performed to date. The objectives of this study were to: (i) assess the genetic diversity, (ii) describe the genetic structure, and (iii) determine the mode of reproduction of B. oryzae and E. rostratum populations from major rice-growing areas of Burkina Faso.

Materials and methods

Sample collection and isolation of strains

Leaves and seeds of cultivated rice (Oryza sativa) showing symptoms of BSR were collected in 2018 and 2019 in different rice fields from Bagré, Bama, Banfora, Banzon and Sourou (Figure 1; Supplementary Tables S1, S2). These sites are on average 365 km apart (minimum distance 59 km and maximum distance 622 km) and represent the main rice-growing areas of Burkina Faso. Diseased leaves and seeds were arranged in Petri dishes containing a double layer of moistened blotting paper (Mathur and Kongsdal, 2003). Incubation was carried out at 25°C for 7 days under alternating 12 h of white light and 12 h of darkness. Using a flame-stretched Pasteur pipette, fungal conidia formed in the lesions were collected and spread on bacto-agar medium. After incubation for 24 h (25°C, 12 h white light and 12 h darkness), a single germinating conidium was transferred with a sterile scalpel onto potato dextrose agar medium (PDA, 39 g/L) and incubated under the same conditions for 7 days.

FIGURE 1
www.frontiersin.org

Figure 1 Schematic location of collection sites for infected rice samples. The number of Bipolaris oryzae and Exserohilum rostratum isolates isolated is indicated by the heptagons and stars respectively. Names in bold refer to the localities where samples were collected and the dashed lines indicate the distance between collection sites.

DNA extraction

Genomic DNA was extracted from fresh mycelium with modifications to the cetyl trimethyl ammonium bromide (CTAB) method described by Doyle and Doyle (1987). A 5-day-old fragment of mycelium was transferred onto Corn Meat Agar (CMA, 18 g/L) medium previously covered with a sterile cellophane disc. The plates were incubated at 25°C (12 h photoperiod) for 7 days, the mycelium was scraped off with a scalpel and approximately 30 mg of the mycelium was transferred to a 2 mL tube containing sterile sand and 2 iron beads. These tubes were placed directly into liquid nitrogen and the mycelium was cold ground using a Retsch mill (Retsch MM 301) for 4 min at 15 Hz. After grinding, each tube received 1 mL of CTAB lysis buffer with 1% polyvinyl pyrrolidone (MW. 360 000), 1% Na2SO3, 0.2 mg/mL RNase and 1 mg/mL proteinase K. The tubes were vortexed (Vortex-Genie 2) and incubated in the oven at 65°C for 2 h. The tubes were then centrifuged for 3 min at 8,000 rpm and approximately 900 µL of the supernatant was transferred to a new 2 mL tube. One mL of chloroform isoamilalcohol (CIAA, 24/1 v/v) was placed in each tube and mixed by inversion before being centrifuged at 10,000 rpm for 10 min. After centrifugation, 700 µL of the supernatant was transferred to 1.5 mL tubes. Six hundred µL of isopropanol cooled to -20°C and 210 µL of 5M NaCl were added and mixed by inversion. The tubes were kept at -20°C overnight. They were then centrifuged at 4°C for 30 min at 13,000 rpm and the supernatant was removed by inverting the tubes. Five hundred µL of 70° ethanol was added to each tube and centrifuged at 4°C for 5 min at 13,000 rpm. Ethanol was removed by inverting the tubes and draining on absorbent paper. The tubes were dried for 1 h in the fume cupboard. DNA was re-suspended in 100 µL of sterile milliQ water. The quality and quantity of the extracted DNA was assessed by electrophoresis in 1% (w/v) agarose gels and spectrophotometry (Nanodrop™ 2000).

Preparation of genotyping-by-sequencing libraries

The preparation of GBS libraries was carried out in collaboration with the genotyping platform of the UMR AGAP/CIRAD according to the method described by Elshire et al. (2011). Genomic DNA of the strains were normalized to a concentration of 10 ng/µL and digested with the restriction enzyme ApeKI (5 U) at 75°C for 2 h. The restriction digestion of DNA was followed by ligation with adapters. The adapters included different barcodes for tagging individual samples and common adapters. The ligation was performed using T4 DNA ligase at 22°C for 30 min and the ligase was inactivated by holding at 65°C for 10 min. Two pools were created. To assess the reproducibility of GBS, each pool contained 6 randomly selected and duplicated isolates. Ligation products were purified with the QIAquick purification kit from QIAgen. The library was enriched by PCR to 50 µL of reaction solution including 25 µL of 2X Taq Master mix NEB, 19 µL of water, 2 µL of DNA and 2 μL each of the forward and reverse primers (10 μM). Amplification was done according to the following program: 98 ˚C for 30 s, followed by 18 cycles of 98 ˚C for 10 s, 65 ˚C for 30 s, 72 ˚C for 30 s, ending at 72 ˚C for 5 min. PCR products were purified with QIAgen’s QIAquick Purification Kit before being evaluated with BioAnalyzer 2100 (Agilent Technologies) for fragment size distribution. The GBS libraries were sequenced on an Illumina HiSeq 4000 producing 150 bp paired-end reads.

Sequence preprocessing and SNP identification

The bioinformatics tool RattleSNP (https://github.com/sravel/RattleSNP) was used to process the data. The workflow performs demultiplexing, sequence mapping to the reference genome and variant calling. The sequences of B. oryzae isolates were mapped onto the genome of strain ATCC 44560 (Condon et al., 2013), and those of E. rostratum onto the genome of BF9006 sequenced by CIRAD (https://www.ncbi.nlm.nih.gov/biosample/SAMN30522240/). Prior to variant calling, strains with a mapping rate of less than 75% to the reference genome, a library size of less than 100,000 reads and a sequencing depth of less than 3 x were considered low quality data and excluded from the dataset. SNPs from strains meeting the criteria (61 isolates for B. oryzae and 151 isolates of E. rostratum) were filtered for missing rate (–max-missing 0.70), minor allele frequency (–maf 0.05), depth (–minDP 2) and site quality (–minQ 30) using VCFtools version 01.16 (Danecek et al., 2011)

Population structure and genetic differentiation analyses

Discriminant Analysis of Principal Components (DAPC) (Jombart et al., 2010) and sNMF (Frichot et al., 2014) analysis were used to investigate the genetic structure of these fungal populations. DAPC was carried out with the adegenet package of the R software. First, the find.clusters function was used to detect the number of clusters in the population (Jombart, 2008; Jombart et al., 2010). The function uses K-means clustering which decomposes the total variance of a variable into two components: inter-group and intra-group. This model maximizes the inter-group variation (Jombart et al., 2010).

sNMF of the LEA R software package was used to estimate the number of discrete genetic clusters (K) and admixture coefficients (Frichot et al., 2014). The optimal number of sNMF clusters was defined as the one corresponding to the lowest cross-entropy value.

To further characterize the genetic structure of the population described by the DAPC and sNMF analyses, a hierarchical analysis of molecular variance (AMOVA) (Excoffier et al., 1992) were performed with the poppr package of R (Kamvar et al., 2014). Pairwise genetic differentiation (FST) was calculated using the ‘genet.dist’ function in the hierfstat package of R (Goudet, 2005), following the method described by Weir and Cockerham (1984). RAxML (Randomized Axelerated Maximum Likelihood) version 8.2.12 was used to assess phylogenetic relationships between isolates using the maximum likelihood method (Stamatakis, 2014).

Genetic and genotypic diversity analyses

Analysis of genotypic diversity and association indices of fully or partially clonal populations is recommended before and after clone correction (Milgroom, 1996; Grünwald and Hoheisel, 2006; Kamvar et al., 2014). To identify clones, the dissimilarity rate between duplicates of the same isolate was calculated using the bitwise.dis function of the R poppr package (Kamvar et al., 2014). For a given species, the maximum dissimilarity between duplicates was measured and used as a threshold to define clones. Two isolates with a dissimilarity rate lower than the threshold were considered as clones.

To estimate genotypic diversity, the poppr package of R (Kamvar et al., 2014), was used to calculate Shannon-Wiener index (H) (Shannon, 1948), Simpson’s complement index of genotypic diversity index (λ) corrected for sample size (Simpson, 1949; Grünwald et al., 2003) and evenness (E5) (Pielou, 1975; Ludwig and Reynolds, 1988; Grünwald et al., 2003). To adjust for the population size, the raw value of the Simpson’s complement index of genotypic diversity index was multiplied by N/(N-1) with N being the number of isolates in the population (Grünwald et al., 2003).

To assess gene diversity, nucleotide diversity (π) values were calculated with VCFtools across the genome with 10 kb windows (Danecek et al., 2011).

Linkage disequilibrium and recombination

Detection of recombination was carried out by calculating linkage disequilibrium (LD) and clonality was tested with the PHI test from genotyping data. Pairwise LD between loci was calculated based on allelic frequency correlations (r2) over a 10 kb window using the PopLDdecay program (Zhang et al., 2019).

A FASTA file was used to generate a NeighborNet in SPLITSTREE4 and the PHI test was performed to test the null hypothesis of clonality (Huson and Bryant, 2006).

Mating type characterization and in vitro sexual reproduction

Identification of mating types was done by amplification of MAT1-1 and MAT1-2 genes with primers specific to B. oryzae (Ahmadpour et al., 2018) and E. rostratum (Kusai et al., 2016). Each PCR reaction using the kit Go-Taq (Promega) had a final volume of 25 μL, containing 1 μL of template genomic DNA (20 ng/µL),1 μL forward primer (10 µM), 1 μL reverse primer (10 µM), 1µL of dNTPs (10 µM each), 5 μL of the buffer 5x, 0.25 µL Taq polymerase (5U/µL) and 14.75 μL sterile water. For the amplifications of E. rostratum MAT genes, the amount of water was 13.75 µL and 1 µL MgCl2 5X was added. PCR cycling was performed as follows: 94°C for 2 min; followed by 35 cycles of 94°C for 1 min, 52°C for 1 min and 72°C for 1 min; plus a final extension at 72°C for 7 min for B. oryzae. For E. rostratum, initial denaturation at 95°C for 2 min; followed by 35 cycles of denaturation at 95 C for 30 s, annealing at 60°C for 30 s, extension at 72°C for 1 min; and final extension at 72°C for 5 min. The null hypothesis of a 1:1 ratio between mating types was tested using the chi-square goodness-of-fit.

The sexual fertility of B. oryzae and E. rostratum isolates was tested in vitro according to the protocol described by Tsuda and Ueyama (1976). A dry, autoclaved rice leaf of about 4 cm was placed in the center of a petri dish on Sach’s culture medium. Pieces of mycelium from the pre-determined MAT1-1 and MAT1-2 strains were placed on either side of the rice leaf. The dishes were sealed with tape to prevent dehydration of the medium and then placed in an oven at 24°C for 30 days.

Results

Genetic structure of B. oryzae and E. rostratum populations

The genetic diversity of isolates of B. oryzae (n=61) and E. rostratum (n=151) collected in Burkina Faso, was estimated using genetic markers (SNPs) obtained by genotyping-by-sequencing. After alignment to the B. oryzae ATCC 44560 reference genome composed of 619 scaffolds and filtering, 28,907 SNPs were identified in the B. oryzae samples (61 isolates). Sequences of E. rostratum isolates (n=151) aligned to the BF9006 isolate reference genome (377 scaffolds) yielded 102,464 SNPs.

To study the B. oryzae and E. rostratum populations from Burkina Faso, the isolates were clustered by DAPC according to their multilocus genotypes. Inference of ancestry was also performed with sNMF.

For B. oryzae, the bayesian information criterion (BIC) of DAPC suggested the best clustering for four genetically distinct populations (K = 4; Supplementary Figure S1, Figure 2A). Similarly, the sNMF analysis indicated an optimum for four clusters (Supplementary Figure 2, Figure 3A). The phylogenetic tree made with RaxML also highlighted four clades (Supplementary Figure S3). The clustering of isolates according to DAPC, sNMF and RaxML was similar. Isolates from cluster 1 and 2 were found in all five geographical collection areas (Table 1, Figure 4), while cluster 4 was only found in Banfora. Isolates from cluster 3 were detected in three areas namely Bagré, Bama and Sourou. The sNMF analysis showed a significant admixture between the four clusters, but particularly between cluster 1 and 2 (Figure 3A). The FST values between the B. oryzae isolates from the different rice-growing areas in Burkina Faso showed little differentiation between sites, ranging from 0.005 to 0.132 (Table 2) and with an average value of 0.08. No correlation was observed between geographic distances of rice growing sites and genetic differentiation index values (Supplementary Figure S4). The AMOVA analysis revealed that the variability was 88.65% within localities and 11.35% among the five rice-growing zones (Table 3).

FIGURE 2
www.frontiersin.org

Figure 2 Discriminant analysis of principal components (DAPC) scatterplot drawn using SNPs in the R package adegenet (K =4). (A) Bipolaris oryzae, (B) Exserohilum rostratum. Each color represents a cluster obtained by DAPC analysis of isolates from different localities in Burkina Faso. For Bipolaris oryzae, Cluster 1 is composed of 27 isolates (Bagre = 4, Bama = 3, Banfora = 5, Banzon = 8 and Sourou =7), Cluster 2 is composed of 13 isolates (Bagre = 4, Bama = 3, Banfora = 1, Banzon = 2 and Sourou =3), Cluster 3 is composed of 16 isolates (Bagre = 1, Bama = 4 and Sourou =11), Cluster 4 is composed of 5 isolates (Banfora = 5). For Exserohilum rostratum, Cluster 1 is composed of 58 isolates (Bagre = 38, Bama = 4, Banfora = 4, Banzon = 2 and Sourou =10), Cluster 2 is composed of 60 isolates (Bagre = 23, Bama = 4, Banfora = 7, Banzon = 2 and Sourou =24), Cluster 3 is composed of 26 isolates (Bagre = 5, Bama = 5, Banfora = 2 and Sourou =14) and Cluster 4 is composed of 7 isolates (Bagre = 3, Banzon = 1 and Sourou = 3).

FIGURE 3
www.frontiersin.org

Figure 3 Population structure analysis of Bipolaris oryzae and Exserohilum rostratum. Ancestry coefficient for each isolate calculated by using sNMF with K=4. (A) B. oryzae, (B) E. rostratum..

TABLE 1
www.frontiersin.org

Table 1 Distribution of Bipolaris oryzae and Exserohilum rostratum clusters obtained by DAPC analysis in five rice-growing areas of Burkina Faso.

FIGURE 4
www.frontiersin.org

Figure 4 Geographic distribution of the four genetic clusters of Bipolaris oryzae in Burkina Faso. The size of sectors is proportional to the proportion each cluster defined by DAPC in each area. Cluster 1 (n = 27 isolates), Cluster 2 (n = 13 isolates), Cluster 3 (n = 16 isolates) and Cluster 4 (n = 5 isolates).

TABLE 2
www.frontiersin.org

Table 2 Pairwise FST estimates between Bipolaris oryzae isolates from five rice-growing zones in Burkina Faso.

TABLE 3
www.frontiersin.org

Table 3 Analysis of molecular variance (AMOVA) of Bipolaris oryzae among and within five rice-growing zones in Burkina Faso.

For E. rostratum, the bayesian information criterion (BIC) of DAPC suggested the best clustering for four genetically distinct populations (K = 4; Supplementary Figure 5, Figure 2B). Similarly, the sNMF analysis indicated an optimum of four clusters (Supplementary Figure S6; Figure 3B). The phylogenetic tree made with RaxML also highlighted four clades (Supplementary Figure S7). The clustering of isolates was identical between methods. Isolates from clusters 1 and 2 were found in all five localities sampled. Isolates from cluster 3 were found in all localities (Figure 5) except Banzon, and isolates cluster 4 was sampled in Bagré, Banzon and Sourou. Very low FST values ranging from 0 to 0.129 and with an average value of 0.022 were observed between populations of different rice production areas in Burkina Faso (Table 4). Distances between the different rice production sites had no impact on population structure (Supplementary Figure S8). The AMOVA analysis revealed that the variability was 94.4% within and 5.6% among the five rice-growing areas (Table 5).

FIGURE 5
www.frontiersin.org

Figure 5 Geographic distribution of the four genetic clusters of Exserohilum rostratum in Burkina Faso. The size of sectors is proportional to the proportion each cluster defined by DAPC in each area. Cluster 1 (n = 58 isolates), Cluster 2 (n = 60 isolates), Cluster 3 (n = 26 isolates) and Cluster 4 (n = 7 isolates).

TABLE 4
www.frontiersin.org

Table 4 Pairwise FST estimates between Exserohilum rostratum isolates from five rice-growing zones in Burkina Faso.

TABLE 5
www.frontiersin.org

Table 5 Analysis of molecular variance (AMOVA) of Exserohilum rostratum among and within five rice-growing zones in Burkina Faso.

Genotypic and genetic diversity

Samples were analyzed for possible clones. The genetic distance between duplicates of B. oryzae isolates varied from 0.035 to 0.087. Those of E. rostratum varied between 0.015 and 0.033. Bipolaris oryzae and E. rostratum isolates were considered as clones if their dissimilarity rates were less than or equal to 0.087 and 0.033 respectively. Based on this criterion, 8 and 11 isolates of B. oryzae and E. rostratum respectively were considered as clones of two or more isolates. Clones represented 13% (8/61) and 7% (11/151) of the isolates in each sample of B. oryzae and E. rostratum, respectively (Table 6).

TABLE 6
www.frontiersin.org

Table 6 Genotypic and gene diversity indices of Bipolaris oryzae and Exserohilum rostratum isolates from five rice-growing areas of Burkina Faso.

Genotypic and gene diversity indexes are presented in Table 6. For B. oryzae, the genotypic diversity (H) between isolates from different rice-growing areas ranged from 2.20 to 3.04 with an average value of 2.45. These indices are influenced by population size and do not allow comparison between different localities. The Simpson index corrected for sample size ranged from 0.999 to 1 (mean value 0.999) indicating a high level of genotypic diversity within each population. The regularity index (E.5) was equal to 1 for all populations.

In terms of genetic diversity, nucleotide diversity values ranged from 1.3 x10-4 to 2.4 x10-4 with an average of 1.9 x10-4. The Banzon population also showed the lowest value of nucleotide diversity (1.3 x10-4). The Sourou population showed the highest diversity (2.4 x10-4). It was followed by the populations of Bama, Bagré and Banfora.

For E. rostratum, the genotypic diversity (H) between isolates from different rice-growing areas ranged from 1.61 to 4.23 with an average value of 2.98. The Simpson index corrected for sample size ranged from 0.999 to 1 (mean value 0.999) indicating high genotypic diversity within each population. The regularity index (E.5) was equal to 1 for all populations. The values of nucleotide diversity varied between 5.7 x10-4 and 9.8 x10-4. The average value was 4.8 x10-4. The isolates from Bama showed the highest value of nucleotide diversity (9.8 x10-4), while the lowest value was observed in Bagré (5.7 x10-4).

Recombination, mating type distribution and in vitro sexual reproduction

In B. oryzae populations from Burkina Faso, Linkage Disequilibrium (LD) pattern showed a rapid decay to less than 1 kb in all populations (Figure 6) suggesting high rate of recombination. The PHI coefficient calculated on the populations with sufficient number of isolates, rejected the hypothesis of clonality (p = 0.0). Finally, a reticulated NeighborNet was observed (Figure 7A), also supporting the hypothesis of recombination in populations of B. oryzae from Burkina Faso. Finally, the mating types MAT1-1 and MAT1-2 were observed in balanced proportions (Table 7).

FIGURE 6
www.frontiersin.org

Figure 6 Estimated linkage disequilibrium decay (LD-decay) of Bipolaris oryzae and Exserohilum rostratum isolates. (A) LD of B. oryzae isolates per locality, (B) LD of B. oryzae isolates per cluster, (C) LD of E. rostratum isolates per locality, (D) LD of E. rostratum isolates per cluster. Localities and genetic clusters with less than 8 isolates were not presented individually (ie Banzon for E. rostratum, cluster 4 for B. oryzae and E. rostratum).

FIGURE 7
www.frontiersin.org

Figure 7 NeighborNet phylogenetic network of Bipolaris oryzae and Exserohilum rostratum isolates from Burkina Faso estimated with SPLITSTREE. (A) Bipolaris oryzae, (B) Exserohilum rostratum. Each color represents a genetic group previously defined by DAPC. For Bipolaris oryzae, Cluster 1 is composed of 27 isolates (Bagre = 4, Bama = 3, Banfora = 5, Banzon = 8 and Sourou =7), Cluster 2 is composed of 13 isolates (Bagre = 4, Bama = 3, Banfora = 1, Banzon = 2 and Sourou =3), Cluster 3 is composed of 16 isolates (Bagre = 1, Bama = 4 and Sourou =11), Cluster 4 is composed of 5 isolates (Banfora = 5). For Exserohilum rostratum, Cluster 1 is composed of 58 isolates (Bagre = 38, Bama = 4, Banfora = 4, Banzon = 2 and Sourou =10), Cluster 2 is composed of 60 isolates (Bagre = 23, Bama = 4, Banfora = 7, Banzon = 2 and Sourou =24), Cluster 3 is composed of 26 isolates (Bagre = 5, Bama = 5, Banfora = 2 and Sourou =14) and Cluster 4 is composed of 7 isolates (Bagre = 3, Banzon = 1 and Sourou = 3).

TABLE 7
www.frontiersin.org

Table 7 Distribution of mating types in B. oryzae and E. rostratum populations in five rice growing areas of Burkina Faso.

In E. rostratum, LD showed a rapid decay to less than 1 kb in all populations (Figure 6) and PHI coefficient rejected the hypothesis of clonality (p = 0.0), supporting the hypothesis of recombination. Low reticulated NeighborNet was observed (Figure 7B). The dominance (79%) of only one mating type (MAT1-2) was observed and the null hypothesis of a 1:1 ratio was rejected by statistical analysis. Ten strains repeatedly amplified markers of the two mating types (Table 7). Potential mixture of strains was discarded by repeating monoconidial isolation and potential contamination was discarded by repeated independent DNA extractions.

To test the sexual fertility of B. oryzae and E. rostratum isolates, crosses between isolates of opposite mating type from the same rice field were performed in vitro. Out of a total of 156 in vitro crosses for B. oryzae, two pairs of isolates (BF9259 and BF9268; BF9294 and BF9296) produced pseudothecia containing ascospores. Concerning E. rostratum, out of 364 crosses made, no pseudothecia were observed.

Discussion

Brown spot of rice (BSR) is an emerging disease in Burkina Faso. Although it was already described, the incidence of the disease increased very recently to reach 80% of farmer’s fields from the western part of the country (Barro et al., 2021). Knowledge of the structure and genetic diversity of the fungal species predominantly responsible for the disease is an important information for disease management and to understand why the disease emerged. This study is the first comprehensive analysis of the genetic diversity of B. oryzae isolated from rice in Burkina Faso and the first study of E. rostratum populations. It is also the first study using genotyping-by-sequencing based analysis to characterize populations of B. oryzae and E. rostratum.

The B. oryzae and E. rostratum populations were sampled in five distinct rice-growing areas of Burkina Faso with an average distance of 365 km between sites (Bagré, Bama, Banfora, Banzon and Sourou). For each fungal species, the five geographic populations did not display a significant genetic differentiation. The absence of genetic differentiation among B. oyzae geographical populations was already observed in Bangladesh (Kamal and Mia, 2009), India (Archana et al., 2014a; Archana et al., 2014b), Iran (Nazari et al., 2015) and Thailand (Chaijuckam et al., 2019). In contrast to this situation, Burgos et al. (2013) established a relationship between the population structure and the geographical origin of isolates from different provinces in the Philippines. Although airborne spores of B. oryzae are important in the initiation of infections (Picco and Rodolfi, 2002; Chakraborty et al., 2020), their spatial spread seems relatively limited (about 5 m) (Chakraborty et al., 2020). It is therefore possible that the dissemination of B. oryzae isolates is mainly carried out by seeds. This hypothesis likely applies to E. rostratum that is frequently sampled from infected rice seeds (Cardona and Gonzáles, 2007; Ouédraogo et al., 2016; Silva et al., 2016; Imrani et al., 2019; Kaboré et al., 2022). In Burkina Faso there are significant exchanges of rice seed between producers across rice production areas. This agricultural practice likely supported the dissemination of B. oryzae and E. rostratum isolates across areas and could account for the lack of relationship between genetic groups and geographic origin of isolates.

The genetic diversity of B. oryzae has been studied using several molecular tools in some countries (Kamal and Mia, 2009; Castell-Miller and Samac, 2012; Archana et al., 2014a; Archana et al., 2014b; Nazari et al., 2015; Ahmadpour et al., 2018; Chaijuckam et al., 2019). Most of them have revealed that these populations had a high level of genetic and genotypic diversities. Bipolaris oryzae populations in the five rice-growing areas of Burkina Faso also displayed a high genotypic diversity, equal frequencies of MAT1-1 and MAT1-2 mating types and genomic evidence of recombination. Balanced frequencies of mating types are expected in populations of bipolar heterothallic species experiencing sexual reproduction because mating type is under the control of a single locus with two pseudo-alleles (idiomorphs). So, a 1:1 segregation of mating types is expected after sexual reproduction. All these genetic and biological characteristics are consistent with sexual reproduction. However, sexual reproduction of B. oryzae was seldom observed in the crosses we performed in vitro and pseudothecia (sexual reproductive structure) have never been observed in the field. Timing sexual reproduction phenomena remains problematic and it is difficult to know if they are still ongoing or if they occurred in the past (Castell-Miller and Samac, 2012). The hypothesis of sexual reproduction of B. oryzae isolates from Burkina Faso will have to be verified by repeating in vitro reproduction tests and actively searching for pseudothecia in rice fields.

In contrast to B. oryzae populations, E. rostratum populations showed a dominant distribution of a single mating type, MAT1-2 (79%). These results are similar to those of Kusai et al. (2016), observed on E. rostratum isolates from rice in Malaysia. The linkage disequilibrium and homoplasy tests did not reject the hypothesis of recombination within the E. rostratum populations of Burkina Faso. However, sexual reproduction of E. rostratum was not observed in the crosses performed in vitro and pseudothecia have never been observed in the field. In addition, the strong bias toward one mating type (79%) suggested that actual E. rostratum populations are predominantly asexual. In addition, we identified 10 isolates of E. rostratum, i.e. 6% of our sample that carried both the MAT1-1 and MAT1-2 idiomorphs. Such isolates with two mating types were also identified within E. turcicum populations from maize in China (Fan et al., 2007). Sequencing the genome of the isolates with both mating types should be used to confirm that they carry both idiomorphs (mero-diploid at these loci).

This study provides new knowledge useful to adapt disease management. First, during this study we have shown that two species are causing BSR in Burkina Faso and that E. rostratum is more frequently encountered than B. oryzae. This information is important to consider for breeding for resistance: screening of new rice varieties should be done with strains of both species. Second, the fact that there is no obvious population structure at the country suggests that there is gene flow on a broad geographic scale. We hypothesize that this is due to transportation of infected seeds. If true, a prophylactic method based on healthy seeds could have an impact. Third, the absence of geographic structure also informs on the strategy of breeding for and deployment of resistant varieties. Since the populations of both pathogenic species are not differentiated at the country scale, a resistant variety effective in a region is expected to be effective in another and selection for resistance in different regions is not mandatory. Conversely, the existence of genetic clusters in the two B. oryzae and E. rostratum probably require to screen for resistance with representative strains of these genetic groups. Finally, knowing that recombination is frequent must draw the attention on the adaptive potential of both species since favorable alleles controlling, for example, virulence to a particular resistance gene or resistance to fungicides, will easily be transmitted.

Data availability statement

The data presented in the study are deposited in the European Variation Archive (EVA) repository, accession number PRJEB56783.

Author contributions

KK: Collection of samples, Data production, Data analysis, Writing and proofreading AK: Collection of samples, Writing and proofreading. HA: Data production, Writing and proofreading. JM: Data production, Writing and proofreading. SG: Data production, Writing and proofreading. LB: Data production, Writing and proofreading. LC: Data production, Writing and proofreading. SR: Data analysis, Writing and proofreading. FC: Data analysis, Writing and proofreading. MB: Collection of samples, Writing and proofreading. CT: Collection of samples, Writing and proofreading. M-HL: Writing and proofreading. DT: Data analysis, Writing and proofreading. All authors contributed to the article and approved the submitted version.

Acknowledgements

We thank the RICE CRP and the Islamic Development Bank for their financial support. This work was facilitated by the “International joint Laboratory LMI PathoBios: Observatory of plant pathogens in West Africa: biodiversity and biosafety” (www.pathobios.com; twitter.com/PathoBios). We are very grateful to the fieldworkers in Burkina Faso notably Kader Guigma, Yacouba Kone, Compaoré Abraham Terry Cedric and Hébié Bakary. We thank the technical support structures (Bagrépole and Autorité de Mise en Valeur de la Vallée du Sourou) and the rice farmers from the study areas for their kind collaboration during sampling rice leaves and seeds. We also thank Pierre Mournet and Julien Frouin for their help in the development of the GBS technique at the UMR AGAP/CIRAD genotyping platform. Part of the experiments were supported by Agropolis Fondation (Labex Agro:ANR 10 LABX 0001 01) under RiPaBIOME project (ID 1702 011).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.1022348/full#supplementary-material

References

Ahmadpour, A., Castell-Miller, C., Javan-Nikkhah, M., Naghavi, M. R., Dehkaei, F. P., Leng, Y., et al. (2018). Population structure, genetic diversity, and sexual state of the rice brown spot pathogen Bipolaris oryzae from three Asian countries. Plant Pathol. 67, 181–192. doi: 10.1111/ppa.12714

CrossRef Full Text | Google Scholar

Archana, B., Kini, K. R., Prakash, H. S. (2014a). Genetic diversity and population structure among isolates of the brown spot fungus, bipolaris oryzae, as revealed by inter-simple sequence repeats (ISSR). Afr. J. Biotechnol. 13, 238–244. doi: 10.5897/AJB2013.12063

CrossRef Full Text | Google Scholar

Archana, B., Kini, K. R., Prakash, H. S. (2014b). Study of genomic fingerprints profile of Indian bipolaris oryzae from rice (Oryza sativa) by RAPD-PCR. Int. J. Life Sci. Biotechnol. Pharma Res. 3, 83–92.

Google Scholar

Aslam, H. M. U., Naveed, K., Hussain, S. I., Shakeel, Q., Ashraf, W., Anwaar, H. A., et al. (2021). First report of brown leaf spot of rice caused by Bipolaris zeicola in Pakistan. Plant Dis. 105, 212. doi: 10.1094/PDIS-04-20-0838-PDN

CrossRef Full Text | Google Scholar

Barnwal, M. K., Kotasthane, A., Magculia, N., Mukherjee, P. K., Savary, S., Sharma, A. K., et al. (2013). A review on crop losses, epidemiology and disease management of rice brown spot to identify research priorities and knowledge gaps. Eur. J. Plant Pathol. 136, 443–457. doi: 10.1007/s10658-013-0195-6

CrossRef Full Text | Google Scholar

Barro, M., Kassankogno, A. I., Wonni, I., Sereme, D., Somda, I., Kaboré, K. H., et al. (2021). Spatiotemporal survey of multiple rice diseases in irrigated areas compared to rainfed lowlands in the Western Burkina Faso. Plant Dis. 105 (12), 3889–3899. doi: 10.1094/PDIS-03-21-0579-RE

PubMed Abstract | CrossRef Full Text | Google Scholar

Boka, A., Bouet, A., Tiendrebeogo, A., Kassankogno, A. I., Ouedraogo, I., Nda, G. N. E., et al. (2018). Pathogenic variability of Bipolaris oryzae causing leaf spot disease of rice in West Africa. Int. J. Phytopathol. 7, 103–110. doi: 10.33687/phytopath.007.03.2643

CrossRef Full Text | Google Scholar

Burgos, M. R. G., Katimbang, M. L. B., Dela Paz, M. A. G., Beligan, G. A., Goodwin, P. H., Ona, I. P., et al. (2013). Genotypic variability and aggressiveness of Bipolaris oryzae in the Philippines. Eur. J. Plant Pathol. 137, 415–429. doi: 10.1007/s10658-013-0256-x

CrossRef Full Text | Google Scholar

Cardona, R., Gonzáles, M. S. (2007). First report of Exserohilum rostratum associated with rice seed in Venezuela. Plant Dis. 91, 226. doi: 10.1094/PDIS-91-2-0226C

PubMed Abstract | CrossRef Full Text | Google Scholar

Castell-Miller, C. V., Samac, D. A. (2012). Population genetic structure, gene flow and recombination of Cochliobolus miyabeanus on cultivated wildrice (Zizania palustris). Plant Pathol. 61, 903–914. doi: 10.1111/j.1365-3059.2011.02581.x

CrossRef Full Text | Google Scholar

Castell-Miller, C. V., Gutierrez-Gonzalez, J. J., Tu, Z., Bushley, K. E., Hainaut, M., Henrissat, B., et al (2016). Genome Assembly of the Fungus Cochliobolus miyabeanus, and Transcriptome Analysis during Early Stages of Infection on American Wildrice (Zizania palustris L.). PLoS ONE 11, 31. doi: 10.1371/journal.pone.0154122

CrossRef Full Text | Google Scholar

Chaijuckam, P., Songkumarn, P., Guerrero, J. (2019). Genetic diversity and aggressiveness of Bipolaris oryzae in north-central Thailand. Appl. Sci. Eng. Prog. 12, 116–125. doi: 10.14416/j.asep.2019.03.004

CrossRef Full Text | Google Scholar

Chakraborty, P., Chakraborty, A., Bhattacharya, S. G. (2020). Dispersal of airborne pathogenic conidia of Bipolaris oryzae inciting brown spot disease of paddy in West Bengal, India. J. Trop. Agric. 58, 199–208.

Google Scholar

Condon, B. J., Leng, Y., Wu, D., Bushley, K. E., Ohm, R. A., Otillar, R., et al. (2013). Comparative genome structure, secondary metabolite, and effector coding capacity across Cochliobolus pathogens. PloS Genet. 9, 29. doi: 10.1371/journal.pgen.1003233

CrossRef Full Text | Google Scholar

Dallagnol, L. J., Rodrigues, F. Á., Mielli, M. V. B., Ma, J. F., Datnoff, L. E. (2009). Defective active silicon uptake affects some components of rice resistance to brown spot. Phytopathology 99, 116–121. doi: 10.1094/PHYTO-99-1-0116

PubMed Abstract | CrossRef Full Text | Google Scholar

Danecek, P., Auton, A., Abecasis, G., Albers, C. A., Banks, E., DePristo, M. A., et al. (2011). The variant call format and VCFtools. Bioinformatics 27, 2156–2158. doi: 10.1093/bioinformatics/btr330

PubMed Abstract | CrossRef Full Text | Google Scholar

Davet, P. (1962). Observations phytosanitaires sur les riz cultivés en côte d’Ivoire et haute Volta pendant la campagne 1962 (Montpellier, France: ORSTOM).

Google Scholar

Dembélé, A. S. (2014) (Burkina Faso). Caractérisation morphologique et pathogénique de quelques isolats de Bipolaris oryzae (Breda de Haan), agent causal de l’helminthosporiose du riz et identification de sources de résistance génétique. [Mémoire d'ingéniorat].

Google Scholar

Diagne, A., Amovin-Assagba, E., Futakachi, K., Wopereis, M. C. S. (2013). “Farmers perceptions of the biophysical constraints to rice production in sub-Saharan Africa, and potential impact of research,” in Realizing africa’s rice promise (Wallingford, Oxfordshire, UK: CABI), 46–68.

Google Scholar

Doyle, J. J., Doyle, J. L. (1987). A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochem. Bull. 19, 11–15.

Google Scholar

Elshire, R. J., Glaubitz, J. C., Sun, Q., Poland, J. A., Kawamoto, K., Buckler, E. S., et al. (2011). A robust, simple genotyping-by-Sequencing (GBS) approach for high diversity species. PloS One 6, e19379. doi: 10.1371/journal.pone.0019379

PubMed Abstract | CrossRef Full Text | Google Scholar

Excoffier, L., Smouse, P. E., Quattro, J. M. (1992). Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics 131, 479–491. doi: 10.1093/genetics/131.2.479

PubMed Abstract | CrossRef Full Text | Google Scholar

Fan, Y., Ma, J., Gui, X., An, X., Sun, S., Dong, J. (2007). Distribution of mating types and genetic diversity induced by sexual recombination in Setosphaeria turcica in northern China. Front. Agric. China 1, 368–376. doi: 10.1007/s11703-007-0062-3

CrossRef Full Text | Google Scholar

Frichot, E., Mathieu, F., Trouillon, T., Bouchard, G., François, O. (2014). Fast and efficient estimation of individual ancestry coefficients. Genetics 196, 973–983. doi: 10.1534/genetics.113.160572

PubMed Abstract | CrossRef Full Text | Google Scholar

Goudet, J. (2005). Hierfstat, a package for r to compute and test hierarchical f-statistics. Mol. Ecol. Notes 5, 184–186. doi: 10.1111/j.1471-8286.2004.00828.x

CrossRef Full Text | Google Scholar

Grünwald, N. J., Goodwin, S. B., Milgroom, M. G., Fry, W. E. (2003). Analysis of genotypic diversity data for populations of microorganisms. Phytopathology® 93, 738–746. doi: 10.1094/PHYTO.2003.93.6.738

CrossRef Full Text | Google Scholar

Grünwald, N. J., Hoheisel, G.-A. (2006). Hierarchical analysis of diversity, selfing, and genetic differentiation in populations of the oomycete Aphanomyces euteiches. Phytopathology® 96, 1134–1141. doi: 10.1094/PHYTO-96-1134

CrossRef Full Text | Google Scholar

Hau, F. C., Rush, M. C. (1982). Preinfectional interaction between Helminthosporium oryzae and resistant and susceptible rice plants. Phytopathology 72, 285–292.

Google Scholar

Huson, D. H., Bryant, D. (2006). Application of phylogenetic networks in evolutionary studies. Mol. Biol. Evol. 23, 254–267. doi: 10.1093/molbev/msj030

PubMed Abstract | CrossRef Full Text | Google Scholar

Imrani, N., Touhami, A. O., Benkirane, R., Douira, A. (2019). Helminthosporium rostrata, foliar parasite of rice in Morocco. Plant Arch. 19, 2045–2018.

Google Scholar

Johnson, D. R., Percich, J. A. (1992). Wild rice domestication, fungal brown spot disease, and the future of commercial production in Minnesota. Plant Dis. 76, 1193–1198. doi: 10.1094/PD-76-1193

CrossRef Full Text | Google Scholar

Jombart, T. (2008). Adegenet: a r package for the multivariate analysis of genetic markers. Bioinformatics 24, 1403–1405. doi: 10.1093/bioinformatics/btn129

PubMed Abstract | CrossRef Full Text | Google Scholar

Jombart, T., Devillard, S., Balloux, F. (2010). Discriminant analysis of principal components: a new method for the analysis of genetically structured populations. BMC Genet. 11, 94. doi: 10.1186/1471-2156-11-94

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaboré, K. B. (2022). Helminthosporiose du riz en afrique de l’Ouest: identification des espèces responsables et diversité génétique de bipolaris oryzae et exserohilum rostratum. [Thèse de doctorat] (France: Université Paris-Saclay).

Google Scholar

Kaboré, K. H., Diagne, D., Milazzo, J., Adreit, H., Lebrun, M.-H., Tharreau, D. (2022). First report of rice brown spot caused by Exserohilum rostratum in Mali. Plant Dis. 106 (6), 1747. doi: 10.1094/PDIS-03-21-0662-PDN

CrossRef Full Text | Google Scholar

Kamal, M. M., Mia, M. A. T. (2009). Diversity and pathogenicity of the rice brown spot pathogen, Bipolaris oryzae (Breda de haan) shoem. in Bangladesh assessed by genetic fingerprint analysis. Bangladesh J. Bot. 38, 119–125. doi: 10.3329/bjb.v38i2.5135

CrossRef Full Text | Google Scholar

Kamvar, Z. N., Tabima, J. F., Grünwald, N. J. (2014). Poppr : an r package for genetic analysis of populations with clonal, partially clonal, and/or sexual reproduction. PeerJ 2, e281. doi: 10.7717/peerj.281

PubMed Abstract | CrossRef Full Text | Google Scholar

Korra, T., Navathe, S., Chand, R. (2022). Pathogenicity and infection behaviour of Exserohilum rostratum on wheat and associated collateral hosts. Research Square doi: 10.21203/rs.3.rs-2027087/v1

CrossRef Full Text | Google Scholar

Kusai, N. A., Azmi, M. M. Z., Zainudin, N. A. I. M., Yusof, M. T., Razak, A. A. (2016). Morphological and molecular characterization, sexual reproduction, and pathogenicity of Setosphaeria rostrata isolates from rice leaf spot. Mycologia 108, 905–914. doi: 10.3852/15-175

PubMed Abstract | CrossRef Full Text | Google Scholar

Ludwig, J. A., Reynolds, J. F. (1988). Statistical ecology: A primer on methods and computing (New York: Wiley-Interscience Pub).

Google Scholar

Ma, Q., Geng, Y., Li, Q., Cheng, C., Zang, R., Guo, Y., et al. (2022). Comparative mitochondrial genome analyses reveal conserved gene arrangement but massive expansion/contraction in two closely related Exserohilum pathogens. Comput. Struct. Biotechnol. J. 20, 1456–1469. doi: 10.1016/j.csbj.2022.03.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Majeed, R. A., Shahid, A. A., Liaqat, G. A., Saleem, K., Asif, M., Shafiq, M., et al. (2016). First report of setosphaeria rostrata causing brown leaf spot of rice in Pakistan. Plant Dis. 100 (10), 2162. doi: 10.1094/PDIS-02-16-0220-PDN

CrossRef Full Text | Google Scholar

Mathur, S. B., Kongsdal, O. (2003). Common laboratory seed health testing methods for detecting fungi, 1st ed. (Bassendorf, Switzerland: International Seed Testing Association).

Google Scholar

Milgroom, M. G. (1996). Recombination and the multilocus structure of fungal populations. Annu. Rev. Phytopathol. 34, 457–477. doi: 10.1146/annurev.phyto.34.1.457

PubMed Abstract | CrossRef Full Text | Google Scholar

Motlagh, M. R. S., Kaviani, B. (2008). Characterization of new Bipolaris spp.: The causal agent of rice brown spot disease in the north of Iran. Int. J. Agric. Biol. 10, 638–642.

Google Scholar

Nazari, S., Javan-Nikkhah, M., Fotouhifar, K.-B., Khosravi, V., Alizadeh, A. (2015). Bipolaris species associated with rice plant: pathogenicity and genetic diversity of Bipolaris oryzae using rep-PCR in mazandaran province of Iran. J. Crop Prot. 4, 497–508.

Google Scholar

Nigatu, G., Hansen, J., Childs, N., Seeley, R. (2017). Sub-Saharan Africa Is projected to be the leader in global rice imports. Amber Waves: The Economics of Food, Farming, Natural Resources, and Rural America 9. doi: 10.22004/ag.econ.266026

CrossRef Full Text | Google Scholar

Nwanze, K. F., Mohapatra, S., Kormawa, P., Keya, S., Bruce-Oliver, S. (2006). Rice development in sub-Saharan Africa. J. Sci. Food Agric. 86, 675–677. doi: 10.1002/jsfa.2415

CrossRef Full Text | Google Scholar

Ouédraogo, I. (2008). Incidence de l’helminthosporiose du riz au Burkina Faso et caractérisation des populations de l’agent pathogène [Bipolaris oryzae (Breda de haan) shoemaker]. [Thèse de doctorat] (Burkina Faso: Université de Ouagadougou).

Google Scholar

Ouédraogo, I., Correll, J. C., Boza, E. J., Cartwright, R. D., Lee, F. N., Sankara, P. (2004). Pathogenic, molecular, and genetic diversity among Bipolaris, Drechslera and Exserohilum species on rice. AAES Res. Ser. 529, 111–119.

Google Scholar

Ouédraogo, I., Wonni, I., Sérémé, D., Kaboré, K. B. (2016). Survey of fungal seed-borne diseases of rice in Burkina Faso. Int. J. Agric. Innov. Res. 5, 476–480.

Google Scholar

Picco, A. M., Rodolfi, M. (2002). Pyricularia grisea and Bipolaris oryzae: a preliminary study on the occurrence of airborne spores in a rice field. Aerobiologia 18, 163–167. doi: 10.1023/A:1020654319130

CrossRef Full Text | Google Scholar

Pielou, E. C. (1975). Ecological diversity (New York: John Wiley & Sons).

Google Scholar

Savary, S., Castilla, N. P., Elazegui, F. A., Teng, P. S. (2005). Multiple effects of two drivers of agricultural change, labour shortage and water scarcity, on rice pest profiles in tropical Asia. Field Crops Res. 91, 263–271. doi: 10.1016/j.fcr.2004.07.017

CrossRef Full Text | Google Scholar

Savary, S., Willocquet, L., Elazegui, F. A., Castilla, N. P., Teng, P. S. (2000). Rice pest constraints in tropical Asia: Quantification of yield losses due to rice pests in a range of production situations. Plant Dis. 84, 357–369. doi: 10.1094/PDIS.2000.84.3.357

PubMed Abstract | CrossRef Full Text | Google Scholar

Savary, S., Willocquet, L., Pethybridge, S. J., Esker, P., McRoberts, N., Nelson, A. (2019). The global burden of pathogens and pests on major food crops. Nat. Ecol. Evol. 3, 430–439. doi: 10.1038/s41559-018-0793-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Shannon, C. E. (1948). A mathematical theory of communication. Bell Syst. Tech. J. 27, 623–656. doi: 10.1002/j.1538-7305.1948.tb00917.x

CrossRef Full Text | Google Scholar

Silva, F. J. A., Maich, S. L. P., Meneses, P. R., Bellé, C., de Barros, D. R., Farias, C. R. J. (2016). First report on Exserohilum rostratum pathogenicity causing brown spot to rice in Brazil. Plant Dis. 100 (12), 2531. doi: 10.1094/PDIS-05-16-0709-PDN

CrossRef Full Text | Google Scholar

Simpson, E. H. (1949). Measurement of diversity. Nature 163, 688. doi: 10.1038/163688a0

CrossRef Full Text | Google Scholar

Stamatakis, A. (2014). RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30, 1312–1313. doi: 10.1093/bioinformatics/btu033

PubMed Abstract | CrossRef Full Text | Google Scholar

Sunder, S., Singh, R., Agarwal, R. (2014). Brown spot of rice: An overview. Indian Phytopathol. 67, 201–215.

Google Scholar

Toher, A. S., Mior Ahmad, Z. A., Wong, M. Y. (2016). First report of Exserohilum rostratum as pathogen of rice brown spot in Malaysia. Plant Dis. 100, 226–226. doi: 10.1094/PDIS-03-15-0276-PDN

CrossRef Full Text | Google Scholar

Tsuda, M., Ueyama, A. (1976). Distribution of two mating types of Cochliobolus miyabeanus in field and a laboratory attempt to produce hybrids between isolates of Helminthosporium oryzae from the USA and Japan. Ann. Phytopath. Soc. Jpn. 42, 7–11. doi: 10.3186/jjphytopath.42.7

CrossRef Full Text | Google Scholar

Tullis, E. C. (1935). Histological studies of rice leaves infected with helminthosporium oryzae. J. Agric. Res. 50, 81–90.

Google Scholar

Ueyama, A., Tsuda, M. (1976). Mating type and sexuality of Cochliobolus miyabeanus, the perfect state of Helminthosporium oryzae. Ann. Phytopathol. Soc Jpn. 42, 1–6. doi: 10.3186/jjphytopath.42.1

CrossRef Full Text | Google Scholar

van Oort, P. A. J. (2018). Mapping abiotic stresses for rice in Africa: Drought, cold, iron toxicity, salinity and sodicity. Field Crops Res. 219, 55–75. doi: 10.1016/j.fcr.2018.01.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Weir, B. S., Cockerham, C. (1984). Estimating f-statistics for the analysis of population structure. Evolution 38, 1358–1370. doi: 10.2307/2408641

PubMed Abstract | CrossRef Full Text | Google Scholar

Zadoks, J. C. (1974). The role of epidemiology in modern phytopathology. Phytopathology 64, 765–906.

Google Scholar

Zhang, C., Dong, S.-S., Xu, J.-Y., He, W.-M., Yang, T.-L. (2019). PopLDdecay: A fast and effective tool for linkage disequilibrium decay analysis based on variant call format files. Bioinformatics 35, 1786–1788. doi: 10.1093/bioinformatics/bty875

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: linkage disequilibrium, mating type, recombination, reproduction, rice disease, fungal pathogen

Citation: Kaboré KH, Kassankogno AI, Adreit H, Milazzo J, Guillou S, Blondin L, Chopin L, Ravel S, Charriat F, Barro M, Tollenaere C, Lebrun M-H and Tharreau D (2022) Genetic diversity and structure of Bipolaris oryzae and Exserohilum rostratum populations causing brown spot of rice in Burkina Faso based on genotyping-by-sequencing. Front. Plant Sci. 13:1022348. doi: 10.3389/fpls.2022.1022348

Received: 18 August 2022; Accepted: 07 November 2022;
Published: 25 November 2022.

Edited by:

Giampiero Valè, Università del Piemonte Orientale, Italy

Reviewed by:

Malkhan Singh Gurjar, Indian Agricultural Research Institute (ICAR), India
Maria Aragona, Council for Agricultural and Economics Research (CREA), Italy

Copyright © 2022 Kaboré, Kassankogno, Adreit, Milazzo, Guillou, Blondin, Chopin, Ravel, Charriat, Barro, Tollenaere, Lebrun and Tharreau. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Kouka Hilaire Kaboré, a291a2FzMkB5YWhvby5mcg==

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.