Co-variation among major classes of LRR-encoding genes in two pairs of plant species

Regulation of plant immunity via small RNA-mediated control of NLR expression

Plants use different receptors to detect potential pathogens: membrane-anchored pattern recognition receptors (PRRs) activated upon perception of pathogen-associated molecular patterns (PAMPs) that elicit pattern-triggered immunity (PTI); and intracellular nucleotide-binding leucine-rich repeat proteins (NLRs) activated by detection of pathogen-derived effectors, activating effector-triggered immunity (ETI). The interconnections between PTI and ETI responses have been increasingly reported. Elevated NLR levels may cause autoimmunity, with symptoms ranging from fitness cost to developmental arrest, sometimes combined with run-away cell death, making accurate control of NLR dosage key for plant survival. Small RNA-mediated gene regulation has emerged as a major mechanism of control of NLR dosage. Twenty-two nucleotide miRNAs with the unique ability to trigger secondary siRNA production from target transcripts are particularly prevalent in NLR regulation. They enhance repression of the primary NLR target, but also bring about repression of NLRs only complementary to secondary siRNAs. We summarize current knowledge on miRNAs and siRNAs in the regulation of NLR expression with an emphasis on 22 nt miRNAs and propose that miRNA and siRNA regulation of NLR levels provides additional links between PTI and NLR defense pathways to increase plant responsiveness against a broad spectrum of pathogens and control an efficient deployment of defenses.

Comprehensive analysis of peptide-coding genes and initial characterization of an LRR-only microprotein in Marchantia polymorpha

In the past two decades, many plant peptides have been found to play crucial roles in various biological events by mediating cell-to-cell communications. However, a large number of small open reading frames (sORFs) or short genes capable of encoding peptides remain uncharacterized. In this study, we examined several candidate genes for peptides conserved between two model plants: Arabidopsis thaliana and Marchantia polymorpha. We examined their expression pattern in M. polymorpha and subcellular localization using a transient assay with Nicotiana benthamiana. We found that one candidate, MpSGF10B, was expressed in meristems, gemma cups, and male reproductive organs called antheridiophores. MpSGF10B has an N-terminal signal peptide followed by two leucine-rich repeat (LRR) domains and was secreted to the extracellular region in N. benthamiana and M. polymorpha. Compared with the wild type, two independent Mpsgf10b mutants had a slightly increased number of antheridiophores. It was revealed in gene ontology enrichment analysis that MpSGF10B was significantly co-expressed with genes related to cell cycle and development. These results suggest that MpSGF10B may be involved in the reproductive development of M. polymorpha. Our research should shed light on the unknown role of LRR-only proteins in land plants.

Large-scale stage-specific regulation of gene expression during host-pathogen interactions in CSP44 bread wheat carrying APR gene Lr48

Genome-wide transcriptome analysis was undertaken in a leaf-rust resistant bread wheat line CSP44 (selected from Australian cv. Condor) carrying the adult plant resistance (APR) gene Lr48. Two pre-adult plant (P-AP) susceptible stages (S48 and S96) and two adult plant (AP) resistant stages (R48 and R96) were used for RNA-seq. At the susceptible P-AP stage (during S48 to S96), expression increased in 2062 genes, and declined in 130 genes; 1775 of 2062 differentially expressed genes (DEGs) also exhibited high expression during early incompatible stage R48. Comparison of S96 with R96 showed that the expression of 80 genes was enhanced and that of 208 genes declined at the AP stage. At the resistant AP stage (during R48 to R96), expression of mere 25 genes increased and that of 126 genes declined. Apparently, the resistance during late adult stage (R96) is caused by regulation of the expression of relatively fewer genes, although at pre-adult stage (S48 to S96), expression of large number of genes increased; expression of majority of these genes kept on increasing during adult stage at R48 also. These and other results of the present study suggest that APR may mimic some kind of systemic acquired resistance (SAR). The host-specific DEGs belonged to 10 different classes including genes involved in defence, transport, epigenetics, photosynthesis, genes encoding some transcription factors etc. The pathogen (Puccinia triticina) specific DEGs (including three genes encoding known biotrophic effectors) seem to help the pathogen in infection/growth through large-scale stage-specific enhanced expression of host's genes. A putative candidate gene for Lr48 containing protein kinase domain (its ortholog in rice encoding OsWAK8) was also identified.

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.

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.

