Skip to main content
  • Research article
  • Open access
  • Published:

The I2 resistance gene homologues in Solanum have complex evolutionary patterns and are targeted by miRNAs

Abstract

Background

Several resistance traits, including the I2 resistance against tomato fusarium wilt, were mapped to the long arm of chromosome 11 of Solanum. However, the structure and evolution of this locus remain poorly understood.

Results

Comparative analysis showed that the structure and evolutionary patterns of the I2 locus vary considerably between potato and tomato. The I2 homologues from different Solanaceae species usually do not have orthologous relationship, due to duplication, deletion and frequent sequence exchanges. At least 154 sequence exchanges were detected among 76 tomato I2 homologues, but sequence exchanges between I2 homologues in potato is less frequent. Previous study showed that I2 homologues in potato were targeted by miR482. However, our data showed that I2 homologues in tomato were targeted by miR6024 rather than miR482. Furthermore, miR6024 triggers phasiRNAs from I2 homologues in tomato. Sequence analysis showed that miR6024 was originated after the divergence of Solanaceae. We hypothesized that miR6024 and miR482 might have facilitated the expansion of the I2 family in Solanaceae species, since they can minimize their potential toxic effects by down-regulating their expression.

Conclusions

The I2 locus represents a most divergent resistance gene cluster in Solanum. Its high divergence was partly due to frequent sequence exchanges between homologues. We propose that the successful expansion of I2 homologues in Solanum was at least partially attributed to miRNA mediated regulation.

Background

Most of the disease resistance genes cloned from plant species encode nucleotide-binding site (NBS) and leucine-rich repeat (LRR) domains. The NBS-LRR encoding genes are often called disease resistance genes (or R-genes), since their main functions are protecting plants from pathogens with only a few exceptions, such as CHS3, a TIR-NBS-LRR-LIM encoding gene, playing role in cold stress and ADR1, a CC-NB-LRR encoding gene, involved in drought tolerance [14]. The R-genes belong to a large gene family, with dozens or hundreds of copies in a genome [57]. The R proteins are composed of a variable N terminus, a conserved central NBS domain and a C-terminus with various number of short LRR motifs [8]. The N terminus usually has a Toll/interleukin-1 receptor (TIR) motif or a coiled-coil (CC) motif [9, 10]. The NBS domains can bind and hydrolyze ATP or GTP, which are involved in signaling leading to pathogen resistance [11]. The LRR domains are most likely to be involved in recognition specificity [12, 13].

R-genes tend to be clustered in plant genomes [6, 9]. Many R-gene sub-families were shown to have been under diversifying selection, particularly at the hyper-variable solvent-exposed residues in the LRR region [14, 15]. Consequently, R-genes represent the most divergent gene families in plant genomes [16]. Besides diversifying selection, frequent sequence exchanges among homologues have also contributed to the high divergence of R-gene families in plants [17, 18]. However, R-genes vary dramatically in frequency of sequence exchanges. Some R-genes, termed Type I R-genes, had frequent sequence exchanges with homologues, resulting in extensive chimeric structures [15, 18]. Frequent sequence exchanges may completely abolished allelic relationship of Type I R-genes from different genotypes. Different genotypes of the same species may harbor a various set of chimeric R-genes (Type I), leading to enormous distinct R-gene sequences in a population [19]. In contrast to the Type I R-genes, some R-genes (termed Type II) do not recombine with paralogues and are highly conserved in different genotypes of a species or closely related species. However, some Type II R-genes were found to be frequently lost in some genotypes, exhibiting presence/absence polymorphism [6, 20]. The mechanism underlying the differentiation of above two distinct evolutionary patterns for R-genes remains unknown.

The resistance gene I2 in tomato (Solanum lycopersicum) encodes resistance against race 2 of the fusarium wilt pathogen Fusarium oxysporum f sp lycopersici [21]. The I2 gene is a member of a gene cluster located on chromosome 11 of tomato. The I2 protein is a typical CC-type R protein. A resistance trait against tomato yellow leaf curve virus (TYLCY) was also mapped in the vicinity of the I2 locus though it remains unclear whether the trait is encoded by a I2 homologue [22, 23]. I2 homologues are found in the syntenic region of potato (S. tuberosum); two resistance genes (R3a and R3b) against potato late blight have been cloned from this region of potato [24, 25]; and at least nine additional R-genes against potato late blight were mapped to this locus, which harbors dozens of I2 homologues [24, 26, 27]. The R3a and R3b from potato exhibit 88% and 80% nucleotide identity with the I2 gene from tomato,respectively. I2 homologues were also found at the corresponding region in pepper (Capsicum annuum), and this locus may represent the most important R-gene locus in Solanaceae [28]. Previous studies of the I2 locus mainly focused on the cloning of functional R-genes, but the structure, evolution and gene content of the entire cluster remain unclear.

R-genes are believed to be physiologically toxic to plant cells and their expression is usually kept at low level when no pathogens exist [29, 30]. Some R-genes are up-regulated when challenged by pathogens [31]. One mechanism to keep R-gene at low level in absence of pathogen is through miRNA mediated gene silencing [32, 33]. Some R-genes could be cleaved by 22-nt miRNAs and successfully produce secondary small interfering RNAs (tasiRNAs) in a phased fashion (phasiRNA). The resistance gene R3a from potato was shown to be cleaved by 22-nt miR482 family and produce phasiRNA [33]. Whether the I2 family from tomato is regulated by miR482 or other miRNAs remains unknown.

In this study, we analyzed the structure and evolution of the I2 locus in tomato, and compared it with that in potato. First, I2 homologues were cloned from several genotypes of cultivated and wild tomato, including genotypes with resistance against TYLCV. Their evolution was investigated in detail through sequence analysis of I2 homologues from multiple genotypes of tomato, potato, pepper and tobacco. The miRNAs targeting I2 homologues in potato and tomato were investigated and compared. The role of miRNAs on R-gene expansion was discussed. The comprehensive study of the structure of the I2 locus may facilitate the cloning and the use of R-genes located in this region, and the evolutionary study may shed light on the mechanism of R-gene expansion in a genome.

Methods

Tomato materials used in this study

Since resistance gene Ty-2 against TYLCV was genetically mapped to the I2 locus, a commercial tomato hybrid cultivar Hongxiaoli with Ty-2 was included in this study. The hybrid was selfed and a F2 population with 736 individuals was generated. A homozygote Ty-2/Ty-2 (flanking markers are homozygous) was obtained from the F2 population to represent the Ty-2 (T) haplotype. Similarly, a susceptible homozygote of ty-2 was identified to represent the susceptible ty-2 (t) haplotype. The F2 population was used to map I2 homologues from the T and t haplotypes. The T haplotype was introgressed into cultivated tomato from S. habrochaites [23, 34]. Four additional S. habrochaites accessions, LA1777, LA1740, LA2158 and LA2860, were also included in this study for I2 homologue analysis.

Inoculation of TYLCV and phenotyping

To phenotype above tomato genotypes on reaction to TYLCV, seedlings of 4–6 leaf stages were inoculated with a TYLCV agrobacterium infectious clone [35]. Inoculated seedlings were kept in a growth chamber under 16 h light/8 h dark cycle for a month before investigating symptom of TYLCV (curve leaf).

PCR amplification, cloning and sequencing of I2homologues

