Efficiency of IRAP and ITS-RFLP marker systems in accessing genetic variation of Pyrenophora graminea

PLANT GENETICS

SHORT COMMUNICATION

Efficiency of IRAP and ITS-RFLP marker systems in accessing genetic variation of Pyrenophora graminea

Imad Zein; Mohammed Jawhar; Mohammed Imad Eddin Arabi

Department of Molecular Biology and Microbiology, AECS, Damascus, Syria

^{Send correspondence to}

ABSTRACT

The usefulness of IRAP (inter-retrotransposon amplified polymorphism) and ITS-RFLP (restriction of PCR-amplified internal transcribed spacers of the rDNA) markers in the analysis of 39 Pyrenophora graminea isolates was determined. Each marker system could discriminate between all of the isolates in detecting polymorphism, albeit with variable efficiency. IRAP and ITS-RFLP produced 85% and 77% polymorphic bands, respectively, with a corresponding mean polymorphic information content (PIC) of 0.38 and 0.36. The IRAP marker index ratio (2.41) was higher than ITS-RFLP (1.50). On one hand, the quality nature of data (QND) was higher for ITS-RFLP (0.169) than IRAP (0.093). However, correlation between both marker similarity matrices was significant (r = 0.34, p < 0.05). These findings suggest their combined use in phylogenetic analysis. To our knowledge, this is the first report of a comparison involving these two advanced DNA marker systems.

Key words:Pyrenophora graminea, IRAP, internal transcribed spacer, rDNA.

Pyrenophora graminea [anamorph Drechslera graminea (Rabenh. ex. Schltdl.) Ito], the seed-borne pathogen responsible for leaf-stripe in barley (Hordeum vulgare L.), inflicts heavy losses on crops (Porta-Puglia et al., 1986; Arabi et al., 2004). Several studies on morphological, physiological and biochemical aspects have already been undertaken (Zriba and Harrabi, 1995; Jawhar et al., 2000).

Traditionally, the classification of P. graminea isolates, besides requiring a certain expertise in taxonomy, may be further complicated by the inherent variation in morphological features among isolates, besides being time consuming, especially in those cases where similar species may be present in one and the same field (Gatti et al., 1992).

Over the years, the methods for detecting and assessing genetic diversity have extended from the analysis of discrete morphological traits to those of biochemical and molecular origin. Two classes of molecular markers which have received much attention in recent studies on genetic diversity in natural populations, are inter-retrotransposon amplified polymorphism (IRAP) (Kalendar et al., 1999; Pasquali et al., 2007), and restriction fragment length polymorphisms (RFLP) of PCR amplified internal transcribed spacer (ITS) regions (ITS-RFLP) (Hsiang and Wu 2000; Nilsson et al., 2008). The usefulness of these two markers types extends to resolve genetic variation among species within a genus or among populations (Redecker et al., 1997; Martin and Rygiewicz 2005; Branco et al., 2007).

Despite the general interest, it is not clear whether these two markers have comparable power for quantifying differentiation among populations. Thus, it would be of interest to determine whether IRAP and ITS-RFLP markers are equally efficient at detecting genetic patterns existent among P. graminea isolates. However, differences would suggest that one marker may be more appreciated for detecting isolation, which has implications for the use of either type of marker for defining demographically independent management units (Mortiz, 1994).

The present study aimed to evaluate the usefulness of both markers in assessing and analyzing the nature and extent of genetic diversity among isolates of P. graminea collected from various regions in Syria.

The 39 monosporic isolates of P. graminea used in the study were identified, cultivated, and maintained as described by Arabi et al. (2002, 2004). They were isolated from leaf-stripe infected barley leaves, originating from various regions in Syria, and selected from among 93 isolates, according to morphological and physiological criteria (virulence). The isolates were grown separately in 9 cm Petri dishes containing potato dextrose agar (PDA, DIFCO, Detroit, MI. USA), and incubated for 10 days at 21 ± 1ºCin the dark to facilitate mycelia growth.

Genomic DNA was extracted from fungal cultures as previously described (Arabi and Jawhar, 2007).

ITS regions and 5.8S rDNAs were amplified for all the isolates using the ITS1 (5' TCCGTAGGTGAACCTGCGG 3') and ITS4 (5'TCCTCCGCTTATTGATATGC 3') primers designed by White et al. (1990). The amplification protocol was as described by Arabi and Jawhar (2007). In separate reactions, 10 µL of PCR reaction were digested for 3 h with six different endonucleases (AluI, EcoR1, BsurI, BamHI, RsaI and HindIII), according to manufacturer's recommendations (MBI Fermentas, York, UK). DNA fragments were size-fractionated by electrophoresis through 1.5% agarose gels. The sizes were determined by comparison with their molecular weight relative to a DNA ladder (Q.BIOgene).

