BACKGROUND
the potato, along with the tomato and pepper, are members of solanaceae family, with the potato being one of the world’s most important food crops  <cit> .
yet the potato is susceptible to a wide range of pathogens, with potato late blight being the most devastating.
potato late blight is caused by phytophthora infestans  <cit> , which greatly effects potato cultivation.
annual potato production is currently over  <dig> million tons and potential production could exceed  <dig> million if potato late blight could be properly controlled  <cit> .
additionally, late blight disease is also responsible for the european potato famine in the 19th century which almost completely destroyed potato crops and led to the starvation of millions of people  <cit> .
thus, the control of such plant diseases is of fundamental importance.
plants have evolved sophisticated systems to recognize pathogenic proteins.
these recognition proteins, usually referred to as resistance  genes  <cit> , work in a highly specific manner according to the gene-for-gene interaction model, in which one plant r gene responds to a pathogen carrying a particular avirulence  gene  <cit> .
this places the host r gene and pathogen avr gene in a co-evolutionary model in which the avr gene evolves different genotypes to avoid recognition while the r genes evolve to recognize the avr genes.
studies have shown that mutations and deletions occur at high frequencies in avr genes, as a response of high rates also noted in r genes  <cit> .
during the past years, many plant r genes have been identified.
most cloned r genes belong to a large gene family and encode proteins with nucleotide-binding sites and leucine-rich repeat  domains  <cit> .
in potatoes, nbs-lrr genes such as r <dig>  rb, r <dig>  r <dig>  rpi-bib <dig>  and rpi-vnt <dig>  have been associated with late blight resistance .
among these genes, rpi-blb <dig> is of particular interest due to its homologs conferring resistance in  <dig> different species and against at least  <dig> different pathogens.
the first member of the rpi-blb <dig> family found to confer resistance was mi- <dig>   <cit> , which provides root knot nematode resistance in tomatoes, with this gene also found to confer aphid resistance  <cit> .
later, mi- <dig>  was found to display resistance to sweet potato whitefly  <cit> , making it the first and the only plant r gene that confers resistance against three distantly related pathogens.
in  <dig>  the late blight resistance gene rpi-blb <dig> identified in potatoes was found to be a mi- <dig>  homolog of mi- <dig>  .
the third r gene cloned from this family, cami, is from hot peppers and has been found to confer root knot nematode resistance, with a 99 % amino acid sequence identity shared with mi- <dig>   <cit> .
three genes, mi- <dig> , rpi-blb <dig>  and cami, derived from different solanaceae organisms confer resistance against various highly divergent pathogens and have originated from the same gene family, making a common ancestor highly possible.
what is it that makes this gene family so special and are there more gene families like this?
to answer these questions, we studied the evolutionary history of this gene family.
in brief, all homologs of these three resistance genes are referred to as rpi-blb <dig> homologs hereafter.
first, homologs of this gene family across the tomato, potato, and pepper genomes were identified, with a varied copy number  and a great nucleotide diversity of this gene family was found within species and between species.
next, the potato lrr regions of this gene family were examined and  <dig> homologs were cloned from  <dig> potato accessions.
a variation in copy number was common between the potato accessions and a high level of nucleotide divergence was found both between accessions and within accessions, suggesting rapid intraspecies and interspecies evolution in this gene family.
meanwhile, strong positive selection and frequent sequence exchanges were also identified in the lrr regions.
overall, our findings indicate the presence of a fast evolving r gene locus with a dramatic variation in copy number, high genetic abundance, and a strong diversifying selection.
the above findings provide further insight enabling the identification of more novel r genes from this special family.
methods
<dig>  genomic source and homolog identification
three fully sequenced genomes were employed in this study, including genome of pepper , genome of tomato  and genome of potato .
the re-annotated nbs-lrr gene sets of potato and tomato were downloaded from previous studies  <cit> .
to identify the homologs of rpi-blb2/mi- <dig> /cami, the coding sequence  of these  <dig> genes were used as query and blast against the whole-genome cdss of the  <dig> genomes and  <dig> re-annotated gene sets.
the e-value was set at 1e- <dig> and other parameters used default settings.
hits with > 60 % coverage and > 60 % identity were deemed homologs.
<dig>  cloning of lrr domain of this gene family from potato accessions
degenerate primers were designed on conserved sites at the edges of the lrr regions in order to gain comprehensive lrr domain sequence information .
pcr was performed in  <dig> potato accessions/genotypes, including  <dig> wild accessions  and  <dig> cultivated accessions .
pcr products were cloned into a pgem-t easy vector and ~ <dig> colonies of each cultivar were then sequenced by abi3100a automated sequencer until no new homolog sequence could be identified.
sequences for all cloned lrr sequences have been deposited in the genbank under accession number from kr <dig> to kr <dig>
<dig>  sequence alignment and phylogenetic analysis
the cdss or cloned lrr sequences were first translated into amino acid sequence and aligned by muscle implanted in mega <dig>  <cit>  and then back translated into nucleotide sequences.
based on the alignments, a maximum likelihood tree was constructed.
first, jmodeltest  <dig> .
<dig>  <cit>  was used to test the best fitted nucleotide substitution model for tree construction, then this model was used by phyml  <dig>   <cit>  with  <dig> bootstrap replicates for maximum likelihood tree construction.
the nonsynonymous  and synonymous  nucleotide substitutions were calculated by dnasp version <dig>  based on the nei-gojobori method  <cit> .
nucleotide diversity  and divergence  were estimated by π and dxy with jukes and cantor correction  <cit>  using dnasp version <dig> .
gene conversion events were analyzed by genecov <dig>  .
RESULTS
<dig>  characterization of rpi-blb <dig> homologs in three species
to construct the evolutionary relationship of the three nbs-lrr genes mi- <dig> , rpi-bib <dig>  and cami, their homologs were first identified in genomes of the  <dig> species .
although the genome size of pepper  was significantly larger than tomato  and potato , the annotated gene number of these  <dig> species was on the same level .
it is reported that the genome expansion of the hot pepper is attributed to accumulation of repetitive sequences  <cit> .
more nbs-encoding genes have been detected in the pepper reference genome than in the tomato or potato reference genomes .
since most nbs genes are duplicated genes, it is difficult to assembly and annotate these genes, and studies have showed that a certain number of nbs genes were missed in potato and tomato reference genomes  <cit> .
jupe et al.
and andolfo et al.
have identified  <dig> and  <dig> nbs genes by renseq in potato and tomato, respectively  <cit> .
as a result, we identified  <dig>   <dig> and  <dig> homologs in the genomes of potato, tomato and pepper, respectively.
after removing incomplete sequences,  <dig>   <dig> and  <dig> homologs from potato, tomato and pepper were used for further analysis.table  <dig> number of rpi-bib <dig> homologs in tomato, pepper and potato genomes
ka/ks
*these data are from re-annotated nbs gene sets by previous studies  <cit>
the nucleotide diversity  and ka/ks were calculated for the homologs from each genome.
the divergence varied a lot between genomes.
the homologs from potato genome were found to be the most conservative with the lowest nucleotide diversity , while the homologs from tomato genome were the most divergent .
such a huge difference was not found when examining ka/ks ratios which ranged from  <dig>  to  <dig>  among the different genomes.
the nucleotide diversity and ka/ks were also determined for the nbs domain and lrr domain separately.
across all of the three genomes, the lrr domain exhibited a higher nucleotide diversity and a higher ka/ks than the nbs domain .
most of the identified homologs were located in tandem on chromosome  <dig> in each genome.
tandem repeat is a major mechanism in nbs-lrr expansion, with previous studies showing that nbs-lrr gene homologues often cluster in tandem on the chromosome  <cit> .
to explore the relation of genomic position among these homologs, a comparative position map of all of the genes on chromosome  <dig> was constructed , extensive conservation of the gene order between these genomes was found among these  <dig> species, suggesting a good synteny in this chromosome.
these homologs all clustered on the same end of the chromosome, further reflecting that tandem repeat events took place in the expansion of this gene family.
additionally, all the homologs on chromosome  <dig> in potato and tomato have the highest similarity with the same gene in pepper, indicating all the homologs in these species share the same ancestor, and duplication events happened several times independently in each species.fig.
<dig> positional relationships among rpi-blb2/mi- <dig> /cami homologs from potato, pepper and tomato reference genomes.
genes are represented by blue arrows and genes with highest similarity between species are linked by green lines
<dig>  phylogenetic analysis of rpi-bib <dig> homologs from tomato, pepper and potato
to further explore the evolutionary history of this gene family, a phylogenetic tree was constructed based on the cloned functional genes and their homologs within the  <dig> reference genomes.
as shown in fig.
<dig>  genes from the same species mostly clustered together on the tree.
this kind of topology further indicates that the expansion of this gene family took place separately in these  <dig> species after they split from their common ancestor.
this tree also shows a closer homology between the tomato and potato than with pepper, coinciding with their genetic relationships on the species level.fig.
<dig> phylogenetic tree derived from rpi-blb2/mi- <dig> /cami homologs.
the maximum-likelihood phylogenetic tree is constructed by phyml  <dig>  based on the cds sequences of each homolog
subfamilies were divided with an identity > 80 % and bootstrap > 60 %, with one major subfamily comprised of  <dig> genes, two families with  <dig> and  <dig> genes and  <dig> minor subfamilies, each comprised of 2– <dig> genes, identified .
all of the  <dig> known r-genes were found in subfamily  <dig>  the major subfamily, which includes  <dig> homologs from tomato,  <dig> genes from potato and  <dig> gene from pepper.
then nucleotide divergence and ka/ks were calculated within each subfamily.
nucleotide divergence varies from  <dig>  to  <dig>  and ka/ks varies from  <dig>  to  <dig> .
to further compare evolutionary patterns among homologs in the largest subfamily, the divergence of genes within tomato, potato and between them was calculated.
the divergence of paralogs within tomato and potato was  <dig>  and  <dig> , while the divergence of homologs between tomato and potato was a little higher than within species .
this shows that in subfamily  <dig>  homologs from tomato are more conserved while homologs from potato are more divergent.table  <dig> nucleotide diversity, ka/ks and the number of gene conversions in each gene subfamily
ka/ks
<dig>  nucleotide variation in lrr regions among potato accessions
typically, in nbs-lrr genes, the lrr domain is responsible for recognition of the avr gene and is more divergent than the nbs domain  <cit> .
to study the evolutionary history of the lrr domain in this gene family,  <dig> homologs from  <dig> wild and  <dig> cultivated potato accessions were cloned  averaging  <dig>  homologs in each accession, ranging from  <dig> to  <dig>  reflecting frequent copy number variations in these accessions.
interestingly, more copy number was detected in cultivated potatoes, with  <dig> accessions having  <dig> copies and  <dig> having  <dig> copies.
it might be because the cultivated potatoes are all tetraploid species, whereas most of the wild species are diploid.table  <dig> statistics of the lrr regions cloned in each potato accession
ka/ks
pp10
ka/ks
to estimate the relative evolutionary rate of the lrr regions, the nucleotide divergence and ka/ks were calculated within each accession.
the average nucleotide divergence of the  <dig> homologs in all of the  <dig> accessions was  <dig> .
the average nucleotide divergence for the wild-type and cultivated accessions were  <dig>  and  <dig> , suggesting that the lrr region in this resistance gene locus is polymorphic.
moreover, the divergence in cultivated accession was significantly higher than that in wild accessions , which implies that the lrr domains in cultivated accession might have been diversified during the process of cultivation.
additionally, a relatively high ka/ks ratio was found both in the wild and cultivated accessions, with  <dig> accessions >  <dig>   <dig> accessions >  <dig>  and the other two accessions ≈  <dig> .
on average, the ka/ks value of the lrr domain was higher than the whole genes from tomato, pepper and potato, indicating that the lrr domains in particular have been positively selected.
when plotting nucleotide divergence against ka/ks ratio, a significant negative relationship was found , suggesting that the nucleotide divergence reduces as the pressure of positive selection increases.
lastly, frameshift mutations  were annotated in each cloned lrr regions and ~30 % of the lrr regions were characterized with frequent frameshift mutations.
the most frequent mutations were detected in the wild species s.demissum, which is a hexaploid, while the other wild accessions or genotypes had fewer frameshift mutations .
in contrast, the mutations detected in the cultivated species distributed relatively evenly.
collectively, these pseuodo-lrr domains might provide a resource or raw material for producing more new lrrs.
<dig>  lrr regions are driven by diversified selection
to investigate the genetic relations and their evolutionary history between lrr homologs from different accessions, a phylogenetic tree was constructed based on cloned lrr regions, together with lrr sequences from potato genome and the functional genes .
the topological structure of the tree showed the cloned lrr regions to be scattered, suggesting that the lrr regions had a high degree of polymorphisms and varied between accessions and within accessions.
moreover, the lrr regions were not distinctly divided into two clear groups, wild and cultivated clades , indicating a frequent introgression during the cultivation of potato.fig.
<dig> phylogenetic tree derived from the  <dig> lrr regions cloned in this study and  <dig> lrr domains from previous studies.
the  <dig> lrr domains were cloned from  <dig> different potato accessions/genotypes,  <dig> of the  <dig> lrr domains were drawn from the potato genomes and the rest  <dig> from previous cloned genes in potato, pepper, and tomato
however, the lrr regions from wild and cultivated accessions were not totally confused or mixed together in the tree, with most small subclades being wild-specific or cultivated-specific.
additionally, in the wild-specific or cultivated-specific subclades, the lrr regions cloned from an accession or genotype were not always clustered, illustrating the complexity of orthologous relationship among the lrr regions.
it also suggests that frequent sequence exchanges may have taken place among accessions of cultivated potatoes or wild-type species.
subfamilies were also divided based on an identity > 80 % and bootstrap > 60 %.
in total,  <dig> subfamilies were identified, each with 2– <dig> members.
the resistance lrrs of tomato and pepper  were clustered in subfamily  <dig>  whereas the lrr of rpi-blb <dig> in potato was isolated.
except the resistance-specific subfamily  <dig> and subfamily  <dig> which contained only genome sequences,  <dig> of  <dig> subfamilies were wild-specific and  <dig> subfamilies are cultivated-specific.
to further investigate evolutionary patterns of subfamilies, the nucleotide diversity, ka/ks and gene conversions were calculated and identified in each subfamily .
to obtain reliable results, subfamilies with few members  were excluded from further analysis.
<dig> subfamilies were kept, including  <dig> subfamilies with ≥  <dig> members.
the nucleotide divergence of subfamily  <dig>  which consists of two functional lrrs, was  <dig>  with a ka/ks >  <dig>  an indicator of positive selection.
additionally, the average number of members in  <dig> wild-specific subfamilies was  <dig> , much higher than the number in cultivated-specific subfamilies .
when comparing wild- and cultivated-specific subfamilies, the wild-specific had higher nucleotide diversity , a higher ka/ks  and an elevated rate of gene conversions events , suggesting that more abundant genetic resources could be re-used in wild species.
moreover, four complex subfamilies with 7– <dig> members had exceptionally high levels of polymorphisms , high ka/ks ratios  and  <dig> detected gene conversion events averaging  <dig>  each subfamily.
for the subfamilies with ≥  <dig> members, the nucleotide divergence of lrr domains in paralogs, orthologs and between wild and cultivated accessions were also explored .
the average nucleotide divergence for the  <dig> large subfamilies was  <dig> , with an average paralog divergence of  <dig> , and ortholog divergence of  <dig> .
the highest divergence was noted between wild and cultivated accessions , which is consistent with their division distribution on the phylogenetic tree.table  <dig> nucleotide diversity, ka/ks and the number of gene conversions in each lrr subfamily
ka/ks
*represents wild-specific or cultivated-specific family
discussion
<dig>  copy number variation between species and between accessions
in the past years, many r genes have been identified and cloned in a variety of species, with the evolutionary histories of some r genes studied in detail  <cit> .
several patterns emerged in these studies.
generally, compared with other functional genes, r genes are always found with a higher rate of evolution.
several r gene families have been identified with a heterozygous evolutionary rate, with some homologs evolving faster and some homologs evolving more slowly  <cit> .
another distinct feature is that r genes tend to have dramatic variations in copy number and cluster together on the chromosome  <cit> .
in our study, the copy number of this r gene family was found to vary not only between species but also within species.
although the hot pepper is very closely related to the tomato and potato, the copy number of rpi-blb <dig> homologs in pepper was about  <dig> fold higher than tomato and potato homologs.
this suggests that after the speciation of the three species, these genes may have undergone rapid copy number variations, contributing to the significant gene expansion in pepper or frequent gene loss in tomato and potato genomes.
between wild and cultivated potato, as shown by the number of lrr regions cloned, the copy numbers also showed differences.
variations were even noted within the same species, such as among the  <dig> accessions of s. tuberosum and  <dig> accessions of s. bulbocastanum that displayed lrrs numbers from  <dig> to  <dig>  copy number variation is a common form of genome diversity that can create new genes, influence phenotypic characters and gene expression.
the variation in copy number of rpi-bib <dig> homologs in different species or accessions may be a strategy to minimize fitness costs in combating pathogens or pests  <cit> .
there are many other examples of variation in r gene copy number.
copy number variation  at the rp <dig> rust resistance locus was observed in different maize cultivars .
the copy number of rgc <dig> genes varied from  <dig> to  <dig> per genome in the seven lettuce genotypes  <cit> .
gene copy number varies between haplotypes or accessions consistent with a birth-and-death model of resistance gene evolution  <cit> .
some nbs-lrr genes are lost and new lineages evolve whilst others are retained.
the main mechanism for copy number variation is unequal crossover.
in our study, the variation identified among accessions within such short divergent times indicates that frequent unequal crossing over had taken place in the evolution of this gene family, and that the chromosome regions where r genes located may be the hot spots for recombination.
<dig>  diversifying selection in lrr regions
till now, most of cloned plant r genes are identified with a nucleotide binding site and a region of leucine rich repeats, also known as nbs-lrr genes.
the nbs-lrr genes are similar to nod-lrr genes in mammals, which are responsible for inflammatory and immune responses  <cit> .
the nbs domain, including atpase and g proteins, plays an important role in the plant immune response through signal transduction, while the lrr domain, a major determinant of recognition specificity, is responsible for avr gene recognition  <cit> .
lrr domains are supposed to evolve much faster and undergo stronger diversify selection than nbs domains.
in our study, we found that in the rpi-blb <dig> gene family, the nbs domains have lower divergence and are under purifying selection, while lrr domains have higher divergence and tend to suffer potential positive selection .
and a great divergence and much diversifying selection were detected in the  <dig> lrr regions we cloned.
the high genetic abundance of lrr regions was not only found among related species, but also within species and even within accessions.
these kinds of evolutionary patterns were in accordance with their functions.
the nbs domain is responsible for signal transduction and thus more conserved, while the lrr domain is responsible for avr gene recognition and thus part of a co-evolutionary relationship with the avr gene.
pathogens evolve diverse avr genes to infect plants and plants evolve diverse lrr domains in order to recognize a broad spectrum of avr genes.
previous studies have indicated that tomato mi- <dig> genes from this gene family have evolved rapidly by gene duplications and frequent sequence exchanges among homologs.
here, numerous gene conversion events haven been detected in the lrr regions, promoting the generation of chimeric genes and novel resistance specificities.
additionally, ~30 % of the examined lrr domains were found to contain frameshift mutations, which is consistent with the evolutionary pattern of nbs-lrrs in arabidopsis populations  <cit> .
these lrr domains provide a potential library without fitness cost for the production of new lrrs through sequence exchange and diversifying selection  <cit> .
most interestingly, in the gene family we studied, lrr domains from one gene family could recognize avr genes from four very different species, including nematode, aphid, whitefly and fungi, and the secret of recognizing multiple pathogens might lie in the variety of lrr domains.
<dig>  prediction of novel r genes
identifying new plant r gene is of fundamental importance in agriculture, especially for gene families like the rpi-blb2/mi- <dig> /cami gene family, which confer resistance in more than one pathogen species.
some other r gene families have also been reported to include r genes of different species carry resistance to different pathogens, such as the rp1/pi <dig> and rp3/pc-b gene families  <cit> .
the maize rust resistance gene rp <dig> and rice blast resistance gene pi <dig> are orthologous and have high sequence similarity.
another two r genes from the rp3/pc-b gene family, the maize gene rp <dig>  which confers resistance against maize common rust, and the sorghum gene pc-b, which provides resistance to root and crown rot.
however, these r genes from grass species confer resistance to the same type of pathogen, fungi.
the rpi-blb2/mi- <dig> /cami gene family in our study is an extremely special one, which includes r genes carrying resistance to four highly divergent organisms, including root-knot nematodes, aphids, whiteflies and the oomycete pathogen, making the presence of more resistance genes highly probable.
furthermore, members of the rpi-blb <dig> gene family, including three functional resistance genes, reside in the short arm of chromosome  <dig> .
actually, this region is a resistance hot spot having many cloned and mapped resistance genes, including the genes coding for nbs-lrr proteins and other types of proteins, such as receptor-like proteins .
in addition to rpi-blb <dig>  mi- <dig> , and cami, the tomato rlp genes cf- <dig> and cf- <dig>  which confer resistance to the leaf mold pathogen cladosporium fulvum  <cit> , and the ol- <dig> and ol- <dig> genes, which provide resistance against the tomato powdery mildew oidium neolycopersici  <cit> , are all physically located within the same region of chromosome  <dig>  additionally, bacterial wilt-resistance gene bw- <dig>  <cit> , a quantitative trait locus  ty- <dig>  which confers resistance to yellow leaf curl virus in tomato  <cit>  and qtls with late blight and blackleg resistance in potato are also mapped to this region  <cit> .
obviously, this complex r gene region is a mine for exploiting more useful resistance resources in solanaceae.
and the recent study found that fast evolving r genes in many grass species confer resistance to blast disease in rice  <cit> .
this gives us two signs.
i) maybe it is a common mechanism that a single r gene family can confer resistance in multi-species.
ii) the r genes in rapidly evolving families are most likely to confer resistance to fast-evolving pathogens.
the rapidly evolving homologs of functional r genes are potential new r genes.
thus studying the evolutionary patterns of cloned functional r genes might helpful to predict novel r genes.
likewise, fast evolving members of this gene family, such as the lrr subfamily  <dig> or  <dig>  exhibit high divergence and frequent sequence exchanges and thus are likely to harbor r genes.
clearly, the rpi-blb <dig> gene family studied in solanaceae is very important both in scientific research and agronomy.
further investigation and comparison of the gene family in solanaceae could provide a rich source of information for studying the evolution of r genes and subsequently enable the identification of potential candidate r genes for agricultural studies.
CONCLUSIONS
by studying the first gene family with resistant to multiple pathogens in plant, we revealed the specialty of this gene family lies in its unusually high level of genetic abundance.
our study provided insights for the evolutionary dynamic of multiple resistance genes and also provided further evidences that maybe it is common r genes like this are of multi-resistance.
additional files
additional file 1:
primers for cloning of lrrs.
additional file 2:
source of the  <dig> potato accessions used in this study.
additional file 3:
relative genomic positions of genes among potato , pepper  and tomato  along chromosome  <dig>
lina zhao and qijun zhang contributed equally to this work.
competing interests
the authors declare that they have no competing interests.
authors’ contributions
zx, ys, and lh conceived the study and designed the experiments.
zl, zq and gr performed the experiments.
zx and lh analyzed the data and wrote the paper.
all authors have read and approved the manuscript for publication.
authors’ information
not applicable.

 the potato, along with the tomato and pepper, are members of solanaceae family, with the potato being one of the world’s most important food crops. yet the potato is susceptible to a wide range of pathogens, with potato late blight being the most devastating. potato late blight is caused by phytophthora infestans, which greatly effects potato cultivation. annual potato production is currently over million tons and potential production could exceed million if potato late blight could be properly controlled. additionally, late blight disease is also responsible for the european potato famine in the 19th century which almost completely destroyed potato crops and led to the starvation of millions of people. thus, the control of such plant diseases is of fundamental importance. plants have evolved sophisticated systems to recognize pathogenic proteins. these recognition proteins, usually referred to as resistance genes, work in a highly specific manner according to the gene-for-gene interaction model, in which one plant r gene responds to a pathogen carrying a particular avirulence gene. this places the host r gene and pathogen avr gene in a co-evolutionary model in which the avr gene evolves different genotypes to avoid recognition while the r genes evolve to recognize the avr genes. studies have shown that mutations and deletions occur at high frequencies in avr genes, as a response of high rates also noted in r genes. during the past years, many plant r genes have been identified. most cloned r genes belong to a large gene family and encode proteins with nucleotide-binding sites and leucine-rich repeat domains. in potatoes, nbs-lrr genes such as r rb, r r rpi-bib and rpi-vnt have been associated with late blight resistance. among these genes, rpi-blb is of particular interest due to its homologs conferring resistance in different species and against at least different pathogens. the first member of the rpi-blb family found to confer resistance was mi-, which provides root knot nematode resistance in tomatoes, with this gene also found to confer aphid resistance. later, mi- was found to display resistance to sweet potato whitefly, making it the first and the only plant r gene that confers resistance against three distantly related pathogens. in the late blight resistance gene rpi-blb identified in potatoes was found to be a mi- homolog of mi-. the third r gene cloned from this family, cami, is from hot peppers and has been found to confer root knot nematode resistance, with a 99 % amino acid sequence identity shared with mi-. three genes, mi-, rpi-blb and cami, derived from different solanaceae organisms confer resistance against various highly divergent pathogens and have originated from the same gene family, making a common ancestor highly possible. what is it that makes this gene family so special and are there more gene families like this? to answer these questions, we studied the evolutionary history of this gene family. in brief, all homologs of these three resistance genes are referred to as rpi-blb homologs hereafter. first, homologs of this gene family across the tomato, potato, and pepper genomes were identified, with a varied copy number and a great nucleotide diversity of this gene family was found within species and between species. next, the potato lrr regions of this gene family were examined and homologs were cloned from potato accessions. a variation in copy number was common between the potato accessions and a high level of nucleotide divergence was found both between accessions and within accessions, suggesting rapid intraspecies and interspecies evolution in this gene family. meanwhile, strong positive selection and frequent sequence exchanges were also identified in the lrr regions. overall, our findings indicate the presence of a fast evolving r gene locus with a dramatic variation in copy number, high genetic abundance, and a strong diversifying selection. the above findings provide further insight enabling the identification of more novel r genes from this special family. methods genomic source and homolog identification three fully sequenced genomes were employed in this study, including genome of pepper, genome of tomato and genome of potato. the re-annotated nbs-lrr gene sets of potato and tomato were downloaded from previous studies. to identify the homologs of rpi-blb2/mi- /cami, the coding sequence of these genes were used as query and blast against the whole-genome cdss of the genomes and re-annotated gene sets. the e-value was set at 1e- and other parameters used default settings. hits with > 60 % coverage and > 60 % identity were deemed homologs. cloning of lrr domain of this gene family from potato accessions degenerate primers were designed on conserved sites at the edges of the lrr regions in order to gain comprehensive lrr domain sequence information. pcr was performed in potato accessions/genotypes, including wild accessions and cultivated accessions. pcr products were cloned into a pgem-t easy vector and ~ colonies of each cultivar were then sequenced by abi3100a automated sequencer until no new homolog sequence could be identified. sequences for all cloned lrr sequences have been deposited in the genbank under accession number from kr to kr sequence alignment and phylogenetic analysis the cdss or cloned lrr sequences were first translated into amino acid sequence and aligned by muscle implanted in mega and then back translated into nucleotide sequences. based on the alignments, a maximum likelihood tree was constructed. first, jmodeltest. was used to test the best fitted nucleotide substitution model for tree construction, then this model was used by phyml with bootstrap replicates for maximum likelihood tree construction. the nonsynonymous and synonymous nucleotide substitutions were calculated by dnasp version based on the nei-gojobori method. nucleotide diversity and divergence were estimated by π and dxy with jukes and cantor correction using dnasp version. gene conversion events were analyzed by genecov.  characterization of rpi-blb homologs in three species to construct the evolutionary relationship of the three nbs-lrr genes mi-, rpi-bib and cami, their homologs were first identified in genomes of the species. although the genome size of pepper was significantly larger than tomato and potato, the annotated gene number of these species was on the same level. it is reported that the genome expansion of the hot pepper is attributed to accumulation of repetitive sequences. more nbs-encoding genes have been detected in the pepper reference genome than in the tomato or potato reference genomes. since most nbs genes are duplicated genes, it is difficult to assembly and annotate these genes, and studies have showed that a certain number of nbs genes were missed in potato and tomato reference genomes. jupe et al. and andolfo et al. have identified and nbs genes by renseq in potato and tomato, respectively. as a result, we identified and homologs in the genomes of potato, tomato and pepper, respectively. after removing incomplete sequences, and homologs from potato, tomato and pepper were used for further analysis.table number of rpi-bib homologs in tomato, pepper and potato genomes ka/ks *these data are from re-annotated nbs gene sets by previous studies the nucleotide diversity and ka/ks were calculated for the homologs from each genome. the divergence varied a lot between genomes. the homologs from potato genome were found to be the most conservative with the lowest nucleotide diversity, while the homologs from tomato genome were the most divergent. such a huge difference was not found when examining ka/ks ratios which ranged from to among the different genomes. the nucleotide diversity and ka/ks were also determined for the nbs domain and lrr domain separately. across all of the three genomes, the lrr domain exhibited a higher nucleotide diversity and a higher ka/ks than the nbs domain. most of the identified homologs were located in tandem on chromosome in each genome. tandem repeat is a major mechanism in nbs-lrr expansion, with previous studies showing that nbs-lrr gene homologues often cluster in tandem on the chromosome. to explore the relation of genomic position among these homologs, a comparative position map of all of the genes on chromosome was constructed, extensive conservation of the gene order between these genomes was found among these species, suggesting a good synteny in this chromosome. these homologs all clustered on the same end of the chromosome, further reflecting that tandem repeat events took place in the expansion of this gene family. additionally, all the homologs on chromosome in potato and tomato have the highest similarity with the same gene in pepper, indicating all the homologs in these species share the same ancestor, and duplication events happened several times independently in each species.fig. positional relationships among rpi-blb2/mi- /cami homologs from potato, pepper and tomato reference genomes. genes are represented by blue arrows and genes with highest similarity between species are linked by green lines phylogenetic analysis of rpi-bib homologs from tomato, pepper and potato to further explore the evolutionary history of this gene family, a phylogenetic tree was constructed based on the cloned functional genes and their homologs within the reference genomes. as shown in fig.  genes from the same species mostly clustered together on the tree. this kind of topology further indicates that the expansion of this gene family took place separately in these species after they split from their common ancestor. this tree also shows a closer homology between the tomato and potato than with pepper, coinciding with their genetic relationships on the species level.fig. phylogenetic tree derived from rpi-blb2/mi- /cami homologs. the maximum-likelihood phylogenetic tree is constructed by phyml based on the cds sequences of each homolog subfamilies were divided with an identity > 80 % and bootstrap > 60 %, with one major subfamily comprised of genes, two families with and genes and minor subfamilies, each comprised of 2– genes, identified. all of the known r-genes were found in subfamily the major subfamily, which includes homologs from tomato, genes from potato and gene from pepper. then nucleotide divergence and ka/ks were calculated within each subfamily. nucleotide divergence varies from to and ka/ks varies from to. to further compare evolutionary patterns among homologs in the largest subfamily, the divergence of genes within tomato, potato and between them was calculated. the divergence of paralogs within tomato and potato was and, while the divergence of homologs between tomato and potato was a little higher than within species. this shows that in subfamily homologs from tomato are more conserved while homologs from potato are more divergent.table nucleotide diversity, ka/ks and the number of gene conversions in each gene subfamily ka/ks nucleotide variation in lrr regions among potato accessions typically, in nbs-lrr genes, the lrr domain is responsible for recognition of the avr gene and is more divergent than the nbs domain. to study the evolutionary history of the lrr domain in this gene family, homologs from wild and cultivated potato accessions were cloned averaging homologs in each accession, ranging from to reflecting frequent copy number variations in these accessions. interestingly, more copy number was detected in cultivated potatoes, with accessions having copies and having copies. it might be because the cultivated potatoes are all tetraploid species, whereas most of the wild species are diploid.table statistics of the lrr regions cloned in each potato accession ka/ks pp10 ka/ks to estimate the relative evolutionary rate of the lrr regions, the nucleotide divergence and ka/ks were calculated within each accession. the average nucleotide divergence of the homologs in all of the accessions was. the average nucleotide divergence for the wild-type and cultivated accessions were and, suggesting that the lrr region in this resistance gene locus is polymorphic. moreover, the divergence in cultivated accession was significantly higher than that in wild accessions, which implies that the lrr domains in cultivated accession might have been diversified during the process of cultivation. additionally, a relatively high ka/ks ratio was found both in the wild and cultivated accessions, with accessions >  accessions >  and the other two accessions ≈ . on average, the ka/ks value of the lrr domain was higher than the whole genes from tomato, pepper and potato, indicating that the lrr domains in particular have been positively selected. when plotting nucleotide divergence against ka/ks ratio, a significant negative relationship was found, suggesting that the nucleotide divergence reduces as the pressure of positive selection increases. lastly, frameshift mutations were annotated in each cloned lrr regions and ~30 % of the lrr regions were characterized with frequent frameshift mutations. the most frequent mutations were detected in the wild species s.demissum, which is a hexaploid, while the other wild accessions or genotypes had fewer frameshift mutations. in contrast, the mutations detected in the cultivated species distributed relatively evenly. collectively, these pseuodo-lrr domains might provide a resource or raw material for producing more new lrrs. lrr regions are driven by diversified selection to investigate the genetic relations and their evolutionary history between lrr homologs from different accessions, a phylogenetic tree was constructed based on cloned lrr regions, together with lrr sequences from potato genome and the functional genes. the topological structure of the tree showed the cloned lrr regions to be scattered, suggesting that the lrr regions had a high degree of polymorphisms and varied between accessions and within accessions. moreover, the lrr regions were not distinctly divided into two clear groups, wild and cultivated clades, indicating a frequent introgression during the cultivation of potato.fig. phylogenetic tree derived from the lrr regions cloned in this study and lrr domains from previous studies. the lrr domains were cloned from different potato accessions/genotypes, of the lrr domains were drawn from the potato genomes and the rest from previous cloned genes in potato, pepper, and tomato however, the lrr regions from wild and cultivated accessions were not totally confused or mixed together in the tree, with most small subclades being wild-specific or cultivated-specific. additionally, in the wild-specific or cultivated-specific subclades, the lrr regions cloned from an accession or genotype were not always clustered, illustrating the complexity of orthologous relationship among the lrr regions. it also suggests that frequent sequence exchanges may have taken place among accessions of cultivated potatoes or wild-type species. subfamilies were also divided based on an identity > 80 % and bootstrap > 60 %. in total, subfamilies were identified, each with 2– members. the resistance lrrs of tomato and pepper were clustered in subfamily whereas the lrr of rpi-blb in potato was isolated. except the resistance-specific subfamily and subfamily which contained only genome sequences, of subfamilies were wild-specific and subfamilies are cultivated-specific. to further investigate evolutionary patterns of subfamilies, the nucleotide diversity, ka/ks and gene conversions were calculated and identified in each subfamily. to obtain reliable results, subfamilies with few members were excluded from further analysis. subfamilies were kept, including subfamilies with ≥  members. the nucleotide divergence of subfamily which consists of two functional lrrs, was with a ka/ks >  an indicator of positive selection. additionally, the average number of members in wild-specific subfamilies was, much higher than the number in cultivated-specific subfamilies. when comparing wild- and cultivated-specific subfamilies, the wild-specific had higher nucleotide diversity, a higher ka/ks and an elevated rate of gene conversions events, suggesting that more abundant genetic resources could be re-used in wild species. moreover, four complex subfamilies with 7– members had exceptionally high levels of polymorphisms, high ka/ks ratios and detected gene conversion events averaging each subfamily. for the subfamilies with ≥  members, the nucleotide divergence of lrr domains in paralogs, orthologs and between wild and cultivated accessions were also explored. the average nucleotide divergence for the large subfamilies was, with an average paralog divergence of, and ortholog divergence of. the highest divergence was noted between wild and cultivated accessions, which is consistent with their division distribution on the phylogenetic tree.table nucleotide diversity, ka/ks and the number of gene conversions in each lrr subfamily ka/ks *represents wild-specific or cultivated-specific family discussion copy number variation between species and between accessions in the past years, many r genes have been identified and cloned in a variety of species, with the evolutionary histories of some r genes studied in detail. several patterns emerged in these studies. generally, compared with other functional genes, r genes are always found with a higher rate of evolution. several r gene families have been identified with a heterozygous evolutionary rate, with some homologs evolving faster and some homologs evolving more slowly. another distinct feature is that r genes tend to have dramatic variations in copy number and cluster together on the chromosome. in our study, the copy number of this r gene family was found to vary not only between species but also within species. although the hot pepper is very closely related to the tomato and potato, the copy number of rpi-blb homologs in pepper was about fold higher than tomato and potato homologs. this suggests that after the speciation of the three species, these genes may have undergone rapid copy number variations, contributing to the significant gene expansion in pepper or frequent gene loss in tomato and potato genomes. between wild and cultivated potato, as shown by the number of lrr regions cloned, the copy numbers also showed differences. variations were even noted within the same species, such as among the accessions of s. tuberosum and accessions of s. bulbocastanum that displayed lrrs numbers from to copy number variation is a common form of genome diversity that can create new genes, influence phenotypic characters and gene expression. the variation in copy number of rpi-bib homologs in different species or accessions may be a strategy to minimize fitness costs in combating pathogens or pests. there are many other examples of variation in r gene copy number. copy number variation at the rp rust resistance locus was observed in different maize cultivars. the copy number of rgc genes varied from to per genome in the seven lettuce genotypes. gene copy number varies between haplotypes or accessions consistent with a birth-and-death model of resistance gene evolution. some nbs-lrr genes are lost and new lineages evolve whilst others are retained. the main mechanism for copy number variation is unequal crossover. in our study, the variation identified among accessions within such short divergent times indicates that frequent unequal crossing over had taken place in the evolution of this gene family, and that the chromosome regions where r genes located may be the hot spots for recombination. diversifying selection in lrr regions till now, most of cloned plant r genes are identified with a nucleotide binding site and a region of leucine rich repeats, also known as nbs-lrr genes. the nbs-lrr genes are similar to nod-lrr genes in mammals, which are responsible for inflammatory and immune responses. the nbs domain, including atpase and g proteins, plays an important role in the plant immune response through signal transduction, while the lrr domain, a major determinant of recognition specificity, is responsible for avr gene recognition. lrr domains are supposed to evolve much faster and undergo stronger diversify selection than nbs domains. in our study, we found that in the rpi-blb gene family, the nbs domains have lower divergence and are under purifying selection, while lrr domains have higher divergence and tend to suffer potential positive selection. and a great divergence and much diversifying selection were detected in the lrr regions we cloned. the high genetic abundance of lrr regions was not only found among related species, but also within species and even within accessions. these kinds of evolutionary patterns were in accordance with their functions. the nbs domain is responsible for signal transduction and thus more conserved, while the lrr domain is responsible for avr gene recognition and thus part of a co-evolutionary relationship with the avr gene. pathogens evolve diverse avr genes to infect plants and plants evolve diverse lrr domains in order to recognize a broad spectrum of avr genes. previous studies have indicated that tomato mi- genes from this gene family have evolved rapidly by gene duplications and frequent sequence exchanges among homologs. here, numerous gene conversion events haven been detected in the lrr regions, promoting the generation of chimeric genes and novel resistance specificities. additionally, ~30 % of the examined lrr domains were found to contain frameshift mutations, which is consistent with the evolutionary pattern of nbs-lrrs in arabidopsis populations. these lrr domains provide a potential library without fitness cost for the production of new lrrs through sequence exchange and diversifying selection. most interestingly, in the gene family we studied, lrr domains from one gene family could recognize avr genes from four very different species, including nematode, aphid, whitefly and fungi, and the secret of recognizing multiple pathogens might lie in the variety of lrr domains. prediction of novel r genes identifying new plant r gene is of fundamental importance in agriculture, especially for gene families like the rpi-blb2/mi- /cami gene family, which confer resistance in more than one pathogen species. some other r gene families have also been reported to include r genes of different species carry resistance to different pathogens, such as the rp1/pi and rp3/pc-b gene families. the maize rust resistance gene rp and rice blast resistance gene pi are orthologous and have high sequence similarity. another two r genes from the rp3/pc-b gene family, the maize gene rp which confers resistance against maize common rust, and the sorghum gene pc-b, which provides resistance to root and crown rot. however, these r genes from grass species confer resistance to the same type of pathogen, fungi. the rpi-blb2/mi- /cami gene family in our study is an extremely special one, which includes r genes carrying resistance to four highly divergent organisms, including root-knot nematodes, aphids, whiteflies and the oomycete pathogen, making the presence of more resistance genes highly probable. furthermore, members of the rpi-blb gene family, including three functional resistance genes, reside in the short arm of chromosome. actually, this region is a resistance hot spot having many cloned and mapped resistance genes, including the genes coding for nbs-lrr proteins and other types of proteins, such as receptor-like proteins. in addition to rpi-blb mi-, and cami, the tomato rlp genes cf- and cf- which confer resistance to the leaf mold pathogen cladosporium fulvum, and the ol- and ol- genes, which provide resistance against the tomato powdery mildew oidium neolycopersici, are all physically located within the same region of chromosome additionally, bacterial wilt-resistance gene bw-, a quantitative trait locus ty- which confers resistance to yellow leaf curl virus in tomato and qtls with late blight and blackleg resistance in potato are also mapped to this region. obviously, this complex r gene region is a mine for exploiting more useful resistance resources in solanaceae. and the recent study found that fast evolving r genes in many grass species confer resistance to blast disease in rice. this gives us two signs. i) maybe it is a common mechanism that a single r gene family can confer resistance in multi-species. ii) the r genes in rapidly evolving families are most likely to confer resistance to fast-evolving pathogens. the rapidly evolving homologs of functional r genes are potential new r genes. thus studying the evolutionary patterns of cloned functional r genes might helpful to predict novel r genes. likewise, fast evolving members of this gene family, such as the lrr subfamily or exhibit high divergence and frequent sequence exchanges and thus are likely to harbor r genes. clearly, the rpi-blb gene family studied in solanaceae is very important both in scientific research and agronomy. further investigation and comparison of the gene family in solanaceae could provide a rich source of information for studying the evolution of r genes and subsequently enable the identification of potential candidate r genes for agricultural studies.  by studying the first gene family with resistant to multiple pathogens in plant, we revealed the specialty of this gene family lies in its unusually high level of genetic abundance. our study provided insights for the evolutionary dynamic of multiple resistance genes and also provided further evidences that maybe it is common r genes like this are of multi-resistance. additional files additional file 1: primers for cloning of lrrs. additional file 2: source of the potato accessions used in this study. additional file 3: relative genomic positions of genes among potato, pepper and tomato along chromosome lina zhao and qijun zhang contributed equally to this work. competing interests the authors declare that they have no competing interests. authors’ contributions zx, ys, and lh conceived the study and designed the experiments. zl, zq and gr performed the experiments. zx and lh analyzed the data and wrote the paper. all authors have read and approved the manuscript for publication. authors’ information not applicable. 
