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

Quantitative trait locus (QTL) analysis and fine-mapping for Fusarium oxysporum disease resistance in Raphanus sativus using GRAS-Di technology

Fusarium wilt is a significant disease in radish, but the genetic mechanisms controlling yellows resistance (YR) are not well understood. This study aimed to identify YR-QTLs and to fine-map one of them using F2:3 populations developed from resistant and susceptible radish parents. In this study, two high-density genetic maps each containing shared co-dominant markers and either female or male dominant markers that spanned 988.6 and 1127.5 cM with average marker densities of 1.40 and 1.53 cM, respectively, were generated using Genotyping by Random Amplicon Sequencing-Direct (GRAS-Di) technology. We identified two YR-QTLs on chromosome R2 and R7, and designated the latter as ForRs1 as the major QTL. Fine mapping narrowed down the ForRs1 locus to a 195 kb region. Among the 16 predicted genes in the delimited region, 4 genes including two receptor-like protein and -kinase genes (RLP/RLK) were identified as prime candidates for ForRs1 based on the nucleotide sequence comparisons between the parents and their predicted functions. This study is the first to use a GRAS-Di for genetic map construction of cruciferous crops and fine map the YR-QTL on the R7 chromosome of radish. These findings will provide groundbreaking insights into radish YR breeding and understanding the genetics of YR mechanism.

Diversity of Late Blight Resistance Genes in the VIR Potato Collection

Late blight (LB) caused by the oomycete Phytophthora infestans (Mont.) de Bary is the greatest threat to potato production worldwide. Current potato breeding for LB resistance heavily depends on the introduction of new genes for resistance to P. infestans (Rpi genes). Such genes have been discovered in highly diverse wild, primitive, and cultivated species of tuber-bearing potatoes (Solanum L. section Petota Dumort.) and introgressed into the elite potato cultivars by hybridization and transgenic complementation. Unfortunately, even the most resistant potato varieties have been overcome by LB due to the arrival of new pathogen strains and their rapid evolution. Therefore, novel sources for germplasm enhancement comprising the broad-spectrum Rpi genes are in high demand with breeders who aim to provide durable LB resistance. The Genbank of the N.I. Vavilov Institute of Plant Genetic Resources (VIR) in St. Petersburg harbors one of the world's largest collections of potato and potato relatives. In this study, LB resistance was evaluated in a core selection representing 20 species of seven Petota series according to the Hawkes (1990) classification: Bulbocastana (Rydb.) Hawkes, Demissa Buk., Longipedicellata Buk., Maglia Bitt., Pinnatisecta (Rydb.) Hawkes, Tuberosa (Rydb.) Hawkes (wild and cultivated species), and Yungasensa Corr. LB resistance was assessed in 96 accessions representing 18 species in the laboratory test with detached leaves using a highly virulent and aggressive isolate of P. infestans. The Petota species notably differed in their LB resistance: S. bulbocastanum Dun., S. demissum Lindl., S. cardiophyllum Lindl., and S. berthaultii Hawkes stood out at a high frequency of resistant accessions (7-9 points on a 9-point scale). Well-established specific SCAR markers of ten Rpi genes-Rpi-R1, Rpi-R2/Rpi-blb3, Rpi-R3a, Rpi-R3b, Rpi-R8, Rpi-blb1/Rpi-sto1, Rpi-blb2, and Rpi-vnt1-were used to mine 117 accessions representing 20 species from seven Petota series. In particular, our evidence confirmed the diverse Rpi gene location in two American continents. The structural homologs of the Rpi-R2, Rpi-R3a, Rpi-R3b, and Rpi-R8 genes were found in the North American species other than S. demissum, the species that was the original source of these genes for early potato breeding, and in some cases, in the South American Tuberosa species. The Rpi-blb1/Rpi-sto1 orthologs from S. bulbocastanum and S. stoloniferum Schlechtd et Bché were restricted to genome B in the Mesoamerican series Bulbocastana, Pinnatisecta, and Longipedicellata. The structural homologs of the Rpi-vnt1 gene that were initially identified in the South American species S. venturii Hawkes and Hjert. were reported, for the first time, in the North American series of Petota species.

Comparative Analysis of the Genetic Diversity of Chilean Cultivated Potato Based on a Molecular Study of Authentic Herbarium Specimens and Present-Day Gene Bank Accessions

At the end of the 1920s, Vavilov organized several potato-collecting missions in South and Central America. Vavilov and his colleagues, Juzepczuk and Bukasov, participated in these expeditions and worked on gathered material, designated two centers of potato varietal riches and diversity-the Peru-Bolivia high-mountain center and the southern coast of Chile. The WIR Herbarium holds authentic specimens of many taxa described by Russian taxonomists. Here, a set of 20 plastid DNA-specific markers was applied for 49 authentic herbarium specimens of Solanum tuberosum L. from the WIR Herbarium to analyze the genetic diversity of the landrace population collected by Juzepczuk in 1928 in southern-central Chile. Two plastid DNA types, T and A, and two chlorotypes were identified in herbarium specimens, with a clear predominance (96%) of chlorotype cpT_III. In addition, we analyzed 46 living Chilean accessions from the VIR field potato gene bank that were collected after the appearance of Phytophthora infestans in Chile. These living accessions were differentiated into four chlorotypes. Finding a D-type cytoplasm in living Chilean accessions that possess two new chlorotypes indicates a replacement of native cultivars and introgression from the wild Mexican species S. demissum that was actively used in breeding as a source of race-specific resistance to late blight.

Genome architecture and tetrasomic inheritance of autotetraploid potato

Potato (Solanum tuberosum) is the most consumed non-cereal food crop. Most commercial potato cultivars are autotetraploids with highly heterozygous genomes, severely hampering genetic analyses and improvement. By leveraging the state-of-the-art sequencing technologies and polyploid graph binning, we achieved a chromosome-scale, haplotype-resolved genome assembly of a cultivated potato, Cooperation-88 (C88). Intra-haplotype comparative analyses revealed extensive sequence and expression differences in this tetraploid genome. We identified haplotype-specific pericentromeres on chromosomes, suggesting a distinct evolutionary trajectory of potato homologous centromeres. Furthermore, we detected double reduction events that are unevenly distributed on haplotypes in 1021 of 1034 selfing progeny, a feature of autopolyploid inheritance. By distinguishing maternal and paternal haplotype sets in C88, we simulated the origin of heterosis in cultivated tetraploid with a survey of 3110 tetra-allelic loci with deleterious mutations, which were masked in the heterozygous condition by two parents. This study provides insights into the genomic architecture of autopolyploids and will guide their breeding.

Genome-wide SNP and InDel analysis of three Philippine mango species inferred from whole-genome sequencing

Background

The Philippines is among the top 10 major exporters of mango worldwide. However, genomic studies of Philippine mangoes remain largely unexplored and lacking. Here, we sequenced the whole genome of the three Philippine mango species, namely, Mangifera odorata (Huani), Mangifera altissima (Paho), and Mangifera indica “Carabao” variety using Illumina HiSeq 2500, to identify and analyze their genome-wide variants (SNPs and InDels).

Results

The high confidence variants were identified by successfully mapping 93–95% of the quality-filtered reads to the Alphonso and Tommy Atkins mango reference genomes. Using these two currently available mango genomes, most variants were observed in M. odorata (4,353,063 and 4,277,287), followed by M. altissima (3,392,763 and 3,449,917), and lastly, M. indica Carabao (2,755,267 and 2,852,480). Approximately 50, 46, and 38% of the variants were unique in the three Philippine mango genomes. The analysis of variant effects and functional annotation across the three mango species revealed 56,982 variants with high-impact effects mapped onto 37,746 genes, of which 25% were found to be novel. The affected mango genes include those with potential economic importance such as 6945 genes for defense/resistance/immune response, 323 genes for fruit development, and 338 genes for anthocyanin production.

Conclusions

To date, this is the first sequencing effort to comprehensively analyze genome-wide variants essential for the development of genome-wide markers specific to these mango species native to the Philippines. This study provides an important genomic resource that can be used for the genetic improvement of mangoes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s43141-022-00326-3.

Resequencing and genome-wide association studies of autotetraploid potato

