The disease resistance gene Dm3 is infrequent in natural populations of Lactuca serriola due to deletions and frequent gene conversions at the RGC2 locus

Creation and judicious application of a wheat resistance gene atlas

Disease-resistance (R) gene cloning in wheat (Triticum aestivum) has been accelerated by the recent surge of genomic resources, facilitated by advances in sequencing technologies and bioinformatics. However, with the challenges of population growth and climate change, it is vital not only to clone and functionally characterize a few handfuls of R genes, but also to do so at a scale that would facilitate the breeding and deployment of crops that can recognize the wide range of pathogen effectors that threaten agroecosystems. Pathogen populations are continually changing, and breeders must have tools and resources available to rapidly respond to those changes if we are to safeguard our daily bread. To meet this challenge, we propose the creation of a wheat R-gene atlas by an international community of researchers and breeders. The atlas would consist of an online directory from which sources of resistance could be identified and deployed to achieve more durable resistance to the major wheat pathogens, such as wheat rusts, blotch diseases, powdery mildew, and wheat blast. We present a costed proposal detailing how the interacting molecular components governing disease resistance could be captured from both the host and the pathogen through biparental mapping, mutational genomics, and whole-genome association genetics. We explore options for the configuration and genotyping of diversity panels of hexaploid and tetraploid wheat, as well as their wild relatives and major pathogens, and discuss how the atlas could inform a dynamic, durable approach to R-gene deployment. Set against the current magnitude of wheat yield losses worldwide, recently estimated at 21%, this endeavor presents one route for bringing R genes from the lab to the field at a considerable speed and quantity.

Proteomics Analysis Indicates Greater Abundance of Proteins Involved in Major Metabolic Pathways in Lactuca sativa cv. Salinas than Lactuca serriola Accession US96UC23

Lettuce (Lactuca sativa), cultivated mainly for its edible leaves and stems, is an important vegetable crop worldwide. Genomes of cultivated lettuce (L. sativa cv. Salinas) and its wild relative L. serriola accession US96UC23 are sequenced, but a clear understanding of the genetic basis for divergence in phenotypes of the two species is lacking. Tandem mass tag (TMT) based mass spectrometry is used to quantitatively compare protein levels between these two species. Four-day old seedlings is transplanted into 500 mL pots filled with soil. Plants are grown for 8 weeks under 250 µmol m-2 sec-1 continuous light, 20 °C and relative humidity between 50-70%. Leaf discs (1 cm diameter) from three individuals per biological replicate are analyzed. A total of 3000 proteins are identified, of which the levels of 650 are significantly different between 'Salinas' and US96UC23. Pathway analysis indicated a higher flux of carbon in 'Salinas' than US96UC23. Many essential metabolic pathways such as tetrapyrrole metabolism and fatty acid biosynthesis are upregulated in 'Salinas' compared with US96UC23. This study provides a reference proteome for researchers interested in understanding lettuce biology and improving traits for cultivation.

Genome-Wide Architecture of Disease Resistance Genes in Lettuce

Genome-wide motif searches identified 1134 genes in the lettuce reference genome of cv. Salinas that are potentially involved in pathogen recognition, of which 385 were predicted to encode nucleotide binding-leucine rich repeat receptor (NLR) proteins. Using a maximum-likelihood approach, we grouped the NLRs into 25 multigene families and 17 singletons. Forty-one percent of these NLR-encoding genes belong to three families, the largest being RGC16 with 62 genes in cv. Salinas. The majority of NLR-encoding genes are located in five major resistance clusters (MRCs) on chromosomes 1, 2, 3, 4, and 8 and cosegregate with multiple disease resistance phenotypes. Most MRCs contain primarily members of a single NLR gene family but a few are more complex. MRC2 spans 73 Mb and contains 61 NLRs of six different gene families that cosegregate with nine disease resistance phenotypes. MRC3, which is 25 Mb, contains 22 RGC21 genes and colocates with Dm13. A library of 33 transgenic RNA interference tester stocks was generated for functional analysis of NLR-encoding genes that cosegregated with disease resistance phenotypes in each of the MRCs. Members of four NLR-encoding families, RGC1, RGC2, RGC21, and RGC12 were shown to be required for 16 disease resistance phenotypes in lettuce. The general composition of MRCs is conserved across different genotypes; however, the specific repertoire of NLR-encoding genes varied particularly of the rapidly evolving Type I genes. These tester stocks are valuable resources for future analyses of additional resistance phenotypes.

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.

