Efficient Discovery of Single-Nucleotide Variations in Cochliobol
Journal of Plant Biochemistry & Physiology

Journal of Plant Biochemistry & Physiology
Open Access

ISSN: 2329-9029

+44 7868 792050

Research Article - (2018) Volume 6, Issue 2

Efficient Discovery of Single-Nucleotide Variations in Cochliobolus sativus Vegetative Compatibility Groups by Ecotilling

Jawhar M1*, Arabi MIE1, MirAli N1 and Till BJ2
1Department of Molecular Biology and Biotechnology, AECS, Damascus, Syria
2Plant Breeding and Genetics Laboratory Vienna, Austria
*Corresponding Author: Jawhar M, Department of Molecular Biology and Biotechnology, AECS, Damascus, Syria, Tel: +2138520 Email:


Cochliobolus sativus, the causal agent of barley spot blotch disease, is able to undergo spontaneous, vegetative cell fusion. In this study, the utility of enzymatic mismatch cleavage for discovery of single-nucleotide variations in Syrian Vegetative Compatibility Groups (VCGs) obtained from pairing complementary nit mutants of C. sativus isolates was investigated for the first time. Gene-specific primers were designed from the whole fungal genome for use in the Ecotilling assays on VCGs isolates. Differences of band patterns among the digested products of different isolates were clearly observed by agarose gel electrophoresis. Interestingly, isolates belonging to the same VCG had some different sets of single nucleotide polymorphisms. Twenty-six putative nucleotide polymorphisms were identified among the C. sativus VCGs. Therefore, the results suggested that this marker type may be informative for confirmation the conventionally determined VCGs by increasing the reliability of heterokaryons group classification.

Keywords: Ecotilling; Nit mutants; Cochliobolus sativus; Barley; Vegetative compatibility group; SNPs


Spot blotch (SB) of barley (Hordeum vulgare L.) caused by the fungus Cochliobolus sativus (Ito and kuribayashi) Drechs. Ex Dastur (anamorph: Bipolaris sorokiniana (Sacc.) Shoem), is a devastating disease that occurs in major barley-growing regions of the world [1,2]. Studies during the last decade have demonstrated that genetic diversity within C. sativus populations is higher than previously thought, and that limited understanding of the nature of such diversity has hampered the progress of understanding SB development and the efficacy of its management [3,4].

Vegetative compatibility groups (VCGs) assay has been used to measure population diversity of fungi and to enable appropriate identification and characterization of individual isolates. They demonstrate the ability of fungal isolates to undergo heterokaryon formation that reflects phenotypic differences or similarity among individuals representing the population of a species [5].

Genetic studies in C. sativus populations yielded contradictory results but, overall, indicated a lack of clear genetic distinction between its VCGs. However, since this fungus is a haploid, strictly asexually reproducing fungus, hyphal fusion is a prerequisite to genetic exchange among different isolates of the pathogen [6,7]. Therefore, C. sativus isolates from different VCGs are thought to be genetically distinct populations which may vary in a number of ecological, physiological, and virulence traits [7,8].

In previous work [9], we found that Syrian C. sativus population grouped into three major VCGs based on complementation tests with nitrate- nonutilizing (nit) mutants, which are impaired in their ability to utilize nitrate from the growing media [10,11]. However, a little is known about their genetic diversity, and this lack of a clear molecular distinction among these VCGs suggests that some gene flow might occur among them. This poses a risk that new virulence phenotypes may develop in the pathogen population.

Ecotilling, a variant of Tilling (Targeting Induced Local Lesions IN Genomes) method, allows high-throughput analyses of natural genetic diversity in populations [12-14]. In this technique, a single-strand specific endonuclease extracted from celery (CEL I) is used to cleave DNA where there are mismatched bases in the DNA heteroduplexes of PCR products. Cleaved fragments observed by gel electrophoresis indicate the presence of nucleotide variations. Because cleavage reveals only heterozygous polymorphisms, a sample pooling strategy involving at least one reference DNA and one test DNA is commonly employed as a means of uncovering homozygous differences between two samples [15].

The present study was undertaken to discover, validate and to genotype the single nucleotide polymorphisms (SNPs) in C. sativus VCGs using Ecotilling assay for better characterize molecular variability among them.

Materials and Methods

Fungal isolates

Members of six C. sativus VCGs isolates were collected from barley leaves grown in four regions of Syria that are separated by over 350 km (Table 1). VCGs characterization for these isolates has been described previously [9]. Five C. sativus isolates (Cs1,11,15,16,22) belonging to the nit1/nit3 (colonies that grew with a thin mycelium on both CDA and as a wild type phyenotype on CDA amended with sodium nitrate) whereas, Cs23 characterized as Nit M (colonies that grew with a thin mycelium on both CDA and CDA media amended with sodium nitrate) according to Puhalla [10].