Potato is the fourth most important food crop in the world. Although with a long history for breeding approaches, genomic information and association between genes and agronomic traits remain largely unknown particularly in autotetraploid potato cultivars, which limit the molecular breeding progression. By resequencing the genome of 108 main cultivar potato accessions with rich genetic diversity and population structure from International Potato Center, with approximate 20-fold coverage, we revealed more than 27 million Single Nucleotide Polymorphisms and ~ 3 million Insertion and Deletions with high quality and accuracy. Domestication analysis and genome-wide association studies (GWAS) identified candidate loci related to photoperiodic flowering time and temperature sensitivity as well as disease resistance, providing informative insights into the selection and domestication of cultivar potato. In addition, GWAS with GWASploy for 25 agronomic traits identified candidate loci by association signals, especially those related to tuber size, small-sized tuber weight and tuber thickness that was also validated by transcriptome analysis. Our study provides a valuable resource that facilitates the elucidation of domestication process as well as the genetic studies and agronomic improvement of autotetraploid potato.

Characterisation of Endogenous Peptides Present in Virgin Olive Oil

The low molecular weight peptide composition of virgin olive oil (VOO) is mostly unknown. We aimed to investigate the composition of the endogenous peptides present in VOO, the protein sources from which those peptides originate and their biological activities. A water-soluble extract containing peptides was obtained from VOO. The peptides were separated by size-exclusion using fast protein liquid chromatography, and the low molecular weight fraction (1600–700 kDa) was analysed by nanoscale liquid chromatography Orbitrap coupled with tandem mass spectrometry and de novo sequencing. Nineteen new peptides were identified by Peaks database algorithm, using the available Olea europaea (cv. Farga) genome database. Eight new peptides were also identified by Peaks de novo sequencing. The protein sources of the peptides detected in the database by Peaks DB were identified by BLAST-P search. Seed storage proteins were among the most frequent sources of VOO peptides. BIOPEP software was used to predict the biological activities of peptides and to simulate (in silico) the proteolytic activity of digestive enzymes on the detected peptide sequences. A selection of synthetic peptides was obtained for investigation of their bioactivities. Peptides VCGEAFGKA, NALLCSNS, CPANGFY, CCYSVY and DCHYFL possessed strong ACE-inhibitory and antioxidant activities in vitro. Antioxidant peptides could play a role in VOO quality.

Insertion of Badnaviral DNA in the Late Blight Resistance Gene (R1a) of Brinjal Eggplant (Solanum melongena)

Endogenous viral elements (EVEs) are widespread in plant genomes. They result from the random integration of viral sequences into host plant genomes by horizontal DNA transfer and have the potential to alter host gene expression. We performed a large-scale search for co-transcripts including caulimovirid and plant sequences in 1,678 plant and 230 algal species and characterized 50 co-transcripts in 45 distinct plant species belonging to lycophytes, ferns, gymnosperms and angiosperms. We found that insertion of badnavirus EVEs along with Ty-1 copia mobile elements occurred into a late blight resistance gene (R1) of brinjal eggplant (Solanum melongena) and wild relatives in genus Solanum and disrupted R1 orthologs. EVEs of two previously unreported badnaviruses were identified in the genome of S. melongena, whereas EVEs from an additional novel badnavirus were identified in the genome of S. aethiopicum, the cultivated scarlet eggplant. Insertion of these viruses in the ancestral lineages of the direct wild relatives of the eggplant would have occurred during the last 3 Myr, further supporting the distinctiveness of the group of the eggplant within the giant genus Solanum.

Functional analysis of RXLR effectors from the New Zealand kauri dieback pathogen Phytophthora agathidicida

New Zealand kauri is an ancient, iconic, gymnosperm tree species that is under threat from a lethal dieback disease caused by the oomycete Phytophthora agathidicida. To gain insight into this pathogen, we determined whether proteinaceous effectors of P. agathidicida interact with the immune system of a model angiosperm, Nicotiana, as previously shown for Phytophthora pathogens of angiosperms. From the P. agathidicida genome, we defined and analysed a set of RXLR effectors, a class of proteins that typically have important roles in suppressing or activating the plant immune system. RXLRs were screened for their ability to activate or suppress the Nicotiana plant immune system using Agrobacterium tumefaciens transient transformation assays. Nine P. agathidicida RXLRs triggered cell death or suppressed plant immunity in Nicotiana, of which three were expressed in kauri. For the most highly expressed, P. agathidicida (Pa) RXLR24, candidate cognate immune receptors associated with cell death were identified in Nicotiana benthamiana using RNA silencing-based approaches. Our results show that RXLRs of a pathogen of gymnosperms can interact with the immune system of an angiosperm species. This study provides an important foundation for studying the molecular basis of plant-pathogen interactions in gymnosperm forest trees, including kauri.

Potato Virus Y Infection Alters Small RNA Metabolism and Immune Response in Tomato

Potato virus Y (PVY) isolate PVY^C-to induces growth reduction and foliar symptoms in tomato, but new vegetation displays symptom recovery at a later stage. In order to investigate the role of micro(mi)RNA and secondary small(s)RNA-regulated mechanisms in tomato defenses against PVY, we performed sRNA sequencing from healthy and PVY^C-to infected tomato plants at 21 and 30 days post-inoculation (dpi). A total of 792 miRNA sequences were obtained, among which were 123 canonical miRNA sequences, many isomiR variants, and 30 novel miRNAs. MiRNAs were mostly overexpressed in infected vs. healthy plants, whereas only a few miRNAs were underexpressed. Increased accumulation of isomiRs was correlated with viral infection. Among miRNA targets, enriched functional categories included resistance (R) gene families, transcription and hormone factors, and RNA silencing genes. Several 22-nt miRNAs were shown to target R genes and trigger the production of 21-nt phased sRNAs (phasiRNAs). Next, 500 phasiRNA-generating loci were identified, and were shown to be mostly active in PVY-infected tissues and at 21 dpi. These data demonstrate that sRNA-regulated host responses, encompassing miRNA alteration, diversification within miRNA families, and phasiRNA accumulation, regulate R and disease-responsive genes. The dynamic regulation of miRNAs and secondary sRNAs over time suggests a functional role of sRNA-mediated defenses in the recovery phenotype.

Structural, Functional, and Genomic Diversity of Plant NLR Proteins: An Evolved Resource for Rational Engineering of Plant Immunity

Plants employ a diverse intracellular system of NLR (nucleotide binding-leucine-rich repeat) innate immune receptors to detect pathogens of all types. These receptors represent valuable agronomic traits that plant breeders rely on to maximize yield in the face of devastating pathogens. Despite their importance, the mechanistic underpinnings of NLR-based disease resistance remain obscure. The rapidly increasing numbers of plant genomes are revealing a diverse array of NLR-type immune receptors. In parallel, mechanistic studies are describing diverse functions for NLR immune receptors. In this review, we intend to broadly describe how the structural, functional, and genomic diversity of plant immune receptors can provide a valuable resource for rational engineering of plant immunity.

Transcriptome and Degradome of microRNAs and Their Targets in Response to Drought Stress in the Plants of a Diploid and Its Autotetraploid Paulownia australis

MicroRNAs (miRNAs) are small, non-coding RNAs that play vital roles in plant growth, development, and stress response. Increasing numbers of studies aimed at discovering miRNAs and analyzing their functions in plants are being reported. In this study, we investigated the effect of drought stress on the expression of miRNAs and their targets in plants of a diploid and derived autotetraploid Paulownia australis. Four small RNA (sRNA) libraries and four degradome libraries were constructed from diploid and autotetraploid P. australis plants treated with either 75% or 25% relative soil water content. A total of 33 conserved and 104 novel miRNAs (processing precision value > 0.1) were identified, and 125 target genes were identified for 36 of the miRNAs by using the degradome sequencing. Among the identified miRNAs, 54 and 68 were differentially expressed in diploid and autotetraploid plants under drought stress (25% relative soil water content), respectively. The expressions of miRNAs and target genes were also validated by quantitative real-time PCR. The results showed that the relative expression trends of the randomly selected miRNAs were similar to the trends predicted by Illumina sequencing. And the correlations between miRNAs and their target genes were also analyzed. Furthermore, the functional analysis showed that most of these miRNAs and target genes were associated with plant development and environmental stress response. This study provided molecular evidence for the possible involvement of certain miRNAs in the drought response and/or tolerance in P. australis, and certain level of differential expression between diploid and autotetraploid plants.

