Illegitimate recombination is a major evolutionary mechanism for initiating size variation in plant resistance genes

Rapid isolation of Yr9 via MutIsoSeq and QTL analysis of durable stripe rust resistance in wheat cultivar Xingzi 9104

The fungus Puccinia striiformis f. sp. tritici (Pst) is the causal agent of wheat stripe rust which constitutes a major limitation to wheat production. Cloning and applying disease-resistant genes are considered as an effective solution. Chinese wheat cultivar Xingzi 9104 (XZ9104) has exhibited durable resistance across multiple environments since its release. Through quantitative trait loci (QTL) analysis, eight QTL were found on chromosome arms 1BS, 1BL, 2AL, 2BL, 3BS, 4BL, 5BL and 7BL. YrXZ identified as 1RS.1BL translocation conferred race-specific all-stage resistance to Pst race CYR23. QYrxz.nwafu-1BL.6 and QYrxz.nwafu-3BS.7 were considered as the adult plant resistance genes Yr29 and Yr30, respectively. Notably, QYrxz.nwafu-2BL.5 accounted for 15.75-47.63% of the phenotypic variation across diverse environments and its pyramiding with Yr29 and Yr30 can confer high level of resistance. Other QTL were environment-dependent with minor effects. To clone the above resistance genes, we created a population of over 2,000 M5 mutants in XZ9104 using ethylmethane sulfonate (EMS) mutagenesis and screened various types of susceptible mutants. Using the MutIsoseq approach with five mutant lines susceptible to race CYR23, we rapid isolated a candidate gene for YrXZ encoding coiled-coil nucleotide-binding site leucine-rich repeat (CC-NBS-LRR) protein. Integrating cytological analysis, gene-based association analysis, transcriptomic profiling and virus-induced gene silencing (VIGS), we confirmed that the causal gene for YrXZ was indeed Yr9. This study demonstrated that multiple QTL with different effects contributed to the durable resistance in XZ9104. Understanding the molecular mechanisms and pathways involved in plant defense can inform future strategies for deploying resistance gene and engineering of genetic resistance against evolving diseases.

Structural variation in the pangenome of wild and domesticated barley

Pangenomes are collections of annotated genome sequences of multiple individuals of a species1. The structural variants uncovered by these datasets are a major asset to genetic analysis in crop plants2. Here we report a pangenome of barley comprising long-read sequence assemblies of 76 wild and domesticated genomes and short-read sequence data of 1,315 genotypes. An expanded catalogue of sequence variation in the crop includes structurally complex loci that are rich in gene copy number variation. To demonstrate the utility of the pangenome, we focus on four loci involved in disease resistance, plant architecture, nutrient release and trichome development. Novel allelic variation at a powdery mildew resistance locus and population-specific copy number gains in a regulator of vegetative branching were found. Expansion of a family of starch-cleaving enzymes in elite malting barleys was linked to shifts in enzymatic activity in micro-malting trials. Deletion of an enhancer motif is likely to change the developmental trajectory of the hairy appendages on barley grains. Our findings indicate that allelic diversity at structurally complex loci may have helped crop plants to adapt to new selective regimes in agricultural ecosystems.

Asymmetric Evolution of Protein Domains in the Leucine-Rich Repeat Receptor-Like Kinase Family of Plant Signaling Proteins

The coding sequences of developmental genes are expected to be deeply conserved, with cis-regulatory change driving the modulation of gene function. In contrast, proteins with roles in defense are expected to evolve rapidly, in molecular arms races with pathogens. However, some gene families include both developmental and defense genes. In these families, does the tempo and mode of evolution differ between genes with divergent functions, despite shared ancestry and structure? The leucine-rich repeat receptor-like kinase (LRR-RLKs) protein family includes members with roles in plant development and defense, thus providing an ideal system for answering this question. LRR-RLKs are receptors that traverse plasma membranes. LRR domains bind extracellular ligands; RLK domains initiate intracellular signaling cascades in response to ligand binding. In LRR-RLKs with roles in defense, LRR domains evolve faster than RLK domains. To determine whether this asymmetry extends to LRR-RLKs that function primarily in development, we assessed evolutionary rates and tested for selection acting on 11 subfamilies of LRR-RLKs, using deeply sampled protein trees. To assess functional evolution, we performed heterologous complementation assays in Arabidopsis thaliana (Arabidopsis). We found that the LRR domains of all tested LRR-RLK proteins evolved faster than their cognate RLK domains. All tested subfamilies of LRR-RLKs had strikingly similar patterns of molecular evolution, despite divergent functions. Heterologous transformation experiments revealed that multiple mechanisms likely contribute to the evolution of LRR-RLK function, including escape from adaptive conflict. Our results indicate specific and distinct evolutionary pressures acting on LRR versus RLK domains, despite diverse organismal roles for LRR-RLK proteins.

The N-terminally truncated helper NLR NRG1C antagonizes immunity mediated by its full-length neighbors NRG1A and NRG1B

Both plants and animals utilize nucleotide-binding leucine-rich repeat immune receptors (NLRs) to perceive the presence of pathogen-derived molecules and induce immune responses. NLR genes are far more abundant and diverse in vascular plants than in animals. Truncated NLRs, which lack one or more of the canonical domains, are also commonly encoded in plant genomes. However, little is known about their functions, especially the N-terminally truncated ones. Here, we show that the Arabidopsis thaliana N-terminally truncated helper NLR (hNLR) gene N REQUIREMENT GENE1 (NRG1C) is highly induced upon pathogen infection and in autoimmune mutants. The immune response and cell death conferred by some Toll/interleukin-1 receptor-type NLRs (TNLs) were compromised in Arabidopsis NRG1C overexpression lines. Detailed genetic analysis revealed that NRG1C antagonizes the immunity mediated by its full-length neighbors NRG1A and NRG1B. Biochemical tests suggested that NRG1C might interfere with the EDS1-SAG101 complex, which functions in immunity signaling together with NRG1A/1B. Interestingly, Brassicaceae NRG1Cs are functionally exchangeable and that the Nicotiana benthamiana N-terminally truncated hNLR NRG2 also antagonizes NRG1 activity. Together, our study uncovers an unexpected negative role of N-terminally truncated hNLRs in immunity in different plant species.

TdPm60 identified in wild emmer wheat is an ortholog of Pm60 and constitutes a strong candidate for PmG16 powdery mildew resistance

We identified TdPm60 alleles from wild emmer wheat (WEW), an ortholog of Pm60 from T. urartu, which constitutes a strong candidate for PmG16 mildew resistance. Deployment of PmG16 in Israeli modern bread wheat cultivar Ruta improved the resistance to several local Bgt isolates. Wild emmer wheat (WEW), the tetraploid progenitor of durum and bread wheat, is a valuable genetic resource for resistance to powdery mildew fungal disease caused by Blumeria graminis f. sp. tritici (Bgt). PmG16 gene, derived from WEW, confers high resistance to most tested Bgt isolates. We mapped PmG16 to a 1.4-cM interval between the flanking markers uhw386 and uhw390 on Chromosome 7AL. Based on gene annotation of WEW reference genome Zavitan_V1, 34 predicted genes were identified within the ~ 3.48-Mb target region. Six genes were annotated as associated with disease resistance, of which TRIDC7AG077150.1 was found to be highly similar to Pm60, previously cloned from Triticum urartu, and resides in the same syntenic region. The functional molecular marker (FMM) for Pm60 (M-Pm60-S1) co-segregated with PmG16, suggesting the Pm60 ortholog from WEW (designated here as TdPm60) as a strong candidate for PmG16. Sequence alignment identified only eight SNPs that differentiate between TdPm60 and TuPm60. Furthermore, TdPm60 was found to be present also in the WEW donor lines of the powdery mildew resistance genes MlIW172 and MlIW72, mapped to the same region of Chromosome 7AL as PmG16, suggesting that TdPm60 constitutes a candidate also for these genes. Furthermore, screening of additional 230 WEW accessions with Pm60 specific markers revealed 58 resistant accessions from the Southern Levant that harbored TdPm60, while none of the susceptible accessions showed the presence of this gene. Deployment of PmG16 in Israeli modern bread wheat cultivar Ruta conferred resistance against several local Bgt isolates.

