Rapidly evolving R genes in diverse grass species confer resistance to rice blast disease

Recombination and repeat-induced point mutation landscapes reveal trade-offs between the sexual and asexual cycles of Magnaporthe oryzae

The fungal disease caused by Magnaporthe oryzae is one of the most devastating diseases that endanger many crops worldwide. Evidence shows that sexual reproduction can be advantageous for fungal diseases as hybridization facilitates host-jumping. However, the pervasive clonal lineages of M. oryzae observed in natural fields contradict this expectation. A better understanding of the roles of recombination and the fungi-specific repeat-induced point mutation (RIP) in shaping its evolutionary trajectory is essential to bridge this knowledge gap. Here we systematically investigate the RIP and recombination landscapes in M. oryzae using a whole genome sequencing data from 252 population samples and 92 cross progenies. Our data reveal that the RIP can robustly capture the population history of M. oryzae, and we provide accurate estimations of the recombination and RIP rates across different M. oryzae clades. Significantly, our results highlight a parent-of-origin bias in both recombination and RIP rates, tightly associating with their sexual potential and variations of effector proteins. This bias suggests a critical trade-off between generating novel allelic combinations in the sexual cycle to facilitate host-jumping and stimulating transposon-associated diversification of effectors in the asexual cycle to facilitate host coevolution. These findings provide unique insights into understanding the evolution of blast fungus.

Genetic architecture of disease resistance and tolerance in Douglas-fir trees

Understanding the genetic basis of how plants defend against pathogens is important to monitor and maintain resilient tree populations. Swiss needle cast (SNC) and Rhabdocline needle cast (RNC) epidemics are responsible for major damage of forest ecosystems in North America. Here we investigate the genetic architecture of tolerance and resistance to needle cast diseases in Douglas-fir (Pseudotsuga menziesii) caused by two fungal pathogens: SNC caused by Nothophaeocryptopus gaeumannii, and RNC caused by Rhabdocline pseudotsugae. We performed case-control genome-wide association analyses and found disease resistance and tolerance in Douglas-fir to be polygenic and under strong selection. We show that stomatal regulation as well as ethylene and jasmonic acid pathways are important for resisting SNC infection, and secondary metabolite pathways play a role in tolerating SNC once the plant is infected. We identify a major transcriptional regulator of plant defense, ERF1, as the top candidate for RNC resistance. Our findings shed light on the highly polygenic architectures underlying fungal disease resistance and tolerance and have important implications for forestry and conservation as the climate changes.

Fine-mapping and evolutionary history of R-BPMV, a dominant resistance gene to Bean pod mottle virus in Phaseolus vulgaris L

KEY MESSAGE

R-BPMV is located within a recently expanded TNL cluster in the Phaseolus genus with suppressed recombination and known for resistance to multiple pathogens including potyviruses controlled by the I gene. Bean pod mottle virus (BPMV) is a comovirus that infects common bean and legumes in general. BPMV is distributed throughout the world and is a major threat on soybean, a closely related species of common bean. In common bean, BAT93 was reported to carry the R-BPMV resistance gene conferring resistance to BPMV and linked with the I resistance gene. To fine map R-BPMV, 182 recombinant inbred lines (RILs) derived from the cross BAT93 × JaloEEP558 were genotyped with polymerase chain reaction (PCR)-based markers developed using genome assemblies from G19833 and BAT93, as well as BAT93 BAC clone sequences. Analysis of RILs carrying key recombination events positioned R-BPMV to a target region containing at least 16 TIR-NB-LRR (TNL) sequences in BAT93. Because the I cluster presents a suppression of recombination and a large number of repeated sequences, none of the 16 TNLs could be excluded as R-BPMV candidate gene. The evolutionary history of the TNLs for the I cluster were reconstructed using microsynteny and phylogenetic analyses within the legume family. A single I TNL was present in Medicago truncatula and lost in soybean, mirroring the absence of complete BPMV resistance in soybean. Amplification of TNLs in the I cluster predates the divergence of the Phaseolus species, in agreement with the emergence of R-BPMV before the separation of the common bean wild centers of diversity. This analysis provides PCR-based markers useful in marker-assisted selection (MAS) and laid the foundation for cloning of R-BPMV resistance gene in order to transfer the resistance into soybean.

Structure-Function Characterisation of Eop1 Effectors from the Erwinia-Pantoea Clade Reveals They May Acetylate Their Defence Target through a Catalytic Dyad

The YopJ group of acetylating effectors from phytopathogens of the genera Pseudomonas and Ralstonia have been widely studied to understand how they modify and suppress their host defence targets. In contrast, studies on a related group of effectors, the Eop1 group, lag far behind. Members of the Eop1 group are widely present in the Erwinia-Pantoea clade of Gram-negative bacteria, which contains phytopathogens, non-pathogens and potential biocontrol agents, implying that they may play an important role in agroecological or pathological adaptations. The lack of research in this group of YopJ effectors has left a significant knowledge gap in their functioning and role. For the first time, we perform a comparative analysis combining AlphaFold modelling, in planta transient expressions and targeted mutational analyses of the Eop1 group effectors from the Erwinia-Pantoea clade, to help elucidate their likely activity and mechanism(s). This integrated study revealed several new findings, including putative binding sites for inositol hexakisphosphate and acetyl coenzyme A and newly postulated target-binding domains, and raises questions about whether these effectors function through a catalytic triad mechanism. The results imply that some Eop1s may use a catalytic dyad acetylation mechanism that we found could be promoted by the electronegative environment around the active site.

Molecular and genetic insights into secondary metabolic regulation underlying insect-pest resistance in legumes

Insect pests pose a major threat to agricultural production, resulting in significant economic losses for countries. A high infestation of insects in any given area can severely reduce crop yield and quality. This review examines the existing resources for managing insect pests and highlights alternative eco-friendly techniques to enhance insect pest resistance in legumes. Recently, the application of plant secondary metabolites has gained popularity in controlling insect attacks. Plant secondary metabolites encompass a wide range of compounds such as alkaloids, flavonoids, and terpenoids, which are often synthesized through intricate biosynthetic pathways. Classical methods of metabolic engineering involve manipulating key enzymes and regulatory genes to enhance or redirect the production of secondary metabolites in plants. Additionally, the role of genetic approaches, such as quantitative trait loci mapping, genome-wide association (GWAS) mapping, and metabolome-based GWAS in insect pest management is discussed, also, the role of precision breeding, such as genome editing technologies and RNA interference for identifying pest resistance and manipulating the genome to develop insect-resistant cultivars are explored, highlighting the positive contribution of plant secondary metabolites engineering-based resistance against insect pests. It is suggested that by understanding the genes responsible for beneficial metabolite compositions, future research might hold immense potential to shed more light on the molecular regulation of secondary metabolite biosynthesis, leading to advancements in insect-resistant traits in crop plants. In the future, the utilization of metabolic engineering and biotechnological methods may serve as an alternative means of producing biologically active, economically valuable, and medically significant compounds found in plant secondary metabolites, thereby addressing the challenge of limited availability.

Structural mechanism of heavy metal-associated integrated domain engineering of paired nucleotide-binding and leucine-rich repeat proteins in rice

Plant nucleotide-binding and leucine-rich repeat (NLR) proteins are immune sensors that detect pathogen effectors and initiate a strong immune response. In many cases, single NLR proteins are sufficient for both effector recognition and signaling activation. These proteins possess a conserved architecture, including a C-terminal leucine-rich repeat (LRR) domain, a central nucleotide-binding (NB) domain, and a variable N-terminal domain. Nevertheless, many paired NLRs linked in a head-to-head configuration have now been identified. The ones carrying integrated domains (IDs) can recognize pathogen effector proteins by various modes; these are known as sensor NLR (sNLR) proteins. Structural and biochemical studies have provided insights into the molecular basis of heavy metal-associated IDs (HMA IDs) from paired NLRs in rice and revealed the co-evolution between pathogens and hosts by combining naturally occurring favorable interactions across diverse interfaces. Focusing on structural and molecular models, here we highlight advances in structure-guided engineering to expand and enhance the response profile of paired NLR-HMA IDs in rice to variants of the rice blast pathogen MAX-effectors (Magnaporthe oryzae AVRs and ToxB-like). These results demonstrate that the HMA IDs-based design of rice materials with broad and enhanced resistance profiles possesses great application potential but also face considerable challenges.

MG1 interacts with a protease inhibitor and confers resistance to rice root-knot nematode