Origin and diversification of leucine-rich repeat receptor-like protein kinase (LRR-RLK) genes in plants

Background

Leucine-rich repeat receptor-like protein kinases (LRR-RLKs) are the largest group of receptor-like kinases in plants and play crucial roles in development and stress responses. The evolutionary relationships among LRR-RLK genes have been investigated in flowering plants; however, no comprehensive studies have been performed for these genes in more ancestral groups. The subfamily classification of LRR-RLK genes in plants, the evolutionary history and driving force for the evolution of each LRR-RLK subfamily remain to be understood.

Results

We identified 119 LRR-RLK genes in the Physcomitrella patens moss genome, 67 LRR-RLK genes in the Selaginella moellendorffii lycophyte genome, and no LRR-RLK genes in five green algae genomes. Furthermore, these LRR-RLK sequences, along with previously reported LRR-RLK sequences from Arabidopsis thaliana and Oryza sativa, were subjected to evolutionary analyses. Phylogenetic analyses revealed that plant LRR-RLKs belong to 19 subfamilies, eighteen of which were established in early land plants, and one of which evolved in flowering plants. More importantly, we found that the basic structures of LRR-RLK genes for most subfamilies are established in early land plants and conserved within subfamilies and across different plant lineages, but divergent among subfamilies. In addition, most members of the same subfamily had common protein motif compositions, whereas members of different subfamilies showed variations in protein motif compositions. The unique gene structure and protein motif compositions of each subfamily differentiate the subfamily classifications and, more importantly, provide evidence for functional divergence among LRR-RLK subfamilies. Maximum likelihood analyses showed that some sites within four subfamilies were under positive selection.

Conclusions

Much of the diversity of plant LRR-RLK genes was established in early land plants. Positive selection contributed to the evolution of a few LRR-RLK subfamilies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0891-5) contains supplementary material, which is available to authorized users.

The Diversification of Plant NBS-LRR Defense Genes Directs the Evolution of MicroRNAs That Target Them

High expression of plant nucleotide binding site leucine-rich repeat (NBS-LRR) defense genes is often lethal to plant cells, a phenotype perhaps associated with fitness costs. Plants implement several mechanisms to control the transcript level of NBS-LRR defense genes. As negative transcriptional regulators, diverse miRNAs target NBS-LRRs in eudicots and gymnosperms. To understand the evolutionary benefits of this miRNA-NBS-LRR regulatory system, we investigated the NBS-LRRs of 70 land plants, coupling this analysis with extensive small RNA data. A tight association between the diversity of NBS-LRRs and miRNAs was found. The miRNAs typically target highly duplicated NBS-LRRs In comparison, families of heterogeneous NBS-LRRs were rarely targeted by miRNAs in Poaceae and Brassicaceae genomes. We observed that duplicated NBS-LRRs from different gene families periodically gave birth to new miRNAs. Most of these newly emerged miRNAs target the same conserved, encoded protein motif of NBS-LRRs, consistent with a model of convergent evolution for these miRNAs. By assessing the interactions between miRNAs and NBS-LRRs, we found nucleotide diversity in the wobble position of the codons in the target site drives the diversification of miRNAs. Taken together, we propose a co-evolutionary model of plant NBS-LRRs and miRNAs hypothesizing how plants balance the benefits and costs of NBS-LRR defense genes.

Small RNAs Add Zing to the Zig-Zag-Zig Model of Plant Defenses

