Diversity and Evolution of Disease Resistance Genes in Barley (Hordeum vulgare L.)

Exon shuffling potentiates a diverse repertoire of brown algal NB-ARC-TPR candidate immune receptor proteins via alternative splicing

Like other organisms, brown algae are subject to diseases caused by bacteria, fungi, and viruses. Brown algal immunity mechanisms are not well characterized; however, there is evidence suggesting that pathogen receptors exist in brown algae. One key protein family likely associated with brown algal innate immunity possesses an NB-ARC domain analogous to innate immune proteins in plants and animals. In this study, we conducted an extensive survey of NB-ARC genes in brown algae and obtained insights into the domain organization and evolutionary history of the encoded proteins. Our data show that brown algae possess an ancient NB-ARC-tetratricopeptide repeat (NB-TPR) domain architecture. We identified an N-terminal effector domain, the four-helix bundle, which was not previously found associated with NB-ARC domains. The phylogenetic tree including NB-ARC domains from all kingdoms of life suggests the three clades of brown algal NB-TPRs are likely monophyletic, whereas their TPRs seem to have distinct origins. One group of TPRs exhibit intense exon shuffling, with various alternative splicing and diversifying selection acting on them, suggesting exon shuffling is an important mechanism for evolving ligand-binding specificities. The reconciliation of gene duplication and loss events of the NB-ARC genes reveals that more independent gene gains than losses have occurred during brown algal evolution, and that tandem duplication has played a major role in the expansion of NB-ARC genes. Our results substantially enhance our understanding of the evolutionary history and exon shuffling mechanisms of the candidate innate immune repertoire of brown algae.

An improved assembly of the "Cascade" hop (Humulus lupulus) genome uncovers signatures of molecular evolution and refines time of divergence estimates for the Cannabaceae family

We present a chromosome-level assembly of the Cascade hop (Humulus lupulus L. var. lupulus) genome. The hop genome is large (2.8 Gb) and complex, and early attempts at assembly were fragmented. Recent advances have made assembly of the hop genome more tractable, transforming the extent of investigation that can occur. The chromosome-level assembly of Cascade was developed by scaffolding the previously reported Cascade assembly generated with PacBio long-read sequencing and polishing with Illumina short-read DNA sequencing. We developed gene models and repeat annotations and used a controlled bi-parental mapping population to identify significant sex-associated markers. We assessed molecular evolution in gene sequences, gene family expansion and contraction, and time of divergence from Cannabis sativa and other closely related plant species using Bayesian inference. We identified the putative sex chromosome in the female genome based on significant sex-associated markers from the bi-parental mapping population. While the estimate of repeat content (~64%) is similar to the estimate for the hemp genome, syntenic blocks in hop contain a greater percentage of LTRs. Hop is enriched for disease resistance-associated genes in syntenic gene blocks and expanded gene families. The Cascade chromosome-level assembly will inform cultivation strategies and serve to deepen our understanding of the hop genomic landscape, benefiting hop researchers and the Cannabaceae genomics community.

A high-quality pseudo-phased genome for Melaleuca quinquenervia shows allelic diversity of NLR-type resistance genes

BACKGROUND

Melaleuca quinquenervia (broad-leaved paperbark) is a coastal wetland tree species that serves as a foundation species in eastern Australia, Indonesia, Papua New Guinea, and New Caledonia. While extensively cultivated for its ornamental value, it has also become invasive in regions like Florida, USA. Long-lived trees face diverse pest and pathogen pressures, and plant stress responses rely on immune receptors encoded by the nucleotide-binding leucine-rich repeat (NLR) gene family. However, the comprehensive annotation of NLR encoding genes has been challenging due to their clustering arrangement on chromosomes and highly repetitive domain structure; expansion of the NLR gene family is driven largely by tandem duplication. Additionally, the allelic diversity of the NLR gene family remains largely unexplored in outcrossing tree species, as many genomes are presented in their haploid, collapsed state.

RESULTS

We assembled a chromosome-level pseudo-phased genome for M. quinquenervia and described the allelic diversity of plant NLRs using the novel FindPlantNLRs pipeline. Analysis reveals variation in the number of NLR genes on each haplotype, distinct clustering patterns, and differences in the types and numbers of novel integrated domains.

