The large repertoire of conifer NLR resistance genes includes drought responsive and highly diversified RNLs

Evolutionary and immune-activating character analyses of NLR genes in algae suggest the ancient origin of plant intracellular immune receptors

Nucleotide-binding leucine-rich repeat (NLR) proteins are crucial intracellular immune receptors in plants, responsible for detecting invading pathogens and initiating defense responses. While previous studies on the evolution and function of NLR genes were mainly limited to land plants, the evolutionary trajectory and immune-activating character of NLR genes in algae remain less explored. In this study, genome-wide NLR gene analysis was conducted on 44 chlorophyte species across seven classes and seven charophyte species across five classes. A few but variable number of NLR genes, ranging from one to 20, were identified in five chlorophytes and three charophytes, whereas no NLR gene was identified from the remaining algal genomes. Compared with land plants, algal genomes possess fewer or usually no NLR genes, implying that the expansion of NLR genes in land plants can be attributed to their adaptation to the more complex terrestrial pathogen environments. Through phylogenetic analysis, domain composition analysis, and conserved motifs profiling of the NBS domain, we detected shared and lineage-specific features between NLR genes in algae and land plants, supporting the common origin and continuous evolution of green plant NLR genes. Immune-activation assays revealed that both TNL and RNL proteins from green algae can elicit hypersensitive responses in Nicotiana benthamiana, indicating the molecular basis for immune activation has emerged in the early evolutionary stage of different types of NLR proteins. In summary, the results from this study suggest that NLR proteins may have taken a role as intracellular immune receptors in the common ancestor of green plants.

Maritime Pine Rootstock Genotype Modulates Gene Expression Associated with Stress Tolerance in Grafted Stems

Climate change-induced hazards, such as drought, threaten forest resilience, particularly in vulnerable regions such as the Mediterranean Basin. Maritime pine (Pinus pinaster Aiton), a model species in Western Europe, plays a crucial role in the Mediterranean forest due to its genetic diversity and ecological plasticity. This study characterizes transcriptional profiles of scion and rootstock stems of four P. pinaster graft combinations grown under well-watered conditions. Our grafting scheme combined drought-sensitive and drought-tolerant genotypes for scions (GAL1056: drought-sensitive scion; and Oria6: drought-tolerant scion) and rootstocks (R1S: drought-sensitive rootstock; and R18T: drought-tolerant rootstock). Transcriptomic analysis revealed expression patterns shaped by genotype provenance and graft combination. The accumulation of differentially expressed genes (DEGs) encoding proteins, involved in defense mechanisms and pathogen recognition, was higher in drought-sensitive scion stems and also increased when grafted onto drought-sensitive rootstocks. DEGs involved in drought tolerance mechanisms were identified in drought-tolerant genotypes as well as in drought-sensitive scions grafted onto drought-tolerant rootstocks, suggesting their establishment prior to drought. These mechanisms were associated with ABA metabolism and signaling. They were also involved in the activation of the ROS-scavenging pathways, which included the regulation of flavonoid and terpenoid metabolisms. Our results reveal DEGs potentially associated with the conifer response to drought and point out differences in drought tolerance strategies. These findings suggest genetic trade-offs between pine growth and defense, which could be relevant in selecting more drought-tolerant Pinus pinaster trees.

Conifers Concentrate Large Numbers of NLR Immune Receptor Genes on One Chromosome

Nucleotide-binding domain and leucine-rich repeat (NLR) immune receptor genes form a major line of defense in plants, acting in both pathogen recognition and resistance machinery activation. NLRs are reported to form large gene clusters in limber pine (Pinus flexilis), but it is unknown how widespread this genomic architecture may be among the extant species of conifers (Pinophyta). We used comparative genomic analyses to assess patterns in the abundance, diversity, and genomic distribution of NLR genes. Chromosome-level whole genome assemblies and high-density linkage maps in the Pinaceae, Cupressaceae, Taxaceae, and other gymnosperms were scanned for NLR genes using existing and customized pipelines. The discovered genes were mapped across chromosomes and linkage groups and analyzed phylogenetically for evolutionary history. Conifer genomes are characterized by dense clusters of NLR genes, highly localized on one chromosome. These clusters are rich in TNL-encoding genes, which seem to have formed through multiple tandem duplication events. In contrast to angiosperms and nonconiferous gymnosperms, genomic clustering of NLR genes is ubiquitous in conifers. NLR-dense genomic regions are likely to influence a large part of the plant's resistance, informing our understanding of adaptation to biotic stress and the development of genetic resources through breeding.

Epigenetic stress memory in gymnosperms

Gymnosperms are long-lived, cone-bearing seed plants that include some of the most ancient extant plant species. These relict land plants have evolved to survive in habitats marked by chronic or episodic stress. Their ability to thrive in these environments is partly due to their phenotypic flexibility, and epigenetic regulation likely plays a crucial part in this plasticity. We review the current knowledge on abiotic and biotic stress memory in gymnosperms and the possible epigenetic mechanisms underlying long-term phenotypic adaptations. We also discuss recent technological improvements and new experimental possibilities that likely will advance our understanding of epigenetic regulation in these ancient and hard-to-study plants.

Comparative analysis, diversification, and functional validation of plant nucleotide-binding site domain genes

