Prevalence of Tobacco mosaic virus in Iran and Evolutionary Analyses of the Coat Protein Gene

Article information

Plant Pathol J. 2013;29(3):260-273
1Department of Plant Pathology, College of Agriculture and Natural Resources, Science and Research Branch, Islamic Azad University, Tehran 14515-775, Iran
2Department of Horticultural Sciences, Seoul Women’s University, Seoul 139-774, Korea
*Corresponding author. Phone) +989124786687, FAX) +982144865464, E-mail) rakhshandehroo_fa@srbiau.ac.ir
Received 2012 September 13; Revised 2013 February 7; Accepted 2013 March 5.

Abstract

The incidence and distribution of Tobacco mosaic virus (TMV) and related tobamoviruses was determined using an enzyme-linked immunosorbent assay on 1,926 symptomatic horticultural crops and 107 asymptomatic weed samples collected from 78 highly infected fields in the major horticultural crop-producing areas in 17 provinces throughout Iran. The results were confirmed by host range studies and reverse transcription-polymerase chain reaction. The overall incidence of infection by these viruses in symptomatic plants was 11.3%. The coat protein (CP) gene sequences of a number of isolates were determined and disclosed to be a high identity (up to 100%) among the Iranian isolates. Phylogenetic analysis of all known TMV CP genes showed three clades on the basis of nucleotide sequences with all Iranian isolates distinctly clustered in clade II. Analysis using the complete CP amino acid sequence showed one clade with two subgroups, IA and IB, with Iranian isolates in both subgroups. The nucleotide diversity within each sub-group was very low, but higher between the two clades. No correlation was found between genetic distance and geographical origin or host species of isolation. Statistical analyses suggested a negative selection and demonstrated the occurrence of gene flow from the isolates in other clades to the Iranian population.

Iran, with an area of 1.65 million square kilometres, contains several different types of climate. Horticultural crop farming is carried out in most of the geographical regions in Iran, except for the deserts and mountains. The climate of Iran is conducive to the growing of many horticultural crops such as broad beans (Vicia faba L.), carrots (Daucus carota subsp. sativus), varieties of cucurbit plants [including cantaloupe (Cucumis melo var. cantalupensis), courgettes (Cucurbita pepo L. cv. Zucchini), cucumber (Cucumis sativus L.), melon (C. melo L.), pumpkin (Cucurbita moschata Duch.), watermelon (Citrullus lanatus var. lanatus) and winter squash (Cucurbita maxima Duch. ex Lam)], eggplant (Solanum melongena L.), garlic (Allium sativum L.), green beans (Phaseolus vulgaris L.), onion (Allium cepa L.), pepper (Capsicum annum L.), potato (Solanum tuberosum subsp. tuberosum), spinach (Spinacia oleraceae L.) and tomato (Solanum lycopersicum L.). Such crop production in the field may be either rain-fed or under irrigation on large scale or small holder commercial farms. Iran is among the main producers of horticulural crops in the world with an approximate production of 1,800,400 tonnes yearly (Faostat, 2009). Cucurbits and tomato are major horticultural crops in Iran (Faostat, 2009). Viral diseases of horticultural crops cause important economic impact worldwide. Among all other viral diseases of such crops, diseases caused by tobamoviruses are often the most destructive and difficult to control (Alonso et al., 1989).

On the basis of genome organization and phylogenetic clustering, tobamoviruses are classified into three sub-groups (Lartey et al., 1996). Tobamoviruses in subgroup 1, including Tobacco mosaic virus (TMV), Tomato mosaic virus (ToMV), Pepper mild mottle virus (PMMV) and Tobacco mild green mosaic virus (TMGMV), are known to infect mostly solanaceous plants. Several tobamoviruses that previously were considered as TMV strains, are now classified as members of new species of the genus Tobamovirus on the basis of differences in host range, serological properties and amino acid sequences of encoded virus proteins (Antignus et al., 2001; Wang et al., 1997). TMV is the type member of the genus Tobamovirus. The TMV virion is a rigid rod (18 nm × 300 nm) containing a single-stranded, positive-sense RNA (Zaitlin, 1999). TMV can be transmitted by mechanical inoculation, grafting, contact between plants and by seed, but not by any known vector (Broadbent, 1965; Hollings and Huttinga, 1976). Tobamoviruses collectively have a very wide host range and cause serious economic impact and significantly yield losses in many crops such as brassicas, cucurbits, solanaceous crops and different ornamental plants such as chrysanthemums (Chrysanthemum indicum L.), impatiens (Impatiens balsamina) and petunia (Petunia × hybrida) (Alexandre et al., 2000; Choi et al., 2009; Kumar et al., 2011; Nassar et al., 2012). Infected plants often show different types of symptoms such as malformations, mosaic and stunting on leaves, flowers, and fruit. TMV has been reported to cause considerable reduction in the tonnage of tomato fruit up to 59% (Cherian and Muniyapppa, 1998). Significant yield losses of about 90% due to TMV infection in pepper also have been reported (Chitra et al., 2002). There are many different strains identified for TMV, some of which are capable of overcoming of their hosts resistance (Padgett and Beachy, 1993). Regardless of the importance of field-grown horticultural crops, there is little information available on the incidence and distribution of TMV and its genetic diversity in cultivated crops in Iran. Thus, in this study we determined the incidence of TMV in plants showing symptoms in commercial fields in the main important provinces for horticultural crop cultivation in Iran, and also determined the genetic diversity of the isolates detected in the surveyed regions. We also detected TMV in some weed species growing in and around crop fields which may be potential reservoir hosts for maintaining TMV in these fields. The information obtained is required as the first step toward the search for control strategies of virus diseases in horticulural crops in Iran.

Materials and Methods

Surveys and sample collection.

Samples of horticultural crop species including bean, cabbage (Brassica oleracea var. capitata), cucumber, eggplant, pepper, potato, pumpkin, radish (Raphanus sativus), tomato, watermelon, and zucchini. were collected from April 2009 to May 2010 from open fields in the following 17 provinces (and specified districts) of Iran: Alburz (Karaj), Boushehr (Boushehr), East Azarbaijan (Maragheh), Fras (Darab, Fasa, Jahrum and Larestan), Golestan (Minoodasht), Hamedan (Tuyserkan), Hormozgan (Hajiabad, Kehorestan, Minab, Roudan), Isfahan (Isfahan), Kerman (Jiroft), Kermanshah (Mahidasht), Khuzestan (Dezful), Mazandaran (Noshahr), Qazvin (Qazvin), Qom (Qom), Sistan and Baluchestan (Zabul and Zahedan), Tehran (Tehran, Varamin, Shahriar), and West Azarbaijan (Urmia) (Table 1). Overall, 78 fields were surveyed and 1926 samples (20–30 samples per field) were collected during the growing season. The fields were randomly selected using a predetermined distance criterion, where distance between the fields ranged from 10 to 30 km. A sample consisted of two symptomatic leaves per plant. Viral-like symptoms observed included dwarfing, mosaic, mottling, yellowing, leaf malformation, necrosis of the leaves and stems and vein clearing. In addition, 107 leaf samples were also collected from different asymptomatic weed species belonging to 12 different botanical families within visited horticultural crop fields in the surveyed provinces for the serological detection of TMV and related tobamoviruses (Table 2). Weed samples were collected from around the bases of crops showing symptoms characteristic of infection by TMV.

Incidence of tobamoviruses on horticultural crops collected from the main provinces of cultivation in Iran

Incidence of tobamovirus infection on selected weed species

Virus testing.

Mature leaf blades were tested in duplicate using double antibody sandwich enzyme-linked immunosorbent assay (DAS–ELISA), with coating antibodies against TMV and related tobamoviruses, as well as alkaline phosphatase conjugates obtained from Agdia (Agdia, Inc., South Bend, IN, USA) (Clark and Adams, 1977).

Host range.