The rice root-knot nematode (Meloidogyne graminicola) is one of the most destructive pests threatening rice (Oryza sativa L.) production in Asia; however, no rice resistance genes have been cloned. Here, we demonstrate that M. GRAMINICOLA-RESISTANCE GENE 1 (MG1), an R gene highly expressed at the site of nematode invasion, determines resistance against the nematode in several rice varieties. Introgressing MG1 into susceptible varieties increases resistance comparable to resistant varieties, for which the leucine-rich repeat domain is critical for recognizing root-knot nematode invasion. We also report transcriptome and cytological changes that are correlated with a rapid and robust response during the incompatible interaction that occurs in resistant rice upon nematode invasion. Furthermore, we identified a putative protease inhibitor that directly interacts with MG1 during MG1-mediated resistance. Our findings provide insight into the molecular basis of nematode resistance as well as valuable resources for developing rice varieties with improved nematode resistance.

An Agrobacterium-Mediated Transient Expression Method for Functional Assay of Genes Promoting Disease in Monocots

Agrobacterium-mediated transient expression (AMTE) has been widely used for high-throughput assays of gene function in diverse plant species. However, its application in monocots is still limited due to low expression efficiency. Here, by using histochemical staining and a quantitative fluorescence assay of β-glucuronidase (GUS) gene expression, we investigated factors affecting the efficiency of AMTE on intact barley plants. We found prominent variation in GUS expression levels across diverse vectors commonly used for stable transformation and that the vector pCBEP produced the highest expression. Additionally, concurrent treatments of plants with one day of high humidity and two days of darkness following agro-infiltration also significantly increased GUS expression efficiency. We thus established an optimized method for efficient AMTE on barley and further demonstrated its efficiency on wheat and rice plants. We showed that this approach could produce enough proteins suitable for split-luciferase assays of protein-protein interactions on barley leaves. Moreover, we incorporated the AMTE protocol into the functional dissection of a complex biological process such as plant disease. Based on our previous research, we used the pCBEP vector to construct a full-length cDNA library of genes upregulated during the early stage of rice blast disease. A subsequent screen of the library by AMTE identified 15 candidate genes (out of ~2000 clones) promoting blast disease on barley plants. Four identified genes encode chloroplast-related proteins: OsNYC3, OsNUDX21, OsMRS2-9, and OsAk2. These genes were induced during rice blast disease; however, constitutive overexpression of these genes conferred enhanced disease susceptibility to Colletotrichum higginsianum in Arabidopsis. These observations highlight the power of the optimized AMTE approach on monocots as an effective tool for facilitating functional assays of genes mediating complex processes such as plant-microbe interactions.

Analysis of a rice blast resistance gene Pita-Fuhui2663 and development of selection marker

Rice blast is a detrimental rice disease caused by the fungus Magnaporthe oryzae. Here, we identified a resistance gene from the rice cultivar Fuhui 2663 which is resistant to the rice blast isolate KJ201. Through isolated population analyses and sequencing approaches, the candidate gene was traced to chromosome 12. With the use of a map-based cloning strategy, the resistance gene was ultimately mapped to an 80-kb resistance locus region containing the Pita gene. Candidate gene prediction and cDNA sequencing indicated that the target resistance gene in Fuhui 2663 was allelic to Pita, thus being referred to as Pita-Fuhui2663 hereafter. Further analysis showed that the Fuhui 2663 protein had one amino acid change: Ala (A) residue 918 in Pita-Fuhui2663 was replaced by Ser (S) in Pita-S, leading to a significant change in the 3D structure of the Pita-S protein. CRISPR/Cas9 knockout experiments confirmed that Pita-Fuhui2663 is responsible for the resistance phenotype of Fuhui 2663. Importantly, Pita-Fuhui2663 did not affect the main agronomic traits of the variety compared to the Pita gene as verified by knockout experiments, indicative of potential applications of Pita-Fuhui2663 in broader breeding programs. Furthermore, a Pita-Fuhui2663-dCAPS molecular marker with good specificity and high efficiency was developed to facilitate rice breeding for resistance to this devastating disease.

Acquisition of a complex root microbiome reshapes the transcriptomes of rice plants

Soil microorganisms can colonize plant roots and assemble in communities engaged in symbiotic relationships with their host. Though the compositional dynamics of root-associated microbiomes have been extensively studied, the host transcriptional response to these communities is poorly understood. Here, we developed an experimental system by which rice plants grown under axenic conditions can acquire a defined endosphere microbiome. Using this setup, we performed a cross-sectional characterization of plant transcriptomes in the presence or absence of a complex microbial community. To account for compositional variation, plants were inoculated with soil-derived microbiomes harvested from three distinct agricultural sites. Soil microbiomes triggered a major shift in the transcriptional profiles of rice plants that included the downregulation of one-third to one-fourth of the families of leucine-rich repeat receptor-like kinases and nucleotide-binding leucine-rich repeat receptors expressed in roots. Though the expression of several genes was consistent across all soil sources, a large fraction of this response was differentially impacted by soil type. These results demonstrate the role of root microbiomes in sculpting the transcriptomes of host plants and highlight the potential involvement of the two main receptor families of the plant immune system in the recruitment and maintenance of an endosphere microbiome.

Development and Characterization of Triticum aestivum-Aegilops longissima 6Sl Recombinants Harboring a Novel Powdery Mildew Resistance Gene Pm6Sl

Powdery mildew of wheat is a foliar disease that is spread worldwide. Cultivation of resistant varieties is the most effective, economical, and environmentally friendly strategy to curb this disease. Powdery mildew resistance genes (Pm) are the primary resources for resistance breeding, and new Pm genes are in constant demand. Previously, we identified Aegilops longissima chromosome 6S^l#3 as a carrier of powdery mildew resistance and designated the resistance gene as Pm6Sl. Here, we reported the design of 24 markers specific to 6S^l#3 on the basis of the full-length cDNA sequences of 6S^l#3 donor Ae. longissma accession TA1910, and the development of wheat-Ae. longissima 6S^l#3 introgression stocks by ph1b-induced homoeologous recombination. Further, 6S^l#3 introgression lines were identified and characterized by integration analysis of powdery mildew responses, in situ hybridization, and molecular markers and Pm6Sl was mapped to a distal interval of 42.80 Mb between markers Ael58410 and Ael57699 in the long arm of 6S^l#3. Two resistant recombinants, R43 (T6BS.6BL-6S^l#3L) and T27 (Ti6AS.6AL-6S^l#3L-6AL), contained segments harboring Pm6Sl with less than 8% of 6S^l#3 genomic length, and two markers were diagnostic for Pm6Sl. This study broadened powdery mildew resistance gene resources for wheat improvement and provided a fundamental basis for fine mapping and cloning of Pm6Sl to further understand its molecular mechanism of disease resistance.

Comparative Association Mapping Reveals Conservation of Major Gene Resistance to White Pine Blister Rust in Southwestern White Pine (Pinus strobiformis) and Limber Pine (P. flexilis)

All native North American white pines are highly susceptible to white pine blister rust (WPBR) caused by Cronartium ribicola. Understanding genomic diversity and molecular mechanisms underlying genetic resistance to WPBR remains one of the great challenges in improvement of white pines. To compare major gene resistance (MGR) present in two species, southwestern white pine (Pinus strobiformis) Cr3 and limber pine (P. flexilis) Cr4, we performed association analyses of Cr3-controlled resistant traits using single nucleotide polymorphism (SNP) assays designed with Cr4-linked polymorphic genes. We found that ∼70% of P. flexilis SNPs were transferable to P. strobiformis. Furthermore, several Cr4-linked SNPs were significantly associated with the Cr3-controlled traits in P. strobiformis families. The most significantly associated SNP (M326511_1126R) almost colocalized with Cr4 on the Pinus consensus linkage group 8, suggesting that Cr3 and Cr4 might be the same R locus, or have localizations very close to each other in the syntenic region of the P. strobiformis and P. flexilis genomes. M326511_1126R was identified as a nonsynonymous SNP, causing amino acid change (Val₃₇₆Ile) in a putative pectin acetylesterase, with coding sequences identical between the two species. Moreover, top Cr3-associated SNPs were further developed as TaqMan genotyping assays, suggesting their usefulness as marker-assisted selection (MAS) tools to distinguish genotypes between quantitative resistance and MGR. This work demonstrates the successful transferability of SNP markers between two closely related white pine species in the hybrid zone, and the possibility for deployment of MAS tools to facilitate long-term WPBR management in P. strobiformis breeding and conservation.

Long Non-Coding RNAs as Emerging Regulators of Pathogen Response in Plants