Patterns and Processes of Diploidization in Land Plants

Most land plants are now known to be ancient polyploids that have rediploidized. Diploidization involves many changes in genome organization that ultimately restore bivalent chromosome pairing and disomic inheritance, and resolve dosage and other issues caused by genome duplication. In this review, we discuss the nature of polyploidy and its impact on chromosome pairing behavior. We also provide an overview of two major and largely independent processes of diploidization: cytological diploidization and genic diploidization/fractionation. Finally, we compare variation in gene fractionation across land plants and highlight the differences in diploidization between plants and animals. Altogether, we demonstrate recent advancements in our understanding of variation in the patterns and processes of diploidization in land plants and provide a road map for future research to unlock the mysteries of diploidization and eukaryotic genome evolution.

Plant NLR diversity: the known unknowns of pan-NLRomes

Plants and pathogens constantly adapt to each other. As a consequence, many members of the plant immune system, and especially the intracellular nucleotide-binding site leucine-rich repeat receptors, also known as NOD-like receptors (NLRs), are highly diversified, both among family members in the same genome, and between individuals in the same species. While this diversity has long been appreciated, its true extent has remained unknown. With pan-genome and pan-NLRome studies becoming more and more comprehensive, our knowledge of NLR sequence diversity is growing rapidly, and pan-NLRomes provide powerful platforms for assigning function to NLRs. These efforts are an important step toward the goal of comprehensively predicting from sequence alone whether an NLR provides disease resistance, and if so, to which pathogens.

The absence of the caffeine synthase gene is involved in the naturally decaffeinated status of Coffea humblotiana, a wild species from Comoro archipelago

Caffeine is the most consumed alkaloid stimulant in the world. It is synthesized through the activity of three known N-methyltransferase proteins. Here we are reporting on the 422-Mb chromosome-level assembly of the Coffea humblotiana genome, a wild and endangered, naturally caffeine-free, species from the Comoro archipelago. We predicted 32,874 genes and anchored 88.7% of the sequence onto the 11 chromosomes. Comparative analyses with the African Robusta coffee genome (C. canephora) revealed an extensive genome conservation, despite an estimated 11 million years of divergence and a broad diversity of genome sizes within the Coffea genus. In this genome, the absence of caffeine is likely due to the absence of the caffeine synthase gene which converts theobromine into caffeine through an illegitimate recombination mechanism. These findings pave the way for further characterization of caffeine-free species in the Coffea genus and will guide research towards naturally-decaffeinated coffee drinks for consumers.

A Truncated Singleton NLR Causes Hybrid Necrosis in Arabidopsis thaliana

Hybrid necrosis in plants arises from conflict between divergent alleles of immunity genes contributed by different parents, resulting in autoimmunity. We investigate a severe hybrid necrosis case in Arabidopsis thaliana, where the hybrid does not develop past the cotyledon stage and dies 3 weeks after sowing. Massive transcriptional changes take place in the hybrid, including the upregulation of most NLR (nucleotide-binding site leucine-rich repeat) disease-resistance genes. This is due to an incompatible interaction between the singleton TIR-NLR gene DANGEROUS MIX 10 (DM10), which was recently relocated from a larger NLR cluster, and an unlinked locus, DANGEROUS MIX 11 (DM11). There are multiple DM10 allelic variants in the global A. thaliana population, several of which have premature stop codons. One of these, which has a truncated LRR-PL (leucine-rich repeat [LRR]-post-LRR) region, corresponds to the DM10 risk allele. The DM10 locus and the adjacent genomic region in the risk allele carriers are highly differentiated from those in the nonrisk carriers in the global A. thaliana population, suggesting that this allele became geographically widespread only relatively recently. The DM11 risk allele is much rarer and found only in two accessions from southwestern Spain-a region from which the DM10 risk haplotype is absent-indicating that the ranges of DM10 and DM11 risk alleles may be nonoverlapping.

What natural variation can teach us about resistance durability

Breeding a crop variety to be resistant to a pathogen usually takes years. This is problematic because pathogens, with short generation times and fluid genomes, adapt quickly to overcome resistance. The triumph of the pathogen is not inevitable, however, as there are numerous examples of durable resistance, particularly in wild plants. Which factors then contribute to such resistance stability over millennia? We review current knowledge of wild and agricultural pathosystems, detailing the importance of genetic, species and spatial heterogeneity in the prevention of pathogen outbreaks. We also highlight challenges associated with increasing resistance diversity in crops, both in light of pathogen (co-)evolution and breeding practices. Historically it has been difficult to incorporate heterogeneity into agriculture due to reduced efficiency in harvesting. Recent advances implementing computer vision and automation in agricultural production may improve our ability to harvest mixed genotype and mixed species plantings, thereby increasing resistance durability.

Integrated Analysis of Large-Scale Omics Data Revealed Relationship Between Tissue Specificity and Evolutionary Dynamics of Small RNAs in Maize (Zea mays)

The evolutionary dynamics and tissue specificity of protein-coding genes are well documented in plants. However, the evolutionary consequences of small RNAs (sRNAs) on tissue-specific functions remain poorly understood. Here, we performed integrated analysis of 195 deeply sequenced sRNA libraries of maize B73, representing more than 10 tissues, and identified a comprehensive list of 419 maize microRNA (miRNA) genes, 271 of which were newly discovered in this study. We further characterized the evolutionary dynamics and tissue specificity of miRNA genes and corresponding miRNA isoforms (isomiRs). Our analysis revealed that tissue specificity of isomiR events tends to be associated with miRNA gene abundance and suggested that the frequencies of isomiR types are affected by the local genomic regions. Moreover, genome duplication (GD) events have dramatic effect on evolutionary dynamics of maize miRNA genes, and the abundance divergence for tissue-specific miRNA genes is associated with GD events. Further study indicated that duplicate miRNA genes with tissue-specific expression patterns, such as miR2275a, a phased siRNA (phasiRNA) trigger, contribute to phenotypic traits in maize. Additionally, our study revealed the expression preference of 21- and 24-nt phasiRNAs in relation to tissue specificity. This large-scale sRNAomic study depicted evolutionary implications of tissue-specific maize sRNAs, which coordinate genome duplication, isomiR modification, phenotypic traits and phasiRNAs differentiation.

Transcriptome-Based Identification and Molecular Evolution of the Cytochrome P450 Genes and Expression Profiling under Dimethoate Treatment in Amur Stickleback (Pungitius sinensis)

Cytochrome P450s (CYPs) are a family of membrane-bound mono-oxygenase proteins, which are involved in cell metabolism and detoxification of various xenobiotic substances. In this study, we identified 58 putative CYP genes in Amur stickleback (Pungitius sinensis) based on the transcriptome sequencing. Conserved motif distribution suggested their functional relevance within each group. Some present recombination events have accelerated the evolution of this gene family. Moreover, a few positive selection sites were identified, which may have accelerated the functional divergence of this family of proteins. Expression patterns of these CYP genes were investigated and indicated that most were affected by dimethoate treatment, suggesting that CYPs were involved in the detoxication of dimethoate. This study will provide a foundation for the further functional investigation of CYP genes in fishes.

RPW8/HR repeats control NLR activation in Arabidopsis thaliana