Specific in planta recognition of two GKLR proteins of the downy mildew Bremia lactucae revealed in a large effector screen in lettuce

Breeding lettuce (Lactuca sativa) for resistance to the downy mildew pathogen Bremia lactucae is mainly achieved by introgression of dominant downy mildew resistance (Dm) genes. New Bremia races quickly render Dm genes ineffective, possibly by mutation of recognized host-translocated effectors or by suppression of effector-triggered immunity. We have previously identified 34 potential RXLR(-like) effector proteins of B. lactucae that were here tested for specific recognition within a collection of 129 B. lactucae-resistant Lactuca lines. Two effectors triggered a hypersensitive response: BLG01 in 52 lines, predominantly L. saligna, and BLG03 in two L. sativa lines containing Dm2 resistance. The N-terminal sequences of BLG01 and BLG03, containing the signal peptide and GKLR variant of the RXLR translocation motif, are not required for in planta recognition but function in effector delivery. The locus responsible for BLG01 recognition maps to the bottom of lettuce chromosome 9, whereas recognition of BLG03 maps in the RGC2 cluster on chromosome 2. Lactuca lines that recognize the BLG effectors are not resistant to Bremia isolate Bl:24 that expresses both BLG genes, suggesting that Bl:24 can suppress the triggered immune responses. In contrast, lettuce segregants displaying Dm2-mediated resistance to Bremia isolate Bl:5 are responsive to BLG03, suggesting that BLG03 is a candidate Avr2 protein.

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.

Mutations in lettuce improvement

Lettuce is a major vegetable in western countries. Mutations generated genetic variations and played an important role in the domestication of the crop. Many traits derived from natural and induced mutations, such as dwarfing, early flowering, male sterility, and chlorophyll deficiency, are useful in physiological and genetic studies. Mutants were also used to develop new lettuce products including miniature and herbicide-tolerant cultivars. Mutant analysis was critical in lettuce genomic studies including identification and cloning of disease-resistance genes. Mutagenesis combined with genomic technology may provide powerful tools for the discovery of novel gene alleles. In addition to radiation and chemical mutagens, unconventional approaches such as tissue or protoplast culture, transposable elements, and space flights have been utilized to generate mutants in lettuce. Since mutation breeding is considered nontransgenic, it is more acceptable to consumers and will be explored more in the future for lettuce improvement.

Natural diversity in flowering responses of Arabidopsis thaliana caused by variation in a tandem gene array

Tandemly arrayed genes that belong to gene families characterize genomes of many organisms. Gene duplication and subsequent relaxation of selection can lead to the establishment of paralogous cluster members that may evolve along different trajectories. Here, we report on the structural variation in MADS AFFECTING FLOWERING 2 (MAF2) gene, one member of the tandemly duplicated cluster of MADS-box-containing transcription factors in Arabidopsis thaliana. The altered gene structure at the MAF2 locus is present as a moderate-frequency polymorphism in Arabidopsis and leads to the extensive diversity in transcript patterns due to alternative splicing. Rearrangements at the MAF2 locus are associated with an early flowering phenotype in BC(5) lines. The lack of suppression of flowering time in a MAF2-insertion line expressing the MAF2-specific artificial miRNA suggests that these MAF2 variants are behaving as loss-of-function alleles. The variation in gene architecture is also associated with segregation distortion, which may have facilitated the spread and the establishment of the corresponding alleles throughout the Eurasian range of the A. thaliana population.