Long non-coding RNAs (lncRNAs) are transcripts without protein-coding potential that contain more than 200 nucleotides that play important roles in plant survival in response to different stresses. They interact with molecules such as DNA, RNA, and protein, and play roles in the regulation of chromatin remodeling, RNA metabolism, and protein modification activities. These lncRNAs regulate the expression of their downstream targets through epigenetic changes, at the level of transcription and post-transcription. Emerging information from computational biology and functional characterization of some of them has revealed their diverse mechanisms of action and possible roles in biological processes such as flowering time, reproductive organ development, as well as biotic and abiotic stress responses. In this review, we have mainly focused on the role of lncRNAs in biotic stress response due to the limited availability of knowledge in this domain. We have discussed the available molecular mechanisms of certain known lncRNAs against specific pathogens. Further, considering that fungal, viral, and bacterial diseases are major factors in the global food crisis, we have highlighted the importance of lncRNAs against pathogen responses and the progress in plant research to develop a better understanding of their functions and molecular mechanisms.

Comparative Genomics and Gene Pool Analysis Reveal the Decrease of Genome Diversity and Gene Number in Rice Blast Fungi by Stable Adaption with Rice

Magnaporthe oryzae caused huge losses in rice and wheat production worldwide. Comparing to long-term co-evolution history with rice, wheat-infecting isolates were new-emerging. To reveal the genetic differences between rice and wheat blast on global genomic scale, 109 whole-genome sequences of M. oryzae from rice, wheat, and other hosts were reanalyzed in this study. We found that the rice lineage had gone through stronger selective sweep and fewer conserved genes than those of Triticum and Lolium lineages, which indicated that rice blast fungi adapted to rice by gene loss and rapid evolution of specific loci. Furthermore, 228 genes associated with host adaptation of M. oryzae were found by presence/absence variation (PAV) analyses. The functional annotation of these genes found that the fine turning of genes gain/loss involved with transport and transcription factor, thiol metabolism, and nucleotide metabolism respectively are major mechanisms for rice adaption. This result implies that genetic base of specific host plant may lead to gene gain/loss variation of pathogens, so as to enhance their adaptability to host. Further characterization of these specific loci and their roles in adaption and evaluation of the fungi may eventually lead to understanding of interaction mechanism and develop new strategies of the disease management.

Genome-wide Identification and Evolutionary Analysis of NBS-LRR Genes From Secale cereale

Secale cereale is an important crop in the Triticeae tribe of the Poaceae family, and it has unique agronomic characteristics and genome properties. It possesses resistance to many diseases and serves as an important resource for the breeding of other Triticeae crops. We performed a genome-wide study on S. cereale to identify the largest group of plant disease resistance genes (R genes), the nucleotide-binding site-leucine-rich repeat receptor (NBS-LRR) genes. In its genome, 582 NBS-LRR genes were identified, including one from the RNL subclass and 581 from the CNL subclass. The NBS-LRR gene number in the S. cereale genome is greater than that in barley and the diploid wheat genomes. S. cereale chromosome 4 contains the largest number of NBS-LRR genes among the seven chromosomes, which is different from the pattern in barley and the genomes B and D of wheat but similar to that in the genome A of wheat. Further synteny analysis suggests that more NBS-LRR genes on chromosome 4 have been inherited from a common ancestor by S. cereale and the wheat genome A than the wheat genomes B and D. Phylogenetic analysis revealed that at least 740 NBS-LRR lineages are present in the common ancestor of S. cereale, Hordeum vulgare and Triticum urartu. However, most of them have only been inherited by one or two species, with only 65 of them preserved in all three species. The S. cereale genome inherited 382 of these ancestral NBS-LRR lineages, but 120 of them have been lost in both H. vulgare and T. urartu. This study provides the full NBS-LRR profile of the S. cereale genome, which is a resource for S. cereale breeding and indicates that S. cereale can be an important material for the molecular breeding of other Triticeae crops.

Identification and Cloning of a CC-NBS-NBS-LRR Gene as a Candidate of Pm40 by Integrated Analysis of Both the Available Transcriptional Data and Published Linkage Mapping

Wheat powdery mildew, caused by the obligate parasite Blumeria graminis f. sp. tritici, severely reduces wheat yields. Identifying durable and effective genes against wheat powdery mildew and further transferring them into wheat cultivars is important for finally controlling this disease in wheat production. Pm40 has been widely used in wheat breeding programs in Southwest China due to the spectrum and potentially durable resistance to powdery mildew. In the present study, a resistance test demonstrated that Pm40 is still effective against the Bgt race E20. We identified and cloned the TraesCS7B01G164000 with a total length of 4883 bp, including three exons and two introns, and encoded a protein carrying the CC-NBS-NBS-LRR domain in the Pm40-linked region flanked by two EST markers, BF478514 and BF291338, by integrating analysis of gene annotation in wheat reference genome and both sequence and expression difference in available transcriptome data. Two missense mutations were detected at positions 68 and 83 in the CC domain. The results of both cosegregation linkage analysis and qRT-PCR also suggested that TraesCS7B01G164000 was a potential candidate gene of Pm40. This study allowed us to move toward the final successfully clone and apply Pm40 in wheat resistance improvement by gene engineering.

Characterization and Evaluation of Transgenic Rice Pyramided with the Pi Genes Pib, Pi25 and Pi54

BACKGROUND

Emergence of new pathogen strains of Magnaporthe oryzae is a major reason for recurrent failure of the resistance mediated by a single resistance gene (Pi) in rice. Stacking various Pi genes in the genome through marker-assisted selection is thus an effective strategy in rice breeding for achieving durable resistance against the pathogen. However, the effect of pyramiding of multiple Pi genes using transgenesis still remains largely unknown.

RESULTS

Three Pi genes Pib, Pi25 and Pi54 were transferred together into two rice varieties, the indica variety Kasalath and the japonica variety Zhenghan 10. Transgenic plants of both Kasalath and Zhenghan 10 expressing the Pi transgenes showed imparted pathogen resistance. All the transgenic lines of both cultivars also exhibited shorter growth periods with flowering 2-4 days early, and shorter plant heights with smaller panicle. Thus, pyramiding of the Pi genes resulted in reduced grain yields in both rice cultivars. However, tiller numbers and grain weight were generally similar between the pyramided lines and corresponding parents. A global analysis of gene expression by RNA-Seq suggested that both enhancement and, to a lesser extent, inhibition of gene transcription occurred in the pyramided plants. A total of 264 and 544 differentially expressed genes (DEGs) were identified in Kasalath and Zhenghan 10, respectively. Analysis of the DEGs suggested that presence of the Pi transgenes did not alter gene expression only related to disease resistance, but also impacted many gene transcriptions in the pathways for plant growth and development, in which several were common for both Kasalath and Zhenghan 10.

CONCLUSION

Pyramiding of the Pi genes Pib, Pi25 and Pi54 via transgenesis is a potentially promising approach for improving rice resistance to the pathogen Magnaporthe oryzae. However, pleiotropic effects of the Pi genes could potentially result in yield loss. These findings support the idea that immunity is often associated with yield penalties. Rational combination of the Pi genes based on the genetic background may be important to balance yield and disease resistance.

A playbook for developing disease-resistant crops through immune receptor identification and transfer

Plants are resistant to most pathogens because of an immune system that perceives invading microbes and activates defense. A large repertoire of innate immune receptors mediates specific direct or indirect recognition of pathogen-derived molecules. Disease is often a consequence of insufficient immune surveillance, and the transfer of immune receptor genes from resistant plants to susceptible crop varieties is an effective strategy for combating disease outbreaks. We discuss approaches for identifying intracellular and cell surface immune receptors, with particular focus on recently developed and emerging methodologies. We also review considerations for the transfer of immune receptor genes into crop species, including additional host factors that may be required for immune receptor function. Together, these concepts lay out a broadly applicable playbook for developing crop varieties with durable disease resistance.

Fine dissection of limber pine resistance to Cronartium ribicola using targeted sequencing of the NLR family

BACKGROUND

Proteins with nucleotide binding site (NBS) and leucine-rich repeat (LRR) domains (NLR) make up one of most important resistance (R) families for plants to resist attacks from various pathogens and pests. The available transcriptomes of limber pine (Pinus flexilis) allow us to characterize NLR genes and related resistance gene analogs (RGAs) in host resistance against Cronartium ribicola, the causal fungal pathogen of white pine blister rust (WPBR) on five-needle pines throughout the world. We previously mapped a limber pine major gene locus (Cr4) that confers complete resistance to C. ribicola on the Pinus consensus linkage group 8 (LG-8). However, genetic distribution of NLR genes as well as their divergence between resistant and susceptible alleles are still unknown.

RESULTS