In many plant species, conflicts between divergent elements of the immune system, especially nucleotide-binding oligomerization domain-like receptors (NLR), can lead to hybrid necrosis. Here, we report deleterious allele-specific interactions between an NLR and a non-NLR gene cluster, resulting in not one, but multiple hybrid necrosis cases in Arabidopsis thaliana. The NLR cluster is RESISTANCE TO PERONOSPORA PARASITICA 7 (RPP7), which can confer strain-specific resistance to oomycetes. The non-NLR cluster is RESISTANCE TO POWDERY MILDEW 8 (RPW8) / HOMOLOG OF RPW8 (HR), which can confer broad-spectrum resistance to both fungi and oomycetes. RPW8/HR proteins contain at the N-terminus a potential transmembrane domain, followed by a specific coiled-coil (CC) domain that is similar to a domain found in pore-forming toxins MLKL and HET-S from mammals and fungi. C-terminal to the CC domain is a variable number of 21- or 14-amino acid repeats, reminiscent of regulatory 21-amino acid repeats in fungal HET-S. The number of repeats in different RPW8/HR proteins along with the sequence of a short C-terminal tail predicts their ability to activate immunity in combination with specific RPP7 partners. Whether a larger or smaller number of repeats is more dangerous depends on the specific RPW8/HR autoimmune risk variant.

Comparative Analysis of the aquaporin Gene Family in 12 Fish Species

Simple Summary

Aquaporins (Aqps) are a group of membrane proteins. In this study, 166 Aqp genes were identified in 12 fish species. Gene organization, motif distribution, recombination, and selection pressure were performed to investigate their evolutionary characteristics. In addition, expression profiles of Aqps were also examined under pathogens infection and organophosphorus pesticide stress. This study will provide a useful reference for further functional study.

Abstract

Aquaporins (Aqps) are a class of water channel proteins that play key roles in many physiological functions and cellular processes. Here, we analyzed 166 putative Aqp genes in 12 fish species and divided them into four groups. Gene organization and motif distribution analyses suggested potentially conserved functions in each group. Several recombination events were identified in some members, which accelerate their divergence in evolution. Furthermore, a few positive selection sites were identified, and mutations at these sites could alter the stability of Aqp proteins. In addition, expression profiles of some Aqp genes under pathogen infection and organophosphorus pesticide stress were also investigated. The result implied that several Aqp genes may affect different immune responses and osmoregulation. This study provides a comparative analysis of the fish Aqp gene family to facilitate further functional analyses.

Genetic Mapping and Molecular Characterization of a Broad-spectrum Phytophthora sojae Resistance Gene in Chinese Soybean

Phytophthora root rot (PRR) causes serious annual soybean yield losses worldwide. The most effective method to prevent PRR involves growing cultivars that possess genes conferring resistance to Phytophthora sojae (Rps). In this study, QTL-sequencing combined with genetic mapping was used to identify RpsX in soybean cultivar Xiu94-11 resistance to all P. sojae isolates tested, exhibiting broad-spectrum PRR resistance. Subsequent analysis revealed RpsX was located in the 242-kb genomic region spanning the RpsQ locus. However, a phylogenetic investigation indicated Xiu94-11 carrying RpsX is distantly related to the cultivars containing RpsQ, implying RpsX and RpsQ have different origins. An examination of candidate genes revealed RpsX and RpsQ share common nonsynonymous SNP and a 144-bp insertion in the Glyma.03g027200 sequence encoding a leucine-rich repeat (LRR) region. Glyma.03g027200 was considered to be the likely candidate gene of RpsQ and RpsX. Sequence analyses confirmed that the 144-bp insertion caused by an unequal exchange resulted in two additional LRR-encoding fragments in the candidate gene. A marker developed based on the 144-bp insertion was used to analyze the genetic population and germplasm, and proved to be useful for identifying the RpsX and RpsQ alleles. This study implies that the number of LRR units in the LRR domain may be important for PRR resistance in soybean.

Evolutionary divergence of the rye Pm17 and Pm8 resistance genes reveals ancient diversity

We have isolated a novel powdery mildew resistance gene in wheat that was originally introgressed from rye. Further analysis revealed evolutionary divergent history of wheat and rye orthologous resistance genes. Wheat production is under constant threat from a number of fungal pathogens, among them is wheat powdery mildew (Blumeria graminis f. sp. tritici). Deployment of resistance genes is the most economical and sustainable method for mildew control. However, domestication and selective breeding have narrowed genetic diversity of modern wheat germplasm, and breeders have relied on wheat relatives for enriching its gene pool through introgression. Translocations where the 1RS chromosome arm was introgressed from rye to wheat have improved yield and resistance against various pathogens. Here, we isolated the Pm17 mildew resistance gene located on the 1RS introgression in wheat cultivar 'Amigo' and found that it is an allele or a close paralog of the Pm8 gene isolated earlier from 'Petkus' rye. Functional validation using transient and stable transformation confirmed the identity of Pm17. Analysis of Pm17 and Pm8 coding regions revealed an overall identity of 82.9% at the protein level, with the LRR domains being most divergent. Our analysis also showed that the two rye genes are much more diverse compared to the variants encoded by the Pm3 gene in wheat, which is orthologous to Pm17/Pm8 as concluded from highly conserved upstream sequences in all these genes. Thus, the evolutionary history of these orthologous loci differs in the cereal species rye and wheat and demonstrates that orthologous resistance genes can take different routes towards functionally active genes. These findings suggest that the isolation of Pm3/Pm8/Pm17 orthologs from other grass species, additional alleles from the rye germplasm as well as possibly synthetic variants will result in novel resistance genes useful in wheat breeding.

Regulation and Evolution of NLR Genes: A Close Interconnection for Plant Immunity

NLR (NOD-like receptor) genes belong to one of the largest gene families in plants. Their role in plants' resistance to pathogens has been clearly described for many members of this gene family, and dysregulation or overexpression of some of these genes has been shown to induce an autoimmunity state that strongly affects plant growth and yield. For this reason, these genes have to be tightly regulated in their expression and activity, and several regulatory mechanisms are described here that tune their gene expression and protein levels. This gene family is subjected to rapid evolution, and to maintain diversity at NLRs, a plethora of genetic mechanisms have been identified as sources of variation. Interestingly, regulation of gene expression and evolution of this gene family are two strictly interconnected aspects. Indeed, some examples have been reported in which mechanisms of gene expression regulation have roles in promotion of the evolution of this gene family. Moreover, co-evolution of the NLR gene family and other gene families devoted to their control has been recently demonstrated, as in the case of miRNAs.

NLR diversity, helpers and integrated domains: making sense of the NLR IDentity

Plant innate immunity relies on genetically predetermined repertoires of immune receptors to detect pathogens and trigger an effective immune response. A large proportion of these receptors are from the Nucletoide Binding Leucine Rich Repeat (NLR) gene family. As plants live longer than most pathogens, maintaining diversity of NLRs and deploying efficient 'pathogen traps' is necessary to withstand the evolutionary battle. In this review, we summarize the sources of diversity in NLR plant immune receptors giving an overview of genomic, regulatory as well as functional studies, including the latest concepts of NLR helpers and NLRs with integrated domains.

Distribution and haplotype diversity of WKS resistance genes in wild emmer wheat natural populations

The wheat stripe rust resistance gene Yr36 ( WKS1 ) with a unique kinase-START domain architecture is highly conserved in wild emmer wheat natural populations. Wild emmer wheat (Triticum dicoccoides) populations have developed various resistance strategies against the stripe rust pathogen Puccinia striiformis f. sp. tritici (Pst). The wild emmer gene, Yr36 (WKS1), which confers partial resistance to a broad spectrum of Pst races, is composed of a kinase and a START lipid-binding domain, a unique gene architecture found only in the Triticeae tribe. The analysis of 435 wild emmer accessions from a broad range of natural habitats revealed that WKS1 and its paralogue WKS2 are present only in the southern distribution range of wild emmer in the Fertile Crescent, supporting the idea that wheat domestication occurred in the northern populations. An analysis of full-length WKS1 sequence from 54 accessions identified 15 different haplotypes and very low-nucleotide diversity (π = 0.00019). The high level of WKS1 sequence conservation among wild emmer populations is in contrast to the high level of diversity previously observed in NB-LRR genes (e.g., Lr10 and Pm3). This phenomenon may reflect the different resistance mechanisms and different evolutionary pathways that shaped these genes, and may shed light on the evolution of genes that confer partial resistance to stripe rust. Only five WKS1 coding sequence haplotypes were revealed among all tested accessions, encoding four different putative WKS1 proteins (designated P0, P1, P2, and P3). Infection tests showed that P0, P1, and P3 haplotypes display a resistance response, while P2 displayed a susceptible response. These results show that the WKS1 proteins (P0, P1, and P3) can be useful to improve wheat resistance to stripe rust.