Leaf samples showing symptoms associated with virus disease and with positive reaction in ELISA that were collected from Alburz, East Azarbaijan, Fars, Golestan and Tehran provinces, which contain the main horticultural crop cultivation regions of Iran, were selected and further evaluated in a host range study. The virus isolates were purified biologically through a single local lesion technique repeated twice on Chenopodium quinoa plants, and then transmitted mechanically to Nicotiana tabacum cv. White Burley for virus propagation according to the published protocol (Adkins et al., 2003). Extracts from the systemically-infected tobacco plants were then used as a source of biologically-purified virus and inoculated mechanically to C. quinoa (family Chenopodiaceae), C. sativus (cucumber; family Cucurbitaceae), Vigna unguiculata cv. Mashhad (cowpea) and Vicia faba (faba bean; family Fabaceae), N. rustica, N. tabacum cv Samsun nn, N. tabacum cv. Virginia,N. tabacum cv. White Burley and S. lycopersicum (tomato), Petunia hybrida (all in the family Solanaceae). (Table 3). Specifically, 1 g of leaf tissue was ground in 2 ml of 0.1 M potassium phosphate buffer, pH 7.5. The extract was rubbed onto leaves of the above herbaceous plants, previously dusted with 300-mesh Carborundum. The experiment was conducted twice using at least three plants for each inoculation. The inoculated plants were then transferred to a greenhouse and maintained at 20 to 30 °C with 12 to 14 h natural lighting for the duration of the test. The final symptom reading was conducted 4 weeks post inoculation, and systemic leaf tissues from each individual plant were collected to confirm the presence of virus by ELISA. The presence of the virus in asymptomatic plants was verified by reverse transcription-polymerase chain reaction (RT-PCR).

Symptom development on indicator plant species after inoculation with Tobacco mosaic virus (TMV) isolates detected in Iran

RT-PCR and sequencing analyses.

Total RNA extraction was performed using LiCl buffer according to the protocol described by Channuntapipat et al. (2001). Total RNAs were isolated from leaves of symptomatic crop plants, asymptomatic weed species and indicator plants (N. rustica) inoculated with selected Iranian TMV isolates including isolates TA, TV, TM, G42 and PU1. First-strand cDNA synthesis was carried out according to the manufacturer’s instructions (MBI, Fermentas, Germany) in a 20 μl reaction containing 10 ng of total RNA, 20 pmol reverse primers (Tob-Uni-1), 50 mM Tris-HCl, pH 8.3, 50 mM KCl, 10 mM of each of the four dNTPs, 20 mM DTT, 20 U Ribonuclease inhibitor and 200 U M-MuLV Reverse transcriptase (MBI, Fermentas, Germany). RT-PCR was carried out with specific degenerate primers designed by Letschert et al. (2002). Degenerate primers Tob-Uni 1 (5′-ATTTAAGTGGASGGAAAAVCACT-3′) and Tob-Uni 2 (5′-GTYGTTGATGAGTTCRTGGA-3′) were used for the genus-specific detection of tobamoviruses that were detected serologically using a specific polyclonal antibody in DAS-ELISA. Specific detection of tobamoviruses species was also performed with a semi-nested PCR method according to Letschert et al. (2002) using Tob-Uni-1 as the downstream primer and TMV-spec (5′-CGGTCAGTGCCGAACAAGAA-3′), ToMV-spec (5′-CGGAAGGCCTAAACCAAAAAG-3′), and TMGMV-spec (5′-AARTAAATAAYAGTGGTAA GAAGGG-3′) forward primers for amplifying the coat protein (CP) gene of TMV, ToMV and TMGMV respectively. The primers were synthesized by MWG Biotech (Germany). The amplified fragment was expected to be ca. 690 bp, 694 bp and 720 bp for TMV, ToMV and TMGMV, respectively, which encompasses 480 bp of the CP coding region plus sequences of the 5′ and 3′ flanking regions. Amplification was carried out in a thermocycler (Eppendorf, Hamburg, Germany) using initial denaturation at 94 °C for 5 min followed by 30 cycles at 94 °C for 60 s, 60 °C for 45 s, 72 °C elongation for 60 s, and a final extension step at 72 °C for 5 min. The PCR mixture (final volume of 25 μl) contained 2.5 μl of cDNA, 2 mM MgCl2, 0.2 mM of each dNTP, 0.4 μM of each primer, and 2.5 U of Taq DNA polymerase (Cinnagene, Iran) in the buffer recommended by the manufacturer. PCR products were checked on gels of 1.5% agarose in 0.5 × TBE buffer (45 mM Tris-borate, 1 mM EDTA, pH 8.0). containing ethidium bromide.

The TMV amplicons (CP gene) for different isolates were obtained from the agarose gels, puried using an Ultra-Clean purification kit as described by the manufacturer (Fermentas Gene JET TM Gel Extraction Kit, Germany), and subjected to dideoxyterminator cycle sequencing by Macrogen (Seoul, South Korea). At least two clones of the CP gene of each TMV isolate were selected and sequenced. The sequences obtained were compared with virus sequences available in the GenBank database using BLAST (Altschul et al., 1997). Phylogenetic analysis of the five representative Iranian TMV isolates was conducted by comparing separately the 480 nucleotides (nt) of the whole CP gene and the 161 nt that included the flanking regions of the CP-coding region with the comparable sequences of other TMV isolates from GenBank. Multiple nucleotide sequence alignments of the Iranian isolates with the previously reported TMV isolates were performed using CLUSTALX 1.8 (Pearson and Lipman, 1988). Phylogenetic relationships were inferred by using the Maximum Likelihood method based on the Kimura 2-parameter model with 1000 bootstraps replicates to estimate node significances (Kimura, 1980), using a CP gene sequence of Odontoglossum ringspot virus (ORSV) as an outgroup. This comparison extended to deduced amino acid sequences. Analyses were performed with the program MEGA 5.05 (Tamura et al., 2011). Using this software and the Tamura-3-parameter nucleotide substitution model (Tamura, 1992), the nucleotide distances and diversities were assessed. Genetic differentiation and the gene flow level between subpopulations also were evaluated by the statistical test Fst(Weir and Cockerham, 1984). Fst ranges from 0 to 1 for undifferentiated to fully differentiated populations, respectively. The statistical tests for estimation of Fst were performed by DnaSP version 5.10.01 (Librado and Rozas, 2009). The ratio of non-synonymous nucleotide diversity to synonymous nucleotide diversity [Pi(a)/Pi(s)] was estimated by the use of DnaSP version 5.10.01 for the Iranian TMV isolates. The value of this ratio is 1, < 1 and > 1 under neutral, negative (purifying selection) and positive (directional selection), respectively.

Results

Field survey and incidence of tobamovirus infection in symptomatic horticultural crops.

Of 1,926 crop leaf samples showing symptoms that were tested by ELISA, tobamoviruses were found in 218 samples (11.3% incidence). The results showed that 53 out of 78 fields surveyed had tobamovirus infections, ranging from 1.5 to 43% of the plants showing symptoms. Apart from Qazvin and Khuzestan provinces, tobamoviruses were distributed among the fifteen surveyed provinces with higher and lower disease incidence in Golestan and Isfahan provinces respectively (Fig. 1, Table 1). Except for East Azarbaijan and Mazandaran, the incidence of symptomatic plants was more frequent in the northern (Alburz, Golestan and West Azarbaijan) or central (Qom and Tehran) provinces compared to the southern (Boushehr, Fars, Hormozgan, Isfahan, Kerman, Khuzestan and Sistan and Baluchestan) and western (Hamedan and Kermanshah) provinces of Iran. (Fig. 1). There were differences in the relative incidence of tobamovirus infection among the horticultural crops tested. The overall picture showed that tobamoviruses were distributed mostly on solanaceous crops (with an average incidence of 12.3%) through the sampled regions compared to cucurbitaceous plants (with an average incidence rate of 8.1%) with the few exceptions of the Larestan, Zabul, Qazvin, Jiroft, Dezful and Minab regions, where no virus was detected in any solanaceous samples tested (Table 1). Tobamovirus infection rates in crop samples showing symptoms, in decreasing order, were tomato (16.9%), eggplant (14.2%), watermelon (12.0%), zucchini (11.0%), pepper (9.7%), potato (8.2%), pumpkin (6.9%), cabbage (6.0%), bean (5%), and cucumber (2.6%) (Table 1). Tobamoviruses were not detected in radish. Symptomatic plants affected by tobamovirus infection were observed in most of the fields surveyed; however, the incidence varied according to the plant species and the geographical regions. Tobamovirus infection was characterized mainly by leaf and fruit mosaics. Plants infected with tobamoviruses showed severe stunting, malformation, yellow spotting on leaves, vein yellowing, mosaic patterns of light and dark green on the leaves and fruits, twisting of leaves, shoestring, interveinal and systemic chlorosis, leaf roll and necrosis.

Fig. 1.