Nucleotide-binding site (NBS) domain genes are one of the superfamily of resistance genes involved in plant responses to pathogens. The current study identified 12,820 NBS-domain-containing genes across 34 species covering from mosses to monocots and dicots. These identified genes are classified into 168 classes with several novel domain architecture patterns encompassing significant diversity among plant species. Several classical (NBS, NBS-LRR, TIR-NBS, TIR-NBS-LRR, etc.) and species-specific structural patterns (TIR-NBS-TIR-Cupin_1-Cupin_1, TIR-NBS-Prenyltransf, Sugar_tr-NBS etc.) were discovered. We observed 603 orthogroups (OGs) with some core (most common orthogroups; OG0, OG1, OG2, etc.) and unique (highly specific to species; OG80, OG82, etc.) OGs with tandem duplications. The expression profiling presented the putative upregulation of OG2, OG6, and OG15 in different tissues under various biotic and abiotic stresses in susceptible and tolerant plants to cotton leaf curl disease (CLCuD). The genetic variation between susceptible (Coker 312) and tolerant (Mac7) Gossypium hirsutum accessions identified several unique variants in NBS genes of Mac7 (6583 variants) and Coker312 (5173 variants). The protein-ligand and proteins-protein interaction showed a strong interaction of some putative NBS proteins with ADP/ATP and different core proteins of the cotton leaf curl disease virus. The silencing of GaNBS (OG2) in resistant cotton through virus-induced gene silencing (VIGS) demonstrated its putative role in virus tittering. The presented study will be further helpful in understanding the plant adaptation mechanism.

Transcriptional Profiling of Early Defense Response to White Pine Blister Rust Infection in Pinus albicaulis (Whitebark Pine)

Pathogen perception generates the activation of signal transduction cascades to host defense. White pine blister rust (WPBR) is caused by Cronartium ribicola J.C. Fisch and affects a number of species of Pinus. One of the most severely affected species is Pinus albicaulis Engelm (whitebark pine). WPBR resistance in the species is a polygenic and complex trait that requires an optimized immune response. We identified early responses in 2-year-old seedlings after four days of fungal inoculation and compared the underlying transcriptomic response with that of healthy non-inoculated individuals. A de novo transcriptome assembly was constructed with 56,796 high quality-annotations derived from the needles of susceptible and resistant individuals in a resistant half-sib family. Differential expression analysis identified 599 differentially expressed transcripts, from which 375 were upregulated and 224 were downregulated in the inoculated seedlings. These included components of the initial phase of active responses to abiotic factors and stress regulators, such as those involved in the first steps of flavonoid biosynthesis. Four days after the inoculation, infected individuals showed an overexpression of chitinases, reactive oxygen species (ROS) regulation signaling, and flavonoid intermediates. Our research sheds light on the first stage of infection and emergence of disease symptoms among whitebark pine seedlings. RNA sequencing (RNA-seq) data encoding hypersensitive response, cell wall modification, oxidative regulation signaling, programmed cell death, and plant innate immunity were differentially expressed during the defense response against C. ribicola.

A genome sequence for the threatened whitebark pine

Whitebark pine (WBP, Pinus albicaulis) is a white pine of subalpine regions in the Western contiguous United States and Canada. WBP has become critically threatened throughout a significant part of its natural range due to mortality from the introduced fungal pathogen white pine blister rust (WPBR, Cronartium ribicola) and additional threats from mountain pine beetle (Dendroctonus ponderosae), wildfire, and maladaptation due to changing climate. Vast acreages of WBP have suffered nearly complete mortality. Genomic technologies can contribute to a faster, more cost-effective approach to the traditional practices of identifying disease-resistant, climate-adapted seed sources for restoration. With deep-coverage Illumina short reads of haploid megagametophyte tissue and Oxford Nanopore long reads of diploid needle tissue, followed by a hybrid, multistep assembly approach, we produced a final assembly containing 27.6 Gb of sequence in 92,740 contigs (N50 537,007 bp) and 34,716 scaffolds (N50 2.0 Gb). Approximately 87.2% (24.0 Gb) of total sequence was placed on the 12 WBP chromosomes. Annotation yielded 25,362 protein-coding genes, and over 77% of the genome was characterized as repeats. WBP has demonstrated the greatest variation in resistance to WPBR among the North American white pines. Candidate genes for quantitative resistance include disease resistance genes known as nucleotide-binding leucine-rich repeat receptors (NLRs). A combination of protein domain alignments and direct genome scanning was employed to fully describe the 3 subclasses of NLRs. Our high-quality reference sequence and annotation provide a marked improvement in NLR identification compared to previous assessments that leveraged de novo-assembled transcriptomes.

High-density genetic linkage mapping in Sitka spruce advances the integration of genomic resources in conifers

In species with large and complex genomes such as conifers, dense linkage maps are a useful resource for supporting genome assembly and laying the genomic groundwork at the structural, populational, and functional levels. However, most of the 600+ extant conifer species still lack extensive genotyping resources, which hampers the development of high-density linkage maps. In this study, we developed a linkage map relying on 21,570 single nucleotide polymorphism (SNP) markers in Sitka spruce (Picea sitchensis [Bong.] Carr.), a long-lived conifer from western North America that is widely planted for productive forestry in the British Isles. We used a single-step mapping approach to efficiently combine RAD-seq and genotyping array SNP data for 528 individuals from 2 full-sib families. As expected for spruce taxa, the saturated map contained 12 linkages groups with a total length of 2,142 cM. The positioning of 5,414 unique gene coding sequences allowed us to compare our map with that of other Pinaceae species, which provided evidence for high levels of synteny and gene order conservation in this family. We then developed an integrated map for P. sitchensis and Picea glauca based on 27,052 markers and 11,609 gene sequences. Altogether, these 2 linkage maps, the accompanying catalog of 286,159 SNPs and the genotyping chip developed, herein, open new perspectives for a variety of fundamental and more applied research objectives, such as for the improvement of spruce genome assemblies, or for marker-assisted sustainable management of genetic resources in Sitka spruce and related species.

Evolution of immunity networks across embryophytes