Arabidopsis thaliana DM2h (R8) within the Landsberg RPP1-like Resistance Locus Underlies Three Different Cases of EDS1-Conditioned Autoimmunity

Plants have a large panel of nucleotide-binding/leucine rich repeat (NLR) immune receptors which monitor host interference by diverse pathogen molecules (effectors) and trigger disease resistance pathways. NLR receptor systems are necessarily under tight control to mitigate the trade-off between induced defenses and growth. Hence, mis-regulated NLRs often cause autoimmunity associated with stunting and, in severe cases, necrosis. Nucleocytoplasmic ENHANCED DISEASE SUSCEPTIBILITY1 (EDS1) is indispensable for effector-triggered and autoimmune responses governed by a family of Toll-Interleukin1-Receptor-related NLR receptors (TNLs). EDS1 operates coincidently or immediately downstream of TNL activation to transcriptionally reprogram cells for defense. We show here that low levels of nuclear-enforced EDS1 are sufficient for pathogen resistance in Arabidopsis thaliana, without causing negative effects. Plants expressing higher nuclear EDS1 amounts have the genetic, phenotypic and transcriptional hallmarks of TNL autoimmunity. In a screen for genetic suppressors of nuclear EDS1 autoimmunity, we map multiple, independent mutations to one gene, DM2h, lying within the polymorphic DANGEROUS MIX2 cluster of TNL RPP1-like genes from A. thaliana accession Landsberg erecta (Ler). The DM2 locus is a known hotspot for deleterious epistatic interactions leading to immune-related incompatibilities between A. thaliana natural accessions. We find that DM2h^Ler underlies two further genetic incompatibilities involving the RPP1-like^Ler locus and EDS1. We conclude that the DM2h^Ler TNL protein and nuclear EDS1 cooperate, directly or indirectly, to drive cells into an immune response at the expense of growth. A further conclusion is that regulating the available EDS1 nuclear pool is fundamental for maintaining homeostatic control of TNL immune pathways.

Plants tune their cellular and developmental programs to different environmental stimuli. Central players in the plant biotic stress response network are intracellular NLR receptors which intercept specific disease-inducing molecules (effectors) produced by pathogenic microbes. Variation in NLR gene repertoires between plant genetic lines is driven by pathogen selection pressure. One evolutionary question is how new, functional NLRs are assembled within a plant genome without mis-activating defense pathways, which can have strong negative effects on growth and fitness. This study focuses on a large, polymorphic sub-class of NLR receptors called TNLs present in dicotyledenous plant lineages. TNL receptors confer immunity to a broad range of pathogens. They also frequently underlie autoimmunity caused by their mis-regulation or deleterious allelic interactions with other genes in crosses between different genetic lines (hybrid incompatibility, HI). TNL pathogen-triggered and autoimmune responses require the conserved nucleocytoplasmic protein EDS1 to transcriptionally reprogram cells for defense. We discover in Arabidopsis thaliana that high levels of nuclear-enriched EDS1 induce transcriptional activation of defenses and growth inhibition without a pathogen effector stimulus. In a mutational screen, we identify one rapidly evolving TNL gene, DM2h^Ler, as a driver of nuclear EDS1 autoimmunity. DM2h^Ler also contributes to two separate cases of EDS1-dependent autoimmunity. Genetic cooperation between DM2h^Ler and EDS1 suggests a functional relationship in the transcriptional feed-forward regulation of defense pathways.

Divergence of the bZIP Gene Family in Strawberry, Peach, and Apple Suggests Multiple Modes of Gene Evolution after Duplication

The basic leucine zipper (bZIP) transcription factors are the most diverse members of dimerizing transcription factors. In the present study, 50, 116, and 47 bZIP genes were identified in Malus domestica (apple), Prunus persica (peach), and Fragaria vesca (strawberry), respectively. Species-specific duplication was the main contributor to the large number of bZIPs observed in apple. After WGD in apple genome, orthologous bZIP genes corresponding to strawberry on duplicated regions in apple genome were retained. However, in peach ancestor, these syntenic regions were quickly lost or deleted. Maybe the positive selection contributed to the expansion of clade S to adapt to the development and environment stresses. In addition, purifying selection was mainly responsible for bZIP sequence-specific DNA binding. The analysis of orthologous pairs between chromosomes indicates that these orthologs derived from one gene duplication located on one of the nine ancient chromosomes in the Rosaceae. The comparative analysis of bZIP genes in three species provides information on the evolutionary fate of bZIP genes in apple and peach after they diverged from strawberry.

Analysis of a plant complex resistance gene locus underlying immune-related hybrid incompatibility and its occurrence in nature

Mechanisms underlying speciation in plants include detrimental (incompatible) genetic interactions between parental alleles that incur a fitness cost in hybrids. We reported on recessive hybrid incompatibility between an Arabidopsis thaliana strain from Poland, Landsberg erecta (Ler), and many Central Asian A. thaliana strains. The incompatible interaction is determined by a polymorphic cluster of Toll/interleukin-1 receptor-nucleotide binding-leucine rich repeat (TNL) RPP1 (Recognition of Peronospora parasitica1)-like genes in Ler and alleles of the receptor-like kinase Strubbelig Receptor Family 3 (SRF3) in Central Asian strains Kas-2 or Kond, causing temperature-dependent autoimmunity and loss of growth and reproductive fitness. Here, we genetically dissected the RPP1-like Ler locus to determine contributions of individual RPP1-like Ler (R1–R8) genes to the incompatibility. In a neutral background, expression of most RPP1-like Ler genes, except R3, has no effect on growth or pathogen resistance. Incompatibility involves increased R3 expression and engineered R3 overexpression in a neutral background induces dwarfism and sterility. However, no individual RPP1-like Ler gene is sufficient for incompatibility between Ler and Kas-2 or Kond, suggesting that co-action of at least two RPP1-like members underlies this epistatic interaction. We find that the RPP1-like Ler haplotype is frequent and occurs with other Ler RPP1-like alleles in a local population in Gorzów Wielkopolski (Poland). Only Gorzów individuals carrying the RPP1-like Ler haplotype are incompatible with Kas-2 and Kond, whereas other RPP1-like alleles in the population are compatible. Therefore, the RPP1-like Ler haplotype has been maintained in genetically different individuals at a single site, allowing exploration of forces shaping the evolution of RPP1-like genes at local and regional population scales.

In plants, naturally evolving disease resistance (R) genes can cause autoimmunity when combined with different genetic backgrounds. This phenomenon, called immune-related hybrid incompatibility (HI), leads to growth inhibition and fitness loss due to inappropriate activation of defense. HI likely reflects different evolutionary paths of immune-related genes in nature. We have examined the genetic architecture of a complex R locus present in a Central European accession (Ler) which underlies HI with Central Asian accessions of Arabidopsis. We show that expression of one gene (R3) within the Ler cluster of eight tandem R genes (R1–R8) controls the balance between growth and defense but that R3 needs at least one other co-acting member within the R locus to condition HI. We traced the R1–R8 haplotype to a local population of Ler relatives in Poland where it also underlies HI with Central Asian accessions. Occurrence of the incompatible haplotype in ∼30% of genetically diverse local individuals, suggests that it has not arisen recently and has been maintained through selection or drift. Co-occurrence in the same population of individuals containing different R genes that do not cause HI provides a basis for determining genetic and environmental forces influencing how plant immunity genes evolve and diversify.

The detection of a de novo allele of the Glu-1Dx gene in wheat-rye hybrid offspring