Map of Iran showing the location of 17 provinces (shaded) where field-grown horticultural crops were surveyed for tobamovirus infection during the 2009 to 2010 growing seasons. Two provinces in which tobamovirus infection was not detected are shown spotted in the map.

Weed plants harboring tobamoviruses in horticultural crop fields.

Of 107 weed leaf samples tested by the DAS-ELISA, tobamoviruses were found in 12 samples. Twenty-six weed species belonging to 12 families were tested for the presence of tobamoviruses. There were seven weed species representing five families (Asteraceae, Chenopodiaceae, Convolvulaceae, Fabaceae and Malvaceae) in which tobamoviruses were detected in one or more locations (Table 2), attesting to the wide natural host range of these viruses. Tobamoviruses were detected serologically in two samples of C. album L., five samples of Chenopodium murale L., and one sample each of Cirsium arvensis L., Convolvulus arvensis L., Malva neglecta Wallr., Sonchus oleraceus L., and Trigonella persica L. (Table 2). Tobamoviruses were distributed mostly on weed species in the Chenopodiaceae family (23.3%). None of the infected weeds showed virus symptoms. Infection of weed plants with tobamoviruses were confirmed following mechanical inoculation to N. tabacum cv. Samsun nn, N. tabacum cv. Virginia, N. tabacum cv. White Burley, and P. hybrida. Tobamoviruses caused necrotic local lesions, mild mosaic, vein clearing and leaf deformation on different inoculated indicator plant species. The presence of tobamoviruses in weed species also was verified by an RT-PCR test.

Host range studies.

The results of host range studies and the symptoms associated with virus infection on indicator plants were generally similar to those of earlier reports (Brunt et al., 1995; Chung et al., 2007; Hollings and Huttinga, 1976; Jung et al., 2002), but with some differences (Table 3). The tobamovirus isolates, which were shown subsequently to be TMV isolates, were grouped into five pathological clusters (represented by isolates TA, TV, TM, G42 and PU1) based on differences in symptoms on host plants and nucleotide sequences of their CP genes. The representative isolates were detected in tomato from Alburz (TA), Tehran (TV) and East Azarbaijan (TM), eggplant from Golestan (G42) and pumpkin (PU1) from Fars province (Table 4). Of the 10 indicator plants tested (eight species from four families, including three varieties of tobacco), eight plants (six species) from three plant families developed symptoms after inoculation with the above Iranian TMV isolates. All five TMV isolates induced necrotic local lesions on inoculated leaves of faba bean and N. rustica (Table 3). Except for TMV-PU1, the other isolates were not infective in cucumber and cowpea cv. Mashhad plants (Table 3). Infections of cucumber with the PU1 isolate remained symptomless. This isolate elicited interveinal chlorosis in cowpea. Infection of cucumber with the PU1 isolate was confirmed by using RT-PCR and by mechanical inoculation onto other indicator plants (data not shown). TMV-PU1 differed from other TMV isolates and known strains of TMV (TMV-U1) in producing leaf curling in C. quinoa, vein clearing in N. tabacum cv. Virginia and necrotic local lesion in P. hybrida. TMV-G42 showed similar responses in the ten indicator plants as TMV-U1.

Accession numbers, host and origin of reported TMV isolates used for phylogenetic comparison of coat protein gene sequence in this study

The TMV isolates TM, TV and TA, all detected first in tomato, caused distinct symptoms on inoculated leaves of C. quinoa, N. tabacum cv. Samsun nn and tomato. While TMV-TM and TMV-U1 caused systemic mosaic symptoms in inoculated N. tabacum cv. Samsun nn, TMV-TV induced necrotic local lesion and TMV-TA induced leaf deformation in the same tobacco cultivar, in addition to systemic mosaic. TMV-TM was also distinct from other tomato isolates and TMV-U1 in not producing leaf curling in tomato but inducing necrotic local lesions in C. quinoa (Table 3).

RT-PCR and phylogenetic relationships based on nucleotide sequence analysis.

To confirm the serological diagnoses, RT-PCR was done using samples from naturally infected crop plants (218 samples that positively reacted with the anti-tobamovirus antiserum in the ELISA test) and selected weed species (12 samples). Except for radish, infection of the different crop and weed samples with TMV was confirmed by RT-PCR and specific primers (Fig. 2 and data not shown). RT-PCR also was successful in detecting TMV in the inoculated indicator plants. A single DNA fragment of the expected size (694 bp) was obtained by RT-PCR with species-specific primers when the infected samples were examined for TMV (Fig. 2B and data not shown). No PCR products were amplified from healthy plants.

Fig. 2.

Electrophoresis analysis of RT-PCR products. The products were obtained using (A) the degenerate primer pair Tob-Uni 1/Tob-Uni 2 for detecting tobamoviruses in general, or species-specific primer pairs for detecting (B) Tobacco mosaic virus (TMV), (C) Tomato mosaic virus (ToMV), or (D) Tobacco mild green mosaic virus (TMGMV). (A), Lanes 1 to 4 represent amplified products from naturally infected bean, cabbage, pumpkin and radish samples, respectively; (B) and (C), Lanes 1 to 3 represent amplified products for naturally infected tomato, eggplant and cucumber samples, respectively; (D), Lane 1 represent amplified products for a naturally infected tomato sample. Lane M, 100-bp DNA ladder (Fermentas, Germany); Lane H, negative controls (PCR conducted using healthy plant material).

Using semi-nested RT-PCR and species-specific primers for ToMV and TMGMV, we were able to ascertain the relative incidence of these two viruses in mixed infection with TMV in infected samples (Table 5). For this purpose, a total of sixty crop samples, for which infection with TMV was demonstrated first by RT-PCR, were tested by species-specific primers for the presence of ToMV or TMGMV. PCR products with the expected size of approximately 690 bp and 720 bp were obtained for ToMV and TMGMV, respectively, in some TMV-infected samples (Figs. 2C and D). Of the 60 tested samples, mixed infections of TMV with ToMV were found in 27 samples (45%) including cucumber, eggplant, pumpkin, tomato, watermelon and zucchini in Fars, Alburz, East and West Azarbaijan, Mazandaran, Golestan and Tehran provinces (Table 5). Mixed infections of TMV with either TMGMV or ToMV were detected only in four tomato samples (6.6%) in Tehran and Alburuz provinces (Table 5). Plants with mixed infections showed more severe symptoms compared to those just singly infected with TMV.

Relative incidence of ToMV and TMGMV present in TMV-infected crop samplesb

A total number of the 116 crop samples with virus-like symptoms including bean (58 samples), cabbage (32 samples) and radish (26 samples) collected from Alburz, East Azarbaijan, Fars and Tehran provinces, which failed to react with polyclonal antisera for tobamoviruses in ELISA, were tested by degenerate primers (Tob-Uni 1/Uni 2) for the presence of other tobamoviruses. A PCR-amplified product of approximately 804 bp, corresponding to the CP gene of tobamoviruses plus sequences of 5′ and 3′ flanking regions, was obtained for all of these samples (Fig. 2A and data not shown). However, no PCR products were amplified when the same samples were tested for ToMV and TMGMV (data not shown), which may indicate the presence of yet other tobamoviruses.

The 480 nt sequences of the entire CP sequence in addition to 3′ and 5′ flanking sequences of the Iranian TMV isolates were obtained and deposited in the GenBank as accession numbers HQ593617, HQ593619 and HQ593621 for tomato isolates TA, TM and TV, HQ593620 for eggplant isolate G42 and HQ593618 for pumpkin isolate PU1, respectively (Table 4). A comparison of the nucleotide sequences for the CP gene of Iranian TMV isolates with those of other TMV isolates in the GenBank confirmed that these detected viruses were TMV. Nucleotide sequences of the complete ToMV and TMGMV CP genes that were obtained from RT-PCR positive samples also were deposited into the GenBank as accession numbers HQ593626 and HQ593622, respectively. Comparison of the nucleotide sequences of these two viruses with those of other isolates in the GenBank confirmed that these detected viruses were ToMV and TMGMV (data not shown). Thus, both ToMV and TMGMV were isolated from tomato in Tehran.