CONCLUSIONS

The high-quality M. quinquenervia genome assembly establishes a new framework for functional and evolutionary studies of this significant tree species. Our findings suggest that maintaining allelic diversity within the NLR gene family is crucial for enabling responses to environmental stress, particularly in long-lived plants.

Unraveling the genetics of polyamine metabolism in barley for senescence-related crop improvement

We explored the polyamine (PA) metabolic pathway genes in barley (Hv) to understand plant development and stress adaptation in Gramineae crops with emphasis on leaf senescence. Bioinformatics and functional genomics tools were utilized for genome-wide identification, comprehensive gene features, evolution, development and stress effects on the expression of the polyamine metabolic pathway gene families (PMGs). Three S-adenosylmethionine decarboxylases (HvSAMDCs), two ornithine decarboxylases (HvODCs), one arginine decarboxylase (HvADC), one spermidine synthase (HvSPDS), two spermine synthases (HvSPMSs), five copper amine oxidases (HvCuAOs) and seven polyamine oxidases (HvPAOs) members of PMGs were identified and characterized in barley. All the HvPMG genes were found to be distributed on all chromosomes of barley. The phylogenetic and comparative assessment revealed that PA metabolic pathway is highly conserved in plants and the prediction of nine H. vulgare miRNAs (hvu-miR) target sites, 18 protein-protein interactions and 961 putative CREs in the promoter region were discerned. Gene expression of HvSAMDC3, HvCuAO7, HvPAO4 and HvSPMS1 was apparent at every developmental stage. SPDS/SPMS gene family was found to be the most responsive to induced leaf senescence. This study provides a reference for the functional investigation of the molecular mechanism(s) that regulate polyamine metabolism in plants as a tool for future breeding decision management systems.

Chromosome-level genome assembly of the aquatic plant Nymphoides indica reveals transposable element bursts and NBS-LRR gene family expansion shedding light on its invasiveness

Nymphoides indica, an aquatic plant, is an invasive species that causes both ecological and economic damage in North America and elsewhere. However, the lack of genomic data of N. indica limits the in-depth analysis of this invasive species. Here, we report a chromosome-level genome assembly of nine pseudochromosomes of N. indica with a total size of ∼ 520 Mb. More than half of the N. indica genome consists of transposable elements (TEs), and a higher density of TEs around genes may play a significant role in response to an ever-changing environment by regulating the nearby gene. Additionally, our analysis revealed that N. indica only experienced a gamma (γ) whole-genome triplication event. Functional enrichment of the N. indica-specific and expanded gene families highlighted genes involved in the responses to hypoxia and plant–pathogen interactions, which may strengthen the ability to adapt to external challenges and improve ecological fitness. Furthermore, we identified 160 members of the nucleotide-binding site and leucine-rich repeat gene family, which may be linked to the defence response. Collectively, the high-quality N. indica genome reported here opens a novel avenue to understand the evolution and rapid invasion of Nymphoides spp.

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.

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.

Distinct Evolutionary Patterns of NBS-Encoding Genes in Three Soapberry Family (Sapindaceae) Species

Nucleotide-binding site (NBS)-type disease resistance genes (R genes) play key roles in plant immune responses and have co-evolved with pathogens over the course of plant lifecycles. Comparative genomic studies tracing the dynamic evolution of NBS-encoding genes have been conducted using many important plant lineages. However, studies on Sapindaceae species have not been performed. In this study, a discrepant number of NBS-encoding genes were identified in the genomes of Xanthoceras sorbifolium (180), Dinnocarpus longan (568), and Acer yangbiense (252). These genes were unevenly distributed and usually clustered as tandem arrays on chromosomes, with few existed as singletons. The phylogenetic analysis revealed that NBS-encoding genes formed three monophyletic clades, RPW8-NBS-LRR (RNL), TIR-NBS-LRR (TNL), and CC-NBS-LRR (CNL), which were distinguished by amino acid motifs. The NBS-encoding genes of the X. sorbifolium, D. longan, and A. yangbiense genomes were derived from 181 ancestral genes (three RNL, 23 TNL, and 155 CNL), which exhibited dynamic and distinct evolutionary patterns due to independent gene duplication/loss events. Specifically, X. sorbifolium exhibited a “first expansion and then contraction” evolutionary pattern, while A. yangbiense and D. longan exhibited a “first expansion followed by contraction and further expansion” evolutionary pattern. However, further expansion in D. longan was stronger than in A. yangbiense after divergence, suggesting that D. longan gained more genes in response to various pathogens. Additionally, the ancient and recent expansion of CNL genes generated the dominance of this subclass in terms of gene numbers, while the low copy number status of RNL genes was attributed to their conserved functions.