Plant small RNAs play important roles in transcriptional and posttranscriptional regulation, with ongoing work demonstrating their functions in diverse pathways. Their roles in defense responses are a topic of active investigation, particularly the rich set of micro (mi)RNAs that target disease resistance genes such as nucleotide binding/leucine-rich repeat (NB-LRR) genes. The miRNA-NB-LRR interactions result in the production of phased, secondary small interfering (phasi)RNAs, and phasiRNAs function in both cis and trans to propagate negative regulatory effects across additional members of the target gene family. Yet, while phasiRNAs have the capacity to trigger targeted decay of specific targets, both in cis and trans, their functional relevance in NB-LRR regulation remains largely a matter of speculation.

Genome-Wide Identification and Characterization of the LRR-RLK Gene Family in Two Vernicia Species

Leucine-rich repeat receptor-like kinases (LRR-RLKs) make up the largest group of RLKs in plants and play important roles in many key biological processes such as pathogen response and signal transduction. To date, most studies on LRR-RLKs have been conducted on model plants. Here, we identified 236 and 230 LRR-RLKs in two industrial oil-producing trees: Vernicia fordii and Vernicia montana, respectively. Sequence alignment analyses showed that the homology of the RLK domain (23.81%) was greater than that of the LRR domain (9.51%) among the Vf/VmLRR-RLKs. The conserved motif of the LRR domain in Vf/VmLRR-RLKs matched well the known plant LRR consensus sequence but differed at the third last amino acid (W or L). Phylogenetic analysis revealed that Vf/VmLRR-RLKs were grouped into 16 subclades. We characterized the expression profiles of Vf/VmLRR-RLKs in various tissue types including root, leaf, petal, and kernel. Further investigation revealed that Vf/VmLRR-RLK orthologous genes mainly showed similar expression patterns in response to tree wilt disease, except 4 pairs of Vf/VmLRR-RLKs that showed opposite expression trends. These results represent an extensive evaluation of LRR-RLKs in two industrial oil trees and will be useful for further functional studies on these proteins.

Genomic and Post-Translational Modification Analysis of Leucine-Rich-Repeat Receptor-Like Kinases in Brassica rapa

Among several receptor-like kinases (RLKs), leucine-rich-repeat receptor-like kinases (LRR-RLKs) are a major group of genes that play crucial roles in growth, development and stress responses in plant systems. Given that they have several functional roles, it is important to investigate their roles in Brassica rapa. In the present study, 303 LRR-RLKs were identified in the genome of B. rapa and comparative phylogenetic analysis of 1213 combined LRR-RLKs of B. rapa, Arabidopsis thaliana, Oryza sativa and Populus trichocarpa helped us to categorize the gene family into 15 subfamilies based on their sequence and structural similarities. The chromosome localizations of 293 genes allowed the prediction of duplicates, and motif conservation and intron/exon patterns showed differences among the B. rapa LRR-RLK (BrLRR-RLK) genes. Additionally, computational function annotation and expression analysis was used to predict their possible functional roles in the plant system. Biochemical results for 11 selected genes showed variations in phosphorylation activity. Interestingly, BrBAK1 showed strong auto-phosphorylation and trans-phosphorylation on its tyrosine and threonine residues compared with AtBAK1 in previous studies. The AtBAK1 receptor kinase is involved in plant growth and development, plant innate immunity, and programmed cell death, and our results suggest that BrBAK1 might also be involved in the same functions. Another interesting result was that BrBAK1, BrBRI1, BrPEPR1 and BrPEPR2 showed activity with both anti-phosphotyrosine and anti-phosphothreonine antibodies, indicating that they might have dual-specificity kinase activity. This study provides comprehensive results for the BrLRR-RLKs, revealing expansion of the gene family through gene duplications, structural similarities and variations among the genes, and potential functional roles according to gene ontology, transcriptome profiling and biochemical analysis.

Evolutionary analysis of RB/Rpi-blb1 locus in the Solanaceae family