Genetics and molecular mapping of resistance to Plasmodiophora brassicae pathotypes 2, 3, 5, 6, and 8 in rutabaga (Brassica napus var. napobrassica)

Clubroot disease, caused by Plasmodiophora brassicae, is a threat to the production of Brassica crops including oilseed B. napus. In Canada, several pathotypes of this pathogen, such as pathotypes 2, 3, 5, 6, and 8, were identified, and resistance to these pathotypes was found in a rutabaga (B. napus var. napobrassica) genotype. In this paper, we report the genetic basis and molecular mapping of this resistance by use of F₂, backcross (BC₁), and doubled haploid (DH) populations generated from crossing of this rutabaga line to a susceptible spring B. napus canola line. The F₁, F₂, and BC₁ populations were evaluated for resistance to pathotype 3, and the DH population was evaluated for resistance to pathotypes 2, 3, 5, 6, and 8. A 3:1 segregation in F₂ and a 1:1 segregation in BC₁ were found for resistance to pathotype 3, and a 1:1 segregation was found in the DH population for resistance to all pathotypes. Molecular mapping by using the DH population identified a genomic region on chromosome A8 carrying resistance to all five pathotypes. This suggests that a single gene or a cluster of genes, located in this genomic region, is involved in the control of resistance to these pathotypes.

A disease resistance locus on potato and tomato chromosome 4 exhibits a conserved multipartite structure displaying different rates of evolution in different lineages

BACKGROUND

In plant genomes, NB-LRR based resistance (R) genes tend to occur in clusters of variable size in a relatively small number of genomic regions. R-gene sequences mostly differentiate by accumulating point mutations and gene conversion events. Potato and tomato chromosome 4 harbours a syntenic R-gene locus (known as the R2 locus in potato) that has mainly been examined in central American/Mexican wild potato species on the basis of its contribution to resistance to late blight, caused by the oomycete pathogen Phytophthora infestans. Evidence to date indicates the occurrence of a fast evolutionary mode characterized by gene conversion events at the locus in these genotypes.

RESULTS

A physical map of the R2 locus was developed for three Solanum tuberosum genotypes and used to identify the tomato syntenic sequence. Functional annotation of the locus revealed the presence of numerous resistance gene homologs (RGHs) belonging to the R2 gene family (R2GHs) organized into a total of 4 discrete physical clusters, three of which were conserved across S. tuberosum and tomato. Phylogenetic analysis showed clear orthology/paralogy relationships between S. tuberosum R2GHs but not in R2GHs cloned from Solanum wild species. This study confirmed that, in contrast to the wild species R2GHs, which have evolved through extensive sequence exchanges between paralogs, gene conversion was not a major force for differentiation in S. tuberosum R2GHs, and orthology/paralogy relationships have been maintained via a slow accumulation of point mutations in these genotypes.

CONCLUSIONS

S. tuberosum and Solanum lycopersicum R2GHs evolved mostly through duplication and deletion events, followed by gradual accumulation of mutations. Conversely, widespread gene conversion is the major evolutionary force that has shaped the locus in Mexican wild potato species. We conclude that different selective forces shaped the evolution of the R2 locus in these lineages and that co-evolution with a pathogen steered selection on different evolutionary paths.

Detection of an inversion in the Ty-2 region between S. lycopersicum and S. habrochaites by a combination of de novo genome assembly and BAC cloning

A chromosomal inversion associated with the tomato Ty - 2 gene for TYLCV resistance is the cause of severe suppression of recombination in a tomato Ty - 2 introgression line. Among tomato and its wild relatives inversions are often observed, which result in suppression of recombination. Such inversions hamper the transfer of important traits from a related species to the crop by introgression breeding. Suppression of recombination was reported for the TYLCV resistance gene, Ty-2, which has been introgressed in cultivated tomato (Solanum lycopersicum) from the wild relative S. habrochaites accession B6013. Ty-2 was mapped to a 300-kb region on the long arm of chromosome 11. The suppression of recombination in the Ty-2 region could be caused by chromosomal rearrangements in S. habrochaites compared with S. lycopersicum. With the aim of visualizing the genome structure of the Ty-2 region, we compared the draft de novo assembly of S. habrochaites accession LYC4 with the sequence of cultivated tomato ('Heinz'). Furthermore, using populations derived from intraspecific crosses of S. habrochaites accessions, the order of markers in the Ty-2 region was studied. Results showed the presence of an inversion of approximately 200 kb in the Ty-2 region when comparing S. lycopersicum and S. habrochaites. By sequencing a BAC clone from the Ty-2 introgression line, one inversion breakpoint was identified. Finally, the obtained results are discussed with respect to introgression breeding and the importance of a priori de novo sequencing of the species involved.

A novel workflow correlating RNA-seq data to Phythophthora infestans resistance levels in wild Solanum species and potato clones

Comparative transcriptomics between species can provide valuable understanding of plant-pathogen interactions. Here, we focus on wild Solanum species and potato clones with varying degree of resistance against Phytophthora infestans, which causes the devastating late blight disease in potato. The transcriptomes of three wild Solanum species native to Southern Sweden, Solanum dulcamara, Solanum nigrum, and Solanum physalifolium were compared to three potato clones, Desiree (cv.), SW93-1015 and Sarpo Mira. Desiree and S. physalifolium are susceptible to P. infestans whereas the other four have different degrees of resistance. By building transcript families based on de novo assembled RNA-seq across species and clones and correlating these to resistance phenotypes, we created a novel workflow to identify families with expanded or depleted number of transcripts in relation to the P. infestans resistance level. Analysis was facilitated by inferring functional annotations based on the family structure and semantic clustering. More transcript families were expanded in the resistant clones and species and the enriched functions of these were associated to expected gene ontology (GO) terms for resistance mechanisms such as hypersensitive response, host programmed cell death and endopeptidase activity. However, a number of unexpected functions and transcripts were also identified, for example transmembrane transport and protein acylation expanded in the susceptible group and a cluster of Zinc knuckle family proteins expanded in the resistant group. Over 400 expressed putative resistance (R-)genes were identified and resistant clones Sarpo Mira and SW93-1015 had ca 25% more expressed putative R-genes than susceptible cultivar Desiree. However, no differences in numbers of susceptibility (S-)gene homologs were seen between species and clones. In addition, we identified P. infestans transcripts including effectors in the early stages of P. infestans-Solanum interactions.

Genome-wide analysis of the gene families of resistance gene analogues in cotton and their response to Verticillium wilt

BACKGROUND

Gossypium raimondii is a Verticillium wilt-resistant cotton species whose genome encodes numerous disease resistance genes that play important roles in the defence against pathogens. However, the characteristics of resistance gene analogues (RGAs) and Verticillium dahliae response loci (VdRLs) have not been investigated on a global scale. In this study, the characteristics of RGA genes were systematically analysed using bioinformatics-driven methods. Moreover, the potential VdRLs involved in the defence response to Verticillium wilt were identified by RNA-seq and correlations with known resistance QTLs.

RESULTS

The G. raimondii genome encodes 1004 RGA genes, and most of these genes cluster in homology groups based on high levels of similarity. Interestingly, nearly half of the RGA genes occurred in 26 RGA-gene-rich clusters (Rgrcs). The homology analysis showed that sequence exchanges and tandem duplications frequently occurred within Rgrcs, and segmental duplications took place among the different Rgrcs. An RNA-seq analysis showed that the RGA genes play roles in cotton defence responses, forming 26 VdRLs inside in the Rgrcs after being inoculated with V. dahliae. A correlation analysis found that 12 VdRLs were adjacent to the known Verticillium wilt resistance QTLs, and that 5 were rich in NB-ARC domain-containing disease resistance genes.

CONCLUSIONS

The cotton genome contains numerous RGA genes, and nearly half of them are located in clusters, which evolved by sequence exchanges, tandem duplications and segmental duplications. In the Rgrcs, 26 loci were induced by the V. dahliae inoculation, and 12 are in the vicinity of known Verticillium wilt resistance QTLs.

