Genomic and epigenomic immunity in common bean: the unusual features of NB-LRR gene family

Fine mapping of a new common bean anthracnose resistance gene (Co-18) to the proximal end of Pv10 in Indian landrace KRC-5

KEY MESSAGE

Mapping and fine mapping of bean anthracnose resistance genes is a continuous process. We report fine mapping of anthracnose resistance gene Co-18 which is the first anthracnose gene mapped to Pv10. The discovery of resistance gene is a major gain in the bean anthracnose pathosystem research. Among the Indian common bean landraces, KRC-5 exhibit high levels of resistance to the bean anthracnose pathogen Colletotrichum lindemuthianum. To precisely map the anthracnose resistance gene, we used a Recombinant Inbred Line (F2:9 RIL) population (KRC-5 × Jawala). The inheritance test revealed that KRC-5 carries a dominant resistance gene temporarily designated as Co-18. We discovered two RAPD markers linked to Co-18 among 287 RAPD markers. These RAPD markers were eventually developed into SCARs (Sc-OPR15 and Sc-OPF6) and flank Co-18 on chromosome Pv10 at a distance of 5.3 and 4.2 cM, respectively. At 4.0-4.1 Mb on Pv10, we detected a SNP (single-nucleotide polymorphism) signal. We synthesized 58 SSRs and 83 InDels from a pool of 135 SSRs and 1134 InDels, respectively. Five SSRs, four InDels, and two SCARs were used to generate the high-density linkage map, which led to the identification of two SSRs (SSR24 and SSR36) that are tightly linked to Co-18. These two SSRs flank the Co-18 to 178 kb genomic region with 13 candidate genes including five NLR (nucleotide-binding and leucine-rich repeat) genes. The closely linked markers SSR24 and SSR36 will be used in cloning and pyramiding of the Co-18 gene with other R genes to develop durable resistant bean varieties.

The pivotal ripening gene SlDML2 participates in regulating disease resistance in tomato

Fruit ripening and disease resistance are two essential biological processes for quality formation and maintenance. DNA methylation, in the form of 5-methylcytosine (5mC), has been elucidated to modulate fruit ripening, but its role in regulating fruit disease resistance remains poorly understood. In this study, we show that mutation of SlDML2, the DNA demethylase gene essential for fruit ripening, affects multiple developmental processes of tomato besides fruit ripening, including seed germination, leaf length and width and flower branching. Intriguingly, loss of SlDML2 function decreased the resistance of tomato fruits against the necrotrophic fungal pathogen Botrytis cinerea. Comparative transcriptomic analysis revealed an obvious transcriptome reprogramming caused by SlDML2 mutation during B. cinerea invasion. Among the thousands of differentially expressed genes, SlβCA3 encoding a β-carbonic anhydrase and SlFAD3 encoding a ω-3 fatty acid desaturase were demonstrated to be transcriptionally activated by SlDML2-mediated DNA demethylation and positively regulate tomato resistance to B. cinerea probably in the same genetic pathway with SlDML2. We further show that the pericarp tissue surrounding B. cinerea infection exhibited a delay in ripening with singnificant decrease in expression of ripening genes that are targeted by SlDML2 and increase in expression of SlβCA3 and SlFAD3. Taken together, our results uncover an essential layer of gene regulation mediated by DNA methylation upon B. cinerea infection and raise the possible that the DNA demethylase gene SlDML2, as a multifunctional gene, participates in modulating the trade-off between fruit ripening and disease resistance.

Genetic analysis of resistance to bean leaf crumple virus identifies a candidate LRR-RLK gene

Bean leaf crumple virus (BLCrV) is a novel begomovirus (family Geminiviridae, genus Begomovirus) infecting common bean (Phaseolus vulgaris L.), threatening bean production in Latin America. Genetic resistance is required to ensure yield stability and reduce the use of insecticides, yet the available resistance sources are limited. In this study, three common bean populations containing a total of 558 genotypes were evaluated in different yield and BLCrV resistance trials under natural infection in the field. A genome-wide association study identified the locus BLC7.1 on chromosome Pv07 at 3.31 Mbp, explaining 8 to 16% of the phenotypic variation for BLCrV resistance. In comparison, whole-genome regression models explained 51 to 78% of the variation and identified the same region on Pv07 to confer resistance. The most significantly associated markers were located within the gene model Phvul.007G040400, which encodes a leucine-rich repeat receptor-like kinase subfamily III member and is likely to be involved in the innate immune response against the virus. The allelic diversity within this gene revealed five different haplotype groups, one of which was significantly associated with BLCrV resistance. As the same genome region was previously reported to be associated with resistance against other geminiviruses affecting common bean, our study highlights the role of previous breeding efforts for virus resistance in the accumulation of positive alleles against newly emerging viruses. In addition, we provide novel diagnostic single-nucleotide polymorphism markers for marker-assisted selection to exploit BLC7.1 for breeding against geminivirus diseases in one of the most important food crops worldwide.