BLAST analysis of the TMV sequences of the CP gene disclosed complete identity (100%) among the Iranian TMV isolates. The maximum and minimum nucleotide sequence identities of the entire CP gene (480 nt) between the Iranian isolates and those deposited previously in GenBank ranged from 93 to 94.5%, whereas when the entire CP gene sequence plus 161 nt flanking regions (641 nt) were compared they showed 90 to 94% sequence identity (data not shown). Phylogenetic analysis based on the nucleotide sequences of the CP gene resulted in classification of TMV isolates to three major clades with high bootstrap support in which isolates in clade I were divided into two subgroups (Fig. 3A). All of the Iranian TMV isolates clustered in clade II, with 94% similarity to TMV-U1 (the type TMV species present in clade I). Isolates Ohio V from USA together with Pet-TW from Taiwan were placed in clade III sharing 86–88% identity to TMV-U1. All the other TMV isolates (Table 4) were clustered in either the larger subgroup IA, sharing 96–100% nt sequence identities with isolate TMV-U1 and comprising isolates from Asia, Africa, Europe and America, or the smaller subgroup IB, composed of only two isolates, one each from China and Brazil (Fig. 3A). The phylogenetic tree based on analysis of the CP amino acid sequences showed a different pattern compared to that of the nucleotide sequences of the CP genes (Fig. 3B). Here, the sequences of the amino acids of the CP gene formed two subgroups radiating from one large clade. Except for the Taiwan and one USA isolate (from clade III of the nucleotide sequence comparisons) in subgroup B, the worldwide isolates were placed in subgroup A, with the highest amino acid identity of 99% to the isolate TMV-U1. The amino acid sequence identity of the complete CP (159 aa) of the Iranian TMV isolates compared to TMV-U1 and other isolates ranged from 99 to 100%. It should be noted that in these analyses we have not used all isolates from the GenBank that were labeled as TMV. This is because the CP of those isolates we did not use (accession numbers AF012927, AF103779,AF318214, AJ310339, and JN204355) were actually ToMV in sequence, as determined by sequence alignments, but the viruses were named TMV.

Fig. 3.

Phylogenetic tree of TMV isolates constructed from the alignment of (A) 480 nucleotide and (B) 159 amino acid sequences of the CP gene of five TMV isolates from Iran, together with the homologous gene of TMV isolates from Brazil, China, Egypt, Germany, India, Japan, Serbia, South Korea, Taiwan, Thailand, USA and Venezuela (Table 4) using the neighbor-joining method based on 1000 replicates. Bootstrap percentages of clades are shown along internal branches of trees derived from bootstrap-resampled data sets. Branch lengths represent genetic distances. The tree was rooted with the CP of Odontoglossum ringspot virus (ORSV), a member of the genus Tobamovirus, which was used as an outgroup species (Table 4). Detailed information for each TMV isolate concerning the country where it was isolated and the accession number is found in the phylogenetic trees.

The pairwise nucleotide distances between isolates within a clade were 0.023, 0.00 and 0.029 for Clades I, II and III, respectively, whereas the distance between two clades was 0.149. The distances of TMV isolates with respect to the ORSV sequence used as an outgroup were between 0.987 and 1.149. This may indicate that TMV could be classified into three genotypes. No correlation was observed between the genetic variation and the geographical origin from which the isolates were obtained. However, in comparing the genetic variation between the isolates in each clade, it seems there is a correlation between genetic variations and the species composing each genotype. Low nucleotide diversity (< 0.017) was observed within each subgroup indicating spatial genetic stability in different TMVs isolated throughout the world. Nucleotide diversity was also high and about 0.046 when any two clades were compared with each other. To estimate the extent of genetic differentiation, we measured the coefficient Fst, which is also an estimate of gene flow (Weir and Hill, 2002). The Fst values comparing clades I and II, as well as clade II and each subgroup (A and B) of clade I, were 0.914, 0.541, and 0.723 respectively, an indication of frequent gene flow. However Fst values were lower than 0.312 between two the subgroups in clade I. The Fst values were lower than 0.1 between Asian and American or between European and either Asian or American isolates, indicating a low intercontinental gene flow between TMV isolates. The ratio of dN/dS was < 1 for CP genomic region of the Iranian TMV isolates, providing evidence of negative selection in this gene. The three statistic tests, Tajima D, Fu and Li’s, were significantly negative (P < 0.05), suggesting negative or purifying selection.

Discussion

In recent years, there have been reports of increasing importance of virus diseases in horticultural crops in Iran (Bananej et al., 2006; Bananej and Vahdat, 2008; Dizadji and Shahraeen, 2009; Massumi et al., 2007; Pourrahim et al., 2007; Soleimani et al., 2011). In these earlier studies, the distribution and incidence of viruses infecting horticultural crops were studied by serological methods, but in none of these cases was TMV studied. Although the presence of TMV has been reported in tomato, tobacco and soybean from Iran (Golnaraghi et al., 2004; Khateri et al., 2008; Massumi et al., 2009), the distribution and incidence of this virus in different horticultural crops was not reported so for. Precise identification of TMV and knowledge about its distribution, genetic variation and alternative host sources in such crops will allow farm advisors and growers to make better management decisions. To achieve this, in the present study we sought to determine the incidence and distribution of TMV on economically important crops in the major horticultural crop-producing areas in 17 provinces throughout Iran. Natural weed hosts of TMV and the genetic variation of TMV also were investigated.

Our results showed that 11.3% of the total samples assessed from different horticultural crops were infected with TMV. This incidence, which is limited to the field plants showing symptoms, can be considered low if compared with that of other countries. For example, the incidence of TMV in pepper fields in Turkey, neighboring the north-eastern region of Iran, was found to be 21.5% (Arli-Sokmen et al., 2005). The infection of horticultural crops by TMV was not equally distributed in Iran, since it ranged from a total absence in some provinces (Khuzestan and Qazvin) to a maximum of 43% in another (Golestan). Widespread prevalence of TMV in most crop-growing areas of Iran would implicate its major impact on crop production throughout the country. To our knowledge, this is the first report on the occurrence of TMV in bean, cabbage, cucumber, pepper, potato, pumpkin, watermelon and zucchini in Iran based on biological reactions, ELISA and RT-PCR. There was a clear-cut correlation between TMV incidence, geographical distribution and the type of vegetable hosts. Although TMV is not considered an economically important bean virus and is not common in bean fields in Iran, we detected TMV infected bean plants in Golestan (Table 1), which may reflect a very high inoculum in this province. Results of our study showed that TMV is present throughout the country, but a particularly extensive spread was recorded in Golestan (43%) and West Azarbaijan (18.7%) provinces (northern regions of Iran), where most of the national tobacco production is concentrated. Previously, TMV had been reported from different cultivars of tobacco with an incidence rate of 28.6% in these regions (Khateri et al., 2008). The relatively high incidence and prevalence of TMV may be related to the presence of tobacco fields near horticultural crops in these provinces. However, this relationship has yet to be confirmed. In the central provinces (Alburz, Qom and Tehran), TMV was detected in 18.8% of the total tested samples. TMV infection was confined to the fields with a history of solanaceous plant production in these regions. It is possible that failure to detect and subsequently eradicate TMV from these farms has led to its persistence in these areas. The virus can survive for many years in dead, dried plant material, which makes it difficult to inactivate (Broadbent, 1965). Plants that remain in the soil may harbor mechanically transmitted virus and act as a source of infection for the next crops, especially where the susceptible crops are grown in the same fields yearly (Hu et al., 1994; Lewandowski et al., 2010).

The results of ELISA tests showed the presence of TMV in various weed species distributed in the surveyed vegetable fields (Table 2). TMV was serologically detected previously in C. murale L. and Chrozophora hierosolymitana from some tomato cultivations in southern regions of Iran (Massumi et al., 2009). However, the natural occurrence of TMV on C. album L, Cirsium arvensis L, Convolvulus arvensis L, M. neglecta Wallr, S. oleraceus L, and T. persica L, based on host range, ELISA, and PCR, is reported for the first time in this study. TMV overwinters in a number of weeds (Shiel and Castello, 1985; Arli-Sokmen et al., 2005). Primary infection of horticultural crops from overwintering sources of the virus can be a contributing factor in the increasing virus incidence in different plantations. Having knowledge of weed reservoirs of TMV is essential for understanding the virus epidemiology and must be understood for successful control through sanitation practices. In our survey, most weed species found to be infected with TMV were symptomless. Previous research also has shown that many TMV-infected herbaceous weed hosts often do not display foliar symptoms (Shiel and Castello, 1985). Nonetheless, further studies comprising more sampling of weed with regard to the symptom status are necessary to improve our knowledge of virus spread in the region.