Land plants (embryophytes), including vascular (tracheophytes) and non-vascular plants (bryophytes), co-evolved with microorganisms since descendants of an algal ancestor colonized terrestrial habitats around 500 million years ago. To cope with microbial pathogen infections, embryophytes evolved a complex immune system for pathogen perception and activation of defenses. With the growing number of sequenced genomes and transcriptome datasets from algae, bryophytes, tracheophytes, and available plant models, comparative analyses are increasing our understanding of the evolution of molecular mechanisms underpinning immune responses in different plant lineages. In this review, recent progress on plant immunity networks is highlighted with emphasis on the identification of key components that shaped immunity against pathogens in bryophytes compared to angiosperms during plant evolution.

Long-insert sequence capture detects high copy numbers in a defence-related beta-glucosidase gene βglu-1 with large variations in white spruce but not Norway spruce

Conifers are long-lived and slow-evolving, thus requiring effective defences against their fast-evolving insect natural enemies. The copy number variation (CNV) of two key acetophenone biosynthesis genes Ugt5/Ugt5b and βglu-1 may provide a plausible mechanism underlying the constitutively variable defence in white spruce (Picea glauca) against its primary defoliator, spruce budworm. This study develops a long-insert sequence capture probe set (Picea_hung_p1.0) for quantifying copy number of βglu-1-like, Ugt5-like genes and single-copy genes on 38 Norway spruce (Picea abies) and 40 P. glauca individuals from eight and nine provenances across Europe and North America respectively. We developed local assemblies (Piabi_c1.0 and Pigla_c.1.0), full-length transcriptomes (PIAB_v1 and PIGL_v1), and gene models to characterise the diversity of βglu-1 and Ugt5 genes. We observed very large copy numbers of βglu-1, with up to 381 copies in a single P. glauca individual. We observed among-provenance CNV of βglu-1 in P. glauca but not P. abies. Ugt5b was predominantly single-copy in both species. This study generates critical hypotheses for testing the emergence and mechanism of extreme CNV, the dosage effect on phenotype, and the varying copy number of genes with the same pathway. We demonstrate new approaches to overcome experimental challenges in genomic research in conifer defences.

Assembly and annotation of the black spruce genome provide insights on spruce phylogeny and evolution of stress response

Black spruce (Picea mariana [Mill.] B.S.P.) is a dominant conifer species in the North American boreal forest that plays important ecological and economic roles. Here, we present the first genome assembly of P. mariana with a reconstructed genome size of 18.3 Gbp and NG50 scaffold length of 36.0 kbp. A total of 66,332 protein-coding sequences were predicted in silico and annotated based on sequence homology. We analyzed the evolutionary relationships between P. mariana and 5 other spruces for which complete nuclear and organelle genome sequences were available. The phylogenetic tree estimated from mitochondrial genome sequences agrees with biogeography; specifically, P. mariana was strongly supported as a sister lineage to P. glauca and 3 other taxa found in western North America, followed by the European Picea abies. We obtained mixed topologies with weaker statistical support in phylogenetic trees estimated from nuclear and chloroplast genome sequences, indicative of ancient reticulate evolution affecting these 2 genomes. Clustering of protein-coding sequences from the 6 Picea taxa and 2 Pinus species resulted in 34,776 orthogroups, 560 of which appeared to be specific to P. mariana. Analysis of these specific orthogroups and dN/dS analysis of positive selection signatures for 497 single-copy orthogroups identified gene functions mostly related to plant development and stress response. The P. mariana genome assembly and annotation provides a valuable resource for forest genetics research and applications in this broadly distributed species, especially in relation to climate adaptation.

A Genome Sequence for the Threatened Whitebark Pine

Whitebark pine (WBP, Pinus albicaulis ) is a white pine of subalpine regions in western contiguous US and Canada. WBP has become critically threatened throughout a significant part of its natural range due to mortality from the introduced fungal pathogen white pine blister rust (WPBR, Cronartium ribicola ) and additional threats from mountain pine beetle ( Dendroctonus ponderosae ), wildfire, and maladaptation due to changing climate. Vast acreages of WBP have suffered nearly complete mortality. Genomic technologies can contribute to a faster, more cost-effective approach to the traditional practices of identifying disease-resistant, climate-adapted seed sources for restoration. With deep-coverage Illumina short-reads of haploid megametophyte tissue and Oxford Nanopore long-reads of diploid needle tissue, followed by a hybrid, multistep assembly approach, we produced a final assembly containing 27.6 Gbp of sequence in 92,740 contigs (N50 537,007 bp) and 34,716 scaffolds (N50 2.0 Gbp). Approximately 87.2% (24.0 Gbp) of total sequence was placed on the twelve WBP chromosomes. Annotation yielded 25,362 protein-coding genes, and over 77% of the genome was characterized as repeats. WBP has demonstrated the greatest variation in resistance to WPBR among the North American white pines. Candidate genes for quantitative resistance include disease resistance genes known as nucleotide-binding leucine-rich-repeat receptors (NLRs). A combination of protein domain alignments and direct genome scanning was employed to fully describe the three subclasses of NLRs (TNL, CNL, RNL). Our high-quality reference sequence and annotation provide a marked improvement in NLR identification compared to previous assessments that leveraged de novo assembled transcriptomes.

Taking the lead: NLR immune receptor N-terminal domains execute plant immune responses

Nucleotide-binding domain and leucine-rich repeat (NLR) proteins are important intracellular immune receptors that activate robust plant immune responses upon detecting pathogens. Canonical NLRs consist of a conserved tripartite architecture that includes a central regulatory nucleotide-binding domain, C-terminal leucine-rich repeats, and variable N-terminal domains that directly participate in immune execution. In flowering plants, the vast majority of NLR N-terminal domains belong to the coiled-coil, Resistance to Powdery Mildew 8, or Toll/interleukin-1 receptor subfamilies, with recent structural and biochemical studies providing detailed mechanistic insights into their functions. In this insight review, we focus on the immune-related biochemistries of known plant NLR N-terminal domains and discuss the evolutionary diversity of atypical NLR domains in nonflowering plants. We further contrast these observations against the known diversity of NLR-related receptors from microbes to metazoans across the tree of life.