A Conserved, Serine-Rich Protein Plays Opposite Roles in N-Mediated Immunity against TMV and N-Triggered Cell Death

Plant nucleotide-binding, leucine-rich, repeat-containing proteins (NLRs) play important roles in plant immunity. NLR expression and function are tightly regulated by multiple mechanisms. In this study, a conserved serine/arginine-rich protein (SR protein) was identified through the yeast one-hybrid screening of a tobacco cDNA library using DNA fragments from the N gene, an NLR that confers immunity to tobacco mosaic virus (TMV). This SR protein showed an interaction with a 3' genomic regulatory sequence (GRS) and has a potential role in regulating the alternative splicing of N. Thus, it was named SR regulator for N, abbreviated SR4N. Further study showed that SR4N plays a positive role in N-mediated cell death but a negative role in N protein accumulation. SR4N also promotes multiple virus replications in co-expression experiments, and this enhancement may not function through RNA silencing suppression, as it did not enhance 35S-GFP expression in co-infiltration experiments. Bioinformatic and molecular studies revealed that SR4N belongs to the SR2Z subtype of the SR protein family, which was conserved in both dicots and monocots, and its roles in repressing viral immunity and triggering cell death were also conserved. Our study revealed new roles for SR2Z family proteins in plant immunity against viruses.

Genome-Wide Transcriptomic Analysis of the Effects of Infection with the Hemibiotrophic Fungus Colletotrichum lindemuthianum on Common Bean

Bean anthracnose caused by the hemibiotrophic fungus Colletotrichum lindemuthianum is one of the most important diseases of common bean (Phaseolus vulgaris) in the world. In the present study, the whole transcriptome of common bean infected with C. lindemuthianum during compatible and incompatible interactions was characterized at 48 and 72 hpi, corresponding to the biotrophy phase of the infection cycle. Our results highlight the prominent role of pathogenesis-related (PR) genes from the PR10/Bet vI family as well as a complex interplay of different plant hormone pathways including Ethylene, Salicylic acid (SA) and Jasmonic acid pathways. Gene Ontology enrichment analysis reveals that infected common bean seedlings responded by down-regulation of photosynthesis, ubiquitination-mediated proteolysis and cell wall modifications. In infected common bean, SA biosynthesis seems to be based on the PAL pathway instead of the ICS pathway, contrarily to what is described in Arabidopsis. Interestingly, ~30 NLR were up-regulated in both contexts. Overall, our results suggest that the difference between the compatible and incompatible reaction is more a question of timing and strength, than a massive difference in differentially expressed genes between these two contexts. Finally, we used RT-qPCR to validate the expression patterns of several genes, and the results showed an excellent agreement with deep sequencing.

Whole-Genome DNA Methylation Analysis in Brassica rapa subsp. perviridis in Response to Albugo candida Infection

DNA methylation is an epigenetic mark associated with several mechanisms in plants including immunity mechanisms. However, little is known about the regulatory role of DNA methylation in the resistance response of Brassica species against fungal diseases. White rust, caused by the fungus Albugo candida, is one of the most widespread and destructive diseases of all the cultivated Brassica species, particularly Brassica rapa L. and Brassica juncea (L.) Czern and Coss. Here, we investigate whole-genome DNA methylation modifications of B. rapa subsp. perviridis in response to white rust. As a result, 233 and 275 differentially methylated regions (DMRs) in the susceptible cultivar "Misugi" and the resistant cultivar "Nanane" were identified, respectively. In both cultivars, more than half of the DMRs were associated with genes (DMR-genes). Gene expression analysis showed that 13 of these genes were also differentially expressed between control and infected samples. Gene ontology enrichment analysis of DMR genes revealed their involvement in various biological processes including defense mechanisms. DMRs were unevenly distributed around genes in susceptible and resistant cultivars. In "Misugi," DMRs tended to be located within genes, while in "Nanane," DMRs tended to be located up and downstream of the genes. However, CG DMRs were predominantly located within genes in both cultivars. Transposable elements also showed association with all three sequence contexts of DMRs but predominantly with CHG and CHH DMRs in both cultivars. Our findings indicate the occurrence of DNA methylation modifications in B. rapa in response to white rust infection and suggest a potential regulatory role of DNA methylation modification in defense mechanisms which could be exploited to improve disease resistance.