to gain more information pertaining to gene resistance within this family, lrr regions were cloned in different wild/cultivated potato accessions. our findings suggest that this family harbors an unusually high level of genetic abundance, making it of particular interest. these genes confer resistance against distantly related pathogen species in different solanaceae hosts. again, frequent copy number variations and a high level of divergence between homolog were observed common among accessions. also, frequent frameshift mutations and abundant gene conversion events were identified in these lrr regions. this revealed a large variation in copy number within this gene family among species and a great divergence was found both between and within species. together with other reported examples, our study also provides evidence that multi-resistance is a common trait in r gene families like this. first, homologs of this gene family were identified in the pepper, tomato and potato genomes. three nbs-lrr genes, rpi-blb mi-, and cami, constitute a very special plant resistance gene family. to characterize this noted resistance, we conducted a series of studies on this gene family.
 three nbs-lrr genes, rpi-blb mi-, and cami, constitute a very special plant resistance gene family. these genes confer resistance against distantly related pathogen species in different solanaceae hosts. to characterize this noted resistance, we conducted a series of studies on this gene family.  first, homologs of this gene family were identified in the pepper, tomato and potato genomes. this revealed a large variation in copy number within this gene family among species and a great divergence was found both between and within species. to gain more information pertaining to gene resistance within this family, lrr regions were cloned in different wild/cultivated potato accessions. again, frequent copy number variations and a high level of divergence between homolog were observed common among accessions. the divergence within species was so high that it reaches the level of divergence between species. also, frequent frameshift mutations and abundant gene conversion events were identified in these lrr regions.  our findings suggest that this family harbors an unusually high level of genetic abundance, making it of particular interest. together with other reported examples, our study also provides evidence that multi-resistance is a common trait in r gene families like this. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords plant r generpi-blb2/mi- /camievolutionary historyissue-copyright-statement© the author 2015 