Challenges of breeding potato cultivars to grow in various environments and to meet different demands

The potato (Solanum tuberosum L.) is cultivated all year round in Japan by using four types of cropping: summer and winter croppings, and double cropping in spring and fall. In each cropping season, growth conditions such as temperature, day length, and growing period, differ drastically; thus, different cultivars adapted to each environment are required. Breeding stations are located in both summer cropping areas and double cropping areas, and cultivars suitable for each cropping system are developed. The required cultivars differ according to cropping type and according to use such as table use, food processing, and starch production. The qualities necessary for each purpose differ and are therefore evaluated accordingly. Improvements in pest and disease resistance and in yield abilities are important as common breeding targets for all purposes. To develop potato cultivars that meet different needs, breeders have continued efforts to improve these traits. In this review, we introduce our approaches to developing new potato cultivars. We also discuss problems predicted in the future and introduce our efforts on broadening genetic diversity.

The evolution of WRKY transcription factors

BACKGROUND

The availability of increasing numbers of sequenced genomes has necessitated a re-evaluation of the evolution of the WRKY transcription factor family. Modern day plants descended from a charophyte green alga that colonized the land between 430 and 470 million years ago. The first charophyte genome sequence from Klebsormidium flaccidum filled a gap in the available genome sequences in the plant kingdom between unicellular green algae that typically have 1-3 WRKY genes and mosses that contain 30-40. WRKY genes have been previously found in non-plant species but their occurrence has been difficult to explain.

RESULTS

Only two WRKY genes are present in the Klebsormidium flaccidum genome and the presence of a Group IIb gene was unexpected because it had previously been thought that Group IIb WRKY genes first appeared in mosses. We found WRKY transcription factor genes outside of the plant lineage in some diplomonads, social amoebae, fungi incertae sedis, and amoebozoa. This patchy distribution suggests that lateral gene transfer is responsible. These lateral gene transfer events appear to pre-date the formation of the WRKY groups in flowering plants. Flowering plants contain proteins with domains typical for both resistance (R) proteins and WRKY transcription factors. R protein-WRKY genes have evolved numerous times in flowering plants, each type being restricted to specific flowering plant lineages. These chimeric proteins contain not only novel combinations of protein domains but also novel combinations and numbers of WRKY domains. Once formed, R protein WRKY genes may combine different components of signalling pathways that may either create new diversity in signalling or accelerate signalling by short circuiting signalling pathways.

CONCLUSIONS

We propose that the evolution of WRKY transcription factors includes early lateral gene transfers to non-plant organisms and the occurrence of algal WRKY genes that have no counterparts in flowering plants. We propose two alternative hypotheses of WRKY gene evolution: The "Group I Hypothesis" sees all WRKY genes evolving from Group I C-terminal WRKY domains. The alternative "IIa + b Separate Hypothesis" sees Groups IIa and IIb evolving directly from a single domain algal gene separate from the Group I-derived lineage.

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

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.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-743) contains supplementary material, which is available to authorized users.

Solanum resistance genes against Phytophthora infestans and their corresponding avirulence genes

Resistance genes against Phytophthora infestans (Rpi genes), the most important potato pathogen, are still highly valued in the breeding of Solanum spp. for enhanced resistance. The Rpi genes hitherto explored are localized most often in clusters, which are similar between the diverse Solanum genomes. Their distribution is not independent of late maturity traits. This review provides a summary of the most recent important revelations on the genomic position and cloning of Rpi genes, and the structure, associations, mode of action and activity spectrum of Rpi and corresponding avirulence (Avr) proteins. Practical implications for research into and application of Rpi genes are deduced and combined with an outlook on approaches to address remaining issues and interesting questions. It is evident that the potential of Rpi genes has not been exploited fully.

Copy number variation in potato - an asexually propagated autotetraploid species

Copy number variation (CNV) has been revealed as a significant contributor to the genetic variation in humans. Although CNV has been reported in several model animal and plant species, the presence of CNV and its biological impact in polyploid species has not yet been documented. We conducted a fluorescence in situ hybridization (FISH)-based CNV survey in potato, a vegetatively propagated autotetraploid species (2n = 4x = 48). We conducted FISH analysis using 18 randomly selected potato bacterial artificial chromosome (BAC) clones in a set of 16 potato cultivars with diverse breeding backgrounds. Six BACs (33%) with insert sizes of 137-145 kb were found to be associated with large CNV events detectable at the cytological level. We demonstrate that the large CNVs associated with two specific BACs (RH102I10 and RH83C08) were widespread among potato cultivars developed in North America and Europe. We measured the transcript abundance of four genes associated with the CNV spanned by BAC RH102I10. All four genes displayed a dosage effect in transcription. Although potato is vegetatively propagated, we observed that female gametes lacking the RH102I10-associated CNV were inferior to those with at least one copy of this CNV, indicating that the RH102I10-associated CNV can impact on the growth and development of the potato plants. Our results show that CNV is highly abundant in the potato genome and may play a significant role in genetic variation of this important food crop.

Frequent loss of lineages and deficient duplications accounted for low copy number of disease resistance genes in Cucurbitaceae

Background

The sequenced genomes of cucumber, melon and watermelon have relatively few R-genes, with 70, 75 and 55 copies only, respectively. The mechanism for low copy number of R-genes in Cucurbitaceae genomes remains unknown.

Results

Manual annotation of R-genes in the sequenced genomes of Cucurbitaceae species showed that approximately half of them are pseudogenes. Comparative analysis of R-genes showed frequent loss of R-gene loci in different Cucurbitaceae species. Phylogenetic analysis, data mining and PCR cloning using degenerate primers indicated that Cucurbitaceae has limited number of R-gene lineages (subfamilies). Comparison between R-genes from Cucurbitaceae and those from poplar and soybean suggested frequent loss of R-gene lineages in Cucurbitaceae. Furthermore, the average number of R-genes per lineage in Cucurbitaceae species is approximately 1/3 that in soybean or poplar. Therefore, both loss of lineages and deficient duplications in extant lineages accounted for the low copy number of R-genes in Cucurbitaceae. No extensive chimeras of R-genes were found in any of the sequenced Cucurbitaceae genomes. Nevertheless, one lineage of R-genes from Trichosanthes kirilowii, a wild Cucurbitaceae species, exhibits chimeric structures caused by gene conversions, and may contain a large number of distinct R-genes in natural populations.

Conclusions

Cucurbitaceae species have limited number of R-gene lineages and each genome harbors relatively few R-genes. The scarcity of R-genes in Cucurbitaceae species was due to frequent loss of R-gene lineages and infrequent duplications in extant lineages. The evolutionary mechanisms for large variation of copy number of R-genes in different plant species were discussed.

Impacts of resistance gene genetics, function, and evolution on a durable future

Studies on resistance gene function and evolution lie at the confluence of structural and molecular biology, genetics, and plant breeding. However, knowledge from these disparate fields has yet to be extensively integrated. This review draws on ideas and information from these different fields to elucidate the influences driving the evolution of different types of resistance genes in plants and the concurrent evolution of virulence in pathogens. It provides an overview of the factors shaping the evolution of recognition, signaling, and response genes in the context of emerging functional information along with a consideration of the new opportunities for durable resistance enabled by high-throughput DNA sequencing technologies.

Evolution of the Rdr1 TNL-cluster in roses and other Rosaceous species

Background

The resistance of plants to pathogens relies on two lines of defense: a basal defense response and a pathogen-specific system, in which resistance (R) genes induce defense reactions after detection of pathogen-associated molecular patterns (PAMPS). In the specific system, a so-called arms race has developed in which the emergence of new races of a pathogen leads to the diversification of plant resistance genes to counteract the pathogens’ effect. The mechanism of resistance gene diversification has been elucidated well for short-lived annual species, but data are mostly lacking for long-lived perennial and clonally propagated plants, such as roses. We analyzed the rose black spot resistance gene, Rdr1, in five members of the Rosaceae: Rosa multiflora, Rosa rugosa, Fragaria vesca (strawberry), Malus x domestica (apple) and Prunus persica (peach), and we present the deduced possible mechanism of R-gene diversification.

Results