Wheat Disease Resistance Genes and Their Diversification Through Integrated Domain Fusions

Plants are in a constant evolutionary arms race with their pathogens. At the molecular level, the plant nucleotide-binding leucine-rich repeat receptors (NLRs) family has coevolved with rapidly evolving pathogen effectors. While many NLRs utilize variable leucine-rich repeats (LRRs) to detect effectors, some have gained integrated domains (IDs) that may be involved in receptor activation or downstream signaling. The major objectives of this project were to identify NLR genes in wheat (Triticum aestivum L.) and assess IDs associated with immune signaling (e.g., kinase and transcription factor domains). We identified 2,151 NLR-like genes in wheat, of which 1,298 formed 547 gene clusters. Among the non-toll/interleukin-1 receptor NLR (non-TNL)-like genes, 1,552 encode LRRs, 802 are coiled-coil (CC) domain-encoding (CC-NBS-LRR or CNL) genes, and three encode resistance to powdery mildew 8 (RPW8) domains (RPW8-NBS-LRR or RNL). The expansion of the NLR gene family in wheat is attributable to its origin by recent polyploidy events. Gene clusters were likely formed by tandem duplications, and wheat NLR phylogenetic relationships were similar to those in barley and Aegilops. We also identified wheat NLR-ID fusion proteins as candidates for NLR functional diversification, often as kinase and transcription factor domains. Comparative analyses of the IDs revealed evolutionary conservation of more than 80% amino acid sequence similarity. Homology assessment indicates that these domains originated as functional non-NLR-encoding genes that were incorporated into NLR-encoding genes through duplication events. We also found that many of the NLR-ID genes encode alternative transcripts that include or exclude IDs, a phenomenon that seems to be conserved among species. To verify this, we have analyzed the alternative transcripts that include or exclude an ID of an NLR-ID from another monocotyledon species, rice (Oryza sativa). This indicates that plants employ alternative splicing to regulate IDs, possibly using them as baits, decoys, and functional signaling components. Genomic and expression data support the hypothesis that wheat uses alternative splicing to include and exclude IDs from NLR proteins.

Revisiting the Origin of Plant NBS-LRR Genes

The NBS-LRR genes are functionally responsible for plant resistance to alien pathogens. Here, we show that NBS-LRR genes originated in the common ancestor of the whole green lineage, and have rapidly diverged into three subclasses with different domain combinations (TNL, CNL, and RNL) before the split of green algae.

Genome-wide identification, characterization, and evolutionary analysis of NBS-encoding resistance genes in barley

In this study, a systematic analysis of Nucleotide-Binding Site (NBS) disease resistance (R) gene family in the barley, Hordeum vulgare L. cv. Bowman, genome was performed. Using multiple computational analyses, we could identify 96 regular NBS-encoding genes and characterize them on the bases of structural diversity, conserved protein signatures, genomic distribution, gene duplications, differential expression, selection pressure, codon usage, regulation by microRNAs and phylogenetic relationships. Depending on the presence or absence of CC and LRR domains; the identified NBS genes were assigned to four distinct groups; NBS-LRR (53.1%), CC-NBS-LRR (14.6%), NBS (26%), and CC-NBS (6.3%). NBS-associated domain analysis revealed the presence of signal peptides, zinc fingers, diverse kinases, and other structural features. Eighty-five of the identified NBS-encoding genes were mapped onto the seven barley chromosomes, revealing that 50% of them were located on chromosomes 7H, 2H, and 3H, with a tendency of NBS genes to be clustered in the distal telomeric regions of the barley chromosomes. Nine gene clusters, representing 22.35% of total mapped barley NBS-encoding genes, were found, suggesting that tandem duplication stands for an important mechanism in the expansion of this gene family in barley. Phylogenetic analysis determined 31 HvNBS orthologs from rice and Brachypodium. 87 out of 96 HvNBSs were supported by expression evidence, exhibiting various and quantitatively uneven expression patterns across distinct tissues, organs, and development stages. Fourteen potential miRNA-R gene target pairs were further identified, providing insight into the regulation of NBS genes expression. These findings offer candidate target genes to engineer disease-resistant barley genotypes, and promote our understanding of the evolution of NBS-encoding genes in Poaceae crops.