The number of solanaceous plants shown to be infected with TMV was relatively high compared to collected cucurbit plant species (Table 1). Solanaceous crops incur the most significant losses from TMV-induced disease (Hollings and Huttinga, 1976). TMV is seed-borne mostly in solanaceous crops and they may spread during transplanting and harvesting after primary infections via infected seeds. Some growers in Iran use their own non-certified seed produced locally and this may contribute to the spread of TMV in these crops. It is noteworthy that TMV, a virus hitherto unreported from cucurbit plants in Iran, prevailed in having high rates of incidence in two cucurbit species (watermelon and zucchini) in the surveyed regions.

Due to the similarity of symptoms induced by TMV with those induced by other tobamoviruses and the possibility of mixed infection of TMV with other tobamoviruses, we tested the TMV infected samples for the presence of TMGMV and ToMV (other tobamoviruses in subgroup I). Mixed infections of tomato with TMV and TMGMV were only observed in Tehran and Alburz. We have previously reported the occurrence of TMGMV in tomato in Tehran (Alishiri et al., 2011). The low incidence of this virus in Iran may suggest that it is a new pathogen in the country. Mixed infection of tomato, eggplant, cucumber, watermelon, zucchini and pumpkin with TMV and ToMV were encountered at Fars, Alburz, East and West Azarbaijan, Mazandaran, Golestan and Tehran as attested by the data in Table 5. Although, mixed infection of solanaceous crops with TMV and ToMV was reported previously in some parts of Iran (Aghamohammadi et al., 2011; Ahoonmanesh et al., 1992; Massumi et al., 2009), the data shown in this study represent the first report on the occurrence of ToMV and its co-infection with TMV in cucurbit plants including cucumber, pumpkin, watermelon and zucchini in Iran. The higher incidence of TMV and ToMV mixed infections compare to TMV and TMGMV mixed infections may indicate that they arose from a common ancestor before colonizing and spreading in Iran.

The host ranges of different TMV isolates collected from different hosts and regions were similar to what has been previously reported for TMV (Chung et al., 2007; Jung et al., 2002; Novikov et al., 2000), but considerable differences were found in the symptoms induced by each TMV isolate on the indicator plants. The five TMV pathology clusters studied in this research showed obvious differences in the symptom development on C. quinoa, cowpea cv. Mashhad, N. tabacum cv. White Burley and tomato. Differences in symptom development on C. quinoa, cowpea and N. tabacum cv. Virginia observed after inoculation with TMV-PU1 (pumpkin isolate) distinguished this isolate from the tomato infecting isolates. TMV-PU1 also was differentiated from other isolates by infecting cucumber plants (Table 3). By contrast, the results also showed that three tomato-infecting TMV isolates (TA, TV and TM), collected from different geographical regions, differed in symptom developments on the indicator plants. Taken together, the data obtained by bioassay suggest clear differences between Iranian isolates of TMV which may be due to the genetic variation among TMV strains present in different hosts and geographical locations and/or climate regions in Iran. Many different strains of TMV are present in nature, which are different in biological properties (Van Regenmortel, 1986). It has been shown that sequences in TMV CP are often responsible for differences in symptom development (Culver, 2002). The work of Aldaoud et al. (1989) has shown that TMV isolates are differentially selected by different plant species and selective pressure plays an important role in TMV genetic divergence. This is in contrast to the results obtained by phylogenetic analysis of Iranian isolates that showed low genetic variations in the CP gene nucleotide sequence, which suggests the negative selection acting on this gene. We found significant evidence for negative selection in the CP gene of a large subset of TMV isolates. In our data set, most parts of the ORF were determined to be under negative selection for the Iranian isolates. For most of these TMV genes, the dN/dS ratio was smaller than 1.0, which was consistent with the negative selection against protein changes. The negative selection in the CP gene would result in a higher protein tolerance against mutation and pressure selection. This may indicate that TMV CP is not flexible in the interaction with other proteins or has little specific interaction with the host or potential vectors. This result also supports the hypothesis that the CP was undergoing intensive evolution. Previously, negative selection was reported for the worldwide ToMV isolates on the basis of the CP gene sequence (Rangel et al., 2011). However, the results indicated a low genetic diversity between isolates within clades and high diversity between the three clades, suggesting effects of selection under different selection pressures. This finding may suggest that different TMV isolates that are best fitted in a given geographical location/ climate have increased in the population. This conclusion has been demonstrated for resistance-breaking isolates of ToMV that overcome the resistance of the Tm-1 gene (Pelham et al., 1970).

Phylogenetic analysis resulted in a tree with three main parts. Iranian TMV-crop isolates obtained from eggplant, pumpkin and tomato, respectively were all clustered together in a distinct branch. No branching pattern based on differences in the host source was found among the TMV-crop isolates studied. All TMV isolates which are distributed in Asia, Africa, Europe and America were located in one clade which could be divided into two subgroups with a branching pattern specific to geographical distribution. These indicate no clear correlation between CP sequence variation and the geographical origins of the virus isolates, in contrast to the results obtained by Jung et al. (2002) who differentiated TMV-potato isolates on the basis of geographical locations/climates. Previous studies elucidated that TMV isolates may be altered by different temperatures (Jones and Dawson, 1978) or in different plant species (Yarwood, 1979). Our results showed that the diversity values are high between the three clades, which may suggest that the populations in each clade could be derived from a single origin. Considering the low number of isolates in clade II, it is possible that the viruses in clade II might represent introductions that have occurred recently. The observed low diversity and the determination that TMV isolates from different host species had the same predominant genotype suggest that host species do not contribute to differentiation of the virus population, as observed with other plant viruses (Garcia-Arenal et al., 2001). No correlation was observed between the genetic distance and the geographical origin or the host species from which the isolates were obtained.

Clade III was composed of TMV isolates from the USA and Taiwan, whereas Clade I was composed of isolates from USA and the rest of the world. Low genetic variation among geographically distant isolates also has been reported for other tobamoviruses (Fraile et al., 1996; Rodriguez-Cerezo et al., 1989). The apparent absence of a natural vector might have prevented population changes and might account for the low genetic diversity observed. However, the possibility of subjection to host selection may not be conceivable for the Iranian isolates, since they were collected from different hosts (Table 2).

In this study, the statistical test Fst was used to assess genetic differentiation and the gene flow level between clades and subgroups. Interestingly, in most cases where Iranian isolates were compared against other geographical groups in subgroup A (for example, Europe and the Americas) the Fst values were low but when pairwise comparisons were repeated using the Iranian isolates in clade II and isolates in other clades, the Fst values were high. This result, together with the phylogenetic analyses, demonstrates the occurrence of gene flow from other isolates to the Iranian population. Our conclusions are consistent with previous research on ToMV (Rangel et al., 2011) showing that limited gene flow occurred between subpopulations from different geographical locations/climates. Overall, a high degree of differentiation was noted between three clades. It may suggest that TMV populations in each clade went through frequent bottlenecks when colonizing new habitats and members of each founder group were from similar population sources. The Fst values between the Iranian–Asian and Iranian–European isolates were 0.678 and 0.914, respectively, indicating that most of the molecular diversity is distributed among these geographical groups.

In conclusion, TMV isolates grouped into three main clades or genotypes worldwide and in combination with the recent work of Rangel et al. (2011), our analysis reveals that high genetic stability of different TMV isolates could be favored by genetic shift after bottlenecks. The genetic structure of TMV populations in different horticultural crops has not been analyzed in other Asian countries and it is not known if the Iranian isolates are present in other countries. Nonetheless, to control TMV infection, sanitary measures should be implemented to minimize TMV dispersal through plant material between countries. These results might provide useful tools to study the interplay between the evolutionary and epidemiological processes acting on TMV population and could help us to design some efficient control strategies for TMV emergence.

Acknowledgements

This work was supported by the Plant Pathology Department of Islamic Azad University (Tehran). PP was support-ed by a grant from the Research Professorship Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (MOEHRD-KRF-2012-0007187) and by Grant No. RDA-PJ007984042013 from the Next Generation BioGreen 21 Program of Rural Development Administration of the Rep. of Korea.

References