We sequenced a 340.4-kb region from R. rugosa orthologous to the Rdr1 locus in R. multiflora. Apart from some deletions and rearrangements, the two loci display a high degree of synteny. Additionally, less pronounced synteny is found with an orthologous locus in strawberry but is absent in peach and apple, where genes from the Rdr1 locus are distributed on two different chromosomes. An analysis of 20 TIR-NBS-LRR (TNL) genes obtained from R. rugosa and R. multiflora revealed illegitimate recombination, gene conversion, unequal crossing over, indels, point mutations and transposable elements as mechanisms of diversification.

A phylogenetic analysis of 53 complete TNL genes from the five Rosaceae species revealed that with the exception of some genes from apple and peach, most of the genes occur in species-specific clusters, indicating that recent TNL gene diversification began prior to the split of Rosa from Fragaria in the Rosoideae and peach from apple in the Spiraeoideae and continued after the split in individual species. Sequence similarity of up to 99% is obtained between two R. multiflora TNL paralogs, indicating a very recent duplication.

Conclusions

The mechanisms by which TNL genes from perennial Rosaceae diversify are mainly similar to those from annual plant species. However, most TNL genes appear to be of recent origin, likely due to recent duplications, supporting the hypothesis that TNL genes in woody perennials are generally younger than those from annuals. This recent origin might facilitate the development of new resistance specificities, compensating for longer generation times in woody perennials.

Dynamic nucleotide-binding site and leucine-rich repeat-encoding genes in the grass family

The proper use of resistance genes (R genes) requires a comprehensive understanding of their genomics and evolution. We analyzed genes encoding nucleotide-binding sites and leucine-rich repeats in the genomes of rice (Oryza sativa), maize (Zea mays), sorghum (Sorghum bicolor), and Brachypodium distachyon. Frequent deletions and translocations of R genes generated prevalent presence/absence polymorphism between different accessions/species. The deletions were caused by unequal crossover, homologous repair, nonhomologous repair, or other unknown mechanisms. R gene loci identified from different genomes were mapped onto the chromosomes of rice cv Nipponbare using comparative genomics, resulting in an integrated map of 495 R loci. Sequence analysis of R genes from the partially sequenced genomes of an African rice cultivar and 10 wild accessions suggested that there are many additional R gene lineages in the AA genome of Oryza. The R genes with chimeric structures (termed type I R genes) are diverse in different rice accessions but only account for 5.8% of all R genes in the Nipponbare genome. In contrast, the vast majority of R genes in the rice genome are type II R genes, which are highly conserved in different accessions. Surprisingly, pseudogene-causing mutations in some type II lineages are often conserved, indicating that their conservations were not due to their functions. Functional R genes cloned from rice so far have more type II R genes than type I R genes, but type I R genes are predicted to contribute considerable diversity in wild species. Type I R genes tend to reduce the microsynteny of their flanking regions significantly more than type II R genes, and their flanking regions have slightly but significantly lower G/C content than those of type II R genes.

Evolution of a complex disease resistance gene cluster in diploid Phaseolus and tetraploid Glycine

We used a comparative genomics approach to investigate the evolution of a complex nucleotide-binding (NB)-leucine-rich repeat (LRR) gene cluster found in soybean (Glycine max) and common bean (Phaseolus vulgaris) that is associated with several disease resistance (R) genes of known function, including Rpg1b (for Resistance to Pseudomonas glycinea1b), an R gene effective against specific races of bacterial blight. Analysis of domains revealed that the amino-terminal coiled-coil (CC) domain, central nucleotide-binding domain (NB-ARC [for APAF1, Resistance genes, and CED4]), and carboxyl-terminal LRR domain have undergone distinct evolutionary paths. Sequence exchanges within the NB-ARC domain were rare. In contrast, interparalogue exchanges involving the CC and LRR domains were common, consistent with both of these regions coevolving with pathogens. Residues under positive selection were overrepresented within the predicted solvent-exposed face of the LRR domain, although several also were detected within the CC and NB-ARC domains. Superimposition of these latter residues onto predicted tertiary structures revealed that the majority are located on the surface, suggestive of a role in interactions with other domains or proteins. Following polyploidy in the Glycine lineage, NB-LRR genes have been preferentially lost from one of the duplicated chromosomes (homeologues found in soybean), and there has been partitioning of NB-LRR clades between the two homeologues. The single orthologous region in common bean contains approximately the same number of paralogues as found in the two soybean homeologues combined. We conclude that while polyploidization in Glycine has not driven a stable increase in family size for NB-LRR genes, it has generated two recombinationally isolated clusters, one of which appears to be in the process of decay.

Identification and localisation of the NB-LRR gene family within the potato genome

Background

The potato genome sequence derived from the Solanum tuberosum Group Phureja clone DM1-3 516 R44 provides unparalleled insight into the genome composition and organisation of this important crop. A key class of genes that comprises the vast majority of plant resistance (R) genes contains a nucleotide-binding and leucine-rich repeat domain, and is collectively known as NB-LRRs.

Results

As part of an effort to accelerate the process of functional R gene isolation, we performed an amino acid motif based search of the annotated potato genome and identified 438 NB-LRR type genes among the ~39,000 potato gene models. Of the predicted genes, 77 contain an N-terminal toll/interleukin 1 receptor (TIR)-like domain, and 107 of the remaining 361 non-TIR genes contain an N-terminal coiled-coil (CC) domain. Physical map positions were established for 370 predicted NB-LRR genes across all 12 potato chromosomes. The majority of NB-LRRs are physically organised within 63 identified clusters, of which 50 are homogeneous in that they contain NB-LRRs derived from a recent common ancestor.

Conclusions

By establishing the phylogenetic and positional relationship of potato NB-LRRs, our analysis offers significant insight into the evolution of potato R genes. Furthermore, the data provide a blueprint for future efforts to identify and more rapidly clone functional NB-LRR genes from Solanum species.

Identification and localisation of the NB-LRR gene family within the potato genome

BACKGROUND

The potato genome sequence derived from the Solanum tuberosum Group Phureja clone DM1-3 516 R44 provides unparalleled insight into the genome composition and organisation of this important crop. A key class of genes that comprises the vast majority of plant resistance (R) genes contains a nucleotide-binding and leucine-rich repeat domain, and is collectively known as NB-LRRs.

RESULTS

As part of an effort to accelerate the process of functional R gene isolation, we performed an amino acid motif based search of the annotated potato genome and identified 438 NB-LRR type genes among the ~39,000 potato gene models. Of the predicted genes, 77 contain an N-terminal toll/interleukin 1 receptor (TIR)-like domain, and 107 of the remaining 361 non-TIR genes contain an N-terminal coiled-coil (CC) domain. Physical map positions were established for 370 predicted NB-LRR genes across all 12 potato chromosomes. The majority of NB-LRRs are physically organised within 63 identified clusters, of which 50 are homogeneous in that they contain NB-LRRs derived from a recent common ancestor.

CONCLUSIONS

By establishing the phylogenetic and positional relationship of potato NB-LRRs, our analysis offers significant insight into the evolution of potato R genes. Furthermore, the data provide a blueprint for future efforts to identify and more rapidly clone functional NB-LRR genes from Solanum species.

A hybrid BAC physical map of potato: a framework for sequencing a heterozygous genome

BACKGROUND

Potato is the world's third most important food crop, yet cultivar improvement and genomic research in general remain difficult because of the heterozygous and tetraploid nature of its genome. The development of physical map resources that can facilitate genomic analyses in potato has so far been very limited. Here we present the methods of construction and the general statistics of the first two genome-wide BAC physical maps of potato, which were made from the heterozygous diploid clone RH89-039-16 (RH).

RESULTS

First, a gel electrophoresis-based physical map was made by AFLP fingerprinting of 64478 BAC clones, which were aligned into 4150 contigs with an estimated total length of 1361 Mb. Screening of BAC pools, followed by the KeyMaps in silico anchoring procedure, identified 1725 AFLP markers in the physical map, and 1252 BAC contigs were anchored the ultradense potato genetic map. A second, sequence-tag-based physical map was constructed from 65919 whole genome profiling (WGP) BAC fingerprints and these were aligned into 3601 BAC contigs spanning 1396 Mb. The 39733 BAC clones that overlap between both physical maps provided anchors to 1127 contigs in the WGP physical map, and reduced the number of contigs to around 2800 in each map separately. Both physical maps were 1.64 times longer than the 850 Mb potato genome. Genome heterozygosity and incomplete merging of BAC contigs are two factors that can explain this map inflation. The contig information of both physical maps was united in a single table that describes hybrid potato physical map.