Genome-wide identification and characterization of NBS-encoding genes in the sweet potato wild ancestor Ipomoea trifida (H.B.K.)

The most predominant type of resistance (R) genes contain nucleotide-binding sites and leucine-rich repeat (NBS-LRR) domains, characterization of which is helpful for plant resistance improvement. However, the NBS genes of Ipomoea trifida (H.B.K.) remain insufficient to date. In this study, a genome-wide analysis of the NBS-encoding gene in I. trifida (H.B.K.) was carried out. A total of 442 NBS encoding genes were identified, amounting to 1.37% of the total genes of I. trifida (H.B.K.). Based on the analysis of the domains, the identified ItfNBS genes were further classified into seven groups: CNL, NL, CN, N, TNL, TN, and RNL. Phylogenetic analysis showed that the I. trifida NBS genes clustered into three independent clades: RNL, TNL, and CNL. Chromosome location analysis revealed that the distribution of ItfNBS genes in chromosomes was uneven, with a number ranging from 3 to 45. Multiple stress-related regulatory elements were detected in the promoters of the NBS-encoding genes, and their expression profiles were obtained. The qRT-PCR analysis revealed that IbNBS10, IbNBS20, IbNBS258, and IbNBS88 responded to stem nematode infection. These results provide critical proof for further characterization and analysis of NBS-encoding genes with important functions.

Phaseolus vulgaris-Colletotrichum lindemuthianum Pathosystem in the Post-Genomic Era: An Update

Phaseolus vulgaris-Colletotrichum lindemuthianum is one among the oldest host and pathogen interface. Researchers have taken painstaking efforts across the world for understanding the dialogue during early and late phases of interaction. Collectively, these efforts resulted in the deluge of information that helped the researchers to underpin the interface. The latest molecular biology techniques furnished novel detection methods for the anthracnose pathogen, refined the understanding of pathogen population dynamics, and provided the insights on co-evolutionary common bean resistance and C. lindemuthianum virulence dynamics. One of the important breakthroughs came when the Phaseolus vulgaris and its corresponding anthracnose pathogen (C. lindemuthianum) genomes were decoded in 2014 and 2017, respectively. Availability of both the genomes yielded a significant genomic information that helped bean communities to fine map the economically important traits and to identify the pathogenicity determinants and effector molecules. The interface is in a continuous development as knowledge of the anthracnose resistance genes, their precise physical locations, and the identification of effector proteins; the fungus arsenals are being routinely updated. Hence, we revisited the interface and tried to provide an overview of host pathogen dialogue in the genomic era. Additionally, we compiled the sporadic information on this pathosystem from India and provided its futuristic road map to shape its research in the world and northern India, the major dry bean area in the country.

Epigenetics: a catalyst of plant immunity against pathogens

The plant immune system protects against pests and diseases. The recognition of stress-related molecular patterns triggers localised immune responses, which are often followed by longer-lasting systemic priming and/or up-regulation of defences. In some cases, this induced resistance (IR) can be transmitted to following generations. Such transgenerational IR is gradually reversed in the absence of stress at a rate that is proportional to the severity of disease experienced in previous generations. This review outlines the mechanisms by which epigenetic responses to pathogen infection shape the plant immune system across expanding time scales. We review the cis- and trans-acting mechanisms by which stress-inducible epigenetic changes at transposable elements (TEs) regulate genome-wide defence gene expression and draw particular attention to one regulatory model that is supported by recent evidence about the function of AGO1 and H2A.Z in transcriptional control of defence genes. Additionally, we explore how stress-induced mobilisation of epigenetically controlled TEs acts as a catalyst of Darwinian evolution by generating (epi)genetic diversity at environmentally responsive genes. This raises questions about the long-term evolutionary consequences of stress-induced diversification of the plant immune system in relation to the long-held dichotomy between Darwinian and Lamarckian evolution.

Genome-Wide Identification of Key Components of RNA Silencing in Two Phaseolus vulgaris Genotypes of Contrasting Origin and Their Expression Analyses in Response to Fungal Infection