The tomato genomic DNA were extracted from mature leaves of each genotypes using CTAB method, and genomic DNA were used as template to PCR amplify I2 homologues [36]. Nearly full-length (~3.5 kb) fragments of I2 homologues were amplified from aforementioned genotypes of tomato using a pair of degenerated primers (Table 1), which were located at +165 and +3,685 of the resistant gene I2, respectively. PCR amplification was in a 25 μl reaction with 1 unit Fast Pfu Taq (TransGen, Beijing, China), treated for 5 mins at 95°C, followed by 32 cycles at 95°C for 30 s, 52°C for 30 s, and a final extension at 72°C for 1.5 min. PCR products were gel purified using Gel Purification Kit (Generay, Shanghai, China) and ligated into vector pZERO5 using TA cloning kit (TransGen, Beijing, China). Individual colonies were sequenced until no new genes were obtained in the last ten colones. If sequences from different clones have higher than 99.7% nucleotide identity, they were considered to be derived from the same gene. The I2 homologues amplified from an accession were named as accession name followed by “I2” then by a number, such as LA1777-I2-1. I2 homologues from the T and t haplotypes were named as T-I2- and t-I2- followed by a number, respectively.

Table 1 Primer sequences used for this study

Sequence analysis

Four I2 homologues [I2 [GenBank: AF118127], I2C-1 [AF004878], I2C-2 [AF004879], I2C-5 [AF408704]], originated from wild tomato species S. pimpinellifilium were downloaded from GenBank [21, 37, 38]. Two partial genes I2C-3 and I2C-4 were excluded from this study. The sequences of R3a [AY849382], R3b [JF900492.1] and sixteen I2 homologues [AY849383-AY849385, EF638450-EF638453, EF638455, EF638456, EF638458, EF638460-EF638465] from potato genotype SH83-92-488 were retrieved from GenBank and named as SH-I2- followed by a number [24, 25]. The I2 homologue [HQ731036] from S. bulbocastanum was named as SB-I2-1, and the one [HQ731037] from S. stoloniferum as SS-I2-1.