This study provides a link between a de novo gene and novel phenotype in wheat-rye hybrids that can be used as a model for induced de novo genetic variation. Wide hybridization can produce de novo DNA variation that may cause novel phenotypes. However, there is still a lack of specific links between changed genes and novel phenotypes in wide hybrids. The well-studied high-molecular-weight glutenin subunit (HMW-GS) genes in tribe Triticeae provide a useful model for addressing this issue. In this study, we investigated the feasibility of a wheat-rye hybridization method for inducing de novo phenotypes using the Glu-1Dx2.2 subunit as an example. We developed three hexaploid wheat lines with normal fertility and a Glu-1Dx2.2 variant, named Glu-1Dx2.2 (v) , derived from three F1 hybrids. The wild-type Glu-1Dx2.2 has two direct repeats of 295 bp length separated by an intervening 101 bp in its central repetitive region. In the mutant Glu-1Dx2.2 (v) , one copy of the repeats and the intervening sequence were deleted, probably through homology-dependent illegitimate recombination (IR). This study provides a direct link between a de novo allele and novel phenotype. Our results indicate that the wheat-rye method may be a useful tool to induce de novo genetic variations that broaden the genetic diversity for wheat improvement.

Genomic variability as a driver of plant-pathogen coevolution?

Pathogens apply one of the strongest selective pressures in plant populations. Understanding plant-pathogen coevolution has therefore been a major research focus for at least sixty years [1]. Recent comparative genomic studies have revealed that the genes involved in plant defense and pathogen virulence are among the most polymorphic in the respective genomes. Which fraction of this diversity influences the host-pathogen interaction? Do coevolutionary dynamics maintain variation? Here we review recent literature on the evolutionary and molecular processes that shape this variation, focusing primarily on gene-for-gene interactions. In summarizing theoretical and empirical studies of the processes that shape this variation in natural plant and pathogen populations, we find a disconnect between the complexity of ecological interactions involving hosts and their myriad microbes, and the models that describe them.

Substitutions of two amino acids in the nucleotide-binding site domain of a resistance protein enhance the hypersensitive response and enlarge the PM3F resistance spectrum in wheat

Proteins with nucleotide-binding site (NBS) and leucine-rich repeat (LRR) domains are major components of the plant immune system. They usually mediate resistance against a subgroup of races of a specific pathogen. For the allelic series of the wheat powdery mildew resistance gene Pm3, alleles with a broad and a narrow resistance spectrum have been described. Here, we show that a broad Pm3 spectrum range correlates with a fast and intense hypersensitive response (HR) in a Nicotiana transient-expression system and this activity can be attributed to two particular amino acids in the ARC2 subdomain of the NBS. The combined substitution of these amino acids in narrow-spectrum PM3 proteins enhances their capacity to induce an HR in Nicotiana benthamiana, and we demonstrate that these substitutions also enlarge the resistance spectrum of the Pm3f allele in wheat. Finally, using Bph14, we show that the region carrying the relevant amino acids also plays a role in the HR regulation of another coiled-coil NBS-LRR resistance protein. These results highlight the importance of an optimized NBS-'molecular switch' for the conversion of initial pathogen perception by the LRR into resistance-protein activation, and we describe a possible approach to extend the effectiveness of resistance genes via minimal targeted modifications in the NBS domain.

The physical map of wheat chromosome 1BS provides insights into its gene space organization and evolution

Background

The wheat genome sequence is an essential tool for advanced genomic research and improvements. The generation of a high-quality wheat genome sequence is challenging due to its complex 17 Gb polyploid genome. To overcome these difficulties, sequencing through the construction of BAC-based physical maps of individual chromosomes is employed by the wheat genomics community. Here, we present the construction of the first comprehensive physical map of chromosome 1BS, and illustrate its unique gene space organization and evolution.

Results

Fingerprinted BAC clones were assembled into 57 long scaffolds, anchored and ordered with 2,438 markers, covering 83% of chromosome 1BS. The BAC-based chromosome 1BS physical map and gene order of the orthologous regions of model grass species were consistent, providing strong support for the reliability of the chromosome 1BS assembly. The gene space for chromosome 1BS spans the entire length of the chromosome arm, with 76% of the genes organized in small gene islands, accompanied by a two-fold increase in gene density from the centromere to the telomere.

Conclusions

This study provides new evidence on common and chromosome-specific features in the organization and evolution of the wheat genome, including a non-uniform distribution of gene density along the centromere-telomere axis, abundance of non-syntenic genes, the degree of colinearity with other grass genomes and a non-uniform size expansion along the centromere-telomere axis compared with other model cereal genomes. The high-quality physical map constructed in this study provides a solid basis for the assembly of a reference sequence of chromosome 1BS and for breeding applications.

Rye Pm8 and wheat Pm3 are orthologous genes and show evolutionary conservation of resistance function against powdery mildew

The improvement of wheat through breeding has relied strongly on the use of genetic material from related wild and domesticated grass species. The 1RS chromosome arm from rye was introgressed into wheat and crossed into many wheat lines, as it improves yield and fungal disease resistance. Pm8 is a powdery mildew resistance gene on 1RS which, after widespread agricultural cultivation, is now widely overcome by adapted mildew races. Here we show by homology-based cloning and subsequent physical and genetic mapping that Pm8 is the rye orthologue of the Pm3 allelic series of mildew resistance genes in wheat. The cloned gene was functionally validated as Pm8 by transient, single-cell expression analysis and stable transformation. Sequence analysis revealed a complex mosaic of ancient haplotypes among Pm3- and Pm8-like genes from different members of the Triticeae. These results show that the two genes have evolved independently after the divergence of the species 7.5 million years ago and kept their function in mildew resistance. During this long time span the co-evolving pathogens have not overcome these genes, which is in strong contrast to the breakdown of Pm8 resistance since its introduction into commercial wheat 70 years ago. Sequence comparison revealed that evolutionary pressure acted on the same subdomains and sequence features of the two orthologous genes. This suggests that they recognize directly or indirectly the same pathogen effectors that have been conserved in the powdery mildews of wheat and rye.

Haplotype variability and identification of new functional alleles at the Rdg2a leaf stripe resistance gene locus

The barley Rdg2a locus confers resistance to the leaf stripe pathogen Pyrenophora graminea and, in the barley genotype Thibaut, it is composed of a gene family with three highly similar paralogs. Only one member of the gene family (called as Rdg2a) encoding for a CC-NB-LRR protein is able to confer resistance to the leaf stripe isolate Dg2. To study the genome evolution and diversity at the Rdg2a locus, sequences spanning the Rdg2a gene were compared in two barley cultivars, Thibaut and Morex, respectively, resistant and susceptible to leaf stripe. An overall high level of sequence conservation interrupted by several rearrangements that included three main deletions was observed in the Morex contig. The main deletion of 13,692 bp was most likely derived from unequal crossing over between Rdg2a paralogs leading to the generation of a chimeric Morex rdg2a gene which was not associated to detectable level of resistance toward leaf stripe. PCR-based analyses of genic and intergenic regions at the Rdg2a locus in 29 H. vulgare lines and one H. vulgare ssp. spontaneum accession indicated large haplotype variability in the cultivated barley gene pool suggesting rapid and recent divergence at this locus. Barley genotypes showing the same haplotype as Thibaut at the Rdg2a locus were selected for a Rdg2a allele mining through allele re-sequencing and two lines with polymorphic nucleotides leading to amino acid changes in the CC-NB and LRR encoding domains, respectively, were identified. Analysis of nucleotide diversity of the Rdg2a alleles revealed that the polymorphic sites were subjected to positive selection. Moreover, strong positively selected sites were located in the LRR encoding domain suggesting that both positive selection and divergence at homologous loci are possibly representing the molecular mechanism for the generation of high diversity at the Rdg2a locus in the barley gene pool.

Divergent evolutionary mechanisms of co-located Tak/Lrk and Glu-D3 loci revealed by comparative analysis of grass genomes