RNA silencing serves key roles in a multitude of cellular processes, including development, stress responses, metabolism, and maintenance of genome integrity. Dicer, Argonaute (AGO), double-stranded RNA binding (DRB) proteins, RNA-dependent RNA polymerase (RDR), and DNA-dependent RNA polymerases known as Pol IV and Pol V form core components to trigger RNA silencing. Common bean (Phaseolus vulgaris) is an important staple crop worldwide. In this study, we aimed to unravel the components of the RNA-guided silencing pathway in this non-model plant, taking advantage of the availability of two genome assemblies of Andean and Meso-American origin. We identified six PvDCLs, thirteen PvAGOs, 10 PvDRBs, 5 PvRDRs, in both genotypes, suggesting no recent gene amplification or deletion after the gene pool separation. In addition, we identified one PvNRPD1 and one PvNRPE1 encoding the largest subunits of Pol IV and Pol V, respectively. These genes were categorized into subgroups based on phylogenetic analyses. Comprehensive analyses of gene structure, genomic localization, and similarity among these genes were performed. Their expression patterns were investigated by means of expression models in different organs using online data and quantitative RT-PCR after pathogen infection. Several of the candidate genes were up-regulated after infection with the fungus Colletotrichum lindemuthianum.

Genome assembly of a Mesoamerican derived variety of lima bean: a foundational cultivar in the Mid-Atlantic USA

Lima bean, Phaseolus lunatus, is closely related to common bean and is high in fiber and protein, with a low glycemic index. Lima bean is widely grown in the state of Delaware, where late summer and early fall weather are conducive to pod production. The same weather conditions also promote diseases such as pod rot and downy mildew, the latter of which has caused previous epidemics. A better understanding of the genes underlying resistance to this and other pathogens is needed to keep this industry thriving in the region. Our current study sought to sequence, assemble, and annotate a commercially available cultivar called Bridgeton, which could then serve as a reference genome, a basis of comparison to other Phaseolus taxa, and a resource for the identification of potential resistance genes. Combined efforts of sequencing, linkage, and comparative analysis resulted in a 623 Mb annotated assembly for lima bean, as well as a better understanding of an evolutionarily dynamic resistance locus in legumes.

Characterization of Nucleotide Binding -Site-Encoding Genes in Sweetpotato, Ipomoea batatas(L.) Lam., and Their Response to Biotic and Abiotic Stresses

Sweetpotato, Ipomoea batatas (L.) Lam., is an important and widely grown crop, yet its production is affected severely by biotic and abiotic stresses. The nucleotide binding site (NBS)-encoding genes have been shown to improve stress tolerance in several plant species. However, the characterization of NBS-encoding genes in sweetpotato is not well-documented to date. In this study, a comprehensive analysis of NBS-encoding genes has been conducted on this species by using bioinformatics and molecular biology methods. A total of 315 NBS-encoding genes were identified, and 260 of them contained all essential conserved domains while 55 genes were truncated. Based on domain architectures, the 260 NBS-encoding genes were grouped into 6 distinct categories. Phylogenetic analysis grouped these genes into 3 classes: TIR, CC (I), and CC (II). Chromosome location analysis revealed that the distribution of NBS-encoding genes in chromosomes was uneven, with a number ranging from 1 to 34. Multiple stress-related regulatory elements were detected in the promoters, and the NBS-encoding genes' expression profiles under biotic and abiotic stresses were obtained. According to the bioinformatics analysis, 9 genes were selected for RT-qPCR analysis. The results revealed that IbNBS75, IbNBS219, and IbNBS256 respond to stem nematode infection; Ib-NBS240, IbNBS90, and IbNBS80 respond to cold stress, while IbNBS208, IbNBS71, and IbNBS159 respond to 30% PEG treatment. We hope these results will provide new insights into the evolution of NBS-encoding genes in the sweetpotato genome and contribute to the molecular breeding of sweetpotato in the future.

Multifaceted Chromatin Structure and Transcription Changes in Plant Stress Response

Sessile plants are exposed throughout their existence to environmental abiotic and biotic stress factors, such as cold, heat, salinity, drought, dehydration, submergence, waterlogging, and pathogen infection. Chromatin organization affects genome stability, and its dynamics are crucial in plant stress responses. Chromatin dynamics are epigenetically regulated and are required for stress-induced transcriptional regulation or reprogramming. Epigenetic regulators facilitate the phenotypic plasticity of development and the survival and reproduction of plants in unfavorable environments, and they are highly diversified, including histone and DNA modifiers, histone variants, chromatin remodelers, and regulatory non-coding RNAs. They contribute to chromatin modifications, remodeling and dynamics, and constitute a multilayered and multifaceted circuitry for sophisticated and robust epigenetic regulation of plant stress responses. However, this complicated epigenetic regulatory circuitry creates challenges for elucidating the common or differential roles of chromatin modifications for transcriptional regulation or reprogramming in different plant stress responses. Particularly, interacting chromatin modifications and heritable stress memories are difficult to identify in the aspect of chromatin-based epigenetic regulation of transcriptional reprogramming and memory. Therefore, this review discusses the recent updates from the three perspectives—stress specificity or dependence of transcriptional reprogramming, the interplay of chromatin modifications, and transcriptional stress memory in plants. This helps solidify our knowledge on chromatin-based transcriptional reprogramming for plant stress response and memory.