Evolutionary Divergence of TNL Disease-Resistant Proteins in Soybean (Glycine max) and Common Bean (Phaseolus vulgaris)

Disease-resistant genes (R genes) encode proteins that are involved in protecting plants from their pathogens and pests. Availability of complete genome sequences from soybean and common bean allowed us to perform a genome-wide identification and analysis of the Toll interleukin-1 receptor-like nucleotide-binding site leucine-rich repeat (TNL) proteins. Hidden Markov model (HMM) profiling of all protein sequences resulted in the identification of 117 and 77 regular TNL genes in soybean and common bean, respectively. We also identified TNL gene homologs with unique domains, and signal peptides as well as nuclear localization signals. The TNL genes in soybean formed 28 clusters located on 10 of the 20 chromosomes, with the majority found on chromosome 3, 6 and 16. Similarly, the TNL genes in common bean formed 14 clusters located on five of the 11 chromosomes, with the majority found on chromosome 10. Phylogenetic analyses of the TNL genes from Arabidopsis, soybean and common bean revealed less divergence within legumes relative to the divergence between legumes and Arabidopsis. Syntenic blocks were found between chromosomes Pv10 and Gm03, Pv07 and Gm10, as well as Pv01 and Gm14. The gene expression data revealed basal level expression and tissue specificity, while analysis of available microRNA data showed 37 predicted microRNA families involved in targeting the identified TNL genes in soybean and common bean.

Genome-Wide Identification of NBS-Encoding Resistance Genes in Sunflower (Helianthus annuus L.)

Nucleotide Binding Site—Leucine-Rich Repeat (NBS-LRR) genes encode disease resistance proteins involved in plants’ defense against their pathogens. Although sunflower is affected by many diseases, only a few molecular details have been uncovered regarding pathogenesis and resistance mechanisms. Recent availability of sunflower whole genome sequences in publicly accessible databases allowed us to accomplish a genome-wide identification of Toll-interleukin-1 receptor-like Nucleotide-binding site Leucine-rich repeat (TNL), Coiled Coil (CC)-NBS-LRR (CNL), Resistance to powdery mildew8 (RPW8)-NBS-LRR (RNL) and NBS-LRR (NL) protein encoding genes. Hidden Markov Model (HMM) profiling of 52,243 putative protein sequences from sunflower resulted in 352 NBS-encoding genes, among which 100 genes belong to CNL group including 64 genes with RX_CC like domain, 77 to TNL, 13 to RNL, and 162 belong to NL group. We also identified signal peptides and nuclear localization signals present in the identified genes and their homologs. We found that NBS genes were located on all chromosomes and formed 75 gene clusters, one-third of which were located on chromosome 13. Phylogenetic analyses between sunflower and Arabidopsis NBS genes revealed a clade-specific nesting pattern in CNLs, with RNLs nested in the CNL-A clade, and species-specific nesting pattern for TNLs. Surprisingly, we found a moderate bootstrap support (BS = 50%) for CNL-A clade being nested within TNL clade making both the CNL and TNL clades paraphyletic. Arabidopsis and sunflower showed 87 syntenic blocks with 1049 high synteny hits between chromosome 5 of Arabidopsis and chromosome 6 of sunflower. Expression data revealed functional divergence of the NBS genes with basal level tissue-specific expression. This study represents the first genome-wide identification of NBS genes in sunflower paving avenues for functional characterization and potential crop improvement.

Disease Resistance Mechanisms in Plants