To identify NLR genes at the Cr4 locus, the present study re-sequenced a total of 480 RGAs using targeted sequencing in a Cr4-segregated seed family. Following a call of single nucleotide polymorphisms (SNPs) and genetic mapping, a total of 541 SNPs from 155 genes were mapped across 12 LGs. Three putative NLR genes were newly mapped in the Cr4 region, including one that co-segregated with Cr4. The tight linkage of NLRs with Cr4-controlled phenotypes was further confirmed by bulked segregation analysis (BSA) using extreme-phenotype genome-wide association study (XP-GWAS) for significance test. Local tandem duplication in the Cr4 region was further supported by syntenic analysis using the sugar pine genome sequence. Significant gene divergences have been observed in the NLR family, revealing that diversifying selection pressures are relatively higher in local duplicated genes. Most genes showed similar expression patterns at low levels, but some were affected by genetic background related to disease resistance. Evidence from fine genetic dissection, evolutionary analysis, and expression profiling suggests that two NLR genes are the most promising candidates for Cr4 against WPBR.

CONCLUSION

This study provides fundamental insights into genetic architecture of the Cr4 locus as well as a set of NLR variants for marker-assisted selection in limber pine breeding. Novel NLR genes were identified at the Cr4 locus and the Cr4 candidates will aid deployment of this R gene in combination with other major/minor genes in the limber pine breeding program.

Molecular Breeding Strategy and Challenges Towards Improvement of Downy Mildew Resistance in Cauliflower (Brassica oleracea var. botrytis L.)

Cauliflower (Brassica oleracea var. botrytis L.) is one of the important, nutritious and healthy vegetable crops grown and consumed worldwide. But its production is constrained by several destructive fungal diseases and most importantly, downy mildew leading to severe yield and quality losses. For sustainable cauliflower production, developing resistant varieties/hybrids with durable resistance against broad-spectrum of pathogens is the best strategy for a long term and reliable solution. Identification of novel resistant resources, knowledge of the genetics of resistance, mapping and cloning of resistance QTLs and identification of candidate genes would facilitate molecular breeding for disease resistance in cauliflower. Advent of next-generation sequencing technologies (NGS) and publishing of draft genome sequence of cauliflower has opened the flood gate for new possibilities to develop enormous amount of genomic resources leading to mapping and cloning of resistance QTLs. In cauliflower, several molecular breeding approaches such as QTL mapping, marker-assisted backcrossing, gene pyramiding have been carried out to develop new resistant cultivars. Marker-assisted selection (MAS) would be beneficial in improving the precision in the selection of improved cultivars against multiple pathogens. This comprehensive review emphasizes the fascinating recent advances made in the application of molecular breeding approach for resistance against an important pathogen; Downy Mildew (Hyaloperonospora parasitica) affecting cauliflower and Brassica oleracea crops and highlights the QTLs identified imparting resistance against this pathogen. We have also emphasized the critical research areas as future perspectives to bridge the gap between availability of genomic resources and its utility in identifying resistance genes/QTLs to breed downy mildew resistant cultivars. Additionally, we have also discussed the challenges and the way forward to realize the full potential of molecular breeding for downy mildew resistance by integrating marker technology with conventional breeding in the post-genomics era. All this information will undoubtedly provide new insights to the researchers in formulating future breeding strategies in cauliflower to develop durable resistant cultivars against the major pathogens in general and downy mildew in particular.

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.

Analysis of intraspecies diversity reveals a subset of highly variable plant immune receptors and predicts their binding sites

The evolution of recognition specificities by the immune system depends on the generation of receptor diversity and on connecting the binding of new antigens with the initiation of downstream signaling. In plant immunity, the innate Nucleotide-Binding Leucine-Rich Repeat (NLR) receptor family enables antigen binding and immune signaling. In this study, we surveyed the NLR complements of 62 ecotypes of Arabidopsis thaliana and 54 lines of Brachypodium distachyon and identified a limited number of NLR subfamilies that show high allelic diversity. We show that the predicted specificity-determining residues cluster on the surfaces of Leucine-Rich Repeat domains, but the locations of the clusters vary among NLR subfamilies. By comparing NLR phylogeny, allelic diversity, and known functions of the Arabidopsis NLRs, we formulate a hypothesis for the emergence of direct and indirect pathogen-sensing receptors and of the autoimmune NLRs. These findings reveal the recurring patterns of evolution of innate immunity and can inform NLR engineering efforts.

Genome-Wide Analysis of NLR Disease Resistance Genes in an Updated Reference Genome of Barley

Barley is one of the top 10 crop plants in the world. During its whole lifespan, barley is frequently infected by various pathogens. In this study, we performed genome-wide analysis of the largest group of plant disease resistance (R) genes, the nucleotide binding site-leucine-rich repeat receptor (NLR) gene, in an updated barley genome. A total of 468 NLR genes were identified from the improved barley genome, including one RNL subclass and 467 CNL subclass genes. Proteins of 43 barley CNL genes were shown to contain 25 different integrated domains, including WRKY and BED. The NLR gene number identified in this study is much larger than previously reported results in earlier versions of barley genomes, and only slightly fewer than that in the diploid wheat Triticum urartu. Barley Chromosome 7 contains the largest number of 112 NLR genes, which equals to seven times of the number of NLR genes on Chromosome 4. The majority of NLR genes (68%) are located in multigene clusters. Phylogenetic analysis revealed that at least 18 ancestral CNL lineages were presented in the common ancestor of barley, T. urartu and Arabidopsis thaliana. Among them fifteen lineages expanded to 533 sub-lineages prior to the divergence of barley and T. urartu. The barley genome inherited 356 of these sub-lineages and duplicated to the 467 CNL genes detected in this study. Overall, our study provides an updated profile of barley NLR genes, which should serve as a fundamental resource for functional gene mining and molecular breeding of barley.

Evidence of Constrained Divergence and Conservatism in Climatic Niches of the Temperate Maples (Acer L.)

Research highlights: The availability of global distribution data and new, fossil-calibrated phylogenies has made it possible to compare the climatic niches of the temperate maple (Acer L.) taxa and assess phylogenetic and continental patterns in niche overlap. Background and Objectives: The maples have radiated from East Asia into two other temperate continental bioregions, North America and Eurasia (Europe and West Asia), over a roughly 60-million-year period. During this time, the Earth’s climate experienced pronounced cooling and drying, culminating in cyclic periods of widespread temperate glaciation in the Pliocene to Pleistocene. The objective of this study is to use newly available data to model the climatic niches of 60% of the temperate maples and assess patterns of niche divergence, constraint, and conservatism in the genus’s radiation out of East Asia. Materials and Methods: I assembled global occurrence data and associated climatic information for 71 maple taxa, including all species endemic to temperate North America and Eurasia and their closely related East Asian congeners. I constructed Maxent niche models for all taxa and compared the climatic niches of 184 taxa pairs and assessed phylogenetic signal in key niche axes for each taxon and in niche overlap at the continental and global scale. Results: Maxent models define a fundamental climatic niche for temperate maples and suggest that drought-intolerant taxa have been lost from the Eurasian maple flora, with little continental difference in temperature optima or breadth. Niche axes and niche overlap show minimal evidence of phylogenetic signal, suggesting adaptive evolution. Pairwise niche comparisons reveal infrequent niche overlap continentally and globally, even among sister pairs, with few taxa pairs sharing ecological niche space, providing evidence for constrained divergence within the genus’s fundamental climatic niche. Evidence of niche conservatism is limited to three somewhat geographically isolated regions of high maple diversity (western North America, the Caucasus, and Japan). Conclusions: Over 60 million years of hemispheric radiation on a cooling and drying planet, the maple genus experienced divergent, though constrained, climatic niche evolution. High climatic niche diversity across spatial and phylogenetic scales along with very limited niche overlap or conservatism suggests that the radiation of the genus has largely been one of adaptive diversification.

A New Catalog of Structural Variants in 1,301 A. thaliana Lines from Africa, Eurasia, and North America Reveals a Signature of Balancing Selection at Defense Response Genes

Genomic variation in the model plant Arabidopsis thaliana has been extensively used to understand evolutionary processes in natural populations, mainly focusing on single-nucleotide polymorphisms. Conversely, structural variation has been largely ignored in spite of its potential to dramatically affect phenotype. Here, we identify 155,440 indels and structural variants ranging in size from 1 bp to 10 kb, including presence/absence variants (PAVs), inversions, and tandem duplications in 1,301 A. thaliana natural accessions from Morocco, Madeira, Europe, Asia, and North America. We show evidence for strong purifying selection on PAVs in genes, in particular for housekeeping genes and homeobox genes, and we find that PAVs are concentrated in defense-related genes (R-genes, secondary metabolites) and F-box genes. This implies the presence of a "core" genome underlying basic cellular processes and a "flexible" genome that includes genes that may be important in spatially or temporally varying selection. Further, we find an excess of intermediate frequency PAVs in defense response genes in nearly all populations studied, consistent with a history of balancing selection on this class of genes. Finally, we find that PAVs in genes involved in the cold requirement for flowering (vernalization) and drought response are strongly associated with temperature at the sites of origin.