Comprehensive genomic resources related to domestication and crop improvement traits in Lima bean

Lima bean (Phaseolus lunatus L.), one of the five domesticated Phaseolus bean crops, shows a wide range of ecological adaptations along its distribution range from Mexico to Argentina. These adaptations make it a promising crop for improving food security under predicted scenarios of climate change in Latin America and elsewhere. In this work, we combine long and short read sequencing technologies with a dense genetic map from a biparental population to obtain the chromosome-level genome assembly for Lima bean. Annotation of 28,326 gene models show high diversity among 1917 genes with conserved domains related to disease resistance. Structural comparison across 22,180 orthologs with common bean reveals high genome synteny and five large intrachromosomal rearrangements. Population genomic analyses show that wild Lima bean is organized into six clusters with mostly non-overlapping distributions and that Mesomerican landraces can be further subdivided into three subclusters. RNA-seq data reveal 4275 differentially expressed genes, which can be related to pod dehiscence and seed development. We expect the resources presented here to serve as a solid basis to achieve a comprehensive view of the degree of convergent evolution of Phaseolus species under domestication and provide tools and information for breeding for climate change resiliency.

Convergent Epigenetic Mechanisms Avoid Constitutive Expression of Immune Receptor Gene Subsets

The gene pool encoding PRR and NLR immune receptors determines the ability of a plant to resist microbial infections. Basal expression of these genes is prevented by diverse mechanisms since their hyperactivity can be harmful. To approach the study of epigenetic control of PRR/NLR genes we here analyzed their expression in mutants carrying abnormal repressive 5-methyl cytosine (5-mC) and histone 3 lysine 9 dimethylation (H3K9me2) marks, due to lack of MET1, CMT3, MOM1, SUVH4/5/6, or DDM1. At optimal growth conditions, none of the mutants showed basal expression of the defense gene marker PR1, but all of them had greater resistance to Pseudomonas syringae pv. tomato than wild type plants, suggesting they are primed to stimulate immune cascades. Consistently, analysis of available transcriptomes indicated that all mutants showed activation of particular PRR/NLR genes under some growth conditions. Under low defense activation, 37 PRR/NLR genes were expressed in these plants, but 29 of them were exclusively activated in specific mutants, indicating that MET1, CMT3, MOM1, SUVH4/5/6, and DDM1 mediate basal repression of different subsets of genes. Some epigenetic marks present at promoters, but not gene bodies, could explain the activation of these genes in the mutants. As expected, suvh4/5/6 and ddm1 activated genes carrying 5-mC and H3K9me2 marks in wild type plants. Surprisingly, all mutants expressed genes harboring promoter H2A.Z/H3K27me3 marks likely affected by the chromatin remodeler PIE1 and the histone demethylase REF6, respectively. Therefore, MET1, CMT3, MOM1, SUVH4/5/6, and DDM1, together with REF6, seemingly contribute to the establishment of chromatin states that prevent constitutive PRR/NLR gene activation, but facilitate their priming by modulating epigenetic marks at their promoters.

North-Western Himalayan Common Beans: Population Structure and Mapping of Quantitative Anthracnose Resistance Through Genome Wide Association Study