Late blight caused by the oomycete Phytophthora infestans is one of the most severe threats to potato production worldwide. Numerous studies suggest that the most effective protective strategy against the disease would be to provide potato cultivars with durable resistance (R) genes. However, little is known about the origin and evolutional history of these durable R-genes in potato. Addressing this might foster better understanding of the dynamics of these genes in nature and provide clues for identifying potential candidate R-genes. Here, a systematic survey was executed at RB/Rpi-blb1 locus, an exclusive broad-spectrum R-gene locus in potato. As indicated by synteny analysis, RB/Rpi-blb1 homologs were identified in all tested genomes, including potato, tomato, pepper, and Nicotiana, suggesting that the RB/Rpi-blb1 locus has an ancient origin. Two evolutionary patterns, similar to those reported on RGC2 in Lactuca, and Pi2/9 in rice, were detected at this locus. Type I RB/Rpi-blb1 homologs have frequent copy number variations and sequence exchanges, obscured orthologous relationships, considerable nucleotide divergence, and high non-synonymous to synonymous substitutions (Ka/Ks) between or within species, suggesting rapid diversification and balancing selection in response to rapid changes in the oomycete pathogen genomes. These characteristics may serve as signatures for cloning of late blight resistance genes.

Positive selection drives neofunctionalization of the UbiA prenyltransferase gene family

Gene duplication provides the key materials for new genes and novel functions. However, the mechanism underlying functional innovation remains unknown. In this study, we revealed the evolutionary pattern of the prenyltransferases of the UbiA gene family in 15 higher plants. Prenyltransferases of the UbiA gene family are involved in many important biological processes of both primary and secondary metabolism. Based on the phylogenetic relationships of the UbiA genes, seven subfamilies are classified. Confirming this classification, genes within each subfamily are characterized by similar exon numbers, exon lengths and patterns of motif combinations. Similar numbers of UbiA genes are found in different species within each subfamily except for Subfamily I, in which a Phaseoleae-specific expansion is detected in clade I-A. Homologous genes in clade I-A evolve rapidly, exchange sequences frequently and experience positive selection. Genes in clade I-A function as flavonoid prenyltransferase synthesis secondary compounds, while other genes from Subfamily I encode homogentisate phytyltransferase, which plays a role in primary metabolism. Thus, our results suggest that the secondary metabolism genes acquire new functions from those of primary metabolism through gene duplication and neofunctionalization driven by positive selection.

Paleo-evolutionary plasticity of plant disease resistance genes

Background

The recent access to a large set of genome sequences, combined with a robust evolutionary scenario of modern monocot (i.e. grasses) and eudicot (i.e. rosids) species from their founder ancestors, offered the opportunity to gain insights into disease resistance genes (R-genes) evolutionary plasticity.

Results

We unravel in the current article (i) a R-genes repertoire consisting in 7883 for monocots and 15758 for eudicots, (ii) a contrasted R-genes conservation with 23.8% for monocots and 6.6% for dicots, (iii) a minimal ancestral founder pool of 384 R-genes for the monocots and 150 R-genes for the eudicots, (iv) a general pattern of organization in clusters accounting for more than 60% of mapped R-genes, (v) a biased deletion of ancestral duplicated R-genes between paralogous blocks possibly compensated by clusterization, (vi) a bias in R-genes clusterization where Leucine-Rich Repeats act as a ‘glue’ for domain association, (vii) a R-genes/miRNAs interome enriched toward duplicated R-genes.

Conclusions

Together, our data may suggest that R-genes family plasticity operated during plant evolution (i) at the structural level through massive duplicates loss counterbalanced by massive clusterization following polyploidization; as well as at (ii) the regulation level through microRNA/R-gene interactions acting as a possible source of functional diploidization of structurally retained R-genes duplicates. Such evolutionary shuffling events leaded to CNVs (i.e. Copy Number Variation) and PAVs (i.e. Presence Absence Variation) between related species operating in the decay of R-genes colinearity between plant species.

Genome-wide identification, characterization and expression analysis of populus leucine-rich repeat receptor-like protein kinase genes

Background