Multispecies comparison of host responses to Fusarium circinatum challenge in tropical pines show consistency in resistance mechanisms

Fusarium circinatum poses a threat to both commercial and natural pine forests. Large variation in host resistance exists between species, with many economically important species being susceptible. Development of resistant genotypes could be expedited and optimised by investigating the molecular mechanisms underlying host resistance and susceptibility as well as increasing the available genetic resources. RNA-seq data, from F. circinatum inoculated and mock-inoculated ca. 6-month-old shoot tissue at 3- and 7-days postinoculation, was generated for three commercially important tropical pines, Pinus oocarpa, Pinus maximinoi and Pinus greggii. De novo transcriptomes were assembled and used to investigate the NLR and PR gene content within available pine references. Host responses to F. circinatum challenge were investigated in P. oocarpa (resistant) and P. greggii (susceptible), in comparison to previously generated expression profiles from Pinus tecunumanii (resistant) and Pinus patula (susceptible). Expression results indicated crosstalk between induced salicylate, jasmonate and ethylene signalling is involved in host resistance and compromised in susceptible hosts. Additionally, higher constitutive expression of sulfur metabolism and flavonoid biosynthesis in resistant hosts suggest involvement of these metabolites in resistance.

Genomic and Epigenomic Mechanisms of the Interaction between Parasitic and Host Plants

Parasitic plants extract nutrients from the other plants to finish their life cycle and reproduce. The control of parasitic weeds is notoriously difficult due to their tight physical association and their close biological relationship to their hosts. Parasitic plants differ in their susceptible host ranges, and the host species differ in their susceptibility to parasitic plants. Current data show that adaptations of parasitic plants to various hosts are largely genetically determined. However, multiple cases of rapid adaptation in genetically homogenous parasitic weed populations to new hosts strongly suggest the involvement of epigenetic mechanisms. Recent progress in genome-wide analyses of gene expression and epigenetic features revealed many new molecular details of the parasitic plants' interactions with their host plants. The experimental data obtained in the last several years show that multiple common features have independently evolved in different lines of the parasitic plants. In this review we discuss the most interesting new details in the interaction between parasitic and host plants.

Chromosome-level assembly of the Phytophthora agathidicida genome reveals adaptation in effector gene families

Phytophthora species are notorious plant pathogens, with some causing devastating tree diseases that threaten the survival of their host species. One such example is Phytophthora agathidicida, the causal agent of kauri dieback - a root and trunk rot disease that kills the ancient, iconic and culturally significant tree species, Agathis australis (New Zealand kauri). A deeper understanding of how Phytophthora pathogens infect their hosts and cause disease is critical for the development of effective treatments. Such an understanding can be gained by interrogating pathogen genomes for effector genes, which are involved in virulence or pathogenicity. Although genome sequencing has become more affordable, the complete assembly of Phytophthora genomes has been problematic, particularly for those with a high abundance of repetitive sequences. Therefore, effector genes located in repetitive regions could be truncated or missed in a fragmented genome assembly. Using a combination of long-read PacBio sequences, chromatin conformation capture (Hi-C) and Illumina short reads, we assembled the P. agathidicida genome into ten complete chromosomes, with a genome size of 57 Mb including 34% repeats. This is the first Phytophthora genome assembled to chromosome level and it reveals a high level of syntenic conservation with the complete genome of Peronospora effusa, the only other completely assembled genome sequence of an oomycete. All P. agathidicida chromosomes have clearly defined centromeres and contain candidate effector genes such as RXLRs and CRNs, but in different proportions, reflecting the presence of gene family clusters. Candidate effector genes are predominantly found in gene-poor, repeat-rich regions of the genome, and in some cases showed a high degree of duplication. Analysis of candidate RXLR effector genes that occur in multicopy gene families indicated half of them were not expressed in planta. Candidate CRN effector gene families showed evidence of transposon-mediated recombination leading to new combinations of protein domains, both within and between chromosomes. Further analysis of this complete genome assembly will help inform new methods of disease control against P. agathidicida and other Phytophthora species, ultimately helping decipher how Phytophthora pathogens have evolved to shape their effector repertoires and how they might adapt in the future.

An Overview of PRR- and NLR-Mediated Immunities: Conserved Signaling Components across the Plant Kingdom That Communicate Both Pathways

Cell-surface-localized pattern recognition receptors (PRRs) and intracellular nucleotide-binding domain and leucine-rich repeat receptors (NLRs) are plant immune proteins that trigger an orchestrated downstream signaling in response to molecules of microbial origin or host plant origin. Historically, PRRs have been associated with pattern-triggered immunity (PTI), whereas NLRs have been involved with effector-triggered immunity (ETI). However, recent studies reveal that such binary distinction is far from being applicable to the real world. Although the perception of plant pathogens and the final mounting response are achieved by different means, central hubs involved in signaling are shared between PTI and ETI, blurring the zig-zag model of plant immunity. In this review, we not only summarize our current understanding of PRR- and NLR-mediated immunities in plants, but also highlight those signaling components that are evolutionarily conserved across the plant kingdom. Altogether, we attempt to offer an overview of how plants mediate and integrate the induction of the defense responses that comprise PTI and ETI, emphasizing the need for more evolutionary molecular plant-microbe interactions (EvoMPMI) studies that will pave the way to a better understanding of the emergence of the core molecular machinery involved in the so-called evolutionary arms race between plants and microbes.

Comparative transcriptomic responses of European and Japanese larches to infection by Phytophthora ramorum

BACKGROUND AND OBJECTIVES