Local Duplication of TIR-NBS-LRR Gene Marks Clubroot Resistance in Brassica napus cv. Tosca

Clubroot, caused by Plasmodiophora brassicae infection, is a disease of growing importance in cruciferous crops, including oilseed rape (Brassica napus). The affected plants exhibit prominent galling of the roots that impairs their capacity for water and nutrient uptake, which leads to growth retardation, wilting, premature ripening, or death. Due to the scarcity of effective means of protection against the pathogen, breeding of resistant varieties remains a crucial component of disease control measures. The key aspect of the breeding process is the identification of genetic factors associated with variable response to the pathogen exposure. Although numerous clubroot resistance loci have been described in Brassica crops, continuous updates on the sources of resistance are necessary. Many of the resistance genes are pathotype-specific, moreover, resistance breakdowns have been reported. In this study, we characterize the clubroot resistance locus in the winter oilseed rape cultivar "Tosca." In a series of greenhouse experiments, we evaluate the disease severity of P. brassicae-challenged "Tosca"-derived population of doubled haploids, which we genotype with Brassica 60 K array and a selection of SSR/SCAR markers. We then construct a genetic map and narrow down the resistance locus to the 0.4 cM fragment on the A03 chromosome, corresponding to the region previously described as Crr3. Using Oxford Nanopore long-read genome resequencing and RNA-seq we review the composition of the locus and describe a duplication of TIR-NBS-LRR gene. Further, we explore the transcriptomic differences of the local genes between the clubroot resistant and susceptible, inoculated and control DH lines. We conclude that the duplicated TNL gene is a promising candidate for the resistance factor. This study provides valuable resources for clubroot resistance breeding programs and lays a foundation for further functional studies on clubroot resistance.

Identification of long noncoding RNAs involved in resistance to downy mildew in Chinese cabbage

Brassica downy mildew, a severe disease caused by Hyaloperonospora brassicae, can cause enormous economic losses in Chinese cabbage (Brassica rapa L. ssp. pekinensis) production. Although some research has been reported recently concerning the underlying resistance to this disease, no studies have identified or characterized long noncoding RNAs involved in this defense response. In this study, using high-throughput RNA sequencing, we analyzed the disease-responding mRNAs and long noncoding RNAs in two resistant lines (T12-19 and 12-85) and one susceptible line (91-112). Clustering and Gene Ontology analysis of differentially expressed genes (DEGs) showed that more DEGs were involved in the defense response in the two resistant lines than in the susceptible line. Different expression patterns and proposed functions of differentially expressed long noncoding RNAs among T12-19, 12-85, and 91-112 indicated that each has a distinct disease response mechanism. There were significantly more cis- and trans-functional long noncoding RNAs in the resistant lines than in the susceptible line, and the genes regulated by these RNAs mostly participated in the disease defense response. Furthermore, we identified a candidate resistance-related long noncoding RNA, MSTRG.19915, which is a long noncoding natural antisense transcript of a MAPK gene, BrMAPK15. Via an agroinfiltration-mediated transient overexpression system and virus-induced gene silencing technology, BrMAPK15 was indicated to have a greater ability to defend against pathogens. MSTRG.19915-silenced seedlings showed enhanced resistance to downy mildew, probably because of the upregulated expression of BrMAPK15. This research identified and characterized long noncoding RNAs involved in resistance to downy mildew, laying a foundation for future in-depth studies of disease resistance mechanisms in Chinese cabbage.

The allelic rice immune receptor Pikh confers extended resistance to strains of the blast fungus through a single polymorphism in the effector binding interface

Arms race co-evolution drives rapid adaptive changes in pathogens and in the immune systems of their hosts. Plant intracellular NLR immune receptors detect effectors delivered by pathogens to promote susceptibility, activating an immune response that halts colonization. As a consequence, pathogen effectors evolve to escape immune recognition and are highly variable. In turn, NLR receptors are one of the most diverse protein families in plants, and this variability underpins differential recognition of effector variants. The molecular mechanisms underlying natural variation in effector recognition by NLRs are starting to be elucidated. The rice NLR pair Pik-1/Pik-2 recognizes AVR-Pik effectors from the blast fungus Magnaporthe oryzae, triggering immune responses that limit rice blast infection. Allelic variation in a heavy metal associated (HMA) domain integrated in the receptor Pik-1 confers differential binding to AVR-Pik variants, determining resistance specificity. Previous mechanistic studies uncovered how a Pik allele, Pikm, has extended recognition to effector variants through a specialized HMA/AVR-Pik binding interface. Here, we reveal the mechanistic basis of extended recognition specificity conferred by another Pik allele, Pikh. A single residue in Pikh-HMA increases binding to AVR-Pik variants, leading to an extended effector response in planta. The crystal structure of Pikh-HMA in complex with an AVR-Pik variant confirmed that Pikh and Pikm use a similar molecular mechanism to extend their pathogen recognition profile. This study shows how different NLR receptor alleles functionally converge to extend recognition specificity to pathogen effectors.

Integrated Analysis of MicroRNA and Target Genes in Brachypodium distachyon Infected by Magnaporthe oryzae by Small RNA and Degradome Sequencing

Rice blast caused by Magnaporthe oryzae is one of the most important diseases that seriously threaten rice production. Brachypodium distachyon is a grass species closely related to grain crops, such as rice, barley, and wheat, and has become a new model plant of Gramineae. In this study, 15 small RNA samples were sequenced to examine the dynamic changes in microRNA (miRNA) expression in B. distachyon infected by M. oryzae at 0, 24, and 48 h after inoculation. We identified 432 conserved miRNAs and 288 predicted candidate miRNAs in B. distachyon. Additionally, there were 7 and 19 differentially expressed miRNAs at 24 and 48 h post-inoculation, respectively. Furthermore, using degradome sequencing, we identified 2,126 genes as targets for 308 miRNAs; using quantitative real-time PCR (qRT-PCR), we validated five miRNA/target regulatory units involved in B. distachyon-M. oryzae interactions. Moreover, using co-transformation technology, we demonstrated that BdNAC21 was negatively regulated by miR164c. This study provides a new approach for identifying resistance genes in B. distachyon by mining the miRNA regulatory network of host-pathogen interactions.

Understanding the Sorghum-Colletotrichum sublineola Interactions for Enhanced Host Resistance

Improving sorghum resistance is a sustainable method to reduce yield losses due to anthracnose, a devastating disease caused by Colletotrichum sublineola. Elucidating the molecular mechanisms of sorghum-C. sublineola interactions would help identify biomarkers for rapid and efficient identification of novel sources for host-plant resistance improvement, understanding the pathogen virulence, and facilitating resistance breeding. Despite concerted efforts to identify resistance sources, the knowledge about sorghum-anthracnose interactions remains scanty. Hence, in this review, we presented an overview of the current knowledge on the mechanisms of sorghum-C. sublineola molecular interactions, sources of resistance for sorghum breeding, quantitative trait loci (QTL), and major (R-) resistance gene sequences as well as defense-related genes associated with anthracnose resistance. We summarized current knowledge about C. sublineola populations and its virulence. Illustration of the sorghum-C. sublineola interaction model based on the current understanding is also provided. We highlighted the importance of genomic resources of both organisms for integrated omics research to unravel the key molecular components underpinning compatible and incompatible sorghum-anthracnose interactions. Furthermore, sorghum-breeding strategy employing rapid sorghum germplasm screening, systems biology, and molecular tools is presented.

Contrasting genetic variation and positive selection followed the divergence of NBS-encoding genes in Asian and European pears

BACKGROUND

The NBS disease-related gene family coordinates the inherent immune system in plants in response to pathogen infections. Previous studies have identified NBS-encoding genes in Pyrus bretschneideri ('Dangshansuli', an Asian pear) and Pyrus communis ('Bartlett', a European pear) genomes, but the patterns of genetic variation and selection pressure on these genes during pear domestication have remained unsolved.

RESULTS