Seed storage and disease resistance proteins are major traits of wheat. The study of their gene organization and evolution has some implications in breeding. In this study, we characterized the hexaploid wheat D-genome BAC clone TaBAC703A9 that contains a low molecular weight glutenin locus (Glu-D3) and a resistance gene analogue cluster. With a gene density of one gene per 4.8 kb, the cluster contains four resistance gene analogues, namely Tak703-1, Lrr703, Tak703, and Lrk703. This structural cluster unit was conserved across nine grass genomes, but divergent evolutionary mechanisms have been involved in shaping the Tak/Lrk loci in the different species. Gene duplication was the major force for the Tak/Lrk evolution in oats, maize, barley, wheat, sorghum, and Brachypodium, while tandem duplication drove the expansion of this locus in japonica rice. Despite the close proximity of the Glu-D3 and the Tak/Lrk loci in wheat, the evolutionary mechanisms that drove their amplification differ. The Glu-D3 region had a lower gene density, and its amplification was driven by retroelements.

Mechanism of random integration of foreign DNA in transgenic mice

Little is known about how foreign DNA is randomly integrated into chromosomes in transgenic animals. In the current study, the insertion sites of 36 transgenic mice were mapped by thermal asymmetric interlaced PCR, and 38 junction sequences were obtained from 30 samples. Analysis of the 38 sequences revealed that 44.7 % of integration events occurred within host gene regions, including 13.2 % (5/38) in exonic regions and 31.6 % (12/38) in intronic regions. The results also revealed that all non-end side integrations of foreign DNA were mediated by short sequence homologies (microhomologies) and that the end side integrations occurred in the presence or absence of microhomologies. In addition, microhomology-mediated mechanisms were also confirmed in four transgenic Arabidopsis thaliana lines. The results indicate that foreign DNA is easily integrated into host gene regions. These results also suggest that the integration of both ends of foreign DNA follows the above-mentioned mechanism in many transgenic/transformed organisms.

Comparative analysis of genome composition in Triticeae reveals strong variation in transposable element dynamics and nucleotide diversity

A 454 sequencing snapshot was utilised to investigate the genome composition and nucleotide diversity of transposable elements (TEs) for several Triticeae taxa, including Triticum aestivum, Hordeum vulgare, Hordeum spontaneum and Secale cereale together with relatives of the A, B and D genome donors of wheat, Triticum urartu (A), Aegilops speltoides (S) and Aegilops tauschii (D). Additional taxa containing the A genome, Triticum monococcum and its wild relative Triticum boeoticum, were also included. The main focus of the analysis was on the genomic composition of TEs as these make up at least 80% of the overall genome content. Although more than 200 TE families were identified in each species, approximately 50% of the overall genome comprised 12-15 TE families. The BARE1 element was the largest contributor to all genomes, contributing more than 10% to the overall genome. We also found that several TE families differ strongly in their abundance between species, indicating that TE families can thrive extremely successfully in one species while going virtually extinct in another. Additionally, the nucleotide diversity of BARE1 populations within individual genomes was measured. Interestingly, the nucleotide diversity in the domesticated barley H. vulgare cv. Barke was found to be twice as high as in its wild progenitor H. spontaneum, suggesting that the domesticated barley gained nucleotide diversity from the addition of different genotypes during the domestication and breeding process. In the rye/wheat lineage, sequence diversity of BARE1 elements was generally higher, suggesting that factors such as geographical distribution and mating systems might play a role in intragenomic TE diversity.

Control of meiotic recombination frequency in plant genomes

Sexual eukaryotes reproduce via the meiotic cell division, where ploidy is halved and homologous chromosomes undergo reciprocal genetic exchange, termed crossover (CO). CO frequency has a profound effect on patterns of genetic variation and species evolution. Relative CO rates vary extensively both within and between plant genomes. Plant genome size varies by over 1000-fold, largely due to differential expansion of repetitive sequences, and increased genome size is associated with reduced CO frequency. Gene versus repeat sequences associate with distinct chromatin modifications, and evidence from plant genomes indicates that this epigenetic information influences CO patterns. This is consistent with data from diverse eukaryotes that demonstrate the importance of chromatin structure for control of meiotic recombination. In this review I will discuss CO frequency patterns in plant genomes and recent advances in understanding recombination distributions.

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.

Genome change in wheat observed through the structure and expression of α/β-gliadin genes

To better understand genome structure and the expression of α/β-gliadin multigenes in hexaploid wheat, bacterial artificial chromosome (BAC) clones containing α/β-gliadin genes from the three loci, Gli-A2, Gli-B2, and Gli-D2, were screened. Based on their restriction fragment patterns, we selected five BAC clones, namely, two clones for Gli-A2, two clones for Gli-B2, and one clone for Gli-D2, to fully sequence. Approximately 200 kb was sequenced for each locus. In total, twelve α/β-gliadin intact genes and four pseudogenes were found, and retrotransposons or other transposons existed in each BAC clone. Dot-plot analysis revealed the pattern of genome segmental duplication within each BAC. We calculated time since duplication of each set of α/β-gliadin genes and insertion of retrotransposons. Duplication of all adjacent genes within the same BAC clone took place before or after allotetrapolyploidization, but duplication of certain genes occurred before diploid differentiation of wheat species. Retrotransposons were also inserted before and after the segmental duplication events. Furthermore, translocation of α/β-gliadin genes from chromosomes 1 to 6 apparently occurred before the diversification of various wheat genomes. Duplication of genome segments containing α/β-gliadin genes and retrotransposons were brought about through unequal crossing-over or saltatory replication and α/β-gliadin genes per se were duplicated without any recombination events. Out of twelve intact α/β-gliadin genes detected from their sequences, nine were expressed, although their patterns of expression were distinct. Since they have similar cis-elements and promoter structures, the mechanisms underlying their distinct gene expression and possible applications are discussed.

Mitotic illegitimate recombination is a mechanism for novel changes in high-molecular-weight glutenin subunits in wheat-rye hybrids

Wide hybrids can have novel traits or changed expression of a quantitative trait that their parents do not have. These phenomena have long been noticed, yet the mechanisms are poorly understood. High-molecular-weight glutenin subunits (HMW-GS) are seed storage proteins encoded by Glu-1 genes that only express in endosperm in wheat and its related species. Novel HMW-GS compositions have been observed in their hybrids. This research elucidated the molecular mechanisms by investigating the causative factors of novel HMW-GS changes in wheat-rye hybrids. HMW-GS compositions in the endosperm and their coding sequences in the leaves of F₁ and F₂ hybrids between wheat landrace Shinchunaga and rye landrace Qinling were investigated. Missing and/or additional novel HMW-GSs were observed in the endosperm of 0.5% of the 2078 F₁ and 22% of 36 F₂ hybrid seeds. The wildtype Glu-1Ax null allele was found to have 42 types of short repeat sequences of 3-60 bp long that appeared 2 to 100 times. It also has an in-frame stop codon in the central repetitive region. Analyzing cloned allele sequences of HMW-GS coding gene Glu-1 revealed that deletions involving the in-frame stop codon had happened, resulting in novel ∼1.8-kb Glu-1Ax alleles in some F₁ and F₂ plants. The cloned mutant Glu-1Ax alleles were expressed in Escherichia coli, and the HMW-GSs produced matched the novel HMW-GSs found in the hybrids. The differential changes between the endosperm and the plant of the same hybrids and the data of E. coli expression of the cloned deletion alleles both suggested that mitotic illegitimate recombination between two copies of a short repeat sequence had resulted in the deletions and thus the changed HMW-GS compositions. Our experiments have provided the first direct evidence to show that mitotic illegitimate recombination is a mechanism that produces novel phenotypes in wide hybrids.

Rapid linkage disequilibrium decay in the Lr10 gene in wild emmer wheat (Triticum dicoccoides) populations

INTRODUCTION

Recombination is a key evolutionary factor enhancing diversity. However, the effect of recombination on diversity in inbreeding species is expected to be low. To estimate this effect, recombination and diversity patterns of Lr10 gene were studied in natural populations of the inbreeder species, wild emmer wheat (Triticum dicoccoides). Wild emmer wheat is the progenitor of most cultivated wheats and it harbors rich genetic resources for disease resistance. Lr10 is a leaf rust resistance gene encoding three domains: a coiled-coil, nucleotide-binding site, and leucine-rich repeat (CC-NBS-LRR).

RESULTS