CONCLUSIONS

The AFLP physical map has already been used by the Potato Genome Sequencing Consortium for sequencing 10% of the heterozygous genome of clone RH on a BAC-by-BAC basis. By layering a new WGP physical map on top of the AFLP physical map, a genetically anchored genome-wide framework of 322434 sequence tags has been created. This reference framework can be used for anchoring and ordering of genomic sequences of clone RH (and other potato genotypes), and opens the possibility to finish sequencing of the RH genome in a more efficient way via high throughput next generation approaches.

Cloning and characterization of r3b; members of the r3 superfamily of late blight resistance genes show sequence and functional divergence

Massive resistance (R) gene stacking is considered to be one of the most promising approaches to provide durable resistance to potato late blight for both conventional and genetically modified breeding strategies. The R3 complex locus on chromosome XI in potato is an example of natural R gene stacking, because it contains two closely linked R genes (R3a and R3b) with distinct resistance specificities to Phytophthora infestans. Here, we report about the positional cloning of R3b. Both transient and stable transformations of susceptible tobacco and potato plants showed that R3b conferred full resistance to incompatible P. infestans isolates. R3b encodes a coiled-coil nucleotide-binding site leucine-rich repeat protein and exhibits 82% nucleotide identity with R3a located in the same R3 cluster. The R3b gene specifically recognizes Avr3b, a newly identified avirulence factor from P. infestans. R3b does not recognize Avr3a, the corresponding avirulence gene for R3a, showing that, despite their high sequence similarity, R3b and R3a have clearly distinct recognition specificities. In addition to the Rpi-mcd1/Rpi-blb3 locus on chromosome IV, the R3 locus on chromosome XI is the second example of an R-gene cluster with multiple genes recognizing different races of P. infestans.

A genome-wide genetic map of NB-LRR disease resistance loci in potato

Like all plants, potato has evolved a surveillance system consisting of a large array of genes encoding for immune receptors that confer resistance to pathogens and pests. The majority of these so-called resistance or R proteins belong to the super-family that harbour a nucleotide binding and a leucine-rich-repeat domain (NB-LRR). Here, sequence information of the conserved NB domain was used to investigate the genome-wide genetic distribution of the NB-LRR resistance gene loci in potato. We analysed the sequences of 288 unique BAC clones selected using filter hybridisation screening of a BAC library of the diploid potato clone RH89-039-16 (S. tuberosum ssp. tuberosum) and a physical map of this BAC library. This resulted in the identification of 738 partial and full-length NB-LRR sequences. Based on homology of these sequences with known resistance genes, 280 and 448 sequences were classified as TIR-NB-LRR (TNL) and CC-NB-LRR (CNL) sequences, respectively. Genetic mapping revealed the presence of 15 TNL and 32 CNL loci. Thirty-six are novel, while three TNL loci and eight CNL loci are syntenic with previously identified functional resistance genes. The genetic map was complemented with 68 universal CAPS markers and 82 disease resistance trait loci described in literature, providing an excellent template for genetic studies and applied research in potato.

Genome sequence and analysis of the tuber crop potato

Potato (Solanum tuberosum L.) is the world's most important non-grain food crop and is central to global food security. It is clonally propagated, highly heterozygous, autotetraploid, and suffers acute inbreeding depression. Here we use a homozygous doubled-monoploid potato clone to sequence and assemble 86% of the 844-megabase genome. We predict 39,031 protein-coding genes and present evidence for at least two genome duplication events indicative of a palaeopolyploid origin. As the first genome sequence of an asterid, the potato genome reveals 2,642 genes specific to this large angiosperm clade. We also sequenced a heterozygous diploid clone and show that gene presence/absence variants and other potentially deleterious mutations occur frequently and are a likely cause of inbreeding depression. Gene family expansion, tissue-specific expression and recruitment of genes to new pathways contributed to the evolution of tuber development. The potato genome sequence provides a platform for genetic improvement of this vital crop.

Profiling of short RNAs during fleshy fruit development reveals stage-specific sRNAome expression patterns

Plants feature a particularly diverse population of short (s)RNAs, the central component of all RNA silencing pathways. Next generation sequencing techniques enable deeper insights into this complex and highly conserved mechanism and allow identification and quantification of sRNAs. We employed deep sequencing to monitor the sRNAome of developing tomato fruits covering the period between closed flowers and ripened fruits by profiling sRNAs at 10 time-points. It is known that microRNAs (miRNAs) play an important role in development but very little information is available about the majority of sRNAs that are not miRNAs. Here we show distinctive patterns of sRNA expression that often coincide with stages of the developmental process such as flowering, early and late fruit maturation. Moreover, thousands of non-miRNA sRNAs are differentially expressed during fruit development and ripening. Some of these differentially expressed sRNAs derived from transposons but many derive from protein coding genes or regions that show homology to protein coding genes, several of which are known to play a role in flower and fruit development. These findings raise the possibility of a regulative role of these sRNAs during fruit onset and maturation in a crop species. We also identified six new miRNAs and experimentally validated two target mRNAs. These two mRNAs are targeted by the same miRNA but do not belong to the same gene family, which is rare for plant miRNAs. Expression pattern and putative function of these targets indicate a possible role in glutamate accumulation, which contributes to establishing the taste of the fruit.

Organization and molecular evolution of a disease-resistance gene cluster in coffee trees

BACKGROUND

Most disease-resistance (R) genes in plants encode NBS-LRR proteins and belong to one of the largest and most variable gene families among plant genomes. However, the specific evolutionary routes of NBS-LRR encoding genes remain elusive. Recently in coffee tree (Coffea arabica), a region spanning the SH3 locus that confers resistance to coffee leaf rust, one of the most serious coffee diseases, was identified and characterized. Using comparative sequence analysis, the purpose of the present study was to gain insight into the genomic organization and evolution of the SH3 locus.

RESULTS

Sequence analysis of the SH3 region in three coffee genomes, Ea and Ca subgenomes from the allotetraploid C. arabica and Cc genome from the diploid C. canephora, revealed the presence of 5, 3 and 4 R genes in Ea, Ca, and Cc genomes, respectively. All these R-gene sequences appeared to be members of a CC-NBS-LRR (CNL) gene family that was only found at the SH3 locus in C. arabica. Furthermore, while homologs were found in several dicot species, comparative genomic analysis failed to find any CNL R-gene in the orthologous regions of other eudicot species. The orthology relationship among the SH3-CNL copies in the three analyzed genomes was determined and the duplication/deletion events that shaped the SH3 locus were traced back. Gene conversion events were detected between paralogs in all three genomes and also between the two sub-genomes of C. arabica. Significant positive selection was detected in the solvent-exposed residues of the SH3-CNL copies.

CONCLUSION

The ancestral SH3-CNL copy was inserted in the SH3 locus after the divergence between Solanales and Rubiales lineages. Moreover, the origin of most of the SH3-CNL copies predates the divergence between Coffea species. The SH3-CNL family appeared to evolve following the birth-and-death model, since duplications and deletions were inferred in the evolution of the SH3 locus. Gene conversion between paralog members, inter-subgenome sequence exchanges and positive selection appear to be the major forces acting on the evolution of SH3-CNL in coffee trees.

Comparative sequence analysis of the potato cyst nematode resistance locus H1 reveals a major lack of co-linearity between three haplotypes in potato (Solanum tuberosum ssp.)