The IRAP method was used for retrotransposon amplification, as described by Kalendar et al. (1999). Primer sequences, as well as retrotransposon type and orientation are shown in Table 1. PCR was carried out using the method described by Jawhar and Arabi (2009). Amplified products were electrophoresed in a 2% agarose gel using a 1 x Tris-borate-EDTA buffer (100 mM Tris-HCl/L, pH 8.3, 83 boric acid/L, 1 mM EDTA/L) at 100 V. Subsequently, the gels were stained with ethidium bromide solution and visualized under ultraviolet illumination. The sizes of the amplified products were determined as mentioned above.

ITS-RFLP and IRAP banding profiles were scored for the presence (1) or absence (0) of bands. The experiments were repeated twice for each isolate and both markers, so as to confirm repeatability and remove monomorphic bands from the analysis. The data were converted to a Jaccards similarity (Jaccard, 1908) coefficient, which was used to construct a dendrogram by the unweighted pair-group method with arithmetic averages (UPGMA) utilising the software package Phylip 3.7 (Felsenstein, 1985). The polymorphism information content (PIC) was calculated for each locus according to Anderson et al. (1993), which provides an estimate of the discriminating power of a locus by taking into account the number of alleles generated by each reaction unit and their frequency distribution in the population. The percent of polymorphic markers (f) was estimated by dividing the number of polymorphic markers by the number of obtained markers. The multiplex ratio (MR) is defined as the number of bands per reaction unit, and the effective multiplex ratio EMR as the product of MR with the fraction of polymorphic markers. Marker utility (MI) for genetic studies was estimated as a marker index according to Powell et al., (1996) and Milbourne et al. (1997).

As PIC values are equal to gene diversity in binary marker systems, the effective number of alleles per marker was calculated as the respective reciprocal of the PIC value. The Mantel test (Mantel, 1967) was applied to ascertain the significance of correlations between pairwise genetic similarities in both marker systems. The probability of calculated correlation was estimated based on 1000 random permutations. These computations were carried out using the Arlequin software package (Excoffier et al., 2005). The quality nature of data (QND) of the marker system and the effective marker index EMI as an overall criterion for the utility of molecular markers were calculated according to Varshney et al. (2007).

Selected IRAP bands were cut with a surgical blade and purified with a QIAgene gel extraction kit according to manufacturers recommendations. Sequencing was carried out on a Genetic Analyzer (ABI 310, Perkin-Elmer, Applied Biosystems, USA). Each sequence was identified by homology search using the Basic Local Alignment Search Tool (BLAST) program (Altschul et al. 1997) against the GenBank nonredundant public sequence database.

PCR amplification with the specific primers ITS1 and ITS4 yielded single DNA fragments present in all isolates with ~ 650 bp in size, which is in agreement with the results obtained by a previous study (Arabi and Jawhar, 2007). Fingerprints generated from the five restriction digestions of the nrDNA ITS region denoted high levels of intraspecific variation within the P. graminea population. A total of 354 scorable DNA bands were scored, 274 of which (77%) being polymorphic, while the number of bands in isolates varied from 3 to 5.

Based on IRAP patterns, 534 bands were obtained, 454 (85%) of which were polymorphic, whereas the number of DNA bands in isolates varied between 4 and 15 (Figure 1). This is sustained by the findings of Taylor et al. (2004), who found a presence of high copy numbers of Pyggy-like sequences in the P. graminea genome by using a primer derived from the LTR (long terminal repeat) of the Pyggy retrotransposon isolated from this fungus. However, the sequence of one IRAP fragment, when using a Sukkula primer, showed similarity 19/27 (70.4%) to the LTR of the Pyggy retrotransposon (AF533703.1). Similarity began from position 34, as position 4 was a G instead of an A, and position 6 an A instead of a G. Furthermore, bases were missing at positions 13 and 24.

In addition, the P.graminea LTR sequence and a rice cDNA clone (Accession No. AK058381) were significantly similar, as attested by 117/147 bp identity (80%). Homology was also evident between P.graminea LTR and an Alternaria alternata LTR (AB025309), with 1556/2096 bp identity (74%). These results indicated the capability of LTR-specific primer to amplify in different target species.

On the other hand, both markers were highly repeatable, although QND was 0.191 for ITS-RFLP markers and only 0.141 for IRAP. PIC values were 0.376 and 0.355 for IRAP and ITS-RFLP, respectively. Furthermore, IRAP markers generated a substantially higher number of markers (7.80) and a superior marker index (2.41) than ITS-RFLP (Table 2).