Lr10 was sequenced from 58 accessions representing 12 diverse habitats in Israel. Diversity analysis revealed a high rate of synonymous and non-synonymous substitutions (d (S) = 0.029, d (N) = 0.018, respectively) in the NBS-LRR domains. Moreover, in contrast to other resistance genes, in Lr10 the CC domain was more diverse than the NBS-LRR domains (d (S) = 0.069 vs. 0.029, d (N) = 0.094 vs. 0.018) and was subjected to positive selection in some of the populations. Seventeen recombination events were detected between haplotypes, especially in the CC domain. Linkage disequilibrium (LD) analysis has shown a rapid decay from r (2) = 0.5 to r (2) = 0.1 within a 2-kb span.

CONCLUSION

These results suggest that recombination is a diversifying force for the R-gene, Lr10, in the selfing species T. dicoccoides. This is the first report of a short-range LD decay in wild emmer wheat.

Specific resistances against Pseudomonas syringae effectors AvrB and AvrRpm1 have evolved differently in common bean (Phaseolus vulgaris), soybean (Glycine max), and Arabidopsis thaliana

*In plants, the evolution of specific resistance is poorly understood. Pseudomonas syringae effectors AvrB and AvrRpm1 are recognized by phylogenetically distinct resistance (R) proteins in Arabidopsis thaliana (Brassicaceae) and soybean (Glycine max, Fabaceae). In soybean, these resistances are encoded by two tightly linked R genes, Rpg1-b and Rpg1-r. To study the evolution of these specific resistances, we investigated AvrB- and AvrRpm1-induced responses in common bean (Phaseolus vulgaris, Fabaceae). *Common bean genotypes of various geographical origins were inoculated with P. syringae strains expressing AvrB or AvrRpm1. A common bean recombinant inbred line (RIL) population was used to map R genes to AvrRpm1. *No common bean genotypes recognized AvrB. By contrast, multiple genotypes responded to AvrRpm1, and two independent R genes conferring AvrRpm1-specific resistance were mapped to the ends of linkage group B11 (Rpsar-1, for resistance to Pseudomonas syringae effector AvrRpm1 number 1) and B8 (Rpsar-2). Rpsar-1 is located in a region syntenic with the soybean Rpg1 cluster. However, mapping of specific Rpg1 homologous genes suggests that AvrRpm1 recognition evolved independently in common bean and soybean. *The conservation of the genomic position of AvrRpm1-specific genes between soybean and common bean suggests a model whereby specific clusters of R genes are predisposed to evolve recognition of the same effector molecules.

The Rvi15 (Vr2) apple scab resistance locus contains three TIR-NBS-LRR genes

Scab caused by the pathogen Venturia inaequalis is considered the most important fungal disease of cultivated apple (Malus x domestica Borkh.). In all, 16 monogenic resistances against scab have been found in different Malus spp. and some of them are currently used in apple breeding for scab-resistant cultivars. However, the self incompatibility and the long generation time of Malus spp. together with the high standards of fruit quality demanded from the fresh market render the breeding of high-quality cultivars in apple a long and expensive task. Therefore, the cloning of disease resistance genes and the use of the cloned genes for the transformation of high-quality apple cultivars could be an approach to solve these drawbacks. We report the construction of a bacterial artificial chromosome (BAC) contig spanning the Rvi15 (Vr2) apple scab resistance locus using two GMAL 2473 BAC libraries. A single BAC clone of the contig was sufficient to span the resistance locus. The BAC clone was completely sequenced, allowing identification of a sequence of 48.6 kb going from the two closest markers (ARGH17 and 77G20RP) bracketing Rvi15 (Vr2). Analysis of the 48.6-kb sequence revealed the presence of three putative genes characterized by a Toll and mammalian interleukin-1 receptor protein nucleotide-binding site leucine-rich repeat structure. All three genes were found to be transcribed.

Sequence analysis of two alleles reveals that intra-and intergenic recombination played a role in the evolution of the radish fertility restorer (Rfo)

BACKGROUND

Land plant genomes contain multiple members of a eukaryote-specific gene family encoding proteins with pentatricopeptide repeat (PPR) motifs. Some PPR proteins were shown to participate in post-transcriptional events involved in organellar gene expression, and this type of function is now thought to be their main biological role. Among PPR genes, restorers of fertility (Rf) of cytoplasmic male sterility systems constitute a peculiar subgroup that is thought to evolve in response to the presence of mitochondrial sterility-inducing genes. Rf genes encoding PPR proteins are associated with very close relatives on complex loci.

RESULTS

We sequenced a non-restoring allele (L7rfo) of the Rfo radish locus whose restoring allele (D81Rfo) was previously described, and compared the two alleles and their PPR genes. We identified a ca 13 kb long fragment, likely originating from another part of the radish genome, inserted into the L7rfo sequence. The L7rfo allele carries two genes (PPR-1 and PPR-2) closely related to the three previously described PPR genes of the restorer D81Rfo allele (PPR-A, PPR-B, and PPR-C). Our results indicate that alleles of the Rfo locus have experienced complex evolutionary events, including recombination and insertion of extra-locus sequences, since they diverged. Our analyses strongly suggest that present coding sequences of Rfo PPR genes result from intragenic recombination. We found that the 10 C-terminal PPR repeats in Rfo PPR gene encoded proteins result from the tandem duplication of a 5 PPR repeat block.

CONCLUSIONS

The Rfo locus appears to experience more complex evolution than its flanking sequences. The Rfo locus and PPR genes therein are likely to evolve as a result of intergenic and intragenic recombination. It is therefore not possible to determine which genes on the two alleles are direct orthologs. Our observations recall some previously reported data on pathogen resistance complex loci.

Relationships between rust resistance genes at the M locus in flax

Genes at the M locus in flax (Linum usitatissimum) that confer resistance to flax rust (Melampsora lini) occur in complex haplotypes containing up to 15 related genes or gene fragments. We have cloned two additional functional resistance genes at this locus, M1 and M3, by transposon tagging and candidate gene approaches, and investigated the genetic relationships between four genes (M, M1, M3 and M4) by recombination analysis. M1 and M3, like M, are members of the nucleotide binding site, leucine-rich repeat (NBS-LRR) family. Comparisons of the predicted M1 and M3 amino acid sequences with M and L6 reveal that: (i) M1 contains four additional LRRs, probably as a result of an unequal crossover event between duplicated regions; (ii) M1 shares large segments of exact identity with M and M3, indicative of intragenic recombination events; and (iii) a large number of amino acid differences are scattered throughout the M, M1 and M3 proteins. Recombination analysis (here and in previous studies) has revealed that M readily recombines with M1, M3 and M4, whereas these three genes fail to recombine despite large family sizes (>5800) in two test-cross families, suggesting that they may occupy allelic positions in the gene cluster. Several restriction fragment length polymorphism markers within or near the M locus were mapped with respect to seven crossover events between M and M1. The results of this and previous studies provide evidence of structural differences between: (i) homoeologous loci in the different genomes of flax; (ii) different haplotypes at the M locus; (iii) different resistance genes in the M group; and (iv) the flanking regions downstream of M locus resistance genes.

Two different CC-NBS-LRR genes are required for Lr10-mediated leaf rust resistance in tetraploid and hexaploid wheat

Comparative study of disease resistance genes in crop plants and their relatives provides insight on resistance gene function, evolution and diversity. Here, we studied the allelic diversity of the Lr10 leaf rust resistance gene, a CC-NBS-LRR coding gene originally isolated from hexaploid wheat, in 20 diploid and tetraploid wheat lines. Besides a gene in the tetraploid wheat variety 'Altar' that is identical to the hexaploid wheat Lr10, two additional, functional resistance alleles showing sequence diversity were identified by virus-induced gene silencing in tetraploid wheat lines. In contrast to most described NBS-LRR proteins, the N-terminal CC domain of LR10 was found to be under strong diversifying selection. A second NBS-LRR gene at the Lr10 locus, RGA2, was shown through silencing to be essential for Lr10 function. Interestingly, RGA2 showed much less sequence diversity than Lr10. These data demonstrate allelic diversity of functional genes at the Lr10 locus in tetraploid wheat, and these new genes can now be analyzed for agronomic relevance. Lr10-based resistance is highly unusual both in its dependence on two, only distantly, related CC-NBS-LRR proteins, as well as in the pattern of diversifying selection in the N-terminal domain. This indicates a new and complex molecular mechanism of pathogen detection and signal transduction.