Using BLASTN method, I2 homologues (>2,500 bp) were retrieved from the sequenced tomato genotype Heinz 1706 [39] and the sequenced potato genotype DM [40]. They were named as HZ-I2- and DM-I2- followed by a number, respectively. Sequences were aligned using program Muscle [41] and manually edited in GeneDoc (http://www.nrbsc.org/gfx/genedoc/). Neighbor-joining (NJ) trees (Kimura two-parameter substitution model) with bootstrap values (1,000 replications) were constructed using Mega 5.0 [42]. Gene conversions among homologues were detected using Geneconv [43] with the default settings and confirmed visually. Dot plot analysis between two sequences was performed using program DOTTER [44].

Mapping of the I2homologues in the T and t haplotypes

A F2 population with 736 individuals was used to map the I2 homologues from the T and t haplotypes. First, primers specific to each I2 homologue were designed (Table 1). A primer pair was considered as specific only when the target gene but none of the other genes in the T and t haplotypes was amplified by this primer combination. The gene-specific primers were used as markers to screen the F2 population and their linkage with CAPS markers M-73000 and M-137 at the I2 locus were analyzed.

Identification of miRNAs targeting I2homologues

To detect miRNAs that potentially target I2 homologues, sequences of all I2 homologues were first used as query to screen for matching small RNAs (sRNA) from database SoMART (http://somart.ist.berkeley.edu/) [45]. Identified sRNAs with more than 100 reads in the database were further confirmed in the psRNAtarget web server (http://plantgrn.noble.org/psRNATarget) using following parameters: maximum expectation = 5.0, and target accessibility - allowed maximum energy to unpair the target site (UPE) = 50. All confirmed matching sRNAs were mapped to tomato genome using bowtie [46]. Then, approximately 800 bp flanking sequences of all confirmed matching sRNAs were extracted, and their fold-back structures were predicted using the RNAfold program with default settings [46, 47]. A folding structure was assumed if there was a central loop and a stem (matching region) of at least 18 bp when folding energy no greater than 18 kcal/mol was used. Then, the fold-back structures were evaluated using MirCheck program with strict parameters: ≤4 mismatches, ≤2 bulged or asymmetrically unpaired nucleotides and ≤2 continuous mismatches in the seed regions to meet the accepted criteria for miRNA annotation [48, 49].

The genomes of S. lycopersicum (ITAG version 2.3), S. tuberosum (version 3_2.1.11), Nicotiana benthemiana (version 0.4.4), N. sylvestris (GenBank Assembly ID: GCA_000393655.1), N. tomentosiformis (GenBank Assembly ID: GCA_000390325.1), Capsicum annuum (version 2.0), Mimulus guttatus (version 2.0), Vitis vinifera (Phytozome v9.0: Vvinifera_145), Carica papaya (Phytozome v9.0: Cpapaya_113) and Arabidopsis lyrata (GEO accession: GSE45676) were chosen to analyze if they contain a certain miRNA. First, a mature miRNA sequence was downloaded from miRBase (release 20), then mapped into these genome sequences using bowtie. If a genome sequence had no more than 4 mismatches with a miRNA sequence, its 800 bp flanking sequences were further investigated for stem-loop structure using program RNAfold and evaluated using program MirCheck as described above.

Three miR6024 hairpin sequences (tomato, potato and tobacco) were downloaded from miRbase (release 20). They were used as query sequence in BLASTN search of ten plant genomes, and significant hits (E value cutoff of 1e−10) were retrieved for further analysis.

Mapping of mRNA cleavage sites in I2homologues

Total RNAs were isolated from leaves of tomato Hongxiaoli plantlets using RNAiso Plus (Takara, Dalian, China) and was purified with RNA cleanup Kit (Kangwei, China). To map the miRNA cleavage site in I2 homologues, modified 5′- RACE was performed using GeneRacer Kit (Invirogen, Carlsbad, CA) as described previously [50]. First, approximately 5 μg total RNA was ligated with the GeneRacer RNA Oligo adapter (250 ng). Then, GeneRacer Oligo dT primer or conserved I2 primer was used to synthesize the first strand cDNA [51]. A degenerate primer (Race-P1) was designed in regions conserved in most tomato I2 homologues to pair with GeneRacer™ 5′ primer to amplify PCR products using the synthesized cDNA as template. Then, GeneRacer™ 5′ nest primer and gene nested primers, Race-P2 and Race-P3, were used for the first round nested PCR and second round nested PCR, respectively (Table 1). The PCR products were gel purified using Gel Purification Kit (Generay, Shanghai, China) and ligated into vector pEASY-T1 using TA cloning kit (TransGen, Beijing, China). Positive clones were sent for sequencing.

Cleavage sites of tomato I2 homologues were also investigated using a degradome database using online program SoMART [45]. First, a randomly chosen tomato I2 homologue was input into programs Silcer Detector and dRNA mapper. The output of these two programs was then analyzed by a third program SMART COMPAR.

Results

The I2locus spans more than 5 Mb on chromosome 11 of tomato and potato

Using BLASTN method, 36 I2 homologues were discovered in the sequenced genome of tomato cultivar Heinz 1706. Nineteen of them are located at the I2 locus, in a 5.4 Mb region on chromosome 11. Of the remaining 17 homologues, five are from other chromosomes and 12 are on chromosome 11 but away from the I2 locus, and they were excluded from further analysis (Additional file 1: Table S1). The positions of I2 homologues and important genetic markers around the I2 locus of tomato Heinz 1706 are shown in Figure 1.

Figure 1
figure 1

Dot plot analysis of the I2 locus between tomato and potato. The vertical line represents the I2 locus in tomato cultivar Heinz 1706, and the horizontal line represents the I2 locus in potato cultivar DM. The positions for nine sub-loci (S-1 to S-9) and important markers are shown in both the horizontal and vertical lines. The numbers in brackets after the sub-locus name are the number of I2 homologues at corresponding positions. No large duplications were found, and duplications were mainly limited to the I2 sequences but not their flanking regions.

The sequenced genome of potato DM has 71 I2 homologues. Of them, 53 are located in the syntenic region of the I2 locus of tomato, spanning a region of at least 7.4 Mb (Additional file 1: Table S1). Fifteen I2 homologues are not on chromosome 11, and three I2 homologues are on chromosome 11 but far away from the syntenic region of the I2 locus. Therefore, the I2 locus was most likely expanded before the speciation of Solanum and may have also experienced amplifications after speciation.

Comparative analysis of the I2/R3locus in tomato and potato

To investigate structure variation between tomato and potato, the I2 locus of tomato Heinz 1706 was compared with that of the potato genotype DM. Good co-linearity was found between tomato and potato, except in the centromeric part of the locus. However, the copy number of I2 homologues as well as their positions varies considerably between the two genomes (Figure 1).

Based on the distribution patterns of I2 homologues, this locus can be further divided into 9 sub-loci (sub-clusters). Dot plot analysis indicated that the duplications resulting in the 9 sub-loci were limited to the I2 sequences but not their flanking regions, since sequences flanking these I2 homologues are unrelated (Additional file 2: Figures S1 and S2). I2 homologues are presented at six sub-loci of both tomato and potato genomes, though their copy number may vary dramatically (Figure 1). At the other three sub-loci (sub-loci 3, 4 and 9), I2 homologues are present in one but absent in the other genome, showing presence/absence divergence between the two species.

To investigate the genetic mechanism for such presence/absence divergence at the three sub-loci, their sequences were used to search the pepper genome. Sequences from the syntenic regions of sub-loci 3 and 4 but not 9 were found in pepper genome. Like in potato, the sub-locus 3 in pepper does not have any I2 homologues (Figure 2A), suggesting that the I2 homologue at sub-locus 3 of tomato was most likely gained after speciation. On the other hand, potato might have lost I2 homologues at sub-locus 4 after speciation since pepper has I2 homologues at this sub-locus (Figure 2B).

Figure 2
figure 2

Duplications and deletions of I2 homologues. Shaded region represents homologous sequences between two species, and the percentage refers to nucleotide identity between the two sequences. A. I2 homologue (black box) is present in tomato but absent in potato and pepper, suggesting that the I2 homologue at sub-locus 3 in tomato was gained after speciation. B. I2 homologues (black vertical lines) are present in tomato and pepper but absent in potato, suggesting that the I2 homologue at sub-locus 4 in potato were lost after speciation.

To study the evolution of I2 homologues from different sub-loci, a distance tree was constructed using approximately 2.5 kb sequences of I2 homologues (their 3′ parts were excluded due to various duplications) from tomato (16), potato (39), pepper (9) and N. sylvestris (18) and N. tomentosiformis (11) (Additional file 2: Figure S3). Homologues HZ-I2-5 and DM-I2-22, which contain many unknown nucleotides, were excluded from further analysis. Based on the tree, the I2 homologues from Solanaceae can be classified into four major groups. One group contains I2 homologues from sub-loci 4 and 5 of pepper; one group contains all I2 homologues from sub-loci 1–7 of tomato and potato; and the other two groups contain I2 homologues from sub-loci 8 and 9 of potato and pepper as well as all I2 homologues from Nicotiana. The phylogenetic tree suggests that all I2 homologues from Nicotiana might be derived from sub-loci 8 and 9 in the common ancestor between Solanum and Nicotiana. The I2 homologues from sub-loci 1 and 2 of tomato are highly similar to those of potato, respectively. However, homologues in other sub-loci showed no obvious orthologous relationship between tomato and potato. The lack of orthologous relationship between I2 homologues in tomato and potato suggest that sequence exchanges may have occurred between homologues from different sub-loci after speciation of Solanum.

Sequencing I2homologues from different genotypes

To better understand the diversity and evolution of the I2 locus, I2 homologues were cloned from six tomato genotypes, including the T and t haplotypes (see MM section), and four wild genotypes of S. habrochaites (LA1777, LA1740, LA2158 and LA2860). Inoculation with the infectious TYLCV clone showed that LA1777 is resistant to TYLCV, consistent with previous conclusion [52], while the other three wild genotypes are susceptible. To clone the I2 homologues from these genotypes, a pair of degenerated primers (I2-F/R) was designed from the conserved regions of I2 homologues (Table 1). Their PCR products were cloned and individual colonies were sequenced. A total of 75, 83, 48, 55, 70 and 70 colonies were sequenced, resulting in 8, 11, 13, 6, 13 and 11 distinct sequences from the T and t haplotypes, S. habrochaites accessions LA1777, LA1740, LA2158 and LA2860, respectively.

Mapping of the I2homologues in the T and t haplotypes

The 8 and 11 I2 homologues obtained from the T and t haplotypes were genetically mapped using a F2 segregating population (736 individuals) derived from the hybrid cultivar Hongxiaoli. First, markers M-73000 and M-137, which flank the I2 locus, were used to screen for recombinants. A total of 44 recombinants were obtained from the 736 F2 individuals. Specific primers were successfully designed for six of the 19 I2 homologues from the T and t haplotypes, and all of them were fine mapped to sub-loci 5–7 of the I2 locus (Figure 3).

Figure 3
figure 3

Distance tree of I2 homologues from different genotypes of tomato. The 76 I2 homologues from tomato form five clades, supported by high bootstrap numbers. Members in clades I, III and V are evenly related homologues, and are extensive chimeras. Members within clades II and IV are highly conserved homologues from different genotypes. The first part of the gene name represents accession/cultivar name, and the numbers in bracket (if any) show gene’s position (sub-locus). Genes from the same genotype spread in the tree. Numbers on nodes are bootstrap values, and values <65 are not shown.

Different evolutionary patterns for I2homologues

The 62 nearly full-length I2 homologues obtained above were combined with the 10 I2 homologues (>2.5 Kb) from Heinz 1706 and four homologues from wild species S. pimpinellifilium for further analysis [Genbank: KJ652840 – KJ652901] (Additional file 1: Table S1). A distance tree of the 76 I2 homologues showed five distinguishable clades (Figure 3). Members from clade II and IV are highly conserved, with all pair-wise nucleotide identities >96.9%. Each genotype/haplotype usually has one representative in clades II and IV, and members within each clade are most likely alleles/orthologues. The exceptions are two I2 homologues from LA1777 in Clade IV, and two I2 homologues from LA1777 and LA1740 in clade II, respectively. The presence of two highly similar sequences is most likely due to the heterozygosity of the I2 locus in these genotypes, which are self-incompatible. Only a few sequence exchanges were detected in clades II and IV, and mainly occurred between alleles (Additional file 1: Table S2). Therefore, clades II and IV have evolved independently from other homologues at the I2 locus, with an evolutionary pattern similar to that of Type II R-genes [15].

In contrast, the other three clades in Figure 3 usually have multiple copies from one genotype and the average pair-wise nucleotide identity within a clade varies from 91.3% to 96.6%. Homologues within a clade are often equally related with each other and do not show obvious allelic/orthologous relationship. High sequence diversity and lack of allelic/orthologous relationship might have been attributed to frequent sequence exchanges among homologues [15]. To test above hypothesis, sequence exchanges were analyzed for all I2 homologues in tomato using software Geneconv, and a total of 154 sequence exchanges were detected. The length of the sequence exchange tracts varied from 77 to 3,022 bp, with an average of 533 bp. Of them, 143 sequence exchanges occurred among members within a clade (98 within clade I, 7 within clade II, 8 within clade III, 1 within clade IV and 29 within clade V) (Additional file 1: Table S2). Clade I has homologues from three sub-loci (5, 6 and 7), and sequence exchanges happened between homologues from different loci. The frequent sequence exchanges and extensive chimeric structures suggest that the I2 homologues in clades I and V have an evolutionary pattern of Type I R-genes [15].

Sequence exchanges occasionally occurred between Type I and Type II I2homologues or between different lineages of Type I genes

Eleven sequence exchanges were detected between genes from different clades. A sequence exchange of 162 bp was detected between LA1740-I2-5 from clade I and three members (LA1740-I2-3, LA2158-I2-5, LA2860-I2-5) from clade II (Additional file 1: Table S2 and Additional file 2: Figure S4). The three members from clade II have identical sequences in the exchange tract. It is unlikely that the three genes had independent sequence exchanges in the same region. It is most likely that a 162 bp sequence of LA1740-I2-3 converted gene LA1740-I2-5, since these two genes are from the same genotype (LA1740). Similarly, gene LA1777-I2-1 from clade II converted gene t-I2-2 from clade I. Sequence exchanges might also occur between different lineages of Type I genes, such as between genes HZ-I2-10 from clade I and t-I2-5 from clade III. Interestingly, the directions of above sequence exchanges were unilaterally from Type II genes to Type I genes or between different lineages of Type I genes, but never from Type I genes to Type II genes.

Rare sequence exchanges among I2homologues in potato

Similar analyses (including phylogenetic analysis and sequence exchange) were applied to I2 homologues in potato. A total of 63 I2 homologues (>2.5 Kb, and one gene with many unknown nucleotides excluded) were obtained from Genbank, mainly from the sequenced genome of DM and the R3a haplotype of the diploid S. tuberosum (Additional file 1: Table S1). Three of them are not from the I2 locus and one has too many missing data, and they were excluded from further sequence analysis. A distance tree was constructed for the remaining 59 potato I2 homologues, and seven tomato homologues including I2, I2C-1 and one from each clade in Figure 3 were included for comparison. Consistent with their locations shown in Figure 1, all I2 homologues from tomato are grouped together with homologues from sub-loci 1–7 of potato (clade I in Figure 4). The homologues from clade I in Figure 3 are quite divergent (mostly < 90% nucleotide identity, with an average of 87.5%) but evenly related. Seven sequence exchanges were detected among the 34 I2 homologues from potato in clade I (Additional file 1: Table S2). Such pattern is in striking contrast to that for the I2 homologues from tomato, in which I2 homologues were differentiated into several groups, including both Type I and Type II genes (see above). No obvious differentiation (i.e. no well supported clades in Figure 4) among potato I2 homologues and occasional sequence exchanges suggest that the I2 homologues in sub-loci 1–7 in potato represent an ancient lineage of Type I genes that were originated at least 7 MYA before the divergence of tomato and potato [53, 54].

Figure 4
figure 4

Distance tree of I2 homologues from potato. Seven representative I2 homologues from tomato are included, which are all grouped into clade I. Genes with name “DM-I2-” are from potato cultivar DM; genes with name “SH-I2-” are from S. tuberosum; genes with name “HZ-I2-” are from tomato Heinz 1706. Numbers on nodes are bootstrap values, and values <65 are not shown.

The homologues from sub-loci 8 and 9 of potato form four clades (II, III, IV and V) in Figure 4. Ten and two sequence exchanges were detected among members in clade II and clade V, respectively. Members within these two clades exhibited 91.8% and 95.3% average nucleotide identity, respectively. Sequence exchanges among members of each clade and high diversity of I2 homologues suggest that clades II and V represent two distinct lineages of Type I I2 homologues [15]. The evolutionary pattern of clade III and IV remains unclear because only two sequences were obtained from this study.

Identification of miRNAs for the I2homologues

The resistance gene R3a from potato was shown to be targeted by members of 22-nt miR482 family [33]. Computational analysis showed that I2 homologues from tomato may also be targeted by miR482 family, which has about 380 reads per million from three tissues (leaf, flower and fruit) of tomato (miRbase release 20). To investigate if I2 homologues from tomato are targeted by other miRNAs, the tomato sRNA libraries (http://somart.ist.berkeley.edu) were first searched using 76 I2 homologues. Consequently, a total of 1,439 distinct sRNAs, including the 22-nt sly-miR6024, were found in the sRNA database when five mismatches were allowed. Sly-miR6024 was previously shown to regulate the expression of R-gene Tm-2 in tomato [33]. To investigate if sly-miR6024 also target the I2 homologues, modified RNA ligase-mediated 5′-RACE was performed [50, 51]. Sequencing the PCR products showed that the mRNAs of at least one I2 homologue from cultivar Hongxiaoli were cleaved at the predicted targeting site of miR6024 (Figure 5 and Additional file 2: Figure S5). However, no cleavage product was detected at the predicted target site of miR482. Therefore, the I2 homologues from tomato are targeted by miR6024 and may not be targeted by miR482. The cleavage of miR6024 on I2 homologues was also confirmed in a degradome database of tomato. Fourteen partial I2 cDNAs were identified in the database. One of partial ones starts from the 485th nucleotide of gene T-I2-3, confirming the cleavage function of miR6024 (Additional file 1: Figure S6A).

Figure 5
figure 5

Cleavage sites of miR6024 and two tasiRNAs in gene T-I2-3 . The 22 nt miR6024 has three mismatches with gene T-I2-3. The blocks of 21 nt (marked by underlines) show the tasiRNAs triggered by miR6024, and the 6th-14th tasiRNAs were omitted as indicated by “//”. The arrow in T-I2-3 matching the miR6024 sequence shows the cleavage site of miR6024, evidenced by RACE-PCR and degradome analysis. The arrows in the third (3′S3) and 15th (3′S15) tasiRNAs are most likely the cleavage sites of tasiRNAs.

The targeting site of sly-miR6024 encodes the 206th – 213th amino acids in the I2 protein, partially overlapping with the conserved P loop. This sequence is highly conserved in all I2 homologues in tomato and potato. We hypothesize that miR6024 regulate expression of most, if not all I2 homologues in tomato, which may facilitate its expansion in a genome.

MiR6024 triggers 21-nt phased siRNA from I2homologues

It was shown that 22-nt miRNAs often trigger the biogenesis of secondary phased siRNA [55, 56]. To investigate if miR6024 can trigger phasiRNA, small RNAs from tomato were analyzed using program SoMART [45]. A total of 847 sRNAs were successfully mapped to a representative I2 homologue (T-I2-3). They were predominantly 21-nt in length and with 5′ U residue, consistent with the features of tasiRNA [33, 57]. Most of them were mapped to the downstream of the cleavage site of miR6042 with a phased pattern (Additional file 2: Figure S6B). The structure of these sRNAs and their phasing with the miR6024 cleavage site indicated that they were tasiRNAs triggered by the 22-nt miRNA6024.

RACE-PCR and degradome database were used to investigate the potential regulating effects of tasiRNAs. Three potential cleavage sites of phasiRNA were identified. A degraded mRNA of gene T-I2-6 starting at 579th nucleotide and a degraded mRNA of gene T-I2-3 starting at 540th nucleotide, were obtained using RACE-PCR strategy (Additional file 2: Figure S6C). These two degraded mRNA were most likely the cleavage products of tasiRNAs (3′S5 and 3′S3, respectively) since their cleavage points are located in the middle of a tasiRNA triggered by the miRNA6024 (Figure 5). In addition, a degradome RNA from tomato degradome database was most likely generated by phasiRNA 3′S15 (Figure 5 and Additional file 2: Figure S6D).

MiR6024 was originated after the divergence of Solanaceae

To gain insight into the evolution of miR6024 in plants, four members of the miR6024 family were downloaded from miRbase (release 20) including sly-miR6024 from tomato, nta-miR6024 from tobacco, stu-miR6024-3p and stu-miR6024-5p from potato. In addition, whole genome sequences of ten plant species were chosen to identify miR6024 genes using bioinformatic approach (see MM section). The miR6024 sequences were used to BLAST search the ten genomes, and significant hits were found in eight of the ten plant species except in C. papaya and A. lyrata. However, the flanking sequences of the significant hits in M. guttatus and V. vinifera could not form hairpin structure, and therefore these two species do not have miR6024. All the six Solanaceae species have the miR6024 sequence and its flanking sequences (miR6024 gene) can form a predicted hairpin structure. The miR6024 gene has no similarity with I2 homologues except that the miR6024/miR6024* can match the target site in I2 homologues. Therefore, the miR6024 was not originated from duplication and inversion of I2 sequences [33]. The 22-nt mature miR6024 was confirmed in a pepper sRNA database (Dr. Li, F., unpublished data). However, point mutations were observed in the miR6024* region of N. sylvestris and N. benthemiana, and it remains unclear if these mutations affect the biogenesis of miR6024 in the two species (Figure 6). Since miR6024 is present in distantly related species in Solanaceae, we hypothesize that the miR6024 was originated in the common ancestor of the Solanaceae family.

Figure 6
figure 6

Alignment of the miR6024 genes in Solanaceae. The bottom line is consensus sequence. Underlines in the consensus sequences represent miR6024, while dash lines under the consensus sequence represent miR6024*. Dots mean nucleotides identical to that of consensus sequence. Dash lines in sequences represent deletions. Sequence name HZ represents tomato Heinz 1706, DM for potato cultivar DM, CA for C. annuum, NTA for N. tobaccum, NT for N. tomentosiformis, NB for N. benthemiana and NS for N. sylvestris.

Discussion

Divergence of I2locus among different Solanaceae species

The I2 locus is a hotspot for R-genes, with several qualitative resistance traits (such as I2, Ty-2 and Sm from tomato, R3, R6 and R7 from potato, and L from pepper) and quantitative resistance traits mapped to this locus. A good understanding of the structure and evolution of this locus will facilitate the cloning and efficient use of R-genes from this region.

The I2 locus contains multiple homologues and spans several megabases on the long arm of chromosome 11 in Solanum. The I2 homologues in this region form nine well-separated sub-cluster (loci). The sub-loci and their expansion were caused by duplication of individual I2 homologues, since there was no evidence of large (>10 kb) duplications in this region. The duplication frequency of I2 homologues might have varied considerably since different species show a large variation in I2 copy number in different sub-loci. Our data showed that deletions also accounted for part of the I2 locus divergence between different species (Figure 2).

I2 homologues from the centromeric part of the I2 locus were grouped into one clade in the phylogenetic tree of all I2 homologues from Solanum, while homologues from sub-loci 8 and 9 formed four clades (Figure 4). The five clades are equally related with each other, suggesting that they were differentiated at similar time. After differentiation, the homologues from the centromeric part of the I2 locus duplicated and generated the sub-loci 1–7.

The structure of I2 locus from the Ty-2 haplotype and S. habrochaitesLA1777

The resistance gene Ty-2 against TYLCV in cultivated tomato was introgressed from S. habrochiates accession“B6013” [34]. The recombination at the I2 locus of the Ty-2 haplotype was shown to be highly suppressed [58]. Recombination suppression is possibly caused by inversion (or other dramatic chromosome change) in this region that may prevent pairing between homologous chromosomes. It remains unclear if such chromosome change also exists in any genotypes susceptible to TYLCV. If such susceptible genotypes are identified, they can be crossed with the Ty-2 haplotype and the resistance trait can be fine mapped.

Four accessions of S. habrochiates were included in this study for analysis of I2 homologues. The number of I2 homologues and their sequences seem to vary considerably between different accessions, showing large diversity at this locus. One of the four accessions, accession LA1777, was resistant against strain TYLCV-Cyprus [59] and tolerance to TYLCV strains from Sardinia and Senegal [52, 60, 61]. Different copy numbers of I2 homologues in the T haplotype (8) and accession LA1777 (13), and no highly similar pairs of genes in the two haplotypes (Figure 3), suggest that they have different I2 locus, consistent with a previous report that the TYLCV resistance in LA1777 is not encoded by Ty-2 but by several recessive loci [61].

Sequence exchanges among I2homologues

All I2 homologues from sub-loci 1–7 in potato were evenly related with each other, similar to the feature of Type I R-genes. However, unlike Type I R-genes, these homologues showed infrequent sequence exchanges, and they are not extensive chimeras though occasional sequence exchanges were detected among them. The presence of a large number of divergent but evenly related homologues and infrequent sequence exchanges suggested that the I2 homologues from sub-loci 1–7 of the I2 locus in potato are Type I R-genes. They are mingled with tomato I2 homologues from sub-loci 1–7 in phylogenetic tree (Figure 4), suggesting that this lineage of Type I genes in potato had existed before the divergence of tomato and potato 7 MYA [53, 54].

The I2 homologues in sub-loci 1–7 in tomato, on the other hand, were diverged into five well-supported clades (Figure 3). Genes physically linked are not necessarily closely related. For example, the genes located in sub-locus 5 were grouped into two different clades (Figure 3). Members from clade I in Figure 3 showed extensive chimeric structures resulted from frequent sequence exchanges; they are evenly related with each other, with most pair-wise nucleotide identity of 85.6-99.8%. Above features are typical of Type I R-genes [15]. In contrast, genes within clades II and IV are highly conserved and each genotype has only one representative. Interestingly, rare sequence exchanges were also found between members from different clades. We conclude that the I2 homologues from sub-loci 1–7 of tomato has started to differentiate. Some groups have evolutionary patterns of Type I while others have evolutionary patterns of Type II R-genes. In contrast, the sub-loci 1–7 from potato did not have such differentiation and maintained as one lineage of Type I genes. Such variation on the organization and evolution of the I2 homologues between tomato and potato, which diverged from their last common ancestor approximately 7 MYA [53, 54], showed that the evolutionary pattern of R-genes may change dramatically within a relatively short evolutionary time.

I2 homologues from sub-locus 9 of potato were grouped into two different clades (Figure 4). Furthermore, the homologues from these two clades are interwoven at sub-locus 9. Members of different lineages of the R1 resistance gene are also interwoven in the R1 resistance gene cluster in potato similar to the organization of the R1 resistance gene family in potato [62]. The interweaving organization of members from different clades suggests that physical proximity is not the main factor for the differentiation of Type I and Type II R-genes [63].

MiRNAs targeting I2homologues

R-genes, though critical for the survival of plants, may pose threats to plant fitness if accumulated to a high copy number [29]. To reduce their fitness cost, R-genes are often kept at low level of expression. One mechanism for down-regulating the expression of R-genes is through miRNAs, which were recently shown to cleave R-genes and generate tasiRNAs [33]. In addition to the miR482 family targeting I2 homologues in potato [33], we identified and confirmed that miR6024 targets I2 family in tomato. Search of the tomato degradome confirmed the cleavage function of miR6024 and also suggest the cleavage function of tasiRNAs triggered by miRNAs (see Results section). The miR482 family has high expression level in three tissues of potato (leaf, flower and stolon) (five members, with an average of 35,760 reads per million), while miR6024 family has a much lower expression in tomato (leaf, flower and fruit) (2,540 reads per million). MiR6024 was found in distantly related genus (Nicotiana and Solanum) in Solanaceae, and therefore it was most likely originated in the common ancestor of Solanaceae. MiR482 was shown to be an ancient miRNA [64]. However, no cleavage of I2 homologues by miR482 was confirmed in tomato though the miR482 does exist in tomato genome. Therefore, the collective regulation by miR482, miR6024 as well as tasiRNAs triggered by them may control the expression level of I2 homologues in Solanaceae, and the down-regulation of the I2 homologues by these miRNAs may make the expansion of I2 homologues less costly in fitness. Though silenced by miRNAs, these R-genes can exercise their full resistance function when challenged by pathogens such as viruses and fungi, which can suppress the silencing machinary of host plants [65, 66].

Conclusions

Comparative analysis of the I2 locus in tomato and potato showed that the I2 locus contains a large number of homologues with high divergence. Its evolutionary patterns varied considerably between tomato and potato. The I2 family was targeted by at least two miRNAs, which may play important roles in down-regulation and evolution of this resistance gene family.

Availability of supporting data

Sequence data described this article can be found in the GenBank data libraries under accession numbers KJ652840 - KJ652901.

References

  1. Dangl J, Jones JDG: Plant pathogens and integrated defence responses to infection. Nature. 2001, 411 (6839): 826-833. 10.1038/35081161.

    Article  CAS  PubMed  Google Scholar 

  2. Chini A, Grant JJ, Seki M, Shinozaki K, Loake GJ: Drought tolerance established by enhanced expression of the CC-NBS-LRR gene, ADR1, requires salicylic acid, EDS1 and ABI1. Plant J. 2004, 38 (5): 810-822. 10.1111/j.1365-313X.2004.02086.x.

    Article  CAS  PubMed  Google Scholar 

  3. Holt BF, Belkhadir Y, Dangl JL: Antagonistic control of disease resistance protein stability in the plant immune system. Science. 2005, 309 (5736): 929-932. 10.1126/science.1109977.

    Article  CAS  PubMed  Google Scholar 

  4. Yang H, Shi Y, Liu J, Guo L, Zhang X, Yang S: A mutant CHS3 protein with TIR-NB-LRR-LIM domains modulates growth, cell death and freezing tolerance in a temperature-dependent manner in Arabidopsis. Plant J. 2010, 63 (2): 283-296. 10.1111/j.1365-313X.2010.04241.x.

    Article  CAS  PubMed  Google Scholar 

  5. Tuskan GA, Difazio S, Jansson S, Bohlmann J, Grigoriev I, Hellsten U, Putnam N, Ralph S, Rombauts S, Salamov A, Schein J, Sterck L, Aerts A, Bhalerao RR, Bhalerao RP, Blaudez D, Boerjan W, Brun A, Brunner A, Busov V, Campbell M, Carlson J, Chalot M, Chapman J, Chen GL, Cooper D, Coutinho PM, Couturier J, Covert S, Cronk Q: The genome of black cottonwood, populus trichocarpa. Science. 2006, 313 (5793): 1596-1604. 10.1126/science.1128691.

    Article  CAS  PubMed  Google Scholar 

  6. Luo S, Zhang Y, Hu Q, Chen J, Li K, Lu C, Liu H, Wang W, Kuang H: Dynamic nucleotide-binding site and leucine-rich repeat-encoding genes in the grass family. Plant Physiol. 2012, 159 (1): 197-210. 10.1104/pp.111.192062.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  7. Ling HQ, Zhao S, Liu D, Wang J, Sun H, Zhang C, Fan H, Li D, Dong L, Tao Y, Gao C, Wu HL, Li YW, Cui Y, Guo XS, Zheng SS, Wang B, Yu K, Liang QS, Yang WL, Lou XY, Chen J, Feng MJ, Jian JB, Zhang XF, Luo GB, Jiang Y, Liu JJ, Wang ZB, Sha YH, et al: Draft genome of the wheat A-genome progenitor Triticum urartu. Nature. 2013, 496 (7443): 87-90. 10.1038/nature11997.

    Article  CAS  PubMed  Google Scholar 

  8. Cannon SB, Zhu H, Baumgarten AM, Spangler R, May G, Cook DR, Young ND: Diversity, distribution, and ancient taxonomic relationships within the TIR and non-TIR NBS-LRR resistance gene subfamilies. J Mol Evol. 2002, 54 (4): 548-562. 10.1007/s00239-001-0057-2.

    Article  CAS  PubMed  Google Scholar 

  9. Meyers BC, Kozik A, Griego A, Kuang H, Michelmore RW: Genome-wide analysis of NBS-LRR-encoding genes in arabidopsis. Plant Cell. 2003, 15 (4): 809-834. 10.1105/tpc.009308.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  10. McHale L, Tan X, Koehl P, Michelmore RW: Plant NBS-LRR proteins: adaptable guards. Genome Biol. 2006, 7 (4): 212-10.1186/gb-2006-7-4-212.

    Article  PubMed Central  PubMed  Google Scholar 

  11. Tameling WI, Elzinga SD, Darmin PS, Vossen JH, Takken FL, Haring MA, Cornelissen BJ: The tomato R gene products I-2 and Mi-1 are functional ATP binding proteins with ATPase activity. Plant Cell. 2002, 14 (11): 2929-2939. 10.1105/tpc.005793.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  12. Kobe B, Deisenhofer J: A structural basis of the interactions between leucine-rich repeats and protein ligands. Nature. 1995, 374 (6518): 183-186. 10.1038/374183a0.

    Article  CAS  PubMed  Google Scholar 

  13. Leister RT, Katagiri F: A resistance gene product of the nucleotide binding site - leucine rich repeats class can form a complex with bacterial avirulence proteins in vivo. Plant J. 2000, 22 (4): 345-354. 10.1046/j.1365-313x.2000.00744.x.

    Article  CAS  PubMed  Google Scholar 

  14. Meyers BC, Shen KA, Rohani P, Gaut BS, Michelmore RW: Receptor-like genes in the major resistance locus of lettuce are subject to divergent selection. Plant Cell. 1998, 10 (11): 1833-1846. 10.1105/tpc.10.11.1833.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  15. Kuang H, Woo SS, Meyers BC, Nevo E, Michelmore RW: Multiple genetic processes result in heterogeneous rates of evolution within the major cluster disease resistance genes in lettuce. Plant Cell. 2004, 16 (11): 2870-2894. 10.1105/tpc.104.025502.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  16. Meyers BC, Kaushik S, Nandety RS: Evolving disease resistance genes. Curr Opin Plant Biol. 2005, 8 (2): 129-134. 10.1016/j.pbi.2005.01.002.

    Article  CAS  PubMed  Google Scholar 

  17. Mondragon-Palomino M, Gaut BS: Gene conversion and the evolution of three leucine-rich repeat gene families in Arabidopsis thaliana. Mol Biol Evol. 2005, 22 (12): 2444-2456. 10.1093/molbev/msi241.

    Article  CAS  PubMed  Google Scholar 

  18. Luo S, Peng J, Li K, Wang M, Kuang H: Contrasting evolutionary patterns of the Rp1 resistance gene family in different species of poaceae. Mol Biol Evol. 2011, 28 (1): 313-325. 10.1093/molbev/msq216.

    Article  CAS  PubMed  Google Scholar 

  19. Kuang H, Ochoa OE, Nevo E, Michelmore RW: The disease resistance gene Dm3 is infrequent in natural populations of Lactuca serriola due to deletions and frequent gene conversions at the RGC2 locus. Plant J. 2006, 47 (1): 38-48. 10.1111/j.1365-313X.2006.02755.x.

    Article  CAS  PubMed  Google Scholar 

  20. Shen J, Araki H, Chen L, Chen JQ, Tian D: Unique evolutionary mechanism in r-genes under the presence/absence polymorphism in Arabidopsis thaliana. Genetics. 2006, 172 (2): 1243-1250.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  21. Simons G, Groenendijk J, Wijbrandi J, Reijans M, Groenen J, Diergaarde P, Van der Lee T, Bleeker M, Onstenk J, de Both M, Haring M, Mes J, Cornelissen B, Zabeau M, Vos P: Dissection of the fusarium I2 gene cluster in tomato reveals six homologs and one active gene copy. Plant Cell. 1998, 10 (6): 1055-1068. 10.1105/tpc.10.6.1055.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  22. Stall RE, Walter JM: Selection and inheritance of resistance in Tomato to isolates of races 1 and 2 of the Fusarium wilt organism. Phytopathology. 1965, 55 (11): 1213-1215.

    Google Scholar 

  23. Ji Y, Scott JW, Schuster DJ: Toward fine mapping of the tomato yellow leaf curl virus resistance gene Ty-2 on chomosome 11 of tomato. HortSci. 2009, 44 (4): 614-618.

    Google Scholar 

  24. Huang S, van der Vossen EA, Kuang H, Vleeshouwers VG, Zhang N, Borm TJ, van Eck HJ, Baker B, Jacobsen E, Visser RG: Comparative genomics enabled the isolation of the R3a late blight resistance gene in potato. Plant J. 2005, 42 (2): 251-261. 10.1111/j.1365-313X.2005.02365.x.

    Article  CAS  PubMed  Google Scholar 

  25. Li G, Huang S, Guo X, Li Y, Yang Y, Guo Z, Kuang H, Rietman H, Bergervoet M, Vleeshouwers VG, van der Vossen EA, Qu D, Visser RG, Jacobsen E, Vossen JH: Cloning and characterization of R3b; members of the R3 superfamily of late blight resistance genes show sequence and functional divergence. Mol Plant Microbe Interact. 2011, 24 (10): 1132-1142. 10.1094/MPMI-11-10-0276.

    Article  CAS  PubMed  Google Scholar 

  26. EI-Kharbotly A, Leonards-Schippers C, Huigen DJ, Jacobsen E, Pereira A, Stiekema WJ, Salamini F, Gebhardt C: Segregation analysis and RFLP mapping of the R1 and R3 alleles conferring race-specific resistance to Phytophthora infestans in progeny of dihaploid potato parents. Mol Gen Genet. 1994, 242 (6): 749-754.

    Google Scholar 

  27. El-Kharbotly A, Palomino-Sanchez C, Salamini F, Jacobsen E, Gebhardt C: R6 and R7 alleles of potato conferring race-specific resistance to Phytophthora infestans (Mont.) de Bary identified genetic loci clustering with the R3 locus on chromosome XI. Theor Appl Genet. 1996, 92 (7): 880-884. 10.1007/BF00221901.

    Article  CAS  PubMed  Google Scholar 

  28. Grube RC, Radwanski ER, Jahn M: Comparative genetics of disease resistance within the solanaceae. Genetics. 2000, 155 (2): 873-887.

    CAS  PubMed Central  PubMed  Google Scholar 

  29. Tian D, Traw MB, Chen JQ, Kreitman M, Bergelson J: Fitness costs of R-gene-mediated resistance in Arabidopsis thaliana. Nature. 2003, 423 (6935): 74-77. 10.1038/nature01588.

    Article  CAS  PubMed  Google Scholar 

  30. Heidel AJ, Clarke JD, Antonovics J, Dong X: Fitness costs of mutations affecting the systemic acquired resistance pathway in Arabidopsis thaliana. Genetics. 2004, 168 (4): 2197-2206. 10.1534/genetics.104.032193.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  31. Wang ZX, Yamanouchi U, Katayose Y, Sasaki T, Yano M: Expression of the Pib rice-blast-resistance gene family is up-regulated by environmental conditions favouring infection and by chemical signals that trigger secondary plant defences. Plant Mol Biol. 2001, 47 (5): 653-661. 10.1023/A:1012457113700.

    Article  CAS  PubMed  Google Scholar 

  32. Zhai J, Jeong DH, De Paoli E, Park S, Rosen BD, Li Y, Gonzalez AJ, Yan Z, Kitto SL, Grusak MA, Jachson SA, Stacey G, Cook DR, Green PJ, Sherrier DJ, Meyers BC: MicroRNAs as master regulators of the plant NB-LRR defense gene family via the production of phased, trans-acting siRNAs. Genes Dev. 2011, 25 (23): 2540-2553. 10.1101/gad.177527.111.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  33. Li F, Pignatta D, Bendix C, Brunkard JO, Cohn MM, Tung J, Sun H, Kumar P, Baker B: MicroRNA regulation of plant innate immune receptors. Proc Natl Acad Sci U S A. 2012, 109 (5): 1790-1795. 10.1073/pnas.1118282109.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  34. Kalloo MKB: Transfer of tomato leaf curl virus resistance from Lycopersicon hirsutum f. glabratum to L. esculentum. Plant Breed. 1990, 105 (2): 156-159. 10.1111/j.1439-0523.1990.tb00469.x.

    Article  Google Scholar 

  35. Xie Y, Jiang T, Zhou X: Agroinoculation Shows Tobacco leaf curl Yunnan virus is a Monopartite Begomovirus. Eur J Plant Pathol. 2006, 115 (4): 369-375. 10.1007/s10658-006-9021-8.

    Article  Google Scholar 

  36. Murray MG, Thompson WF: Rapid isolation of high molecular weight plant DNA. Nucleic Acids Res. 1980, 8 (19): 4321-4325. 10.1093/nar/8.19.4321.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  37. Ori N, Eshed Y, Paran I, Presting G, Aviv D, Tanksley S, Fluhr DZAR: The I2C family from the wilt disease resistance locus i2 belongs to the nucleotide binding, leucine-rich repeat superfamily of plant resistance genes. Plant Cell. 1997, 9: 521-532. 10.1105/tpc.9.4.521.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  38. Sela-Buurlage MB, Budai-Hadrian O, Pan Q, Carmel-Goren L, Vunsch R, Zamir D, Fluhr R: Genome-wide dissection of Fusarium resistance in tomato reveals multiple complex loci. Mol Genet Genomics. 2001, 265 (6): 1104-1111. 10.1007/s004380100509.

    Article  CAS  PubMed  Google Scholar 

  39. Tomato Genome C: The tomato genome sequence provides insights into fleshy fruit evolution. Nature. 2012, 485 (7400): 635-641. 10.1038/nature11119.

    Article  Google Scholar 

  40. Potato Genome Sequencing Consortium: Genome sequence and analysis of the tuber crop potato. Nature. 2011, 475 (7355): 189-195. 10.1038/nature10158.

    Article  Google Scholar 

  41. Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32 (5): 1792-1797. 10.1093/nar/gkh340.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  42. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28 (10): 2731-2739. 10.1093/molbev/msr121.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  43. Sawyer S: Statistical tests for detecting gene conversion. Mol Biol Evol. 1989, 6 (5): 526-538.

    CAS  PubMed  Google Scholar 

  44. Sonnhammer EL, Durbin R: A dot-matrix program with dynamic threshold control suited for genomic DNA and protein sequence analysis. Gene. 1995, 167 (1–2): GC1-10.

    CAS  PubMed  Google Scholar 

  45. Li F, Orban R, Baker B: SoMART: a web server for plant miRNA, tasiRNA and target gene analysis. Plant J. 2012, 70 (5): 891-901. 10.1111/j.1365-313X.2012.04922.x.

    Article  CAS  PubMed  Google Scholar 

  46. Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10 (3): R25-10.1186/gb-2009-10-3-r25.

    Article  PubMed Central  PubMed  Google Scholar 

  47. Hofacker IL: Vienna RNA secondary structure server. Nucleic Acids Res. 2003, 31 (13): 3429-3431. 10.1093/nar/gkg599.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  48. Rhoades MW, Reinhart BJ, Lim LP, Burge CB, Bartel B, Bartel DP: Prediction of plant MicroRNA targets. Cell. 2002, 110 (4): 513-520. 10.1016/S0092-8674(02)00863-2.

    Article  CAS  PubMed  Google Scholar 

  49. Meyers BC, Axtell MJ, Bartel B, Bartel DP, Baulcombe D, Bowman JL, Cao X, Carrington JC, Chen X, Green PJ, Jones SG, Jacobsen SE, Mallory AC, Martienssen RA, Poethig RS, Qi YJ, Vaucheret H, Voinnet O, Watanabe Y, Weigel D, Zhu JK: Criteria for annotation of plant MicroRNAs. Plant Cell. 2008, 20 (12): 3186-3190. 10.1105/tpc.108.064311.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  50. Llave C, Xie Z, Kasschau KD, Carrington JC: Cleavage of Scarecrow-like mRNATargets directed by a class of Arabidopsis miRNA. Science. 2002, 297 (5589): 2053-2056. 10.1126/science.1076311.

    Article  CAS  PubMed  Google Scholar 

  51. Chi M, Bhagwat B, Lane WD, Tang G, Su Y, Sun R, Oomah BD, Wiersma PA, Xiang Y: Reduced polyphenol oxidase gene expression and enzymatic browning in potato (Solanum tuberosum L.) with artificial microRNAs. BMC Plant Biol. 2014, 14: 62-10.1186/1471-2229-14-62.

    Article  PubMed Central  PubMed  Google Scholar 

  52. Zakay Y, Navot N, Zeidan M, Kedar N, Rabinowitch H, Czosnek HZ, Zamir D: Screening Lycopersicon accessions for resistance to tomato yellow leaf curl virus: presence of viral DNA and symptom development. Plant Dis. 1991, 75 (3): 279-281. 10.1094/PD-75-0279.

    Article  Google Scholar 

  53. Chaw SM, Zharkikh A, Sung HM, Lau TC, Li WH: Molecular phylogeny of extant gymnosperms and seed plant evolution: analysis of nuclear 18S rRNA sequences. Mol Biol Evol. 1997, 14 (1): 56-68. 10.1093/oxfordjournals.molbev.a025702.

    Article  CAS  PubMed  Google Scholar 

  54. Wikstrom NSV, Chase MW: Evolution of the angiosperms: calibrating the family tree. Proc Biol Sci. 2001, 268 (1482): 2211-2220. 10.1098/rspb.2001.1782.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  55. Chen HM, Chen LT, Patel K, Li YH, Baulcombe DC, Wu SH: 22-Nucleotide RNAs trigger secondary siRNA biogenesis in plants. Proc Natl Acad Sci U S A. 2010, 107 (34): 15269-15274. 10.1073/pnas.1001738107.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  56. Wang Y, Itaya A, Zhong X, Wu Y, Zhang J, van der Knaap E, Olmstead R, Qi Y, Ding B: Function and evolution of a MicroRNA that regulates a Ca2+-ATPase and triggers the formation of phased small interfering rnas in tomato reproductive growth. Plant Cell. 2011, 23 (9): 3185-3203. 10.1105/tpc.111.088013.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  57. Allen E, Xie Z, Gustafson AM, Carrington JC: microRNA-Directed phasing during trans-acting siRNA biogenesis in plants. Cell. 2005, 121 (2): 207-221. 10.1016/j.cell.2005.04.004.

    Article  CAS  PubMed  Google Scholar 

  58. Xiaohui Yang MC, Hutton SF, Scott JW, Yanmei Guo XW, Md Harunur R, Dora S, de Jong RGFV H, Yuling B, Yongchen D: Fine mapping of the tomato yellow leaf curl virus resistance gene Ty-2 on chromosome 11 of tomato. Mol Breed. 2014, doi10.1007/s11032-014-0072-9

    Google Scholar 

  59. Ioannou N: Yellow leaf curl and other virus diseases of tomato in Cyprus. Plant Pathol. 1985, 34 (3): 428-434. 10.1111/j.1365-3059.1985.tb01383.x.

    Article  Google Scholar 

  60. Fargette DLM, Harrison BD: Serological studies on the accumulation and localisation of three tomato leaf curl geminiviruses in resistant and susceptible Lycopersicon species and tomato cultivars. Ann Appl Biol. 1996, 128 (2): 317-328. 10.1111/j.1744-7348.1996.tb07326.x.

    Article  Google Scholar 

  61. Vidavsky F, Czosnek H: Tomato breeding lines resistant and tolerant to tomato yellow leaf curl virus issued from Lycopersicon hirsutum. Phytopathology. 1998, 88 (9): 910-914. 10.1094/PHYTO.1998.88.9.910.

    Article  CAS  PubMed  Google Scholar 

  62. Kuang H, Wei F, Marano MR, Wirtz U, Wang X, Liu J, Shum WP, Zaborsky J, Tallon LJ, Rensink W, Lobst S, Zhang PF, Torngvist CE, Tek A, Bamberg J, Helgeson J, Fry W, You F, Luo MC, Jiang JM, Buell R, Baker B: The R1 resistance gene cluster contains three groups of independently evolving, type I R1 homologues and shows substantial structural variation among haplotypes of Solanum demissum. Plant J. 2005, 44 (1): 37-51. 10.1111/j.1365-313X.2005.02506.x.

    Article  CAS  PubMed  Google Scholar 

  63. Kuang H, Caldwell KS, Meyers BC, Michelmore RW: Frequent sequence exchanges between homologs of RPP8 in Arabidopsis are not necessarily associated with genomic proximity. Plant J. 2008, 54 (1): 69-80. 10.1111/j.1365-313X.2008.03408.x.

    Article  CAS  PubMed  Google Scholar 

  64. Shivaprasad PV, Chen HM, Patel K, Bond DM, Santos BA, Baulcombe DC: A microRNA superfamily regulates nucleotide binding site-leucine-rich repeats and other mRNAs. Plant Cell. 2012, 24 (3): 859-874. 10.1105/tpc.111.095380.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  65. Li F, Ding SW: Virus counterdefense: diverse strategies for evading the RNA-silencing immunity. Annu Rev Microbiol. 2006, 60: 503-531. 10.1146/annurev.micro.60.080805.142205.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  66. Weiberg A, Wang M, Lin FM, Zhao H, Zhang Z, Kaloshian I, Huang HD, Jin H: Fungal small RNAs suppress plant immunity by hijacking host RNA interference pathways. Science. 2013, 342 (6154): 118-123. 10.1126/science.1239705.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

Download references

Acknowledgments

This research was supported by National Natural Science Foundation of China [grant no. 31272030, 31300299 and 31221062].

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Jiongjiong Chen.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

CW carried out the molecular and genetic studies, participated in analysis and interpretation of the data and drafted the manuscript. FL helped in the miRNA study and revised the manuscript. HK helped in the evolution study. JC designed the study and revised the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12864_2014_6427_MOESM1_ESM.xls

Additional file 1: Table S1: I2 homolgoues used in this study. Table S1-1. I2 homologues indentified in the sequenced tomato genome (version 2.40). Table S1-2. I2 homologues identified in the sequenced potato genome (version 3_2.1.11). Table S1-3. I2 homologues from Solanum used in this study. Table S2. Sequence exchanges between I2 homologues. Table S2-1. Sequence exchanges detected between tomato I2 homologues from different clades in Figure 3. Table S2-2. Sequence exchanges detected between potato I2 homologues from different clades in Figure 4. (XLS 80 KB)

12864_2014_6427_MOESM2_ESM.docx

Additional file 2: Figure S1: Dot plot analysis of the I2 locus in tomato. Figure S2. Dot plot analysis of the I2 locus in potato. Figure S3. Distance tree of I2 homologues from tomato, potato, pepper and tobacco. Figure S4. A sequence exchange tract of 162 bp between a Type I and a Type II I2 homologues. Figure S5. The cleavage site of miR6024 obtained through sequencing RACE-PCR products. Figure S6. A. Alignment of T-I2-3 sequence with a cDNA sequence (03330341_71_2) found from a tomato degradome database. B. The abundance of I2-derived sRNAs. C. Cleavage sites of two tasiRNAs triggered by miR6024. The bottom lines of each panel represent 5′ sequences obtained using RACE-PCR. The line in gene T-I2-3 shows the 3rd tasiRNA, and the line in gene T-I2-6 shows the 5th tasiRNA. D. A degraded RNA retrieved from a degradome database. (DOCX 3 MB)

Authors’ original submitted files for images

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wei, C., Kuang, H., Li, F. et al. The I2 resistance gene homologues in Solanum have complex evolutionary patterns and are targeted by miRNAs. BMC Genomics 15, 743 (2014). https://doi.org/10.1186/1471-2164-15-743

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-15-743

Keywords