Adkins S, Kamenova I, Achor D, Lewandowski DJ. 2003;Biological and molecular characterization of a novel tobamovirus with a unique host range. Plant Dis 87:1190–1196.
Aldaoud R, Dawson WO, Jones GE. 1989;Rapid, random evolution of the genetic structure of replicating tobacco mosaic virus populations. Intervirology 30:227–233.
Alexandre MAV, Soares RM, Rivas EB, Duarte LML, Chagas CM, Saunal H, Van Regenmortel MHV, Richtzehain LJ. 2000;Characterization of a strain of Tobacco mosaic virus from Petunia. J. Phytopathol 148:601–607.
Alishiri A, Rakhshandehroo F, Zamanizadeh HR. 2011;First report of Tobacco mild green mosaic virus infecting tomato in Iran. New Dis. Rep 23:30.
Aghamohammadi V, Rakhshandehroo F, Shams-Bakhsh M. 2011;First report of Tomato mosaic virus in eggplant in Iran. J. Plant Path 93:S4.69.
Ahoonmanesh A, Bahar M, Ghobadi C. 1992;Tomato mosaic virus in Iran. Iran. J. Plant Pathol 28:1–4.
Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ. 1997;Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 25:3389–3402.
Alonso E, Garcia-Luque I, Avila-Rincon MJ, Wicke B, Serra MT, Diaz-Ruiz JR. 1989;A tobamovirus causing heavy losses in protected pepper crops in Spain. J. Phytopathol 125:67–76.
Antignus Y, Wang Y, Pearlsman M, Lachman O, Lavi N, Gal-on A. 2001;Biological and molecular characterization of a new cucurbit-infecting Tobamovirus. Phytopathology 91:565–571.
Arli-Sokmen M, Mennan H, Sevik MA, Ecevit O. 2005;Occurrence of viruses in field-grown pepper crops and some of their reservoir weed hosts in Samsun, Turkey. Phytoparasitica 33:347–358.
Bananej K, Desbiez C, Wipf-Scheibel A, Vahdat A, Kheyr-Pour A, Ahoonmanesh A, Lecoq H. 2006;First report of Cucurbit aphid-borne yellows virus in Iran causing yellows on four cucurbit crops. Plant Dis 90:526.
Banane K, Vahdat A. 2008;Identification, distribution and incidence of viruses in field-grown cucurbit crops of Iran. Phytopath. Medit 47:247–257.
Broadbent LH. 1965;The epidemiology of tomato mosaic. Seed-transmission of TMV. Ann. Appl. Biol 56:177–205.
Brunt AA, Crabtree K, Dallwitz MJ, Gibbs AJ, Watson L. 1995. Viruses of plants: Description and lists from the VIDE database CAB International. Wallingford, UK:
Channuntapipat C, Sedgley M, Collins G. 2001;Sequences of the cDNAs and genomic DNAs encoding the S1, S7, S8, and Sf alleles from almond, Prunus dulcis. Theor. Appl. Genet 103:1115–1122.
Cherian S, Muniyapppa V. 1998;ELISA based survey and host range of tomato mosaic tobamovirus. Ind. J. Virol 14:65–69.
Chitra TR, Prakash HS, Albrechtsen SE, Shetty HS, Mathur SB. 2002;Indexing of leaf and seed samples of tomato and bell pepper for tobamoviruses. Ind. Phytopathol 55:84–86.
Choi SK, Yoon JY, Chung BN. 2009;Genome analysis and characterization of a tobacco mosaic virus isolate infecting balsam (Impatiens balsamina). Arch. Virol 154:881–885.
Chung BN, Kim JS, Cho JD, Cheong SR, Jeong MI. 2007;Tobacco mosaic virus detected in vegetatively propagated petunia hybrid ‘Surfinia’. Plant Pathol. J 23:34–36.
Clark MF, Adams AN. 1977;Characteristics of the micro-plate method of enzyme-linked immunosorbent assay for the detection of plant viruses. J. Gen. Virol 34:475–483.
Culver JN. 2002;Tobacco mosaic virus assembly and disassembly: Determinants in pathogenicity and resistance. Annu. Rev. Phytopathol 40:287–308.
Dizadji A, Shahraeen N. 2009;Occurrence, distribution and seasonal changes of viruses infecting common bean in north-western Iran. Arch. Phytopath. Plant Protec 44:1647–1654.
F.A.O.. 2009. F.A.O. Statistical Databases. Preliminary 2009 production data. http://faostat.fao.org/site/339/default.aspx.
Fraile A, Malpica JM, Aranda MA, Rodriguez-Cerezo E, Garcia-Arenal F. 1996;Genetic diversity in tobacco mild green mosaic tobamovirus infecting the wild plant Nicotiana glauca. Virology 223:148–155.
Fu YX, Li WH. 1993;Statistical tests of neutrality of mutations. Genetics 133:693–709.
Garcia-Arenal F, Fraile A, Malpica JM. 2001;Variability and genetic structure of plant virus populations. Ann. Rev. Phytopathol 39:157–186.
Goelet P, Lomonossoff GP, Butler PJG, Akam ME, Gait MJ, Karn J. 1982;Nucleotide sequence of tobacco mosaic virus RNA. Proc. Natl. Acad. Sci. USA 9:5818–5822.
Golnaraghi AR, Shahraeen N, Pourrahim R, Farzadfar Sh, Ghasemi A. 2004;Occurrence and relative incidence of viruses infecting soybeans in Iran. Plant Dis 88:1069–1074.
Hollings M, Huttinga H. 1976. Tomato mosaic virus. CMI/ AAB Descriptions of plant viruses, no. 156 Commonwealth Mycology Institute/Association for Applied Biology. Kew:
Hu JS, Ferreira S, Xu MQ, Lu M, Iha M, Pflum M, Wang M. 1994;Transmission, movement and inactivation of Cymbidium mosaic and Odontoglossum ringspot viruses. Plant Dis 78:633–636.
Jones GE, Dawson WO. 1978;Stability of mutations conferring temperature sensitivity on tobacco mosaic virus. Intervirology 9:149–155.
Jung HW, Yun WS, Hahm YI, Kim KH. 2002;Characterization of Tobacco mosaic virus isolated from potato showing yellow leaf mosaic and stunting symptoms in Korea. Plant Dis 86:112–117.
Khateri H, Moarrefzadeh N, Mosahebi G, Koohi-Habibi M. 2008;Virus diseases in the tobacco fields of Gilan and West Azarbaijan provinces of Iran. Comm. Appl. Biol. Sci 73:307–310.
Kimura M. 1980;A simple method for estimating evolutionary rate of base substitutions through comparative studies of nucleotide sequences. J. Mol. Evol 16:111–120.
Kumar S, Udaya Shankar AC, Nayaka SC, Lund OS, Prakash HS. 2011;Detection of Tobacco mosaic virus and Tomato mosaic virus in pepper and tomato by multiplex RT-PCR. Lett. Appl. Microbiol 53:359–363.
Lartey RT, Voss TC, Melcher UK. 1996;Tobamovirus evolution: gene overlaps, recombination, and taxonomic implications. Mol. Biol. Evol 13:1327–1338.
Letschert B, Adam G, Lesemann DE, Willingmann P, Heinze P. 2002;Detection and differentiation of serologically cross-reacting tobamoviruses of economical importance by RT-PCR and RT-PCR-RFLP. J. Virol. Method 106:1–10.
Lewandowski DJ, Hayes AJ, Adkins S. 2010;Surprising results from a search for effective disinfectants for Tobacco mosaic virus contaminated tools. Plant Dis 94:542–550.
Librado P, Rozas J. 2009;DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25:1451–1452.
Massumi H, Samei A, Hosseini Pour A, Shaabanian M, Rahimian H. 2007;Occurrence, distribution, and relative incidence of seven viruses infecting greenhouse-grown cucurbits in Iran. Plant Dis 91:159–163.
Massumi H, Shaabanian M, Hosseini Pour A, Heydarnejad J, Rahimian H. 2009;Incidence of viruses infecting tomato and their natural hosts in the southeast and central regions of Iran. Plant Dis 93:67–72.
Moreira SR, Eiras M, Chaves ALR, Galleti SR, Colariccio A. 2003;Caracterização de uma nova estirpe do Tomato mosaic virus isolada de tomateiro no Estado de São Paulo. Fitopatol. Brasil 28:602–607.
Nassar EA, El-Dougdoug KA, Osman ME, Dawoud RA, Kinawy AH. 2012;Characterization and elimination of a TMV isolate infecting Chrysanthemum plants in Egypt. Int. J. Virol 8:14–26.
Nei M, Kumar S. 2000. Molecular Evolution and Phylogenetics Oxford University Press. New York:
Novikov VK, Belenovich EV, Dobrov EN, Zavriev SK. 2000;Kazakh strains of tobacco mosaic virus: two strains with potentially destabilizing amino acid substitutions in the coat protein. Physiol. Mol. Plant Path 56:71–77.
Padgett HS, Beachy RN. 1993;Analysis of a tobacco mosaic virus strain capable of overcoming N gene-mediated resistance. Plant Cell 5:577–586.
Pearson WR, Lipman DJ. 1988;Improved tools for biological sequence comparison. Proc. Natl. Acad. Sci. USA 85:2444–2448.
Pelham J, Fletcher JT, Hawkins JH. 1970;The establishment of a new strain of tobacco mosaic virus resulting from the use of resistant varieties of tomato. Ann. Appl Biol 65:293–297.
Pourrahim R, Farzadfar Sh, Golnaraghi AR, Ahoonmanesh A. 2007;Incidence and distribution of important viral pathogens in some Iranian potato fields. Plant Dis 91:609–615.
Rangel EA, Alfaro-Fernandez A, Font-San-Ambrosio MI, Luis-Arteaga M, Rubio L. 2011;Genetic variability and evolutionary analyses of the coat protein gene of Tomato mosaic virus. Virus Genes 43:435–438.
Rodriguez-Cerezo E, Moya A, Garcı!a-Arenal F. 1989;Variability and evolution of the plant RNA virus Pepper mild mottle virus. J. Virol 63:2198–2203.
Saitou N, Nei M. 1987;The neighbor-joining method: A new method for reconstructing phylogenetic trees. Mol. Biol. Evol 4:406–425.
Shiel PJ, Castello JD. 1985;Detection of tobacco mosaic and tobacco ringspot viruses in herbaceous and woody plants near virus-infected white ash trees in central New York. Plant Dis 69:791–795.
Soleimani P, Mosahebi G, Koohi-Habibi M. 2011;Identification of some viruses causing mosaic on lettuce and characterization of Lettuce mosaic virus from Tehran Province in Iran. Afric. J. Agric. Res 6:3029–3035.
Tajima F. 1989;Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 123:585–595.
Tamura K. 1992;Estimation of the number of nucleotide substitutions when there are strong transition-transversion and G+Ccontent biases. Mol. Biol. Evol 9:678–687.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. 2011;MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol. Biol. Evol 28:2731–2739.
Van Regenmortel MHV. 1986. Tobacco mosaic virus. Antigenic structure. The Plant Viruses 2In : van Regenmortel MHV, Fraenkel-Conrat H, eds. Plenum Press. New York:
Wang H, Culver JN, Stubbs G. 1997;Structure of Ribgrass mosaic virus at 2.9 A resolution: evolution and taxonomy of tobamoviruses. J. Mol. Biol 269:769–779.
Weir BS, Cockerham CC. 1984;Estimating F-statistics for the analysis of population structure. Evolution 38:1358–1370.
Weir BS, Hill WG. 2002;Estimating F-statistics. Ann. Rev. Genet 36:721–750.
Yarwood CE. 1979;Host passage effect with plant viruses. Adv. Virus Res 25:169–190.
Zaitlin M. 1999;Elucidation of the genome organization of tobacco mosaic virus. Philos. Trans. R Soc. Lond. B Biol. Sci 354:587–591.