Common bean (Phaseolus vulgaris L.) is an important legume crop of north-western (NW) Himalayan region and the major disease that causes catastrophic loss to the crop is anthracnose, which is caused by Colletotrichum lindemuthianum. The pathogen is highly diverse and most of the commercial cultivars are susceptible to different races prevalent in the region. The lack of information on the genomic regions associated with anthracnose resistance in NW Himalayan common bean population prompted us to dissect Quantitative Resistance Loci (QRLs) against major anthracnose races. In this study, 188 common bean landraces collected from NW region were screened against five important anthracnose races and 113 bean genotypes showed resistance to one or multiple races. Genotyping by sequencing (GBS) was performed on a panel of 192 bean lines (4 controls plus 188 Indian beans) and 22,589 SNPs were obtained that are evenly distributed. Population structure analysis of 192 bean genotypes categorized 188 Indian beans into two major clusters representing Andean and Mesoamerican gene pools with obvious admixtures. Many QRLs associated with anthracnose resistance to Indian C. lindemuthianum virulences (race 3, 87, and 503) are located at Pv04 within the gene models that encode typical resistance gene signatures. The QRLs associated with race 73 are located on Pv08 and overlaps with Co-4 anthracnose resistance gene. A SNP located at distal end of Pv11 in a gene model Phvul.011G202300 which encodes a LRR with a typical NB-ARC domain showed association with race 73 resistance. Common bean genomic regions located at Pv03, Pv09, and Pv11 showed association with resistance to anthracnose race 2047. The present study showed presence of many novel bean genomic regions associated with anthracnose resistance. The presence of Co-4 and Co-2 genes in our material is encouraging for breeding durable anthracnose resistant cultivars for the region.

Common bean resistance to Xanthomonas is associated with upregulation of the salicylic acid pathway and downregulation of photosynthesis

BACKGROUND

Common bacterial blight (CBB) caused by Xanthomonas phaseoli pv. phaseoli and Xanthomonas citri pv. fuscans is one of the major threats to common bean crops (Phaseolus vulgaris L.). Resistance to CBB is particularly complex as 26 quantitative resistance loci to CBB have been described so far. To date, transcriptomic studies after CBB infection have been very scarce and the molecular mechanisms underlying susceptibility or resistance are largely unknown.

RESULTS

We sequenced and annotated the genomes of two common bean genotypes being either resistant (BAT93) or susceptible (JaloEEP558) to CBB. Reciprocal BLASTp analysis led to a list of 20,787 homologs between these genotypes and the common bean reference genome (G19833), which provides a solid dataset for further comparative analyses. RNA-Seq after inoculation with X. phaseoli pv. phaseoli showed that the susceptible genotype initiated a more intense and diverse biological response than the resistant genotype. Resistance was linked to upregulation of the salicylic acid pathway and downregulation of photosynthesis and sugar metabolism, while susceptibility was linked to downregulation of resistance genes and upregulation of the ethylene pathway and of genes involved in cell wall modification.

CONCLUSIONS

This study helps better understanding the mechanisms occurring during the early colonization phase of common bean by Xanthomonas and unveils new actors potentially important for resistance and susceptibility to CBB. We discuss the potential link between the pathways induced during bean colonization and genes induced by transcription activator-like effectors (TALEs), as illustrated in other Xanthomonas pathovars.

Common bean resistance to Xanthomonas is associated with upregulation of the salicylic acid pathway and downregulation of photosynthesis

Background

Common bacterial blight (CBB) caused by Xanthomonas phaseoli pv. phaseoli and Xanthomonas citri pv. fuscans is one of the major threats to common bean crops (Phaseolus vulgaris L.). Resistance to CBB is particularly complex as 26 quantitative resistance loci to CBB have been described so far. To date, transcriptomic studies after CBB infection have been very scarce and the molecular mechanisms underlying susceptibility or resistance are largely unknown.

Results

We sequenced and annotated the genomes of two common bean genotypes being either resistant (BAT93) or susceptible (JaloEEP558) to CBB. Reciprocal BLASTp analysis led to a list of 20,787 homologs between these genotypes and the common bean reference genome (G19833), which provides a solid dataset for further comparative analyses. RNA-Seq after inoculation with X. phaseoli pv. phaseoli showed that the susceptible genotype initiated a more intense and diverse biological response than the resistant genotype. Resistance was linked to upregulation of the salicylic acid pathway and downregulation of photosynthesis and sugar metabolism, while susceptibility was linked to downregulation of resistance genes and upregulation of the ethylene pathway and of genes involved in cell wall modification.

Conclusions

This study helps better understanding the mechanisms occurring during the early colonization phase of common bean by Xanthomonas and unveils new actors potentially important for resistance and susceptibility to CBB. We discuss the potential link between the pathways induced during bean colonization and genes induced by transcription activator-like effectors (TALEs), as illustrated in other Xanthomonas pathovars.

Impact of DNA Demethylases on the DNA Methylation and Transcription of Arabidopsis NLR Genes