In this study, 338 and 412 NBS-encoding genes were identified from Asian and European pear genomes. This difference between the two pear species was the result of proximal duplications. About 15.79% orthologous gene pairs had Ka/Ks ratio more than one, indicating two pear species undergo strong positive selection after the divergence of Asian and European pear. We identified 21 and 15 NBS-encoding genes under fire blight and black spot disease-related QTL, respectively, suggesting their importance in disease resistance. Domestication caused decreased nucleotide diversity across NBS genes in Asian cultivars (cultivated 6.23E-03; wild 6.47E-03), but opposite trend (cultivated 6.48E-03; wild 5.91E-03) appeared in European pears. Many NBS-encoding coding regions showed Ka/Ks ratio of greater than 1, indicating the role of positive selection in shaping diversity of NBS-encoding genes in pear. Furthermore, we detected 295 and 122 significantly different SNPs between wild and domesticated accessions in Asian and European pear populations. Two NBS genes (Pbr025269.1 and Pbr019876.1) with significantly different SNPs showed >5x upregulation between wild and cultivated pear accessions, and > 2x upregulation in Pyrus calleryana after inoculation with Alternaria alternata. We propose that positively selected and significantly different SNPs of an NBS-encoding gene (Pbr025269.1) regulate gene expression differences in the wild and cultivated groups, which may affect resistance in pear against A. alternata.

CONCLUSION

Proximal duplication mainly led to the different number of NBS-encoding genes in P. bretschneideri and P. communis genomes. The patterns of genetic diversity and positive selection pressure differed between Asian and European pear populations, most likely due to their independent domestication events. This analysis helps us understand the evolution, diversity, and selection pressure in the NBS-encoding gene family in Asian and European populations, and provides opportunities to study mechanisms of disease resistance in pear.

Gene presence-absence variation associates with quantitative Verticillium longisporum disease resistance in Brassica napus

Although copy number variation (CNV) and presence-absence variation (PAV) have been discovered in selected gene families in most crop species, the global prevalence of these polymorphisms in most complex genomes is still unclear and their influence on quantitatively inherited agronomic traits is still largely unknown. Here we analyze the association of gene PAV with resistance of oilseed rape (Brassica napus) against the important fungal pathogen Verticillium longisporum, as an example for a complex, quantitative disease resistance in the strongly rearranged genome of a recent allopolyploid crop species. Using Single Nucleotide absence Polymorphism (SNaP) markers to efficiently trace PAV in breeding populations, we significantly increased the resolution of loci influencing V. longisporum resistance in biparental and multi-parental mapping populations. Gene PAV, assayed by resequencing mapping parents, was observed in 23-51% of the genes within confidence intervals of quantitative trait loci (QTL) for V. longisporum resistance, and high-priority candidate genes identified within QTL were all affected by PAV. The results demonstrate the prominent role of gene PAV in determining agronomic traits, suggesting that this important class of polymorphism should be exploited more systematically in future plant breeding.

Botryosphaeriaceae Fungi as Causal Agents of Dieback and Corky Bark in Rambutan and Longan

Fungi in the Botryosphaeriaceae family cause dieback, fruit rots, and stem cankers in many tropical fruit trees. To identify which species of Botryosphaeriaceae were present in tropical fruit in Puerto Rico and the symptoms they cause in rambutan and longan, a disease survey was conducted throughout the island from 2008 to 2016. Diseased organs of rambutan, longan, mango, and tangerine were collected and 39 isolates belonging to the Botryosphaeriaceae family were isolated and identified. Phylogenetic analysis of three nuclear genes identified nine species: six Lasiodiplodia spp. and three Neofusicoccum spp. All 39 isolates were inoculated on healthy 1-year-old rambutan and longan seedlings to confirm their pathogenicity. Dieback on both rambutan and longan was observed at 14 days after inoculation (DAI). Fourteen isolates from seven Botryosphaeriaceae species (Lasiodiplodia brasiliensis, L. hormozganensis, L. iraniensis, L. pseudotheobromae, L. theobromae, Neofusicoccum batangarum, and N. parvum) caused dieback in rambutan. Five of these pathogenic isolates were collected from rambutan, four from longan, two from mango, and three from tangerine. Ten isolates of four Lasiodiplodia spp. caused dieback in longan: L. hormozganensis, L. iraniensis, L. pseudotheobromae and L. theobromae,. Of these, three were collected from rambutan, three from longan, one from mango, and three from tangerine. Abundant development of pycnidia on branches, called corky bark, were observed on both rambutan and longan 60 DAI. Cross-inoculations showed that pathogenicity is wide in spectrum, indicating that different planting alternatives should be considered for better crop management.

Evolutionary balance between LRR domain loss and young NBS-LRR genes production governs disease resistance in Arachis hypogaea cv. Tifrunner

BACKGROUND

Cultivated peanut (Arachis hypogaea L.) is an important oil and protein crop, but it has low disease resistance; therefore, it is important to reveal the number, sequence features, function, and evolution of genes that confer resistance. Nucleotide-binding site-leucine-rich repeats (NBS-LRRs) are resistance genes that are involved in response to various pathogens.

RESULTS

We identified 713 full-length NBS-LRRs in A. hypogaea cv. Tifrunner. Genetic exchange events occurred on NBS-LRRs in A. hypogaea cv. Tifrunner, which were detected in the same subgenomes and also found in different subgenomes. Relaxed selection acted on NBS-LRR proteins and LRR domains in A. hypogaea cv. Tifrunner. Using quantitative trait loci (QTL), we found that NBS-LRRs were involved in response to late leaf spot, tomato spotted wilt virus, and bacterial wilt in A. duranensis (2 NBS-LRRs), A. ipaensis (39 NBS-LRRs), and A. hypogaea cv. Tifrunner (113 NBS-LRRs). In A. hypogaea cv. Tifrunner, 113 NBS-LRRs were classified as 75 young and 38 old NBS-LRRs, indicating that young NBS-LRRs were involved in response to disease after tetraploidization. However, compared to A. duranensis and A. ipaensis, fewer LRR domains were found in A. hypogaea cv. Tifrunner NBS-LRR proteins, partly explaining the lower disease resistance of the cultivated peanut.

CONCLUSIONS

Although relaxed selection acted on NBS-LRR proteins and LRR domains, LRR domains were preferentially lost in A. hypogaea cv. Tifrunner compared to A. duranensis and A. ipaensis. The QTL results suggested that young NBS-LRRs were important for resistance against diseases in A. hypogaea cv. Tifrunner. Our results provid insight into the greater susceptibility of A. hypogaea cv. Tifrunner to disease compared to A. duranensis and A. ipaensis.

Large-scale identification and functional analysis of NLR genes in blast resistance in the Tetep rice genome sequence

Tetep is a rice cultivar known for broad-spectrum resistance to blast, a devastating fungal disease. The molecular basis for its broad-spectrum resistance is still poorly understood. Is it because Tetep has many more NLR genes than other cultivars? Or does Tetep possess multiple major NLR genes that can individually confer broad-spectrum resistance to blast? Moreover, are there many interacting NLR pairs in the Tetep genome? We sequenced its genome, obtained a high-quality assembly, and annotated 455 nucleotide-binding site leucine-rich repeat (NLR) genes. We cloned and tested 219 NLR genes as transgenes in 2 susceptible cultivars using 5 to 12 diversified pathogen strains; in many cases, fewer than 12 strains were successfully cultured for testing. Ninety cloned NLRs showed resistance to 1 or more pathogen strains and each strain was recognized by multiple NLRs. However, few NLRs showed resistance to >6 strains, so multiple NLRs are apparently required for Tetep's broad-spectrum resistance to blast. This was further supported by the pedigree analyses, which suggested a correlation between resistance and the number of Tetep-derived NLRs. In developing a method to identify NLR pairs each of which functions as a unit, we found that >20% of the NLRs in the Tetep and 3 other rice genomes are paired. Finally, we designed an extensive set of molecular markers for rapidly introducing clustered and paired NLRs in the Tetep genome for breeding new resistant cultivars. This study increased our understanding of the genetic basis of broad-spectrum blast resistance in rice.

Comparative Genomics Analysis in Grass Species Reveals Two Distinct Evolutionary Strategies Adopted by R Genes