Leucine-rich repeat receptor-like kinases (LRR-RLKs) comprise the largest group within the receptor-like kinase (RLK) superfamily in plants. This gene family plays critical and diverse roles in plant growth, development and stress response. Although the LRR-RLK families in Arabidopsis and rice have been previously analyzed, no comprehensive studies have been performed on this gene family in tree species.

Results

In this work, 379 LRR-RLK genes were retrieved from the Populus trichocarpa genome and further grouped into 14 subfamilies based on their structural and sequence similarities. Approximately 82% (312 out of 379) of the PtLRR-RLK genes are located in segmental duplication blocks indicating the role of duplication process in the expansion of this gene family. The conservation and variation in motif composition and intron/exon arrangement among PtLRR-RLK subfamilies were analyzed to provide additional support for their phylogenetic relationship and more importantly to indicate the potential divergence in their functions. Expression profiling of PtLRR-RLKs showed that they were differentially expressed in different organs and tissues and some PtLRR-RLKs were specifically expressed in meristem tissues, which indicated their potential involvement in tissue development and differentiation. For most AtLRR-RLKs with defined functions, Populus homologues exhibiting similar expression patterns could be identified, which might indicate the functional conservation during evolution. Among 12 types of environmental cues analyzed by the genome-wide microarray data, PtLRR-RLKs showed specific responses to shoot organogenesis, wounding, low ammonium feeding, hypoxia and seasonal dormancy, but not to drought, re-watering after drought, flooding, AlCl₃ treatment and bacteria or fungi treatments.

Conclusions

This study provides the first comprehensive genomic analysis of the Populus LRR-RLK gene family. Segmental duplication contributes significantly to the expansion of this gene family. Populus and Arabidopsis LRR-RLK homologues not only share similar genetic structures but also exhibit comparable expression patterns which point to the possible functional conservation of these LRR-RLKs in two model systems. Transcriptome profiling provides the first insight into the functional divergence among PtLRR-RLK gene subfamilies and suggests that they might take important roles in growth and adaptation of tree species.

Tackling drought stress: receptor-like kinases present new approaches

Global climate change and a growing population require tackling the reduction in arable land and improving biomass production and seed yield per area under varying conditions. One of these conditions is suboptimal water availability. Here, we review some of the classical approaches to dealing with plant response to drought stress and we evaluate how research on RECEPTOR-LIKE KINASES (RLKs) can contribute to improving plant performance under drought stress. RLKs are considered as key regulators of plant architecture and growth behavior, but they also function in defense and stress responses. The available literature and analyses of available transcript profiling data indeed suggest that RLKs can play an important role in optimizing plant responses to drought stress. In addition, RLK pathways are ideal targets for nontransgenic approaches, such as synthetic molecules, providing a novel strategy to manipulate their activity and supporting translational studies from model species, such as Arabidopsis thaliana, to economically useful crops.

Adaptive evolution of Xa21 homologs in Gramineae

The XA21 protein has broad spectrum resistance against Xanthomonas oryzae pv. oryzae. Although Xa21-mediated immunity is well characterized, little is known about the origin and evolutionary history of this gene in grasses. Therefore, we analyzed all Xa21 gene homologs in eight whole-genome sequenced rice lines, as well as in four gramineous genomes, rice, Brachypodium, sorghum and maize; using Arabidopsis Xa21 homologs as outgroups, 17, 7, 7 and 3 Xa21 homologs were detected in these four grasses, respectively. Synteny and phylogenetic analysis showed that frequent gene translocation, duplication and/or loss, have occurred at Xa21 homologous loci, suggesting that they have undergone or are undergoing rapid generation of copy number variations. Within the rice species, the high level of nucleotide diversity between Xa21-like orthologs showed a strong association with the presence/absence haplotypes, suggesting that the genetic structure of rice lines plays an important role in the variations between these Xa21-like orthologs. Strongly positive selection was detected in the core region of the leucine-rich repeat domains of the Xa21 subclade among the rice lines, indicating that the rapid gene diversification of Xa21 homologs may be a strategy for a given species to adapt to the changing spectrum of species-specific pathogens.