Active DNA demethylation is an important epigenetic process that plays a key role in maintaining normal gene expression. In plants, active DNA demethylation is mediated by DNA demethylases, including ROS1, DME, DML2, and DML3. In this study, the available bisulfite sequencing and mRNA sequencing data from ros1 and rdd mutants were analyzed to reveal how the active DNA demethylation process shapes the DNA methylation patterns of Arabidopsis nucleotide-binding leucine-rich repeat (NLR) genes, a class of important plant disease resistance genes. We demonstrate that the CG methylation levels of three NLR genes (AT5G49140, AT5G35450, and AT5G36930) are increased in the ros1 mutants relative to the wild-type plants, whereas the CG methylation level of AT2G17050 is decreased. We also observed increased CG methylation levels of AT4G11170 and AT5G47260 and decreased CG methylation levels of AT5G38350 in rdd mutants. We further found that the expression of three NLR genes (AT1G12280, AT1G61180, and AT4G19520) was activated in both ros1 and rdd mutants, whereas the expression of another three NLR genes (AT1G58602, AT1G59620, and AT1G62630) was repressed in these two mutants. Quantitative reverse transcriptase–polymerase chain reaction detection showed that the expression levels of AT1G58602.1, AT4G19520.3, AT4G19520.4, and AT4G19520.5 were decreased in the ros1 mutant; AT3G50950.1 and AT3G50950.2 in the rdd mutant were also decreased in expression compared to Col-0, whereas AT1G57630.1, AT1G58602.2, and AT5G45510.1 were upregulated in the rdd mutant relative to Col-0. These results indicate that some NLR genes are regulated by DNA demethylases. Our study demonstrates that each DNA demethylase (ROS1, DML2, and DML3) exerts a specific effect on the DNA methylation of the NLR genes, and active DNA demethylation is part of the regulation of DNA methylation and transcriptional activity of some Arabidopsis NLR genes.

Gene duplication and transposition of mobile elements drive evolution of the Rpv3 resistance locus in grapevine

A wild grape haplotype (Rpv3-1) confers resistance to Plasmopara viticola. We mapped the causal factor for resistance to an interval containing a TIR-NB-LRR (TNL) gene pair that originated 1.6-2.6 million years ago by a tandem segmental duplication. Transient coexpression of the TNL pair in Vitis vinifera leaves activated pathogen-induced necrosis and reduced sporulation compared with control leaves. Even though transcripts of the TNL pair from the wild haplotype appear to be partially subject to nonsense-mediated mRNA decay, mature mRNA levels in a homozygous resistant genotype were individually higher than the mRNA trace levels observed for the orthologous single-copy TNL in sensitive genotypes. Allelic expression imbalance in a resistant heterozygote confirmed that cis-acting regulatory variation promotes expression in the wild haplotype. The movement of transposable elements had a major impact on the generation of haplotype diversity, altering the DNA context around similar TNL coding sequences and the GC-content in their proximal 5'-intergenic regions. The wild and domesticated haplotypes also diverged in conserved single-copy intergenic DNA, but the highest divergence was observed in intraspecific and not in interspecific comparisons. In this case, introgression breeding did not transgress the genetic boundaries of the domesticated species, because haplotypes present in modern varieties sometimes predate speciation events between wild and cultivated species.

DNA Methylation: Toward Crop Disease Resistance Improvement

Crop diseases, in conjunction with climate change, are a major threat to global crop production. DNA methylation is an epigenetic mark and is involved in plants' biological processes, including development, stress adaptation, and genome evolution. By providing a new source of variation, DNA methylation introduces novel direction to both scientists and breeders with its potential in disease resistance enhancement. Here, we discuss the impact of pathogen-induced DNA methylation modifications on a host's transcriptome reprogramming and genome stability, as part of the plant's defense mechanisms. We also highlight the knowledge gaps that need to be investigated for understanding the entire role of DNA methylation in plant pathogen interactions. This will ultimately assist breeders toward improving resistance and decreasing yield losses.

An evolutionarily conserved non-synonymous SNP in a leucine-rich repeat domain determines anthracnose resistance in watermelon

A non-synonymous SNP of CC-NBS-LRR was firstly mapped to confer resistance to anthracnose in watermelon. Newly proposed LRR domain harboring the SNP is evolutionary conserved in the Cucurbitaceae and Fabaceae. Anthracnose disease caused by Colletotrichum devastates many plants. Despite the importance of the disease, the mechanisms of resistance against it are poorly understood. Here, we identified a non-synonymous single-nucleotide polymorphism (SNP) located in a leucine-rich repeat domain as a marker for resistance to anthracnose race 1 in watermelon, using a combination of genetic analyses. We validated this SNP in segregating populations and 59 watermelon accessions using high-resolution melting assays and Sanger sequencing. We demonstrated that the resulting arginine-to-lysine substitution is particularly conserved among the Cucurbitaceae and Fabaceae. We identified a conserved motif, IxxLPxSxxxLYNLQTLxL, found in 1007 orthologues/paralogues from 89 plant species, and discovered that residue 18 of this motif could determine resistance to disease caused by external invaders. This study provides a step forward in understanding anthracnose resistance in watermelon, as well as functional and evolutionary insight into leucine-rich repeat proteins.