Resistance genes play an important role in the defense of plants against the invasion of pathogens. In Setaria italica and closely related grass species, R genes have been identified through genetic mapping and genome-wide homologous/domain searching. However, there has been to date no systematic analysis of the evolutionary features of R genes across all sequenced grass genomes. Here, we determined and comprehensively compared R genes in all 12 assembled grass genomes and an outgroup species (Arabidopsis thaliana) through synteny and selection analyses of multiple genomes. We found that the two groups of nucleotide binding site (NBS) domains containing R genes—R tandem duplications (TD) and R singletons—adopted different strategies and showed different features in their evolution. Based on K_a/K_s analysis between syntenic R loci pairs of TDs or singletons, we conclude that R singletons are under stronger purifying selection to be conserved among different grass species than R TDs, while R genes located at TD arrays have evolved much faster through diversifying selection. Furthermore, using the variome datasets of S. italica populations, we scanned for selection signals on genes and observed that a part of R singleton genes have been under purifying selection in populations of S. italica, which is consistent with the pattern observed in syntenic R singletons among different grass species. Additionally, we checked the synteny relationships of reported R genes in grass species and found that the functionally mapped R genes for novel resistance traits are prone to appear in TDs and are heavily divergent from their syntenic orthologs in other grass species, such the black streak R gene Rxo1 in Z. mays and the blast R gene Pi37 in O. sativa. These findings indicate that the R genes from TDs adopted tandem duplications to evolve faster and accumulate more mutations to facilitate functional innovation to cope with variable threats from a fluctuating environment, while R singletons provide a way for R genes to maintain sequence stability and retain conservation of function.

Origin and evolution of the plant immune system

Contents Summary 70 I. Introduction 70 II. Ancient associations between plants and microbes 72 III. Evolutionary dynamics of plant-pathogen interactions 74 IV. Evolutionary signature of plant-pathogen interactions 74 V. Origin and evolution of RLK proteins 75 VI. Origin and evolution of NLR proteins 77 VII. Origin and evolution of SA signaling 78 VIII. Origin and evolution of RNA-based defense 79 IX. Perspectives 79 Acknowledgements 80 References 80 SUMMARY: Microbes have engaged in antagonistic associations with plants for hundreds of millions of years. Plants, in turn, have evolved diverse immune strategies to combat microbial pathogens. The conflicts between plants and pathogens result in everchanging coevolutionary cycles known as 'Red Queen' dynamics. These ancient and ongoing plant-pathogen interactions have shaped the evolution of both plant and pathogen genomes. With the recent explosion of plant genome-scale data, comparative analyses provide novel insights into the coevolutionary dynamics of plants and pathogens. Here, we discuss the ancient associations between plants and microbes as well as the evolutionary principles underlying plant-pathogen interactions. We synthesize and review the current knowledge on the origin and evolution of key components of the plant immune system. We also highlight the importance of studying algae and nonflowering land plants in understanding the evolution of the plant immune system.

Connecting genome structural variation with complex traits in crop plants

Structural genome variation is a major determinant of useful trait diversity. We describe how genome analysis methods are enabling discovery of trait-associated structural variants and their potential impact on breeding. As our understanding of complex crop genomes continues to grow, there is growing evidence that structural genome variation plays a major role in determining traits important for breeding and agriculture. Identifying the extent and impact of structural variants in crop genomes is becoming increasingly feasible with ongoing advances in the sophistication of genome sequencing technologies, particularly as it becomes easier to generate accurate long sequence reads on a genome-wide scale. In this article, we discuss the origins of structural genome variation in crops from ancient and recent genome duplication and polyploidization events and review high-throughput methods to assay such variants in crop populations in order to find associations with phenotypic traits. There is increasing evidence from such studies that gene presence-absence and copy number variation resulting from segmental chromosome exchanges may be at the heart of adaptive variation of crops to counter abiotic and biotic stress factors. We present examples from major crops that demonstrate the potential of pangenomic diversity as a key resource for future plant breeding for resilience and sustainability.

NBS-Encoding Genes in Brassica napus Evolved Rapidly After Allopolyploidization and Co-localize With Known Disease Resistance Loci

Genes containing nucleotide-binding sites (NBS) play an important role in pathogen resistance in plants. However, the evolutionary fate of NBS-encoding genes after formation of allotetraploid Brassica napus (A_nA_nC_nC_n, 2n = 38) is still unknown. We performed a genome-wide comparison of putatively functional NBS-encoding genes in B. napus and its progenitor species Brassica rapa (A_rA_r, 2n = 20) and Brassica oleracea (C_oC_o, 2n = 18), identifying 464, 202, and 146 putatively functional NBS-encoding genes respectively, with genes unevenly distributed in several clusters. The A_n-subgenome of B. napus possessed similar numbers of NBS-encoding genes (191 genes) to the A_r genome of B. rapa (202 genes) and similar clustering patterns. However, the C_n genome of B. napus had many more genes (273) than the B. oleracea C_o genome (146), with different clustering trends. Only 97 NBS-encoding genes (66.4%) in B. oleracea were homologous with NBS-encoding genes in B. napus, while 176 NBS-encoding genes (87.1%) were homologous between B. rapa and B. napus. These results suggest a greater diversification of NBS-encoding genes in the C genome may have occurred after formation of B. napus. Although most NBS-encoding genes in B. napus appeared to derive from the progenitors, the birth and death of several NBS-encoding genes was also putatively mediated by non-homologous recombination. The Ka/Ks values of most homologous pairs between B. napus and the progenitor species were less than 1, suggesting purifying selection during B. napus evolution. The majority of NBS-encoding genes (60% in all species) showed higher expression levels in root tissue (out of root, leaf, stem, seed and flower tissue types). Comparative analysis of NBS-encoding genes with mapped resistance QTL against three major diseases of B. napus (blackleg, clubroot and Sclerotinia stem rot) found 204 NBS-encoding genes in B. napus located within 71 resistance QTL intervals. The majority of NBS-encoding genes were co-located with resistance QTLs against a single disease, while 47 genes were co-located with QTLs against two diseases and 3 genes were co-located with QTLs against all three. Our results revealed significant variation as well as interesting evolutionary trajectories of NBS-encoding genes in the different Brassica subgenomes, while co-localization of NBS-encoding genes and resistance QTL may facilitate resistance breeding in oilseed rape.

Genome-wide identification of LRR-containing sequences and the response of these sequences to nematode infection in Arachis duranensis

BACKGROUND

Leucine-rich repeat (LRR)-containing genes are involved in responses to various diseases. Recently, RNA-seq data from A. duranensis after nematode (Meloidogyne arenaria) infection were released. However, the number of LRR-containing genes present in A. duranensis and the response of LRR-containing genes to nematode infection are poorly understood.

RESULTS

In this study, we found 509 amino acid sequences containing nine types of LRR domains in A. duranensis. The inferred phylogenetic relationships revealed that the nine types of LRR domains had two originations. The inferred selective pressure was mainly consistent with LRR domains undergoing purifying selection. Twenty-one LRR-containing genes were associated with possible resistance to nematode infection after 3, 6, and 9 days. Among them, Aradu.T5WNW, Aradu.JM17V, and Aradu.MKP1A were up-regulate at these three time points, while Aradu.QD5DS and Aradu.M0ENQ were up-regulated 6 and 9 days after nematode infection. The expression of the above mentioned five genes was significantly and negatively correlated with the number of LRR8 domain, indicating that fewer LRR8 domains are associated with the promotion of LRR-containing genes that resist nematode infection. Patterns of co-expression and cis-acting elements indicated that WRKY possibly regulate the responses of LRR-containing genes to nematode infection and that expansin genes may work together with LRR-containing genes in response to nematode infection.

CONCLUSIONS

We identified the number and type of LRR-containing genes in A. duranensis. The LRR-containing genes that were found appear to be involved in responses to nematode infection. The number of LRR8 domains was negatively correlated with expression after nematode infection. The WRKY transcription factor may regulate resistance to nematode infection based on LRR-containing genes. Our results could improve the understanding of resistance to nematodes and molecular breeding in peanuts.

Expression of a maize NBS gene ZmNBS42 enhances disease resistance in Arabidopsis

Expression of the ZmNBS42 in Arabidopsis plants conferred resistance to bacterial pathogens, providing potential resistance enhancement of maize in further genetic breeding. Nucleotide-binding site (NBS) domain proteins play critical roles in disease resistance. In this study, we isolate a novel NBS gene ZmNBS42 from maize and systematically investigate its function on disease resistance. We find that the expression levels of ZmNBS42 in maize leaf were strikingly increased in response to Bipolaris maydis inoculation and SA treatment. The spatial expression pattern analysis reveals that, during development, ZmNBS42 is ubiquitously highly expressed in maize root, leaf, stem, internode and seed, but lowly expressed in pericarp and embryo. To better understand the roles of ZmNBS42, we overexpressed ZmNBS42 in heterologous systems. Transient overexpression of ZmNBS42 in the leaves of Nicotiana benthamiana induces a hypersensitive response. ZmNBS42 overexpression (ZmNBS42-OE) Arabidopsis plants produced more SA content than Col-0 plants, and increased the expression levels of some defense-responsive genes compared to Col-0 plants. Moreover, the ZmNBS42-OE Arabidopsis plants displayed enhanced resistance against Pseudomonas syringae pathovar tomato DC3000 (Pst DC3000). These results together suggest that ZmNBS42 can serve as an important regulator in disease resistance, thus better understanding of ZmNBS42 would benefit the resistance enhancement in maize breeding programs.