The UPGMA dendrogram generated from IRAP and ITS-RFLP data demonstrated that isolates clustered into five groups for both markers by a similarity index of 0.341 for IRAP and 0.460 for ITS-RFLP (Figure 2). The correlation between IRAP and ITS-RFLP similarity matrices was moderate but significant (r = 0.34, p < 0.05).

The usefulness of a molecular marker technique depends upon both the polymorphic information content (PIC) of the markers and the number of markers generated by each primer (Varshney et al., 2007). Even though both IRAP and ITS-RFLP markers exhibited comparable PIC values, owing to the higher number of markers per assay, the MI of IRAP was 2.41 higher than for ITS-RFLP (Table 2). Bernardes et al. (2007) compared the performance of REMAP, a retrotransposon based maker technique, and ISSR with the fungus Magnaporthe grisea and has found an MI of 1.54 and 4.25, respectively. These marker systems are similar to IRAP and ITS-RFLP, in that they are PCRbased, anonymous and dominantly inherited. They also depend on repeated patterns in the genome to provide annealing sites for universal primers.

The results showed that band quality could have benefited from the additional restriction step following amplification in the ITS-RFLP protocol, thereby leading to a clearly defined band pattern. The documentation capabilities of band information produced by both marker assays in gene bank systems were comparable, as they were equally evaluated. This was placed on par with the widely used AFLP markers (Varshney et al., 2007), which are far inferior to single locus markers such as SSRs or SNPs with unique primers for each locus. On considering both quantitative and qualitative attributes, IRAP turned out to be superior to ITS-RFLP, as depicted by a more effective marker index (0.338 and 0.255, respectively).

Researchers have examined the existence of correlations between various molecular marker techniques in Fusarium oxysporum f. sp. lentis. Belabid et al. (2004) reported similar genetic relationships through RAPD and AFLP analysis. In barley, Russell et al. (1997) found that RFLP and AFLP, but not SSR, were correlated. In the present study, we found significant and moderate correlation between IRAP and ITS-RFLP in P. graminea pathogen according to the Mantel test, which confirmed in the partially conserved dendrogram topologies inferred from each of the similarity matrices (Figure 2). The moderate correspondence between these markers could possibly be attributed to different amplification targets in the P. graminea genome.

To our knowledge this is the first comparative report on the two advanced IRAP and ITS-RFLP genetic marker systems. The present study emphasized that, besides their effective employment, both of these DNA markers may furnish comparable results in assays of genetic differentiation among P. graminea isolates. Furthermore, due to the specific advantages of each marker, the combination of both marker systems can give us greater confidence that the delineated patterns are real, through drawing on results from multiple genetic systems (Allendorf and Seeb, 2000).

Acknowledgments

The authors wish to thank the Director General of AECS for his support and Dr. N. MirAli the Head of the Biotechnology Department for his critical review of the manuscript.

Received: July 30, 2009; Accepted: December 17, 2009.

Associate Editor: Everaldo Gonçalves de Barros

License information: 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 work is properly cited.