Phytophthora ramorum severely affects both European larch (EL) and Japanese larch (JL) trees as indicated by high levels of mortality particularly in the UK. Field observations suggested that EL is less severely affected and so may be less susceptible to P. ramorum than JL; however, controlled inoculations have produced inconsistent or non-statistically significant differences. The present study aimed to compare RNA transcript accumulation profiles in EL and JL in response to inoculation with P. ramorum to improve our understanding of their defence responses.

METHODOLOGY

RNA-sequencing was carried out on bark tissues following the inoculation with P. ramorum of potted saplings in both EL and JL carried out under controlled environment conditions, with sampling at 1, 3, 10, and 25 days post inoculation in infected and control plants.

RESULTS

All of the inoculated trees rapidly developed lesions but no statistically significant differences were found in lesion lengths between EL and JL. RNA-Sequencing comparing control and inoculate saplings identified key differences in differentially expressed genes (DEGs) between the two larch species. European larch had rapid induction of defence genes within 24 hours of infection followed by sustained expression until 25 days after inoculation. Results in JL were more varied; upregulation was stronger but more transient and represented fewer defence pathways. Gene enrichment analyses highlighted differences in jasmonate signalling and regulation including NPR1 upregulation in EL only, and specific aspects of secondary metabolism. Some DEGs were represented by multiple responsive copies including lipoxygenase, chalcone synthase and nucleotide-binding, leucine-rich-repeat genes.

CONCLUSION

The variations between EL and JL in responsive DEGs of interest as potentially related to differences seen in the field and should be considered in the selection of trees for planting and future breeding.

Genome-Wide Identification, Characterization, and Comparative Analysis of NLR Resistance Genes in Coffea spp

The largest family of disease resistance genes in plants are nucleotide-binding site leucine-rich repeat genes (NLRs). The products of these genes are responsible for recognizing avirulence proteins (Avr) of phytopathogens and triggering specific defense responses. Identifying NLRs in plant genomes with standard gene annotation software is challenging due to their multidomain nature, sequence diversity, and clustered genomic distribution. We present the results of a genome-wide scan and comparative analysis of NLR loci in three coffee species (Coffea canephora, Coffea eugenioides and their interspecific hybrid Coffea arabica). A total of 1311 non-redundant NLR loci were identified in C. arabica, 927 in C. canephora, and 1079 in C. eugenioides, of which 809, 562, and 695 are complete loci, respectively. The NLR-Annotator tool used in this study showed extremely high sensitivities and specificities (over 99%) and increased the detection of putative NLRs in the reference coffee genomes. The NLRs loci in coffee are distributed among all chromosomes and are organized mostly in clusters. The C. arabica genome presented a smaller number of NLR loci when compared to the sum of the parental genomes (C. canephora, and C. eugenioides). There are orthologous NLRs (orthogroups) shared between coffee, tomato, potato, and reference NLRs and those that are shared only among coffee species, which provides clues about the functionality and evolutionary history of these orthogroups. Phylogenetic analysis demonstrated orthologous NLRs shared between C. arabica and the parental genomes and those that were possibly lost. The NLR family members in coffee are subdivided into two main groups: TIR-NLR (TNL) and non-TNL. The non-TNLs seem to represent a repertoire of resistance genes that are important in coffee. These results will support functional studies and contribute to a more precise use of these genes for breeding disease-resistant coffee cultivars.

Evolutionary footprint of plant immunity

There are pieces of evidence from genomic footprints and fossil records indicating that plants have co-evolved with microbes after terrestrialization for more than 407 million years. Therefore, to truly comprehend plant evolution, we need to understand the co-evolutionary process and history between plants and microbes. Recent developments in genomes and transcriptomes of a vast number of plant species as well as microbes have greatly expanded our knowledge of the evolution of the plant immune system. In this review, we summarize recent advances in the co-evolution between plants and microbes with emphasis on the plant side and point out future research needed for understanding plant-microbial co-evolution. Knowledge of the evolution and variation of the plant immune system will better equip us on designing crops with boosted performance in agricultural fields.

A new biochemistry connecting pathogen detection to induced defense in plants

Plant cell surface and intracellular immune receptors recognizing pathogen attack utilize the same defense machineries to mobilize resistance. New genetic, protein structural and biochemical information on receptor activation and signaling is transforming understanding of how their shared defense network operates. We discuss the biochemical activities of two classes of intracellular nucleotide-binding/leucine-rich repeat (NLR) receptor - one forming a Ca²⁺ channel, the other an NADase enzyme - which define engagement of enhanced disease susceptibility 1 (EDS1)-family heterodimers and cofunctioning helper NLRs (RNLs) to connect receptor systems and amplify defenses. Toll-interleukin-1 receptor (TIR) domain NLR receptors and TIR-domain proteins, with a capacity to produce NAD+-derived small molecules, require EDS1 dimers and RNLs for defense induction. The TIR-driven EDS1/RNL modules emerge as central elements in Ca²⁺ -based immunity signaling initiated by receptors outside and inside host cells.

Transcriptome Analyses of Near Isogenic Lines Reveal Putative Drought Tolerance Controlling Genes in Wheat

Drought stress, especially at the grain-filling stage, is a major constraint for wheat production. Drought tolerance is a complex trait controlled by a large array of genes and pathways. This study conducted gene expression profiling on two pairs of near-isogenic lines (NILs) for an important qDSI.4B.1 QTL conferring drought tolerance on the short arm of chromosome 4B in wheat. Analysis showed 1,614 genome-wide differentially expressed genes (DEGs) between the tolerant and susceptible isolines in both NIL pairs. Six common DEGs were found between NIL1 and NIL2 at both 7 and 14 days after stress induction, with two of them having single nucleotide polymorphism (SNP) variants. These six genes that were confirmed by quantitative real-time PCR (qRT-PCR) expression analysis are considered candidate genes for drought tolerance mediated by qDSI.4B.1 QTL with their main contributions to gene regulation, cell elongation, protein quality control, secondary metabolism, and hormone signaling. These six candidate genes and the highest number of DEGs and variants (SNPs/indels) were located between 49 and 137 Mbp of 4BS, making this interval the most probable location for the qDSI.4B.1 locus. Additionally, 765 and 84 DEGs were detected as responsive genes to drought stress in tolerant and susceptible isolines, respectively. According to gene ontology (GO), protein phosphorylation, oxidation reduction, and regulation of transcription were top biological processes involved in the drought response and tolerance. These results provide insights into stress responses regulated by the 4BS locus and have identified candidate genes and genetic markers that can be used for fine mapping of the qDSI.4B.1 locus and, ultimately, in wheat breeding programs for drought tolerance.