Article information Continued

Fig. 1.

Map of Iran showing the location of 17 provinces (shaded) where field-grown horticultural crops were surveyed for tobamovirus infection during the 2009 to 2010 growing seasons. Two provinces in which tobamovirus infection was not detected are shown spotted in the map.

Fig. 2.

Electrophoresis analysis of RT-PCR products. The products were obtained using (A) the degenerate primer pair Tob-Uni 1/Tob-Uni 2 for detecting tobamoviruses in general, or species-specific primer pairs for detecting (B) Tobacco mosaic virus (TMV), (C) Tomato mosaic virus (ToMV), or (D) Tobacco mild green mosaic virus (TMGMV). (A), Lanes 1 to 4 represent amplified products from naturally infected bean, cabbage, pumpkin and radish samples, respectively; (B) and (C), Lanes 1 to 3 represent amplified products for naturally infected tomato, eggplant and cucumber samples, respectively; (D), Lane 1 represent amplified products for a naturally infected tomato sample. Lane M, 100-bp DNA ladder (Fermentas, Germany); Lane H, negative controls (PCR conducted using healthy plant material).

Fig. 3.

Phylogenetic tree of TMV isolates constructed from the alignment of (A) 480 nucleotide and (B) 159 amino acid sequences of the CP gene of five TMV isolates from Iran, together with the homologous gene of TMV isolates from Brazil, China, Egypt, Germany, India, Japan, Serbia, South Korea, Taiwan, Thailand, USA and Venezuela (Table 4) using the neighbor-joining method based on 1000 replicates. Bootstrap percentages of clades are shown along internal branches of trees derived from bootstrap-resampled data sets. Branch lengths represent genetic distances. The tree was rooted with the CP of Odontoglossum ringspot virus (ORSV), a member of the genus Tobamovirus, which was used as an outgroup species (Table 4). Detailed information for each TMV isolate concerning the country where it was isolated and the accession number is found in the phylogenetic trees.

Table 1.

Incidence of tobamoviruses on horticultural crops collected from the main provinces of cultivation in Iran

Sampled horticultural cropsa
Province Region Infected Samplesb Solanaceous Cucurbitaceous Legumes Brassicas




Tomato Potato Eggplant Pepper Cucumber Watermelon Zucchini Pumpkin Bean Radish Cabbage
Fars Larestan 0/69 0/32 0/3 0/9 0/1 0/6 ND 0/3 0/2 0/8 ND 0/5
Fasa 11/92 1/4 9/57 0/2 1/5 ND 0/5 0/4 0/6 0/2 0/4 0/3
Jahrum 3/63 3/42 ND ND 0/14 ND 0/4 ND ND 0/3 ND ND
Darab 36/215 (11.38%)c 36/100 0/17 0/15 ND 0/30 0/1 0/8 ND 0/30 0/2 0/12
Isfahan Isfahan 1/65 (1.5%) 1/47 ND 0/3 ND 0/5 0/2 0/2 0/3 0/3 ND ND
Sistan and Baluchestan Zabul 0/29 0/10 ND ND ND 0/3 0/3 0/3 0/3 0/2 0/4 0/1
Zahedan 3/28 (5.2%) 0/15 ND ND 2/5 1/1 ND 0/2 0/5 ND ND ND
Alburz Karaj 12/72 (16.6%) 3/10 1/6 3/9 2/9 ND 0/6 1/9 2/7 0/5 0/8 0/3
East Azarbaijan Maragheh 12/128 (9.37%) 6/64 ND 1/9 ND 3/40 ND 1/1 1/1 0/4 0/5 0/4
Mazandaran Noshahr 2/62 (3.22%) ND ND 0/13 1/20 ND ND ND 1/12 0/6 0/7 0/4
Golestan minoodasht 34/79 (43%) 14/23 ND 8/16 ND 0/2 2/10 4/18 ND 6/8 ND 0/2
Boushehr Boushehr 1/57 (1.75%) 0/30 ND ND 1/5 ND 0/4 0/3 0/4 ND 0/3 0/8
West Azarbaijan Urmia 20/107 (18.7%) 5/15 ND 1/7 2/9 0/32 9/34 0/1 ND ND ND 3/9
Qazvin Qazvin 0/35 (0%) 0/4 0/8 0/2 ND 0/5 ND 0/4 ND 0/3 0/4 0/5
Tehran 39/110 39/77 ND 0/15 0/4 ND 0/6 ND 0/5 ND 0/3 ND
Tehran Varamin 21/158 4/35 ND 8/45 0/2 0/19 0/3 7/42 0/4 0/1 ND 2/7
Shahriar 2/33 (20.6%) 2/17 ND ND ND ND 0/4 ND ND 0/5 0/4 0/3
Qom Qom 7/57 (12.3%) 1/11 ND 4/13 1/13 0/9 0/3 ND 1/6 ND ND 0/2
Kerman Jiroft 2/65 (3%) 0/22 0/6 ND 0/5 0/9 1/8 ND 1/9 0/6 ND ND
Khuzestan Dezful 0/68 (0%) 0/54 ND ND ND 0/3 ND 0/2 ND ND 0/3 0/6
Hormozgan Roudan 3/83 0/25 ND 2/19 0/10 1/9 0/5 ND ND 0/8 ND 0/7
Minab 0/29 ND ND 0/3 0/5 0/2 0/4 ND 0/7 0/2 0/6 ND
Kehorestan 2/70 1/40 ND 0/1 1/8 0/4 ND 0/7 0/3 0/7 ND ND
Hajiabad 3/66 (3.2%) ND 1/40 ND 1/10 0/3 1/2 0/4 0/3 0/4 ND ND
Hamedan Tuyserkan 3/40 (7.5%) 2/8 1/9 ND ND 0/5 ND ND 0/2 0/8 0/6 0/2
Kermanshah Mahidasht 1/46 (2.1%) 0/14 ND 0/8 1/8 ND 0/4 0/5 0/4 0/3 ND ND