The H1 locus confers resistance to the potato cyst nematode Globodera rostochiensis pathotypes 1 and 4. It is positioned at the distal end of chromosome V of the diploid Solanum tuberosum genotype SH83-92-488 (SH) on an introgression segment derived from S. tuberosum ssp. andigena. Markers from a high-resolution genetic map of the H1 locus (Bakker et al. in Theor Appl Genet 109:146-152, 2004) were used to screen a BAC library to construct a physical map covering a 341-kb region of the resistant haplotype coming from SH. For comparison, physical maps were also generated of the two haplotypes from the diploid susceptible genotype RH89-039-16 (S. tuberosum ssp. tuberosum/S. phureja), spanning syntenic regions of 700 and 319 kb. Gene predictions on the genomic segments resulted in the identification of a large cluster consisting of variable numbers of the CC-NB-LRR type of R genes for each haplotype. Furthermore, the regions were interspersed with numerous transposable elements and genes coding for an extensin-like protein and an amino acid transporter. Comparative analysis revealed a major lack of gene order conservation in the sequences of the three closely related haplotypes. Our data provide insight in the evolutionary mechanisms shaping the H1 locus and will facilitate the map-based cloning of the H1 resistance gene.

Mining disease-resistance genes in roses: functional and molecular characterization of the rdr1 locus

The interaction of roses with the leaf spot pathogen Diplocarpon rosae (the cause of black spot on roses) is an interesting pathosystem because it involves a long-lived woody perennial, with life history traits very different from most model plants, and a hemibiotrophic pathogen with moderate levels of gene flow. Here we present data on the molecular structure of the first monogenic dominant resistance gene from roses, Rdr1, directed against one isolate of D. rosae. Complete sequencing of the locus carrying the Rdr1 gene resulted in a sequence of 265,477 bp with a cluster of nine highly related TIR-NBS-LRR (TNL) candidate genes. After sequencing revealed candidate genes for Rdr1, we implemented a gene expression analysis and selected five genes out of the nine TNLs. We then silenced the whole TNL gene family using RNAi (Rdr1-RNAi) constructed from the most conserved sequence region and demonstrated a loss of resistance in the normally resistant genotype. To identify the functional TNL gene, we further screened the five TNL candidate genes with a transient leaf infiltration assay. The transient expression assay indicated a single TNL gene (muRdr1H), partially restoring resistance in the susceptible genotype. Rdr1 was found to localize within the muRdr1 gene family; the genes within this locus contain characteristic motifs of active TNL genes and belong to a young cluster of R genes. The transient leaf assay can be used to further analyze the rose black spot interaction and its evolution, extending the analyses to additional R genes and to additional pathogenic types of the pathogen.

Understanding and exploiting late blight resistance in the age of effectors

Potato (Solanum tuberosum) is the world's third-largest food crop. It severely suffers from late blight, a devastating disease caused by Phytophthora infestans. This oomycete pathogen secretes host-translocated RXLR effectors that include avirulence (AVR) proteins, which are targeted by resistance (R) proteins from wild Solanum species. Most Solanum R genes appear to have coevolved with P. infestans at its center of origin in central Mexico. Various R and Avr genes were recently cloned, and here we catalog characterized R-AVR pairs. We describe the mechanisms that P. infestans employs for evading R protein recognition and discuss partial resistance and partial virulence phenotypes in the context of our knowledge of effector diversity and activity. Genome-wide catalogs of P. infestans effectors are available, enabling effectoromics approaches that accelerate R gene cloning and specificity profiling. Engineering R genes with expanded pathogen recognition has also become possible. Importantly, monitoring effector allelic diversity in pathogen populations can assist in R gene deployment in agriculture.

Mechanisms of haplotype divergence at the RGA08 nucleotide-binding leucine-rich repeat gene locus in wild banana (Musa balbisiana)

BACKGROUND

Comparative sequence analysis of complex loci such as resistance gene analog clusters allows estimating the degree of sequence conservation and mechanisms of divergence at the intraspecies level. In banana (Musa sp.), two diploid wild species Musa acuminata (A genome) and Musa balbisiana (B genome) contribute to the polyploid genome of many cultivars. The M. balbisiana species is associated with vigour and tolerance to pests and disease and little is known on the genome structure and haplotype diversity within this species. Here, we compare two genomic sequences of 253 and 223 kb corresponding to two haplotypes of the RGA08 resistance gene analog locus in M. balbisiana "Pisang Klutuk Wulung" (PKW).

RESULTS

Sequence comparison revealed two regions of contrasting features. The first is a highly colinear gene-rich region where the two haplotypes diverge only by single nucleotide polymorphisms and two repetitive element insertions. The second corresponds to a large cluster of RGA08 genes, with 13 and 18 predicted RGA genes and pseudogenes spread over 131 and 152 kb respectively on each haplotype. The RGA08 cluster is enriched in repetitive element insertions, in duplicated non-coding intergenic sequences including low complexity regions and shows structural variations between haplotypes. Although some allelic relationships are retained, a large diversity of RGA08 genes occurs in this single M. balbisiana genotype, with several RGA08 paralogs specific to each haplotype. The RGA08 gene family has evolved by mechanisms of unequal recombination, intragenic sequence exchange and diversifying selection. An unequal recombination event taking place between duplicated non-coding intergenic sequences resulted in a different RGA08 gene content between haplotypes pointing out the role of such duplicated regions in the evolution of RGA clusters. Based on the synonymous substitution rate in coding sequences, we estimated a 1 million year divergence time for these M. balbisiana haplotypes.

CONCLUSIONS

A large RGA08 gene cluster identified in wild banana corresponds to a highly variable genomic region between haplotypes surrounded by conserved flanking regions. High level of sequence identity (70 to 99%) of the genic and intergenic regions suggests a recent and rapid evolution of this cluster in M. balbisiana.

Comparison of the chromosome maps around a resistance hot spot on chromosome 5 of potato and tomato using BAC-FISH painting

Potato chromosome 5 harbours numerous genes for important qualitative and quantitative traits, such as resistance to the root cyst nematode Globodera pallida and the late blight fungus, Phytophthora infestans. The genes make up part of a "hot spot" for resistances to various pathogens covering a genetic map length of 3 cM between markers GP21 and GP179. We established the physical size and position of this region on chromosome 5 in potato and tomato using fluorescence in situ hybridization (FISH) on pachytene chromosomes. Five potato bacterial artificial chromosome (BAC) clones with the genetically anchored markers GP21, R1-contig (proximal end), CosA, GP179, and StPto were selected, labeled with different fluorophores, and hybridized in a five-colour FISH experiment. Our results showed the location of the BAC clones in the middle of the long arm of chromosome 5 in both potato and tomato. Based on chromosome measurements, we estimate the physical size of the GP21-GP179 interval at 0.85 Mb and 1.2 Mb in potato and tomato, respectively. The GP21-GP179 interval is part of a genome segment known to have inverted map positions between potato and tomato.

Salicylic acid is involved in the Nb-mediated defense responses to Potato virus X in Solanum tuberosum

To evaluate the role of salicylic acid (SA) in Nb-mediated hypersensitive resistance to Potato virus X (PVX) avirulent strain ROTH1 in Solanum tuberosum, we have constructed SA-deficient transgenic potato plant lines by overexpressing the bacterial enzyme salicylate hydroxylase (NahG), which degrades SA. Evaluation of these transgenic lines revealed hydrogen peroxide accumulation and spontaneous lesion formation in an age- and light-dependent manner. In concordance, NahG potato plants were more sensitive to treatment with methyl viologen, a reactive oxygen species-generating compound. In addition, when challenged with PVX ROTH1, NahG transgenic lines showed a decreased disease-resistance response to infection and were unable to induce systemic acquired resistance. However, the avirulent viral effector, the PVX 25-kDa protein, does induce expression of the pathogenesis-related gene PR-1a in NahG potato plants. Taken together, our data indicate that SA is involved in local and systemic defense responses mediated by the Nb gene in Solanum tuberosum. This is the first report to show that basal levels of SA correlate with hypersensitive resistance to PVX.

Phytophthora infestans isolates lacking class I ipiO variants are virulent on Rpi-blb1 potato

A strategy to control the devastating late blight disease is providing potato cultivars with genes that are effective in resistance to a broad spectrum of Phytophthora infestans isolates. Thus far, most late blight resistance (R) genes that were introgressed in potato were quickly defeated. In contrast, the Rpi-blb1 gene originating from Solanum bulbocastanum has performed as an exclusive broad-spectrum R gene for many years. Recently, the RXLR effector family ipiO was identified to contain Avr-blb1. Monitoring the genetic diversity of the ipiO family in a large set of isolates of P. infestans and related species resulted in 16 ipiO variants in three distinct classes. Class I and class II but not class III ipiO variants induce cell death when coinfiltrated with Rpi-blb1 in Nicotiana benthamiana. Class I is highly diverse and is represented in all analyzed P. infestans isolates except two Mexican P. infestans isolates, and these were found virulent on Rpi-blb1 plants. In its C-terminal domain, IPI-O contains a W motif that is essential for triggering Rpi-blb1-mediated cell death and is under positive selection. This study shows that profiling the variation of Avr-blb1 within a P. infestans population is instrumental for predicting the effectiveness of Rpi-blb1-mediated resistance in potato.