Investigating Host and Parasitic Plant Interaction by Tissue-Specific Gene Analyses on Tomato and Cuscuta campestris Interface at Three Haustorial Developmental Stages

Parasitic weeds cause billions of dollars in agricultural losses each year worldwide. Cuscuta campestris (C. campestris), one of the most widespread and destructive parasitic plants in the United States, severely reduces yield in tomato plants. Reducing the spread of parasitic weeds requires understanding the interaction between parasites and hosts. Several studies have identified factors needed for parasitic plant germination and haustorium induction, and genes involved in host defense responses. However, knowledge of the mechanisms underlying the interactions between host and parasitic plants, specifically at the interface between the two organisms, is relatively limited. A detailed investigation of the crosstalk between the host and parasite at the tissue-specific level would enable development of effective parasite control strategies. To focus on the haustorial interface, we used laser-capture microdissection (LCM) with RNA-seq on early, intermediate and mature haustorial stages. In addition, the tomato host tissue that immediately surround the haustoria was collected to obtain tissue- resolution RNA-Seq profiles for C. campestris and tomato at the parasitism interface. After conducting RNA-Seq analysis and constructing gene coexpression networks (GCNs), we identified CcHB7, CcPMEI, and CcERF1 as putative key regulators involved in C. campestris haustorium organogenesis, and three potential regulators, SlPR1, SlCuRe1-like, and SlNLR, in tomatoes that are involved in perceiving signals from the parasite. We used host-induced gene silencing (HIGS) transgenic tomatoes to knock-down the candidate genes in C. campestris and produced CRISPR transgenic tomatoes to knock out candidate genes in tomatoes. The interactions of C. campestris with these transgenic lines were tested and compared with that in wild-type tomatoes. The results of this study reveal the tissue-resolution gene regulatory mechanisms at the parasitic plant-host interface and provide the potential of developing a parasite-resistant system in tomatoes.

NLR diversity and candidate fusiform rust resistance genes in loblolly pine

Resistance to fusiform rust disease in loblolly pine (Pinus taeda) is a classic gene-for-gene system. Early resistance gene mapping in the P. taeda family 10-5 identified RAPD markers for a major fusiform rust resistance gene, Fr1. More recently, single nucleotide polymorphism (SNP) markers associated with resistance were mapped to a full-length gene model in the loblolly pine genome encoding for a nucleotide-binding site leucine-rich repeat (NLR) protein. NLR genes are one of the most abundant gene families in plant genomes and are involved in effector-triggered immunity. Inter- and intraspecies studies of NLR gene diversity and expression have resulted in improved disease resistance. To characterize NLR gene diversity and discover potential resistance genes, we assembled de novo transcriptomes from 92 loblolly genotypes from across the natural range of the species. In these transcriptomes, we identified novel NLR transcripts that are not present in the loblolly pine reference genome and found significant geographic diversity of NLR genes providing evidence of gene family evolution. We designed capture probes for these NLRs to identify and map SNPs that stably cosegregate with resistance to the SC20-21 isolate of Cronartium quercuum f.sp. fusiforme (Cqf) in half-sib progeny of the 10-5 family. We identified 10 SNPs and 2 quantitative trait loci associated with resistance to SC20-21 Cqf. The geographic diversity of NLR genes provides evidence of NLR gene family evolution in loblolly pine. The SNPs associated with rust resistance provide a resource to enhance breeding and deployment of resistant pine seedlings.

Structural basis of NLR activation and innate immune signalling in plants

Animals and plants have NLRs (nucleotide-binding leucine-rich repeat receptors) that recognize the presence of pathogens and initiate innate immune responses. In plants, there are three types of NLRs distinguished by their N-terminal domain: the CC (coiled-coil) domain NLRs, the TIR (Toll/interleukin-1 receptor) domain NLRs and the RPW8 (resistance to powdery mildew 8)-like coiled-coil domain NLRs. CC-NLRs (CNLs) and TIR-NLRs (TNLs) generally act as sensors of effectors secreted by pathogens, while RPW8-NLRs (RNLs) signal downstream of many sensor NLRs and are called helper NLRs. Recent studies have revealed three dimensional structures of a CNL (ZAR1) including its inactive, intermediate and active oligomeric state, as well as TNLs (RPP1 and ROQ1) in their active oligomeric states. Furthermore, accumulating evidence suggests that members of the family of lipase-like EDS1 (enhanced disease susceptibility 1) proteins, which are uniquely found in seed plants, play a key role in providing a link between sensor NLRs and helper NLRs during innate immune responses. Here, we summarize the implications of the plant NLR structures that provide insights into distinct mechanisms of action by the different sensor NLRs and discuss plant NLR-mediated innate immune signalling pathways involving the EDS1 family proteins and RNLs.

RefPlantNLR is a comprehensive collection of experimentally validated plant disease resistance proteins from the NLR family