• Allendorf FW and Seeb LW (2000) Concordance of genetic divergence among sockeye salmon populations at allozyme, nuclear DNA and mitochondrial DNA. Evolution 54:640-651.
• Anderson JA, Churchill GA, Autrique JG, Tanksley SD and Sorells ME (1993) Optimizing parental selection for genetic linkage maps. Genome 36:181-186.
• Arabi MIE and Jawhar M (2007) Heterogeneity in Pyrenophora graminea as revealed by ITS-RFLP. J Plant Pathol 89:391-395.
• Arabi MIE, MirAli N, Jawhar M and Al-Safadi B (2002) Differentiation of Drechslera graminea isolates by cultural characters and SDS-PAGE. J Plant Pathol 84:153-156.
• Arabi MIE, Jawhar M, Al-Safadi B and MirAli N (2004) Yield response of barley to leaf stripe (Pyrenophora graminea)under experimental conditions in southern Syria. J Phytopathol 152:519-523.
• Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W and Lipman DJ (1997) Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Nucleic Acids Res 25:3389-3402.
• Belabid L, Baum M, Forts Z, Bouznad Z and Eujayl I (2004) Pathogenic and genetic characterization of Algerian isolates of Fusarium oxysporum f. sp. lentis by RAPD and AFLP analyses. Afr J Biotechnol 3:25-31.
• Bernardes A, Mistura CC, Carvalho FIF, Oliveira CA, Chadha S and Gopalakrishna T (2007) Comparative assessment of REMAP and ISSR marker assays for genetic polymorphism studies in Magnaporthe grisea Curr Sci 93:688-692.
• Branco CJS, Vieira EA, Gaspar-Malone G, Kopp MM and Malone E (2007) IRAP and REMAP assessments of genetic similarity in rice. J Appl Genet 48:107-113.
• Excoffier L, Laval G and Schneider S (2005) Arlequin 3.0: An integrated software package for population genetics data analysis. Evol Bioinfo Online 1:47-50.
• Felsenstein J (1985) Confidence limits on phylogenies: An approach using the bootstrap. Evolution 39:783-791.
• Gatti A, Rizza F, Delogu G, Terzi V, Porta-Puglia A and Vannacci G (1992) Physiological and biochemical variability in a population of Drechslera graminea J Genet Breed 46:179-186.
• Hsiang T and Wu C (2000) Genetic relationships of pathogenic Typhula species assessed by RAPD, ITS-RFLP and ITS sequencing. Mycol Res 104:16-22.
• Jaccard P (1908) Nouvelles recherches sur la distribution florale. Bull SocVaud Sci Nat 44:223-270.
• Jawhar M, Sangawn RS and Arabi MIE (2000) Identification of Drechslera graminea isolates by cultural characters and RAPD analysis. Cereal Res Comm 28:89-93.
• Jawhar M and Arabi MIE (2009) Heterogeneity in the internal transcribed spacers of the ribosomal DNA in Cochliobolus sativus as revealed by IRAP. J Plant Pathol 91:123-126.
• Kalendar R, Grob T, Regina M, Suoniem A and Schulman A (1999) IRAP and REMAP: Two new retrotransposon-based DNA fingerprinting techniques. Theor Appl Genet 98:704-711.
• Mantel N (1967) The detection of disease clustering and a generalized regression approach. Cancer Res 27:209-220.
• Martin KJ and Rygiewicz PT (2005) Fungal-specific PCR primers developed for analysis of the ITS region of environmental DNA extracts. BMC Microbiol 5:21-28.
• Milbourne D, Meyer R, Bradshaw JE, Baird E, Bonar N, Provan J, Powell W and Waugh R (1997) Comparison of PCR based marker systems for the analysis of genetic relationships in cultivated potato. Mol Breed 3:127-136.
• Mortiz C (1994) Defining "evolutionary significant units" for conservation. Trends Ecol Evol 9:373-376.
• Nilsson RH, Kristiansson E, Ryberg M, Hallenberg N and Larsson KH (2008) Intraspecific ITS variability in the kingdom fungi as expressed in the international sequence databases and its implications for molecular species identification. Evol Bioinfo 4:193-201.
• Pasquali M, Dematheis F, Gullino ML and Garibaldi A (2007) Identification of race 1 of Fusarium oxysporum f. sp. lactucae on lettuce by inter-retrotransposon sequencecharacterized amplified region technique. Phytopathology 97:987-996.
• Porta-Puglia A, Delogu G and Vannacci G (1986) Pyrenophora graminea on winter barley seed: Effect on disease incidence and yield losses. Phytopathology 117:26-33.
• Powell W, Morgante M, Andre C, Hanafey M, Vogel J, Tingey SV and Rafalski A (1996) The utility of RFLP, RAPD, AFLP and SSRP (microsatellite) markers for germplasm analysis. Mol Breed 2:225-238.
• Redecker D, Thierfelder H, Walker C and Werner D (1997) Restriction analysis of PCR-amplified internal transcribed spacers of ribosomal DNA as a tool for species identification in different genera of the order Glomales. Appl Envir Microbiol 63:1756-1761.
• Russell RJ, Fuller JD, Macaulay M, Hats BG, Jahoor A, Powell W and Waugh R (1997) Direct comparison of levels of genetic variation among barley accessions detected by RFLPs, AFLPs, SSRs and RAPDs. Theor Appl Genet 95:714-722.
• Taylor EJA, Konstantinova P, Leigh F, Bates JA and Lee D (2004) Gypsy-like retrotransposons in Pyrenophera species: An abundant and informative class of molecular markers. Genome 47:519-525.
• Varshney RK, Chabane K, Hendre PS, Aggarwal RK and Graner A (2007) Comparative assessment of EST-SSR, EST-SNP and AFLP markers for evaluation of genetic diversity and conservation of genetic resources using wild, cultivated and elite barleys. Plant Sci 173:638-649.
• White TJ, Bruns T, Lee S and Taylor J (1990) Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In: Innis MA, Gelfand DH, Sninsky JJ and White TJ (eds), PCR Protocol, a Guide to Methods and Applications. Academic Press, San Diego, pp 315-322.
• Zriba W and Harrabi M (1995) Cultural and pathogenic variability in Pyrenophora graminea isolates. Rachis 14:99.

Publication Dates

• Publication in this collection
  21 May 2010
• Date of issue
  2010

History

• Received
  30 July 2009
• Accepted
  17 Dec 2009