In planta expression screens of Phytophthora infestans RXLR effectors reveal diverse phenotypes, including activation of the Solanum bulbocastanum disease resistance protein Rpi-blb2

The Irish potato famine pathogen Phytophthora infestans is predicted to secrete hundreds of effector proteins. To address the challenge of assigning biological functions to computationally predicted effector genes, we combined allele mining with high-throughput in planta expression. We developed a library of 62 infection-ready P. infestans RXLR effector clones, obtained using primer pairs corresponding to 32 genes and assigned activities to several of these genes. This approach revealed that 16 of the 62 examined effectors cause phenotypes when expressed inside plant cells. Besides the well-studied AVR3a effector, two additional effectors, PexRD8 and PexRD36(45-1), suppressed the hypersensitive cell death triggered by the elicitin INF1, another secreted protein of P. infestans. One effector, PexRD2, promoted cell death in Nicotiana benthamiana and other solanaceous plants. Finally, two families of effectors induced hypersensitive cell death specifically in the presence of the Solanum bulbocastanum late blight resistance genes Rpi-blb1 and Rpi-blb2, thereby exhibiting the activities expected for Avrblb1 and Avrblb2. The AVRblb2 family was then studied in more detail and found to be highly variable and under diversifying selection in P. infestans. Structure-function experiments indicated that a 34-amino acid region in the C-terminal half of AVRblb2 is sufficient for triggering Rpi-blb2 hypersensitivity and that a single positively selected AVRblb2 residue is critical for recognition by Rpi-blb2.

Major-effect QTLs for stem and foliage resistance to late blight in the wild potato relatives Solanum sparsipilum and S. spegazzinii are mapped to chromosome X

To find out new resistance sources to late blight in the wild germplasm for potato breeding, we examined the polygenic resistance of Solanum sparsipilum and S. spegazzinii by a quantitative trait locus (QTL) analysis. We performed stem and foliage tests under controlled conditions in two diploid mapping progenies. Four traits were selected for QTL detection. A total of 30 QTLs were mapped, with a large-effect QTL region on chromosome X detected in both potato relatives. The mapping of literature-derived markers highlighted colinearities with published late blight QTLs or R-genes. Results showed (a) the resistance potential of S. sparsipilum and S. spegazzinii for late blight control, and (b) the efficacy of the stem test as a complement to the foliage test to break down the complex late blight resistance into elementary components. The relationships of late blight resistance QTLs with R-genes and maturity QTLs are discussed.

The zig-zag-zig in oomycete-plant interactions

In addition to a range of preformed barriers, plants defend themselves against microbial invasion by detecting conserved, secreted molecules, called pathogen-associated molecular patterns (PAMPs). PAMP-triggered immunity (PTI) is the first inducible layer of plant defence that microbial pathogens must navigate by the delivery of effector proteins that act to suppress or otherwise manipulate key components of resistance. Effectors may themselves be targeted by a further layer of defence, effector-triggered immunity (ETI), as their presence inside or outside host cells may be detected by resistance proteins. This 'zig-zag-zig' of tightly co-evolving molecular interactions determines the outcome of attempted infection. In this article, we consider the complex molecular interplay between plants and plant pathogenic oomycetes, drawing on recent literature to illustrate what is known about oomycete PAMPs and elicitors of defence responses, the effectors they utilize to suppress PTI, and the phenomenal molecular 'battle' between effector and resistance (R) genes that dictates the establishment or evasion of ETI.

Mapping of STS markers obtained from oat resistance gene analog sequences

Two previously isolated resistance gene analogs (RGAs) of oat have been located as RFLPs in the reference map of Avena byzantina ‘Kanota’ × Avena sativa ‘Ogle’ in regions either homologous or homoeologous to loci for resistance to Puccinia coronata , the causal agent of crown rust. In this study, the RGAs were mapped in two recombinant inbred line (RIL) populations that segregate for crown rust resistance: the diploid Avena strigosa × Avena wiestii RIL population (Asw), which has been used for mapping the complex locus PcA, and the hexaploid MN841801-1 × Noble-2 RIL population (MN), in which QTLs have been located. To obtain single-locus markers, RGAs were converted to sequence tagged site (STS) markers using a procedure involving extension of the original RGA sequence lengths by PCR genome walking, amplification and cloning of the parental fragments, and identification of single nucleotide polymorphisms. The procedure successfully obtained STSs from different members of the L7M2 family of sequences, the initial NBS of which have nucleotide similarities of >83%. However, for RGA III2.18, the parental lines were not polymorphic for the STSs assayed. A sequence characterized amplified region (SCAR) marker with features of an RGA had been previously identified for gene Pc94. This marker was also mapped in the above RIL populations. Markers based on RGA L7M2 co-localized with markers defining the QTL Prq1a in linkage group MN3, and were located 15.2 cM from PcA in linkage group AswAC. The SCAR marker for Pc94 was also located in the QTL Prq1a but at 39.5 cM from PcA in AswAC, indicating that the NBS-LRR sequence represented by this marker is not related to PcA. L7M2 was also excluded as a member of the PcA cluster, although it could be an appropriate marker for the Prq1a cluster if chromosome rearrangements are postulated.

A high-resolution map of the Grp1 locus on chromosome V of potato harbouring broad-spectrum resistance to the cyst nematode species Globodera pallida and Globodera rostochiensis

The Grp1 locus confers broad-spectrum resistance to the potato cyst nematode species Globodera pallida and Globodera rostochiensis and is located in the GP21-GP179 interval on the short arm of chromosome V of potato. A high-resolution map has been developed using the diploid mapping population RHAM026, comprising 1,536 genotypes. The flanking markers GP21 and GP179 have been used to screen the 1,536 genotypes for recombination events. Interval mapping of the resistances to G. pallida Pa2 and G. rostochiensis Ro5 resulted in two nearly identical LOD graphs with the highest LOD score just north of marker TG432. Detailed analysis of the 44 recombinant genotypes showed that G. pallida and G. rostochiensis resistance could not be separated and map to the same location between marker SPUD838 and TG432. It is suggested that the quantitative resistance to both nematode species at the Grp1 locus is mediated by one or more tightly linked R genes that might belong to the NBS-LRR class.

The frequent transcriptional readthrough of the tobacco Tnt1 retrotransposon and its possible implications for the control of resistance genes

Retrotransposons are a major component of eukaryote genomes, being especially abundant in plant genomes. They are frequently found inserted in gene-rich regions and have greatly contributed to the evolution of gene coding capacity and regulation. Retrotransposon insertions can influence the expression of neighboring genes in many ways, such as modifying their promoter or terminator sequences and altering their epigenetic control. Plant retrotransposons are highly regulated and their expression is usually associated with stress situations. While the control of transcription of some plant retrotransposons has been analyzed in some detail, little is known about the transcriptional termination of these elements. Here we show that the transcripts of the tobacco retrotransposon Tnt1 display a high variability of polyadenylation sites, only a fraction of them terminating at the major termination site. We also report on the ability of Tnt1 to extend its transcription into flanking genomic sequences and we analyze a particular case in which Tnt1 transcripts include sequences of an oppositely oriented resistance-like gene. The expression of this gene and the neighboring Tnt1 copy generate transcripts overlapping in more that 800 nucleotides, which could anneal and form dsRNAs and enter into silencing regulatory pathways. Resistance gene loci are usually composed of tandem arrays of resistance-like genes, a number of which contain mutations, including retrotransposon insertions, and are considered as to be pseudogenes. Given that plant retrotransposons are usually regulated by stress, the convergent expression of these resistance-like pseudogenes and the interleaving inducible retrotransposons may contribute to the control of plant responses to stress.