Plants have developed a complex defense system against diverse pests and pathogens. Once pathogens overcome mechanical barriers to infection, plant receptors initiate signaling pathways driving the expression of defense response genes. Plant immune systems rely on their ability to recognize enemy molecules, carry out signal transduction, and respond defensively through pathways involving many genes and their products. Pathogens actively attempt to evade and interfere with response pathways, selecting for a decentralized, multicomponent immune system. Recent advances in molecular techniques have greatly expanded our understanding of plant immunity, largely driven by potential application to agricultural systems. Here, we review the major plant immune system components, state of the art knowledge, and future direction of research on plant⁻pathogen interactions. In our review, we will discuss how the decentralization of plant immune systems have provided both increased evolutionary opportunity for pathogen resistance, as well as additional mechanisms for pathogen inhibition of such defense responses. We conclude that the rapid advances in bioinformatics and molecular biology are driving an explosion of information that will advance agricultural production and illustrate how complex molecular interactions evolve.

Identification of QTL conferring resistance to stripe rust (Puccinia striiformis f. sp. hordei) and leaf rust (Puccinia hordei) in barley using nested association mapping (NAM)

The biotrophic rust fungi Puccinia hordei and Puccinia striiformis are important barley pathogens with the potential to cause high yield losses through an epidemic spread. The identification of QTL conferring resistance to these pathogens is the basis for targeted breeding approaches aiming to improve stripe rust and leaf rust resistance of modern cultivars. Exploiting the allelic richness of wild barley accessions proved to be a valuable tool to broaden the genetic base of resistance of barley cultivars. In this study, SNP-based nested association mapping (NAM) was performed to map stripe rust and leaf rust resistance QTL in the barley NAM population HEB-25, comprising 1,420 lines derived from BC₁S₃ generation. By scoring the percentage of infected leaf area, followed by calculation of the area under the disease progress curve and the average ordinate during a two-year field trial, a large variability of resistance across and within HEB-25 families was observed. NAM based on 5,715 informative SNPs resulted in the identification of twelve and eleven robust QTL for resistance against stripe rust and leaf rust, respectively. Out of these, eight QTL for stripe rust and two QTL for leaf rust are considered novel showing no overlap with previously reported resistance QTL. Overall, resistance to both pathogens in HEB-25 is most likely due to the accumulation of numerous small effect loci. In addition, the NAM results indicate that the 25 wild donor QTL alleles present in HEB-25 strongly differ in regard to their individual effect on rust resistance. In future, the NAM concept will allow to select and combine individual wild barley alleles from different HEB parents to increase rust resistance in barley. The HEB-25 results will support to unravel the genetic basis of rust resistance in barley, and to improve resistance against stripe rust and leaf rust of modern barley cultivars.

Pathogen-inducible Ta-Lr34res expression in heterologous barley confers disease resistance without negative pleiotropic effects

Plant diseases are a serious threat to crop production. The informed use of naturally occurring disease resistance in plant breeding can greatly contribute to sustainably reduce yield losses caused by plant pathogens. The Ta-Lr34res gene encodes an ABC transporter protein and confers partial, durable, and broad spectrum resistance against several fungal pathogens in wheat. Transgenic barley lines expressing Ta-Lr34res showed enhanced resistance against powdery mildew and leaf rust of barley. While Ta-Lr34res is only active at adult stage in wheat, Ta-Lr34res was found to be highly expressed already at the seedling stage in transgenic barley resulting in severe negative effects on growth. Here, we expressed Ta-Lr34res under the control of the pathogen-inducible Hv-Ger4c promoter in barley. Sixteen independent barley transformants showed strong resistance against leaf rust and powdery mildew. Infection assays and growth parameter measurements were performed under standard glasshouse and near-field conditions using a convertible glasshouse. Two Hv-Ger4c::Ta-Lr34res transgenic events were analysed in detail. Plants of one transformation event had similar grain production compared to wild-type under glasshouse and near-field conditions. Our results showed that negative effects caused by constitutive high expression of Ta-Lr34res driven by the endogenous wheat promoter in barley can be eliminated by inducible expression without compromising disease resistance. These data demonstrate that Ta-Lr34res is agronomically useful in barley. We conclude that the generation of a large number of transformants in different barley cultivars followed by early field testing will allow identifying barley lines suitable for breeding.