Conserved globulin gene across eight grass genomes identify fundamental units of the loci encoding seed storage proteins

The wheat high molecular weight (HMW) glutenins are important seed storage proteins that determine bread-making quality in hexaploid wheat (Triticum aestivum). In this study, detailed comparative sequence analyses of large orthologous HMW glutenin genomic regions from eight grass species, representing a wide evolutionary history of grass genomes, reveal a number of lineage-specific sequence changes. These lineage-specific changes, which resulted in duplications, insertions, and deletions of genes, are the major forces disrupting gene colinearity among grass genomes. Our results indicate that the presence of the HMW glutenin gene in Triticeae genomes was caused by lineage-specific duplication of a globulin gene. This tandem duplication event is shared by Brachypodium and Triticeae genomes, but is absent in rice, maize, and sorghum, suggesting the duplication occurred after Brachypodium and Triticeae genomes diverged from the other grasses ~35 Ma ago. Aside from their physical location in tandem, the sequence similarity, expression pattern, and conserved cis-acting elements responsible for endosperm-specific expression further support the paralogous relationship between the HMW glutenin and globulin genes. While the duplicated copy in Brachypodium has apparently become nonfunctional, the duplicated copy in wheat has evolved to become the HMW glutenin gene by gaining a central prolamin repetitive domain.

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.

Unlocking wheat genetic resources for the molecular identification of previously undescribed functional alleles at the Pm3 resistance locus

The continuous improvement of crop plants is essential for agriculture in the coming decades and relies on the use of genetic variability through breeding. However, domestication and modern breeding have reduced diversity in the crop germplasm. Global gene banks conserve diversity, but these resources remain underexplored owing to a lack of efficient strategies to isolate important alleles. Here we describe a large-scale allele-mining project at the molecular level. We first selected a set of 1,320 bread wheat landraces from a database of 16,089 accessions, using the focused identification of germplasm strategy. On the basis of a hierarchical selection procedure on this set, we then isolated 7 resistance alleles of the powdery mildew resistance gene Pm3, doubling the known functional allelic diversity at this locus. This targeted approach for molecular utilization of gene bank accessions reveals landraces as a rich resource of new functional alleles. This strategy can be implemented for other studies on the molecular diversity of agriculturally important genes, as well as for molecular breeding.

The fractionated orthology of Bs2 and Rx/Gpa2 supports shared synteny of disease resistance in the Solanaceae

Comparative genomics provides a powerful tool for the identification of genes that encode traits shared between crop plants and model organisms. Pathogen resistance conferred by plant R genes of the nucleotide-binding-leucine-rich-repeat (NB-LRR) class is one such trait with great agricultural importance that occupies a critical position in understanding fundamental processes of pathogen detection and coevolution. The proposed rapid rearrangement of R genes in genome evolution would make comparative approaches tenuous. Here, we test the hypothesis that orthology is predictive of R-gene genomic location in the Solanaceae using the pepper R gene Bs2. Homologs of Bs2 were compared in terms of sequence and gene and protein architecture. Comparative mapping demonstrated that Bs2 shared macrosynteny with R genes that best fit criteria determined to be its orthologs. Analysis of the genomic sequence encompassing solanaceous R genes revealed the magnitude of transposon insertions and local duplications that resulted in the expansion of the Bs2 intron to 27 kb and the frequently detected duplications of the 5'-end of R genes. However, these duplications did not impact protein expression or function in transient assays. Taken together, our results support a conservation of synteny for NB-LRR genes and further show that their distribution in the genome has been consistent with global rearrangements.

Genome sequence comparison of Col and Ler lines reveals the dynamic nature of Arabidopsis chromosomes

Large differences in plant genome sizes are mainly due to numerous events of insertions or deletions (indels). The balance between these events determines the evolutionary direction of genome changes. To address the question of what phenomena trigger these alterations, we compared the genomic sequences of two Arabidopsis thaliana lines, Columbia (Col) and Landsberg erecta (Ler). Based on the resulting alignments large indels (>100 bp) within these two genomes were analysed. There are ∼8500 large indels accounting for the differences between the two genomes. The genetic basis of their origin was distinguished as three main categories: unequal recombination (Urec)-derived, illegitimate recombination (Illrec)-derived and transposable elements (TE)-derived. A detailed study of their distribution and size variation along chromosomes, together with a correlation analyses, allowed us to demonstrate the impact of particular recombination-based mechanisms on the plant genome evolution. The results show that unequal recombination is not efficient in the removal of TEs within the pericentromeric regions. Moreover, we discovered an unexpectedly high influence of large indels on gene evolution pointing out significant differences between the various gene families. For the first time, we present convincing evidence that somatic events do play an important role in plant genome evolution.

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.

Integration of new CAPS and dCAPS-RGA markers into a composite chickpea genetic map and their association with disease resistance

A composite linkage map was constructed based on two interspecific recombinant inbred line populations derived from crosses between Cicer arietinum (ILC72 and ICCL81001) and Cicer reticulatum (Cr5-10 or Cr5-9). These mapping populations segregate for resistance to ascochyta blight (caused by Ascochyta rabiei), fusarium wilt (caused by Fusarium oxysporum f. sp. ciceris) and rust (caused by Uromyces ciceris-arietini). The presence of single nucleotide polymorphisms in ten resistance gene analogs (RGAs) previously isolated and characterized was exploited. Six out of the ten RGAs were novel sequences. In addition, classes RGA05, RGA06, RGA07, RGA08, RGA09 and RGA10 were considerate putatively functional since they matched with several legume expressed sequences tags (ESTs) obtained under infection conditions. Seven RGA PCR-based markers (5 CAPS and 2 dCAPS) were developed and successfully genotyped in the two progenies. Six of them have been mapped in different linkage groups where major quantitative trait loci conferring resistance to ascochyta blight and fusarium wilt have been reported. Genomic locations of RGAs were compared with those of known Cicer R-genes and previously mapped RGAs. Association was detected between RGA05 and genes controlling resistance to fusarium wilt caused by races 0 and 5.

Analysis of intraspecies diversity in wheat and barley genomes identifies breakpoints of ancient haplotypes and provides insight into the structure of diploid and hexaploid triticeae gene pools

A large number of wheat (Triticum aestivum) and barley (Hordeum vulgare) varieties have evolved in agricultural ecosystems since domestication. Because of the large, repetitive genomes of these Triticeae crops, sequence information is limited and molecular differences between modern varieties are poorly understood. To study intraspecies genomic diversity, we compared large genomic sequences at the Lr34 locus of the wheat varieties Chinese Spring, Renan, and Glenlea, and diploid wheat Aegilops tauschii. Additionally, we compared the barley loci Vrs1 and Rym4 of the varieties Morex, Cebada Capa, and Haruna Nijo. Molecular dating showed that the wheat D genome haplotypes diverged only a few thousand years ago, while some barley and Ae. tauschii haplotypes diverged more than 500,000 years ago. This suggests gene flow from wild barley relatives after domestication, whereas this was rare or absent in the D genome of hexaploid wheat. In some segments, the compared haplotypes were very similar to each other, but for two varieties each at the Rym4 and Lr34 loci, sequence conservation showed a breakpoint that separates a highly conserved from a less conserved segment. We interpret this as recombination breakpoints of two ancient haplotypes, indicating that the Triticeae genomes are a heterogeneous and variable mosaic of haplotype fragments. Analysis of insertions and deletions showed that large events caused by transposable element insertions, illegitimate recombination, or unequal crossing over were relatively rare. Most insertions and deletions were small and caused by template slippage in short homopolymers of only a few base pairs in size. Such frequent polymorphisms could be exploited for future molecular marker development.