Comparative large-scale analysis of interactions between several crop species and the effector repertoires from multiple pathovars of Pseudomonas and Ralstonia

Bacterial plant pathogens manipulate their hosts by injection of numerous effector proteins into host cells via type III secretion systems. Recognition of these effectors by the host plant leads to the induction of a defense reaction that often culminates in a hypersensitive response manifested as cell death. Genes encoding effector proteins can be exchanged between different strains of bacteria via horizontal transfer, and often individual strains are capable of infecting multiple hosts. Host plant species express diverse repertoires of resistance proteins that mediate direct or indirect recognition of bacterial effectors. As a result, plants and their bacterial pathogens should be considered as two extensive coevolving groups rather than as individual host species coevolving with single pathovars. To dissect the complexity of this coevolution, we cloned 171 effector-encoding genes from several pathovars of Pseudomonas and Ralstonia. We used Agrobacterium tumefaciens-mediated transient assays to test the ability of each effector to induce a necrotic phenotype on 59 plant genotypes belonging to four plant families, including numerous diverse accessions of lettuce (Lactuca sativa) and tomato (Solanum lycopersicum). Known defense-inducing effectors (avirulence factors) and their homologs commonly induced extensive necrosis in many different plant species. Nonhost species reacted to multiple effector proteins from an individual pathovar more frequently and more intensely than host species. Both homologous and sequence-unrelated effectors could elicit necrosis in a similar spectrum of plants, suggesting common effector targets or targeting of the same pathways in the plant cell.

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.

Independent evolution of functional Pm3 resistance genes in wild tetraploid wheat and domesticated bread wheat

The Pm3 alleles of cultivated bread wheat confer gene for gene resistance to the powdery mildew fungus. They represent a particular case of plant disease resistance gene evolution, because of their recent origin and possible evolution after the formation of hexaploid wheat. The Pm3 locus is conserved in tetraploid wheat, thereby allowing the comparative evolutionary study of the same resistance locus in a domesticated species and in one of its wild ancestors. We have identified 61 Pm3 allelic sequences from wild and domesticated tetraploid wheat subspecies. The Pm3 sequences corresponded to 24 different haplotypes. They showed low sequence diversity, differing by only a few polymorphic sequence blocks that were further reshuffled between alleles by gene conversion and recombination. Polymorphic sequence blocks are different from the blocks found in functional Pm3 alleles of hexaploid wheat, indicating an independent evolution of the Pm3 loci in the two species. A new functional gene was identified in a wild wheat accession from Syria. This gene, Pm3k, conferred intermediate race-specific resistance to powdery mildew, and consists of a mosaic of gene segments derived from non-functional alleles. This demonstrates that Pm3-based resistance is not very frequent in wild tetraploid wheat, and that the evolution of functional resistance genes occurred independently in wild tetraploid and bread wheat. The Pm3 sequence variability and geographic distribution indicated that diversity was higher in wild emmer wheat from the Levant area, compared with the accessions from Turkey. Further screens for Pm3 functional genes in wild wheat should therefore focus on accessions from the Levant region.

Ecological genomics of natural plant populations: the Israeli perspective