Reference datasets are critical in computational biology. They help define canonical biological features and are essential for benchmarking studies. Here, we describe a comprehensive reference dataset of experimentally validated plant nucleotide-binding leucine-rich repeat (NLR) immune receptors. RefPlantNLR consists of 481 NLRs from 31 genera belonging to 11 orders of flowering plants. This reference dataset has several applications. We used RefPlantNLR to determine the canonical features of functionally validated plant NLRs and to benchmark 5 NLR annotation tools. This revealed that although NLR annotation tools tend to retrieve the majority of NLRs, they frequently produce domain architectures that are inconsistent with the RefPlantNLR annotation. Guided by this analysis, we developed a new pipeline, NLRtracker, which extracts and annotates NLRs from protein or transcript files based on the core features found in the RefPlantNLR dataset. The RefPlantNLR dataset should also prove useful for guiding comparative analyses of NLRs across the wide spectrum of plant diversity and identifying understudied taxa. We hope that the RefPlantNLR resource will contribute to moving the field beyond a uniform view of NLR structure and function.

Major QTL confer race-nonspecific resistance in the co-evolved Cronartium quercuum f. sp. fusiforme-Pinus taeda pathosystem

Fusiform rust disease, caused by the endemic fungus Cronartium quercuum f. sp. fusiforme, is the most damaging disease affecting economically important pine species in the southeast United States. Unlike the major epidemics of agricultural crops, the co-evolved pine-rust pathosystem is characterized by steady-state dynamics and high levels of genetic diversity within environments. This poses a unique challenge and opportunity for the deployment of large-effect resistance genes. We used trait dissection to study the genetic architecture of disease resistance in two P. taeda parents that showed high resistance across multiple environments. Two mapping populations (full-sib families), each with ~1000 progeny, were challenged with a complex inoculum consisting of 150 pathogen isolates. High-density linkage mapping revealed three major-effect QTL distributed on two linkage groups. All three QTL were validated using a population of 2057 cloned pine genotypes in a 6-year-old multi-environmental field trial. As a complement to the QTL mapping approach, bulked segregant RNAseq analysis revealed a small number of candidate nucleotide binding leucine-rich repeat genes harboring SNP associated with disease resistance. The results of this study show that in P. taeda, a small number of major QTL can provide effective resistance against genetically diverse mixtures of an endemic pathogen. These QTL vary in their impact on disease liability and exhibit additivity in combination.

EDS1 signalling: At the nexus of intracellular and surface receptor immunity

The conserved lipase-like protein EDS1 transduces signals from pathogen-activated intracellular nucleotide-binding leucine-rich repeat (NLR) receptors to transcriptional defences and host cell death. In this pivotal NLR signalling role, EDS1 works as a heterodimer with each of its partners, SAG101 and PAD4. Different properties of EDS1-SAG101 and EDS1-PAD4 complexes and functional relationships to sensor and helper NLRs have emerged. EDS1-SAG101 dimers confer effector-triggered immunity mediated by intracellular TNL receptors. In contrast, EDS1-PAD4 dimers have a broader role promoting basal immune responses that can be initiated inside cells by TNL- or CNL-type NLRs, and at the cell surface by LRR-receptor proteins. Characterizing the essential elements of these two EDS1 modules will help to connect intracellular and surface receptor signalling networks in the plant immune system.

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.

Plant "helper" immune receptors are Ca2+-permeable nonselective cation channels

Plant nucleotide-binding leucine-rich repeat receptors (NLRs) regulate immunity and cell death. In Arabidopsis, a subfamily of "helper" NLRs is required by many "sensor" NLRs. Active NRG1.1 oligomerized, was enriched in plasma membrane puncta, and conferred cytoplasmic calcium ion (Ca2+) influx in plant and human cells. NRG1.1-dependent Ca2+ influx and cell death were sensitive to Ca2+ channel blockers and were suppressed by mutations affecting oligomerization or plasma membrane enrichment. Ca2+ influx and cell death mediated by NRG1.1 and ACTIVATED DISEASE RESISTANCE 1 (ADR1), another helper NLR, required conserved negatively charged N-terminal residues. Whole-cell voltage-clamp recordings demonstrated that Arabidopsis helper NLRs form Ca2+-permeable cation channels to directly regulate cytoplasmic Ca2+ levels and consequent cell death. Thus, helper NLRs transduce cell death signals directly.

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.

Molecular evidence of the avocado defense response to Fusarium kuroshium infection: a deep transcriptome analysis using RNA-Seq

Fusarium kuroshium is a novel member of the Ambrosia Fusarium Clade (AFC) that has been recognized as one of the symbionts of the invasive Kuroshio shot hole borer, an Asian ambrosia beetle. This complex is considered the causal agent of Fusarium dieback, a disease that has severely threatened natural forests, landscape trees, and avocado orchards in the last 8 years. Despite the interest in this species, the molecular responses of both the host and F. kuroshium during the infection process and disease establishment remain unknown. In this work, we established an in vitro pathosystem using Hass avocado stems inoculated with F. kuroshium to investigate differential gene expression at 1, 4, 7 and 14 days post-inoculation. RNA-seq technology allowed us to obtain data from both the plant and the fungus, and the sequences obtained from both organisms were analyzed independently. The pathosystem established was able to mimic Fusarium dieback symptoms, such as carbohydrate exudation, necrosis, and vascular tissue discoloration. The results provide interesting evidence regarding the genes that may play roles in the avocado defense response to Fusarium dieback disease. The avocado data set comprised a coding sequence collection of 51,379 UniGenes, from which 2,403 (4.67%) were identified as differentially expressed. The global expression analysis showed that F. kuroshium responsive UniGenes can be clustered into six groups according to their expression profiles. The biologically relevant functional categories that were identified included photosynthesis as well as responses to stress, hormones, abscisic acid, and water deprivation. Additionally, processes such as oxidation-reduction, organization and biogenesis of the cell wall and polysaccharide metabolism were detected. Moreover, we identified orthologues of nucleotide-binding leucine-rich receptors, and their possible action mode was analyzed. In F. kuroshium, we identified 57 differentially expressed genes. Interestingly, the alcohol metabolic process biological category had the highest number of upregulated genes, and the enzyme group in this category may play an important role in the mechanisms of secondary metabolite detoxification. Hydrolytic enzymes, such as endoglucanases and a pectate lyase, were also identified, as well as some proteases. In conclusion, our research was conducted mainly to explain how the vascular tissue of a recognized host of the ambrosia complex responds during F. kuroshium infection since Fusarium dieback is an ambrosia beetle-vectored disease and many variables facilitate its establishment.