Genetic signatures of plant resistance genes with known function within and between species

Plant disease resistance (R) genes have undergone significant evolutionary divergence to cope with rapid changes in pathogens. These highly variable evolutionary patterns may have contributed to diversity in R gene protein families or structures. Here, the evolutionary patterns of 76 identified R genes and their homologs were investigated within and between plant species. Results demonstrated that nucleotide binding sites and leucine-rich-repeat genes located in loci with complex evolutionary histories tended to evolve rapidly, have high variation in copy numbers, exhibit high levels of nucleotide variation and frequent gene conversion events, and also exhibit high non-synonymous to synonymous substitution ratios in LRR regions. However, non-NBS-LRR R genes are relatively well conserved with constrained variation and are more likely to participate in the basic defense system of hosts. In addition, both conserved and highly divergent evolutionary patterns were observed for the same R genes and were consistent with inter- and intra-specific distributions of some R genes. These results thus indicate either continuous or altered evolutionary patterns between and within species. The present investigation is the first attempt to investigate evolutionary patterns among all clearly functional R genes. The results reported here thus provide a foundation for future plant disease studies.

Comparative root transcriptome of wild Arachis reveals NBS-LRR genes related to nematode resistance

BACKGROUND

The Root-Knot Nematode (RKN), Meloidogyne arenaria, significantly reduces peanut grain quality and yield worldwide. Whilst the cultivated species has low levels of resistance to RKN and other pests and diseases, peanut wild relatives (Arachis spp.) show rich genetic diversity and harbor high levels of resistance to many pathogens and environmental constraints. Comparative transcriptome analysis can be applied to identify candidate resistance genes.

RESULTS

Transcriptome analysis during the early stages of RKN infection of two peanut wild relatives, the highly RKN resistant Arachis stenosperma and the moderately susceptible A. duranensis, revealed genes related to plant immunity with contrasting expression profiles. These included genes involved in hormone signaling and secondary metabolites production and also members of the NBS-LRR class of plant disease resistance (R) genes. From 345 NBS-LRRs identified in A.duranensis reference genome, 52 were differentially expressed between inoculated and control samples, with the majority occurring in physical clusters unevenly distributed on eight chromosomes with preferential tandem duplication. The majority of these NBS-LRR genes showed contrasting expression behaviour between A. duranensis and A. stenosperma, particularly at 6 days after nematode inoculation, coinciding with the onset of the Hypersensitive Response in the resistant species. The physical clustering of some of these NBS-LRR genes correlated with their expression patterns in the contrasting genotypes. Four NBS-LRR genes exclusively expressed in A. stenosperma are located within clusters on chromosome Aradu. A09, which harbors a QTL for RKN resistance, suggesting a functional role for their physical arrangement and their potential involvement in this defense response.

CONCLUSION

The identification of functional novel R genes in wild Arachis species responsible for triggering effective defense cascades can contribute to the crop genetic improvement and enhance peanut resilience to RKN.

Identifying a large number of high-yield genes in rice by pedigree analysis, whole-genome sequencing, and CRISPR-Cas9 gene knockout

Repeated artificial selection of a complex trait facilitates the identification of genes underlying the trait, especially if multiple selected descendant lines are available. Here we developed a pedigree-based approach to identify genes underlying the Green Revolution (GR) phenotype. From a pedigree analysis, we selected 30 cultivars including the "miracle rice" IR8, a GR landmark, its ancestors and descendants, and also other related cultivars for identifying high-yield genes. Through sequencing of these genomes, we identified 28 ancestral chromosomal blocks that were maintained in all the high-yield cultivars under study. In these blocks, we identified six genes of known function, including the GR gene sd1, and 123 loci with genes of unknown function. We randomly selected 57 genes from the 123 loci for knockout or knockdown studies and found that a high proportion of these genes are essential or have phenotypic effects related to rice production. Notably, knockout lines have significant changes in plant height (P < 0.003), a key GR trait, compared with wild-type lines. Some gene knockouts or knockdowns were especially interesting. For example, knockout of Os10g0555100, a putative glucosyltransferase gene, showed both reduced growth and altered panicle architecture. In addition, we found that in some retained chromosome blocks several GR-related genes were clustered, although they have unrelated sequences, suggesting clustering of genes with similar functions. In conclusion, we have identified many high-yield genes in rice. Our method provides a powerful means to identify genes associated with a specific trait.

Long noncoding RNAs involve in resistance to Verticillium dahliae, a fungal disease in cotton

Long noncoding RNAs (lncRNAs) have several known functions in plant development, but their possible roles in responding to plant disease remain largely unresolved. In this study, we described a comprehensive disease-responding lncRNA profiles in defence against a cotton fungal disease Verticillium dahliae. We further revealed the conserved and specific characters of disease-responding process between two cotton species. Conservatively for two cotton species, we found the expression dominance of induced lncRNAs in the Dt subgenome, indicating a biased induction pattern in the co-existing subgenomes of allotetraploid cotton. Comparative analysis of lncRNA expression and their proposed functions in resistant Gossypium barbadense cv. '7124' versus susceptible Gossypium hirsutum cv. 'YZ1' revealed their distinct disease response mechanisms. Species-specific (LS) lncRNAs containing more SNPs displayed a fiercer inducing level postinfection than the species-conserved (core) lncRNAs. Gene Ontology enrichment of LS lncRNAs and core lncRNAs indicates distinct roles in the process of biotic stimulus. Further functional analysis showed that two core lncRNAs, GhlncNAT-ANX2- and GhlncNAT-RLP7-silenced seedlings, displayed an enhanced resistance towards V. dahliae and Botrytis cinerea, possibly associated with the increased expression of LOX1 and LOX2. This study represents the first characterization of lncRNAs involved in resistance to fungal disease and provides new clues to elucidate cotton disease response mechanism.

Development of Molecular Markers Linked to Powdery Mildew Resistance Gene Pm4b by Combining SNP Discovery from Transcriptome Sequencing Data with Bulked Segregant Analysis (BSR-Seq) in Wheat

Powdery mildew resistance gene Pm4b, originating from Triticum persicum, is effective against the prevalent Blumeria graminis f. sp. tritici (Bgt) isolates from certain regions of wheat production in China. The lack of tightly linked molecular markers with the target gene prevents the precise identification of Pm4b during the application of molecular marker-assisted selection (MAS). The strategy that combines the RNA-Seq technique and the bulked segregant analysis (BSR-Seq) was applied in an F2:3 mapping population (237 families) derived from a pair of isogenic lines VPM1/7∗Bainong 3217 F4 (carrying Pm4b) and Bainong 3217 to develop more closely linked molecular markers. RNA-Seq analysis of the two phenotypically contrasting RNA bulks prepared from the representative F2:3 families generated 20,745,939 and 25,867,480 high-quality read pairs, and 82.8 and 80.2% of them were uniquely mapped to the wheat whole genome draft assembly for the resistant and susceptible RNA bulks, respectively. Variant calling identified 283,866 raw single nucleotide polymorphisms (SNPs) and InDels between the two bulks. The SNPs that were closely associated with the powdery mildew resistance were concentrated on chromosome 2AL. Among the 84 variants that were potentially associated with the disease resistance trait, 46 variants were enriched in an about 25 Mb region at the distal end of chromosome arm 2AL. Four Pm4b-linked SNP markers were developed from these variants. Based on the sequences of Chinese Spring where these polymorphic SNPs were located, 98 SSR primer pairs were designed to develop distal markers flanking the Pm4b gene. Three SSR markers, Xics13, Xics43, and Xics76, were incorporated in the new genetic linkage map, which located Pm4b in a 3.0 cM genetic interval spanning a 6.7 Mb physical genomic region. This region had a collinear relationship with Brachypodium distachyon chromosome 5, rice chromosome 4, and sorghum chromosome 6. Seven genes associated with disease resistance were predicted in this collinear genomic region, which included C2 domain protein, peroxidase activity protein, protein kinases of PKc_like super family, Mlo family protein, and catalytic domain of the serine/threonine kinases (STKc_IRAK like super family). The markers developed in the present study facilitate identification of Pm4b during its MAS practice.