Primer name Forward primer sequence (5'T3') Reverse primer sequence (5'T3')

Table 1: Primer sequences used for PCR amplification.

DNA extraction

DNA was extracted according to standard protocols [16], resuspended in TE buffer (10 mM Tris–HCl, pH 8.0; 1mM EDTA) and stored at -20°C.

Ecotilling protocol

Primers were designed from whole genome sequence information from JGI (htpp://genomejgi.doe.gove) for C. sativus [17,18] using parameters previously described for Tilling [19,20]. PCR primers specific for Carp1 and GlutSynth genes were designed based on the genome sequence of these genes in C. sativus ND90 (Table 1). PCR amplification with genomic DNA at varying concentrations, primer pair mismatches cleavage using celery juice extract, gel electrophoresis and image analysis followed previously described protocols [19]. The PCR reaction was performed in a 20 μl reaction volume containing 1 U AmpliTaq® DNA polymerase, 60 μM dNTPs, 0.75 ng of genomic DNA, and 1.0 μM primers in 1X PCR Buffer (1.5 mM MgCl2, 50 mM KCl, 10 mM Tris HCl, pH 8.3 at RT). After an initial denaturing step at 94°C for 4 min, the DNA was amplified through 20 cycles consisting of denaturing at 94°C for 5 s, annealing at 65°C for 1 min decreasing by 0.5°C per cycle, and extension at 72°C for 1 min. The samples were then subjected to an additional 30 cycles consisting of denaturing at 94°C for 5 s, annealing at 55°C for 1 min, and extension at 72°C for 1 min, with a final extension of 5 min at 72°C. Denaturing and reannealing step is included at the end of the PCR reaction (99°C for 10 min, followed by 70 cycles of 70°C for 20s decreasing by 0.3°C per cycle) to allow the formation of heteroduplexes if nucleotide polymorphisms are present in the pool.

CELI digestion of the heteroduplexes and agarose gel electrophoresis detection

Celery juice extract (CEL I enzyme) was extracted from celery following the method described by Till et al. [21]. Heteroduplexes (20 μl) formed after PCR amplification, denaturation and annealing, were mixed with 20 μl nuclease mixture, containing 6 μl CJE Buffer (10 mM Hepes, pH 7.5, 10 mM MgSO4, 10 mM KCl, 0.002% Triton X-100, and 0.2 μg/ml BSA); 2.4 μl crude celery juice extract (CJE) preparation [19], and 11.6 μl ultrapure water. The digestion was then performed at 45°C for 45 min followed by incubation on ice to stop the reaction. Subsequently, 4 μl of 10 × gel loading buffer was added to each sample. The digested products were separated on 2% agarose gels with 1 × TAE buffer at 110 V for 50 min, stained in an ethidium bromide buffer (0.5 μg/ml) for 15 min, and then visualized by Gel Doc XR (BioRad Laboratories, Inc.). The experiments were repeated twice and resulted similar results.

Sample analysis

The purified products were sequenced with ABI 310 (Perklin-elmer, Applied Biosystems, USA) using the same primers used for PCR amplification. Geneious Pro Software v 8.0 (Biomatters Ltd; [22] was used for detecting SNPs between VCGs.

Results and Discussion

In this work, a VCGs enzymatic mismatch cleavage strategy was developed for Ecotilling in haploid C. sativus fungus using selfextracted nuclease and agarose gel electrophoresis (Figure 1).


Figure 1: Scheme for Ecotilling protocol for analysis of VCG C. sativus single spore culture.

The assay was able to successfully distinguish C. sativus isolates representing three major VCGs previously described by Arabi and Jawhar [9] using nit complementation tests. SNPs were detected among isolates belonging to the same VCG typically, and twenty-six putative nucleotide polymorphisms were identified among the C. sativus VCGs (Table 2).

Isolate no. VCG Region Mutant phenotype Gene Location Nucleotide change Primer Amino acid change
Cs22 A Hasake nit1/nit3 Carp1 1197 C/G Forward -
        " 1224 T/A Forward -
        Glutsynth 1383 G/A Forward -
        " 1407 A/C Forward -
Cs16 A ICARDAA nit1/nit3 Carp1 561 T/G Reverse -
        " 1172 T/G Forward -
        Glutsynth 1353 G/A Forward -
        " 1389 A/G Forward -
Cs11 B Daraa nit1/nit3 Carp1 96 C/T Forward -
        " 135 C/T Forward -
        Glutsynth 323 T/C Forward -
        " 327 G/C Forward -
Cs23 B Allepo nitM Carp1 374 G/C Forward -
        " 395 C/T Reverse I/F
        " 1188 C/A Forward -
        Glutsynth 312 A/G Forward  
        " 1250 A/G Forward -
Cs1 C Daraa nit1/nit3 Carp1 451 C/T Reverse L/F
        " 479 C/T Reverse S/F
        " 485 C/T Reverse A/F
        Glutsynth 650 C/G Reverse -
Cs15 C   nit1/nit3 Carp1 488 C/T Reverse L/F
        " 1022 A/T Reverse -
        " 1149 A/C Forward -
        Glutsynth 660 G/T Reverse V/F
        " 1308 G/C Reverse -

Table 2: Vegetative compatibility (VCG A,B and C), nit mutants, and Ecotilling polymorphisms between the C. sativus isolates and refernce isolate. Mutant phyenotype generated as described by Puhalla [11].

Differences of band patterns among the digested products of different isolates were clearly observed by agarose gel electrophoresis. In our experiments, since C. staivus is a haploid organism, individual isolates did not formed heteroduplexed DNA at the site of a mismatch, therefore, for the success of the method based on using pooling strategy for each VCG consisted of two isolates [13], and each VCG produced unique Ecotilling fingerprints (Figure 2).


Figure 2: Ecotilling in C. sativus VCGs A, B and C using Glutsynth primers. M: molecular weight marker (HinfI; MBI Fermentas, York, UK).

The differences in SNPs detected for isolates within a VCG could be the result of genetic changes in clonal lineages. Leslie [5] noted that isolates from the same VCG are not always members of the same clone. Also, the existence of bridging isolates able to form stable heterokaryons between isolates in different VCGs would eventually facilitate some gene flow between them [5]. Therefore, the evolutionary potential within C. sativus may be intimately related to VCGs [3]. This is relevant to disease control in terms of potential for emergence of new virulence or pathogenicity traits in natural populations of C. sativus . Additionally, by comparing with our previous effort to characterize C. sativus populations collected from different locations in Syria [3], using RAPD and ITS-RFLP markers, the present results suggest that isolates assigned to the same fingerprint group using these DNA markers would also belong to the same VCG.

On the other hand, detecting polymorphisms among C. sativus VCGs are not surprising given the scale of this study and what is known about their different origin over several years [23]. The VCGs isolates tested here were first selected based on morphological and physiological criteria [3,24], and so a high degree of variation may be expected. On the other hand, studying the different VCGs of C. sativus confirms that this is an attractive model organism for both phenotype and nucleotide sequence. Nucleotide diversity [25,26] and diversity at selected genes in C. sativus [27] were reported.


This work demonstrated the utility of enzymatic mismatch cleavage for Ecotilling in VCGs of the haploid organism. The assay was able to successfully distinguish C. sativus isolates representing three Syrian major VCGs. Interestingly, isolates belonging to the same VCG had some different sets of single nucleotide polymorphisms. The data provide an independent confirmation of the genetic structure and diversity of C. sativus populations where all the isolates were initially grouped by VCG analysis [3]. It is now possible to more accurately predict that C. sativus isolates belonging to the same Ecotilling profile are also vegetatively compatible. However, even though this study was carried with a limited number of isolates, these results suggest that the conventional VCG assay together with Ecotilling can give a reliable measure for studying heterokaryons formation in fungi. Moreover, this method compares favorably to the most common alternative–full sequencing, since Ecotilling greatly reduces the number of sequences that needs to be determined and is several fold cheap than full sequencing [19].


The authors would like to thank the Director General of AECS and the Head of molecular biology and biotechnology department for their support. They would also like to thank the Food and Agriculture Organization of the United Nations and the International Atomic Energy Agency through their Joint FAO/IAEA programme of nuclear techniques in food and agriculture for supporting the laboratory experiments. Thanks also extended to Mr. E. Al-Shehadah for his great help for achieving the Geneious analysis.

Conflict of Interest

The authors declare that the study involved no conflict of interest from any party.


  1. Mathre DE (1997) Compendium of barley diseases. American Phytopathological Society Press, St. Paul, MN, USA.
  2. Sharma RC, Duveiller E (2006) Spot blotch continues to cause substantial grain yield reductions under resource limited farming conditions. J Phytopathol 154:  482-488.
  3. Arabi MIE, Jawhar M (2007) Molecular and pathogenic variation identified among isolates of Cochliobolus sativus. Austr Plant Pathology 36: 17-21.
  4. Leisova-Svobodova L, Minarikova V, Kucera L, Pereyra SA (2012) Structure of the Cochliobolus sativus population variability. Plant Pathol 61: 709-718.
  5. Leslie FJ (1993) Fungal vegetative compatibility. Annu Rev Phytopathol 31: 127-150.
  6. Tinline RD, Dixon JG (1958) Cochliobolus sativus I. Perithecial development and the inheritance of spore colour and mating type.  Mycologia 50: 697-706.
  7. Kumar J, Schafer P, Huckelhoven R, Langen G, Baltruschat H, et al. (2002) Bipolarissorokiniana, a cereal pathogen of global concern:  cytological and molecular approaches towards better control. Mol Plant Pathol 3: 185-195.
  8. Pandey SP, Sharma S, Chand R, Shahi P, Joshi AK (2008) Clonal variability and its relevance in generation of new pathotypes in the spot blotch pathogen, Bipolarissorokiniana. CurrMicrob 56: 33-41.
  9. Arabi MIE, Jawhar M (2007b) Vegetative compatibility and pathogenicity in Cochliobolus sativus. Austr Plant Pathol 36: 173-174.
  10. Puhalla JE (1985) Classification of strains of Fusarium oxysporum on the basis of vegetative compatibility. Can J Bot 63: 179-183.
  11. Strausbaugh CA, Eujayl IA, Martin FN (2016) Pathogenicity, vegetative compatibility and genetic diversity of Verticillium dahliae isolates from sugar beet. Can J Plant Pathol 38: 492- 505.
  12. Bajaj D, Srivastava R,  Nath M,  Tripathi S, Bharadwaj C, et al. (2016) Ecotilling based association mapping efficiently delineates functionally relevant natural allelic variants of candidate genes governing agronomic traits in chickpea. Fron Plant Sci 7: 450
  13. Jawhar M, Till BJ, Albaterni A, Skiheita A, Arabi MIE, et al. (2017) Optimizing Ecotilling technique for tracking field Cochliobolus sativus population. J Plant Pathol 99: 339-345.
  14. Jawhar M, Till BJ, Arabi MIE, MirAli N (2017) An efficient throughput protocol for validation of pathotyping by enzymatic mismatch cleavage. J Plant BiochemPhysiol 5: 4.
  15. Hofinger JB, Huynh OA, Jankowicz-Cieslak J, Muller A, Otto I, et al. (2013) Validation of doubled haploid plants by enzymatic mismatch cleavage. Plant Methods 9: 43.
  16. Leach J, Finkelstein DB, Rambosek JA (1986) Rapid miniprep of DNA from filamentous fungi. Fungal Gen Newsl 33: 32-33.
  17. Ohm RA, Feau N, Henrissat B, Schoch CL, Horwitz BA, et al. (2012) Diverse lifestyles and strategies of plant pathogenesis encoded in the genomes of eighteen Dothideomycetes fungi. PLoS Pathogens 8: 1003037.
  18. Condon BJ, Leng Y, Wu D, Bushley KE, Ohm RA, et al. (2013) Comparative genome structure, secondary metabolite, and effector coding capacity across Cochliobolus pathogens. PLoS Gen 9
  19. Till BJ, Zerr T, Comai L, HenikoffS (2006) A protocol for Tilling and Ecotilling in plants and animals. Nature Prot 1: 2465-2477.
  20. Till BJ, Reynolds SH, Greene EA, Codomo CA, Enns LC,et al. (2003) Large-scale discovery of induced point mutations with high-throughput TILLING. Genome Res 13: 524-530.
  21. Till BJ, Burtner C, Comai L, Henikoff S (2004) Mismatch cleavage by single strand specific nucleases. Nucl Acids Res 32: 2632-2641.
  22. Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, et al. (2012) Geneious Basic:  an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 28: 1647-1649.
  23. Arabi MIE, Jawhar M (2004) Identification of Cochliobolus sativus (spot blotch) isolates expressing differential virulence on barley genotypes in Syria. J Phytopathol 152: 461-464
  24. Arabi MIE, Jawhar M (2003) Pathotypes of Cochliobolus sativus (spot blotch) on barley in Syria. J Plant Pathol 85: 193-196.
  25. Zhong S, Steffenson BJ (2001) Virulence and molecular diversity in Cochliobolus sativus. Phytopathology 91: 469-476.
  26. Gurung S, Mahto BN, Gyawali S, Adhikari TB (2013) Phenotypic and molecular diversity of Cochliobolus sativus populations from wheat. Plant Dis 97: 62-73.
  27. Yueqiang L (2011) Functional genomics of Cochliobolus sativus. Phd. Thesis. North Dakota State University, USA, pp: 1-24.
Citation: Jawhar M, Arabi MIE, MirAli N, Till BJ (2018) Efficient Discovery of Single-Nucleotide Variations in Cochliobolus sativus Vegetative Compatibility Groups by Ecotilling. J Plant Biochem Physiol 6: 211.

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