Total 218/1926 (11.3)d 118/699 (16.9)e 12/146 (8.2) 27/189 (14.2) 13/133 (9.7) 5/187 (2.6) 13/108 (12) 13/118 (11) 6/86 (6.9) 6/118 (5) 0/59 (0) 5/83 (6)

The highest incidence for tobamovirus infection in each surveyed host species is shown in bold and was highest for tomato.

a

Data obtained by DAS-ELISA for tobamovirus. Number of infected samples over number of samples collected for each crop. ND = not determined.

b

Nunber of infected samples over number of samples collected for each province.

c

Percentage of tobamovirus infection in each province.

d

Average percentage of tobamovirus infection for all provinces.

e

Average percentage of tobamovirus infection in each crop species.

Table 2.

Incidence of tobamovirus infection on selected weed species

Family species Common name Incidencea
Chenopodiaceae (23.3%)
  Chenopodium murale L. Nettle-Leaf goosefoot 5/17b
  C. album L. Lamb’s quarters 2/7
  C. amaranticolor L. Goosefoot 0/3
  Salsola kali L. Saltwort 0/3
Malvaceae (9%)
  Malva neglecta Wallr. Common mallow 1/11
Portulacaceae (0)
  Portulaca oleraceae L. Common Purslane 0/5
Solanaceae (0)
  Solanum nigrum L. Nightshade 0/4
  Physalis alkekengi L. Bladder cherry 0/3
Convolvulaceae (11.1%)
  Convolvulus arvensis L. Field Bindweed 1/6
  Calystegia sepium L. Hedge Bindweed 0/3
Orobanchaceae (0)
  Orobanche aegyptica Pers. Broomrape 0/3
Polygonaceae (0)
  Rumex vesicarius L. Bladderdock 0/3
  Polygonum aviculare L. Knotweed 0/2
Asteraceae (20%)
  Sonchus oleraceus L Sowthistle 1/3
  Cirsium arvensis L. Canada thistle 1/5
  Lactuca scariola L. Prickly Lettuce 0/2
Amaranthaceae (0)
  Amaranthus retroflexus L. Redroot pigweed 0/3
  Amaranthus hybridus L. Slim amaranthus 0/3
  Spinacia oleracea L. Garden spinach 0/3
Brassicaceae (0)
  Raohanus raphanistrum L. Wild radish 0/3
  Sinaps arvensis L. Wild mustard 0/3
Fabaceae (11.1%)
  Glycyrrhiza glabra L. Licorice root 0/3
  Trigonella persica L. Fenugreek 1/3
  Ononis spinosa L. Spiny rest harrow 0/3
Boraginaceae (0)
  Heliotropium europium L. European heliotrope 0/3

Total 12/107 (11.2%)c
a

Incidence rate of tobamovirus infection on different weed families in the surveyed provinces from Table 1.

b

No. of infected samples/collected samples, determined by DAS-ELISA detection of tobamoviruses infecting different weed species collected from surveyed vegetable fields.

c

Average incidence of tobamovirus infection on selected weed species.

Table 3.

Symptom development on indicator plant species after inoculation with Tobacco mosaic virus (TMV) isolates detected in Iran

Leaf symptoms caused by TMV pathological clustersa

Indicator plant U1b TA TV G42 PU1 TM
Solanaceae
  Nicotiana rustica NL/−c NL/− NL/− NL/− NL/− NL/−
  N. tabacum cv. Samsun nn −/M −/M, Def NL/M −/M −/M −/M
  N. tabacum cv. Virginia NL, M/− NL, M/− NL, M/− NL, M/− NL, M/VC NL, M/−
  N. tabacum cv. White burley −/VC −/VC −/VC −/VC −/VC −/VC
  Solanum lycopersicum M/LC M/LC −/LC M/LC M/− M/−
  Petunia hybrida −/M −/M −/M −/M NL/M −/M
Chenopodiaceae
  Chenopodium quinoa CL/− CL/− CL/− CL/− CL/LC NL/−
Cucurbitaceae
  Cucumis sativus −/− −/− −/− −/− −/−d −/−
Fabaceae
  Vigna unguiculata cv. Mashhad −/− −/− −/− −/− IVC/− −/−
  Vicia faba NL/− NL/− NL/− NL/− NL/− NL/−
a

Abbreviations: CL, cholorotic local lesions; Def, deformation; IVC, interveinal chlorosis; LC, leaf curling; NL, necrotic local lesions; VC, vein clearing; M, mosaic; –, no symptoms.

b

TMV common strain U1 from the Cornell University, Plant Virus collection.

c

Local symptoms/systemic symptoms.

d

The inoculated plants remain symptomless but virus was detectable by ELISA and PCR.

Table 4.

Accession numbers, host and origin of reported TMV isolates used for phylogenetic comparison of coat protein gene sequence in this study

Isolate or Strain Accession no. Host Origin Nt identitya
TA HQ593617 Tomato Iran 94%
PU1 HQ593618 Pumpkin Iran 94%
TM HQ593619 Tomato Iran 94%
G42 HQ593620 Eggplant Iran 94%
TV HQ593621 Tomato Iran 94%
TMV-SoyIn JQ895560 Soybean India 99%
TMV-P AY029262 Petunia Brazil 96%
Chuxiong-1 HE818417 Tobacco China 95%
B935A Z96945 Faba bean China 99%
TMV-Egypt AY686725 Tomato Egypt 98%
Fulvum AF546184 Tobacco Germany 99%
Vulgare AJ429078 Tomato Germany 99%
K98 AB669002 Sweet pepper Japan 99%
699-07 GQ340671 Tobacco Serbia 98%
TMV-Pet DQ981481 Petunia South Korea 98%
IM AB369276 Impatiens South Korea 99%
Pet-TW EF392659 Petunia Taiwan 86%
TMV-pepper AY633749 Chili pepper Thailand 98%
TMV-Pet7 GQ370523 Petunia U.S.A 99%
OhioV FR878069 Petunia U.S.A 87%
Vulgare DQ211985 Canavalia Venezuela 99%
TMV-U1 V01408 Tobacco Germany 100%
ORSVb AB693990 Orchid Indonesia -
a

Nucleotide sequence identity of the CP gene to strain U1

b

Odontoglossum ringspot virus (ORSV), a member of the genus Tobamovirus, which was used as an out-group species.

Table 5.

Relative incidence of ToMV and TMGMV present in TMV-infected crop samplesb

Province Virus Samplesb Horticultural crops
Tomato Eggplant Cucumber Watermelon Zucchini Pumpkin
Fars ToMV
TMGMV
5
ND
2
ND
0
ND
0
ND
ND
ND
ND
ND
ND
ND
Alburz ToMV
TMGMV
5
5
2
3
0
0
0
0
0
0
0
0
0
0
East Azarbaijan ToMV
TMGMV
5
ND
0
ND
0
ND
2
ND
1
ND
ND
ND
ND
ND
Mazandaran ToMV
TMGMV
5
5
1
0
1
0
ND
0
ND
0
ND
0
ND
0
Golestan ToMV
TMGMV
5
5
2
0
2
0
ND
0
ND
0
1
0
1
0
West Azarbaijan ToMV
TMGMV
5
5
2
0
0
0
0
0
3
0
0
0
0
0
Tehran ToMV
TMGMV
5
5
3
2
2
0
0
0
0
0
1
0
1
0
Total ToMV 60 27
TMGMV 5
a

Results obtained by semi-nested RT-PCR from 60 selected TMV-infected crop samples.

b

Number of vegetable samples infected with the specified tobamovirus. ND = not determined.