BABA-Induced DNA Methylome Adjustment to Intergenerational Defense Priming in Potato to Phytophthora infestans

We provide evidence that alterations in DNA methylation patterns contribute to the regulation of stress-responsive gene expression for an intergenerational resistance of β-aminobutyric acid (BABA)-primed potato to Phytophthora infestans. Plants exposed to BABA rapidly modified their methylation capacity toward genome-wide DNA hypermethylation. De novo induced DNA methylation (5-mC) correlated with the up-regulation of Chromomethylase 3 (CMT3), Domains rearranged methyltransferase 2 (DRM2), and Repressor of silencing 1 (ROS1) genes in potato. BABA transiently activated DNA hypermethylation in the promoter region of the R3a resistance gene triggering its downregulation in the absence of the oomycete pathogen. However, in the successive stages of priming, an excessive DNA methylation state changed into demethylation with the active involvement of potato DNA glycosylases. Interestingly, the 5-mC-mediated changes were transmitted into the next generation in the form of intergenerational stress memory. Descendants of the primed potato, which derived from tubers or seeds carrying the less methylated R3a promoter, showed a higher transcription of R3a that associated with an augmented intergenerational resistance to virulent P. infestans when compared to the inoculated progeny of unprimed plants. Furthermore, our study revealed that enhanced transcription of some SA-dependent genes (NPR1, StWRKY1, and PR1) was not directly linked with DNA methylation changes in the promoter region of these genes, but was a consequence of methylation-dependent alterations in the transcriptional network.

Molecular mechanisms that limit the costs of NLR-mediated resistance in plants

Crop diseases cause significant yield losses, and the use of resistant cultivars can effectively mitigate these losses and control many plant diseases. Most plant resistance (R) genes encode immune receptors composed of nucleotide-binding and leucine-rich repeat (NLR) domains. These proteins mediate the specific recognition of pathogen avirulence effectors to induce defence responses. However, NLR-triggered immunity can be associated with a reduction in growth and yield, so-called 'fitness costs'. Recent data have shown that plants use an elaborate interplay of different mechanisms to control NLR gene transcript levels, as well as NLR protein abundance and activity, to avoid the associated cost of resistance in the absence of a pathogen. In this review, we discuss the different levels of NLR regulation (transcriptional, post-transcriptional and at the protein level). We address the apparent need for plants to maintain diverse modes of regulation. A recent model suggesting an equilibrium 'ON/OFF state' of NLR proteins, in the absence of a pathogen, provides the context for our discussion.

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.

Common Bean Subtelomeres Are Hot Spots of Recombination and Favor Resistance Gene Evolution

Subtelomeres of most eukaryotes contain fast-evolving genes usually involved in adaptive processes. In common bean (Phaseolus vulgaris), the Co-2 anthracnose resistance (R) locus corresponds to a cluster of nucleotide-binding-site leucine-rich-repeat (NL) encoding sequences, the prevalent class of plant R genes. To study the recent evolution of this R gene cluster, we used a combination of sequence, genetic and cytogenetic comparative analyses between common bean genotypes from two distinct gene pools (Andean and Mesoamerican) that diverged 0.165 million years ago. Co-2 is a large subtelomeric cluster on chromosome 11 comprising from 32 (Mesoamerican) to 52 (Andean) NL sequences embedded within khipu satellite repeats. Since the recent split between Andean and Mesoamerican gene pools, the Co-2 cluster has experienced numerous gene-pool specific NL losses, leading to distinct NL repertoires. The high proportion of solo-LTR retrotransposons indicates that the Co-2 cluster is located in a hot spot of unequal intra-strand homologous recombination. Furthermore, we observe large segmental duplications involving both Non-Homologous End Joining and Homologous Recombination double-strand break repair pathways. Finally, the identification of a Mesoamerican-specific subtelomeric sequence reveals frequent interchromosomal recombinations between common bean subtelomeres. Altogether, our results highlight that common bean subtelomeres are hot spots of recombination and favor the rapid evolution of R genes. We propose that chromosome ends could act as R gene incubators in many plant genomes.