The genomic era revolutionized evolutionary population biology. The ecological genomics of the wild progenitors of wheat and barley reviewed here was central in the research program of the Institute of Evolution, University of Haifa, since 1975 ( http://evolution.haifa.ac.il ). We explored the following questions: (1) How much of the genomic and phenomic diversity of wild progenitors of cultivars (wild emmer wheat, Triticum dicoccoides, the progenitor of most wheat, plus wild relatives of the Aegilops species; wild barley, Hordeum spontaneum, the progenitor of cultivated barley; wild oat, Avena sterilis, the progenitor of cultivated oats; and wild lettuce species, Lactuca, the progenitor and relatives of cultivated lettuce) are adaptive and processed by natural selection at both coding and noncoding genomic regions? (2) What is the origin and evolution of genomic adaptation and speciation processes and their regulation by mutation, recombination, and transposons under spatiotemporal variables and stressful macrogeographic and microgeographic environments? (3) How much genetic resources are harbored in the wild progenitors for crop improvement? We advanced ecological genetics into ecological genomics and analyzed (regionally across Israel and the entire Near East Fertile Crescent and locally at microsites, focusing on the "Evolution Canyon" model) hundreds of populations and thousands of genotypes for protein (allozyme) and deoxyribonucleic acid (DNA) (coding and noncoding) diversity, partly combined with phenotypic diversity. The environmental stresses analyzed included abiotic (climatic and microclimatic, edaphic) and biotic (pathogens, demographic) stresses. Recently, we introduced genetic maps, cloning, and transformation of candidate genes. Our results indicate abundant genotypic and phenotypic diversity in natural plant populations. The organization and evolution of molecular and organismal diversity in plant populations, at all genomic regions and geographical scales, are nonrandom and are positively correlated with, and partly predictable by, abiotic and biotic environmental heterogeneity and stress. Biodiversity evolution, even in small isolated populations, is primarily driven by natural selection including diversifying, balancing, cyclical, and purifying selection regimes interacting with, but, ultimately, overriding the effects of mutation, migration, and stochasticity. The progenitors of cultivated plants harbor rich genetic resources and are the best hope for crop improvement by both classical and modern biotechnological methods. Future studies should focus on the interplay between structural and functional genome organization focusing on gene regulation.

Frequent sequence exchanges between homologs of RPP8 in Arabidopsis are not necessarily associated with genomic proximity

Disease resistance (R) genes are often clustered in plant genomes and may exhibit heterogeneous rates of evolution. Some (type I R genes) have evolved rapidly through frequent sequence exchanges, while others (type II R genes) have evolved independently and tend to be conserved in different genotypes or related species. The RPP8 resistance gene in Arabidopsis thaliana is located at a complex locus that also harbors the sequence-related resistance genes HRT and RCY1 in different ecotypes. We sequenced 98 homologs of RPP8 from A. thaliana, Arabidopsis arenosa and Arabidopsis lyrata. Three lineages of type II and one lineage of type I RPP8 homologs were identified. Two of the three lineages of type II genes are each represented by a single-copy locus on either chromosomes I or V. Chromosome V contains two small clusters of RPP8 paralogs. One cluster contains both type I and type II genes and the other comprises only type I genes. These multi-copy loci have expanded and contracted through unequal crossovers, which have generated chimeric genes as well as variations in copy number. Sequence exchanges, most likely gene conversions, were detected between RPP8 homologs that are spatially separated by 2.2 Mb and 12 cM. The sequence exchanges between type I homologs within a locus have been more frequent than sequence exchanges between homologs from two different loci, indicating the influence of chromosomal position on the evolution of these R genes. However, physical distance was not the only factor determining the frequency of sequence exchange, because some closely linked paralogs exhibited little sequence exchange. At least two distinct lineages of type II RPP8 homologs were identified in different species, with obvious allelic/orthologous relationships within each lineage. Therefore, the differentiation of type I and type II RPP8 homologs seems to have occurred before speciation of A. thaliana, A. arenosa and A. lyrata.

Evolution and genetic population structure of prickly lettuce (Lactuca serriola) and its RGC2 resistance gene cluster

Genetic structure and diversity of natural populations of prickly lettuce (Lactuca serriola) were studied using AFLP markers and then compared with the diversity of the RGC2 disease resistance gene cluster. Screening of 696 accessions from 41 populations using 319 AFLP markers showed that eastern Turkish and Armenian populations were the most diverse populations and might be located in the origin and center of diversity of L. serriola. Screening 709 accessions using the microsatellite MSATE6 that is located in the coding region of most RGC2 homologs detected 366 different haplotypes. Again, the eastern Turkish and Armenian populations had the highest diversities at the RGC2 cluster. The diversities at the RGC2 cluster in different populations were significantly correlated with their genomewide diversities. There was significant variation of copy number of RGC2 homologs in different populations, ranging from 12 to 22 copies per genome. The nucleotide diversities of two conserved lineages (type II) of RGC2 genes (K and L) were not correlated with diversities calculated using the MSATE6 or AFLP data. We hypothesize that the high genomewide diversity and diversity of the RGC2 cluster in eastern Turkish and Armenian populations resulted from high abiotic and biotic stresses in the regions of origin of L. serriola.

Silencing of the major family of NBS-LRR-encoding genes in lettuce results in the loss of multiple resistance specificities

The RGC2 gene cluster in lettuce (Lactuca sativa) is one of the largest known families of genes encoding nucleotide binding site-leucine-rich repeat (NBS-LRR) proteins. One of its members, RGC2B, encodes Dm3 which determines resistance to downy mildew caused by the oomycete Bremia lactucae carrying the cognate avirulence gene, Avr3. We developed an efficient strategy for analysis of this large family of low expressed genes using post-transcriptional gene silencing (PTGS). We transformed lettuce cv. Diana (carrying Dm3) using chimeric gene constructs designed to simultaneously silence RGC2B and the GUS reporter gene via the production of interfering hairpin RNA (ihpRNA). Transient assays of GUS expression in leaves accurately predicted silencing of both genes and were subsequently used to assay silencing in transgenic T(1) plants and their offspring. Levels of mRNA were reduced not only for RGC2B but also for all seven diverse RGC2 family members tested. We then used the same strategy to show that the resistance specificity encoded by the genetically defined Dm18 locus in lettuce cv. Mariska is the result of two resistance specificities, only one of which was silenced by ihpRNA derived from RGC2B. Analysis of progeny from crosses between transgenic, silenced tester stocks and lettuce accessions carrying other resistance genes previously mapped to the RGC2 locus indicated that two additional resistance specificities to B. lactucae, Dm14 and Dm16, as well as resistance to lettuce root aphid (Pemphigus bursarius L.), Ra, are encoded by RGC2 family members.

Introgressed and endogenous Mi-1 gene clusters in tomato differ by complex rearrangements in flanking sequences and show sequence exchange and diversifying selection among homologues

Many plant disease resistance genes (R-genes) encode proteins characterized by the presence of a nucleotide-binding site (NBS) and a leucine-rich repeat (LRR) region and occur in clusters of related genes in plant genomes. One such gene, Mi-1, confers isolate-specific resistance against root-knot nematodes, aphids and whiteflies in cultivated tomato, Solanum lycopersicon. The DNA region carrying Mi-1 and six closely related sequences was introgressed into tomato from Solanum peruvianum in the 1940s. For both susceptible and resistant tomato, Mi-1 homologues are present in two clusters with 3 and 4 copies each on the short arm of chromosome 6. Two homologues from each source are pseudogenes, and one homologue from each source encodes a truncated product. DNA sequence identity among the homologues including the truncated genes, but excluding the pseudogenes, ranges from 92.9 to 96.7%. All the non-pseudogene homologues are transcribed. Comparison of homologues suggests that extensive sequence exchange has occurred. Regions of diversifying selection are present in the ARC2 domain of the NBS region and dispersed throughout the LRR region, suggesting that these regions are possible locations of specificity determinants. Other sequences in the introgressed region have similarity to the Arabidopsis auxin-receptor protein TIR1, a jumonji-like transcription factor and a Na(+)/H(+) antiporter. Analysis of sequences flanking the Mi-1-homologues reveals blocks of homology, but complex differences in arrangement of these blocks when susceptible and resistant genotypes are compared indicating that the region has undergone considerable rearrangement during evolution, perhaps contributing to evolution of specificity.