Integrity of the Post-LRR Domain Is Required for TIR-NB-LRR Function

Plants trigger appropriate defense responses, notably, through intracellular nucleotide-binding (NB) and leucine-rich repeat (LRR)-containing receptors (NLRs) that detect secreted pathogen effector proteins. In NLR resistance genes, the toll/interleukin-1 receptor (TIR)-NB-LRR proteins (TNLs) are an important subfamily, out of which approximately half the members carry a post-LRR (PL) domain of unknown role. We first investigated the requirement of the PL domain for TNL-mediated immune response by mutating the most conserved amino acids across PL domains of Arabidopsis thaliana TNLs. We identified several amino acids in the PL domain of RPS4, required for its ability to trigger a hypersensitive response to AvrRps4 in a Nicotiana tabacum transient assay. Mutating the corresponding amino acids within the PL domain of the tobacco TNL gene N also affected its function. Consequently, our results indicate that the integrity of the PL domain at conserved positions is crucial for at least two unrelated TNLs. We then tested the PL domain specificity for function by swapping PL domains between the paralogs RPS4 and RPS4B. Our results suggest that the PL domain is involved in their TNL pair specificity, 'off state' stability, and NLR complex activation. Considering genetically paired Arabidopsis TNLs, we finally compared the PL and TIR domains of their sensor and executor sequences, respectively. While TIR and PL domains from executors present complete motifs, sensors showed a lack of conservation with degenerated motifs. We here provide a contribution to the functional analysis of the PL domain in order to decipher its role for TNL function.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Evolution of NLR Resistance Genes in Magnoliids: Dramatic Expansions of CNLs and Multiple Losses of TNLs

Magnoliids are the third-largest group of angiosperms and occupy a critical position in angiosperm evolution. In the past years, due to the lack of sequenced genomes, the disease resistance gene (R gene) profile of magnoliids remains poorly understood. By the genome-wide identification of 1,832 NLR genes from seven magnoliid genomes, we built a framework for the evolution of magnoliid R genes. TNL genes were completely absent from five magnoliids, presumably due to immune pathway deficiencies. A total of 74 ancestral R genes (70 CNLs, 3 TNLs, and 1 RNL) were recovered in a common ancestor of magnoliids, from which all current NLR gene repertoires were derived. Tandem duplication served as the major drive for NLR genes expansion in seven magnoliid genomes, as most surveyed angiosperms. Due to recent rapid expansions, most magnoliids exhibited "a first expansion followed by a slight contraction and a further stronger expansion" evolutionary pattern, while both Litsea cubeba and Persea americana showed a two-times-repeated pattern of "expansion followed by contraction." The transcriptome analysis of seven different tissues of Saururus chinensis revealed a low expression of most NLR genes, with some R genes displaying a relatively higher expression in roots and fruits. Overall, our study sheds light on the evolution of NLR genes in magnoliids, compensates for insufficiency in major angiosperm lineages, and provides an important reference for a better understanding of angiosperm NLR genes.

A Reference Genome Sequence for Giant Sequoia

The giant sequoia (Sequoiadendron giganteum) of California are massive, long-lived trees that grow along the U.S. Sierra Nevada mountains. Genomic data are limited in giant sequoia and producing a reference genome sequence has been an important goal to allow marker development for restoration and management. Using deep-coverage Illumina and Oxford Nanopore sequencing, combined with Dovetail chromosome conformation capture libraries, the genome was assembled into eleven chromosome-scale scaffolds containing 8.125 Gbp of sequence. Iso-Seq transcripts, assembled from three distinct tissues, was used as evidence to annotate a total of 41,632 protein-coding genes. The genome was found to contain, distributed unevenly across all 11 chromosomes and in 63 orthogroups, over 900 complete or partial predicted NLR genes, of which 375 are supported by annotation derived from protein evidence and gene modeling. This giant sequoia reference genome sequence represents the first genome sequenced in the Cupressaceae family, and lays a foundation for using genomic tools to aid in giant sequoia conservation and management.

Convergent Loss of an EDS1/PAD4 Signaling Pathway in Several Plant Lineages Reveals Coevolved Components of Plant Immunity and Drought Response

Plant innate immunity relies on nucleotide binding leucine-rich repeat receptors (NLRs) that recognize pathogen-derived molecules and activate downstream signaling pathways. We analyzed the variation in NLR gene copy number and identified plants with a low number of NLR genes relative to sister species. We specifically focused on four plants from two distinct lineages, one monocot lineage (Alismatales) and one eudicot lineage (Lentibulariaceae). In these lineages, the loss of NLR genes coincides with loss of the well-known downstream immune signaling complex ENHANCED DISEASE SUSCEPTIBILITY 1 (EDS1)/PHYTOALEXIN DEFICIENT 4 (PAD4). We expanded our analysis across whole proteomes and found that other characterized immune genes were absent only in Lentibulariaceae and Alismatales. Additionally, we identified genes of unknown function that were convergently lost together with EDS1/PAD4 in five plant species. Gene expression analyses in Arabidopsis (Arabidopsis thaliana) and Oryza sativa revealed that several homologs of the candidates are differentially expressed during pathogen infection, drought, and abscisic acid treatment. Our analysis provides evolutionary evidence for the rewiring of plant immunity in some plant lineages, as well as the coevolution of the EDS1/PAD4 pathway and drought responses.