Tracking ancestral lineages and recent expansions of NBS-LRR genes in angiosperms

The role of microRNAs in NBS-LRR gene expression and its implications for plant immunity and crop development

Plants evolved, over millions of years, complex defense systems against pathogens. Once infected, the interaction between pathogen effector molecules and host receptors triggers plant immune responses, which include apoptosis, systemic immune response, among others. An important protein family responsible for pathogen effector recognition is the nucleotide binding site-leucine repeat rich (NBS-LRR) proteins. The NBS-LRR gene family is the largest disease resistance gene class in plants. These proteins are widely distributed in vascular plants and have a complex multigenic cluster distribution in plant genomes. To counteract the genetic load of such a large gene family on fitness cost, plants evolved a mechanism using post transcriptional gene silencing induced by small RNAs, particularly microRNAs. For the NBS-LRR gene family, the small RNAs involved in this silencing mechanism are mainly the microRNA482/2118 superfamily. This suppression mechanism is relieved upon pathogen infection, thus allowing increased NBS-LRR expression and triggering plant immunity. In this review, we will discuss the biogenesis of microRNAs and secondary RNAs involved in this silencing mechanism, biochemical and structural features of NBS-LRR proteins in response to pathogen effectors and the evolution of microRNA-based silencing mechanism with a focus on the miR482/2118 family. Furthermore, the biotechnological manipulation of microRNA expression, using both transgenic or genome editing approaches to improve cultivated plants will be discussed, with a focus on the miR482/2118 family in soybean.

QTL mapping of downy mildew resistance in foxtail millet by SLAF‑seq and BSR-seq analysis

KEY MESSAGE

Key message Three major QTLs for resistance to downy mildew were located within an 0.78 Mb interval on chromosome 8 in foxtail millet. Downy mildew, a disease caused by Sclerospora graminicola, is a serious problem that jeopardizes the yield and quality of foxtail millet. Breeding resistant varieties represents one of the most economical and effective solutions, yet there is a lack of molecular markers related to the resistance. Here, a mapping population comprising of 158 F6:7 recombinant inbred lines (RILs) was constructed from the crossing of G1 and JG21. Based on the specific locus amplified fragment sequencing results, a high-density linkage map of foxtail millet with 1031 bin markers, spanning 1041.66 cM was constructed. Based on the high-density linkage map and the phenotype data in four environments, a total of nine quantitative trait loci (QTL) associated with resistance to downy mildew were identified. Further BSR-seq confirmed the genomic regions containing the potential candidate genes related to downy mildew resistance. Interestingly, a 0.78-Mb interval between C8M257 and C8M268 on chromosome 8 was highlighted because of its presence in three major QTL, qDM8_1, qDM8_2, and qDM8_4, which contains 10 NBS-LRR genes. Haplotype analysis in RILs and natural population suggest that 9 SNP loci on Seita8G.199800, Seita8G.195900, Seita8G.198300, and Seita.8G199300 genes were significantly correlated with disease resistance. Furthermore, we found that those genes were taxon-specific by collinearity analysis of pearl millet and foxtail millet genomes. The identification of these new resistance QTL and the prediction of resistance genes against downy mildew will be useful in breeding for resistant varieties and the study of genetic mechanisms of downy mildew disease resistance in foxtail millet.

Plant Immunity: At the Crossroads of Pathogen Perception and Defense Response

Plants are challenged by different microbial pathogens that affect their growth and productivity. However, to defend pathogen attack, plants use diverse immune responses, such as pattern-triggered immunity (PTI), effector-triggered immunity (ETI), RNA silencing and autophagy, which are intricate and regulated by diverse signaling cascades. Pattern-recognition receptors (PRRs) and nucleotide-binding leucine-rich repeat (NLR) receptors are the hallmarks of plant innate immunity because they can detect pathogen or related immunogenic signals and trigger series of immune signaling cascades at different cellular compartments. In plants, most commonly, PRRs are receptor-like kinases (RLKs) and receptor-like proteins (RLPs) that function as a first layer of inducible defense. In this review, we provide an update on how plants sense pathogens, microbe-associated molecular patterns (PAMPs or MAMPs), and effectors as a danger signals and activate different immune responses like PTI and ETI. Further, we discuss the role RNA silencing, autophagy, and systemic acquired resistance as a versatile host defense response against pathogens. We also discuss early biochemical signaling events such as calcium (Ca2+), reactive oxygen species (ROS), and hormones that trigger the activation of different plant immune responses. This review also highlights the impact of climate-driven environmental factors on host-pathogen interactions.

A new comprehensive annotation of leucine-rich repeat-containing receptors in rice

Oryza sativa (rice) plays an essential food security role for more than half of the world's population. Obtaining crops with high levels of disease resistance is a major challenge for breeders, especially today, given the urgent need for agriculture to be more sustainable. Plant resistance genes are mainly encoded by three large leucine-rich repeat (LRR)-containing receptor (LRR-CR) families: the LRR-receptor-like kinase (LRR-RLK), LRR-receptor-like protein (LRR-RLP) and nucleotide-binding LRR receptor (NLR). Using lrrprofiler, a pipeline that we developed to annotate and classify these proteins, we compared three publicly available annotations of the rice Nipponbare reference genome. The extended discrepancies that we observed for LRR-CR gene models led us to perform an in-depth manual curation of their annotations while paying special attention to nonsense mutations. We then transferred this manually curated annotation to Kitaake, a cultivar that is closely related to Nipponbare, using an optimized strategy. Here, we discuss the breakthrough achieved by manual curation when comparing genomes and, in addition to 'functional' and 'structural' annotations, we propose that the community adopts this approach, which we call 'comprehensive' annotation. The resulting data are crucial for further studies on the natural variability and evolution of LRR-CR genes in order to promote their use in breeding future resilient varieties.

Cross Kingdom Immunity: The Role of Immune Receptors and Downstream Signaling in Animal and Plant Cell Death

Both plants and animals are endowed with sophisticated innate immune systems to combat microbial attack. In these multicellular eukaryotes, innate immunity implies the presence of cell surface receptors and intracellular receptors able to detect danger signal referred as damage-associated molecular patterns (DAMPs) and pathogen-associated molecular patterns (PAMPs). Membrane-associated pattern recognition receptors (PRRs), such as Toll-like receptors (TLRs), C-type lectin receptors (CLRs), receptor-like kinases (RLKs), and receptor-like proteins (RLPs) are employed by these organisms for sensing different invasion patterns before triggering antimicrobial defenses that can be associated with a form of regulated cell death. Intracellularly, animals nucleotide-binding and oligomerization domain (NOD)-like receptors or plants nucleotide-binding domain (NBD)-containing leucine rich repeats (NLRs) immune receptors likely detect effectors injected into the host cell by the pathogen to hijack the immune signaling cascade. Interestingly, during the co-evolution between the hosts and their invaders, key cross-kingdom cell death-signaling macromolecular NLR-complexes have been selected, such as the inflammasome in mammals and the recently discovered resistosome in plants. In both cases, a regulated cell death located at the site of infection constitutes a very effective mean for blocking the pathogen spread and protecting the whole organism from invasion. This review aims to describe the immune mechanisms in animals and plants, mainly focusing on cell death signaling pathways, in order to highlight recent advances that could be used on one side or the other to identify the missing signaling elements between the perception of the invasion pattern by immune receptors, the induction of defenses or the transmission of danger signals to other cells. Although knowledge of plant immunity is less advanced, these organisms have certain advantages allowing easier identification of signaling events, regulators and executors of cell death, which could then be exploited directly for crop protection purposes or by analogy for medical research.

Genome- Wide Analysis of the Nucleotide Binding Site Leucine-Rich Repeat Genes of Four Orchids Revealed Extremely Low Numbers of Disease Resistance Genes

Orchids are one of the most diverse flowering plant families, yet possibly maintain the smallest number of the nucleotide-binding site-leucine-rich repeat (NBS-LRR) type plant resistance (R) genes among the angiosperms. In this study, a genome-wide search in four orchid taxa identified 186 NBS-LRR genes. Furthermore, 214 NBS-LRR genes were identified from seven orchid transcriptomes. A phylogenetic analysis recovered 30 ancestral lineages (29 CNL and one RNL), far fewer than other angiosperm families. From the genetics aspect, the relatively low number of ancestral R genes is unlikely to explain the low number of R genes in orchids alone, as historical gene loss and scarce gene duplication has continuously occurred, which also contributes to the low number of R genes. Due to recent sharp expansions, Phalaenopsis equestris and Dendrobium catenatum having 52 and 115 genes, respectively, and exhibited an “early shrinking to recent expanding” evolutionary pattern, while Gastrodia elata and Apostasia shenzhenica both exhibit a “consistently shrinking” evolutionary pattern and have retained only five and 14 NBS-LRR genes, respectively. RNL genes remain in extremely low numbers with only one or two copies per genome. Notably, all of the orchid RNL genes belong to the ADR1 lineage. A separate lineage, NRG1, was entirely absent and was likely lost in the common ancestor of all monocots. All of the TNL genes were absent as well, coincident with the RNL NRG1 lineage, which supports the previously proposed notion that a potential functional association between the TNL and RNL NRG1 genes.

High Resolution Mapping of Rph MBR1012 Conferring Resistance to Puccinia hordei in Barley (Hordeum vulgare L.)

Isolation of disease resistance genes in barley was hampered by the large genome size, but has become easy due to the availability of the reference genome sequence. During the last years, many genomic resources, e.g., the Illumina 9K iSelect, the 50K Infinium arrays, the Barley Genome Zipper, POPSEQ, and genotyping by sequencing (GBS), were developed that enable enhanced gene isolation in combination with the barley genome sequence. In the present study, we developed a fine map of the barley leaf rust resistance gene Rph MBR1012. 537 segmental homozygous recombinant inbred lines (RILs) derived from 4775 F2-plants were used to construct a high-resolution mapping population (HRMP). The Barley Genome Zipper, the 9K iSelect chip, the 50K Infinium chip and GBS were used to develop 56 molecular markers located in the target interval of 8 cM. This interval was narrowed down to about 0.07 cM corresponding to 0.44 Mb of the barley reference genome. Eleven low-confidence and 18 high-confidence genes were identified in this interval. Five of these are putative disease resistance genes and were subjected to allele-specific sequencing. In addition, comparison of the genetic map and the reference genome revealed an inversion of 1.34 Mb located distally to the resistance locus. In conclusion, the barley reference sequence and the respective gene annotation delivered detailed information about the physical size of the target interval, the genes located in the target interval and facilitated the efficient development of molecular markers for marker-assisted selection for RphMBR1012.

Complex regulation of microRNAs in roots of competitively-grown isogenic Nicotiana attenuata plants with different capacities to interact with arbuscular mycorrhizal fungi

BACKGROUND

Nicotiana attenuata is an ecological model plant whose 2.57 Gb genome has recently been sequenced and assembled and for which miRNAs and their genomic locations have been identified. To understand how this plant's miRNAs are reconfigured during plant-arbuscular mycorrhizal fungal (AMF) interactions and whether hostplant calcium- and calmodulin dependent protein kinase (CCaMK) expression which regulates the AMF interaction also modulates miRNAs levels and regulation, we performed a large-scale miRNA analysis of this plant-AMF interaction.

RESULTS

Next generation sequencing of miRNAs in roots of empty vector (EV) N. attenuata plants and an isogenic line silenced in CCaMK expression (irCCaMK) impaired in AMF-interactions grown under competitive conditions with and without AMF inoculum revealed a total of 149 unique miRNAs: 67 conserved and 82 novel ones. The majority of the miRNAs had a length of 21 nucleotides. MiRNA abundances were highly variable ranging from 400 to more than 25,000 reads per million. The miRNA profile of irCCaMK plants impaired in AMF colonization was distinct from fully AMF-functional EV plants grown in the same pot. Six conserved miRNAs were present in all conditions and accumulated differentially depending on treatment and genotype; five (miR6153, miR403a-3p, miR7122a, miR167-5p and miR482d, but not miR399a-3p) showed the highest accumulation in AMF inoculated EV plants compared to inoculated irCCaMK plants. Furthermore, the accumulation patterns of sequence variants of selected conserved miRNAs showed a very distinct pattern related to AMF colonization - one variant of miR473-5p specifically accumulated in AMF-inoculated plants. Also abundances of miR403a-3p, miR171a-3p and one of the sequence variants of miR172a-3p increased in AMF-inoculated EV compared to inoculated irCCaMK plants and to non-inoculated EV plants, while miR399a-3p was most strongly enriched in AMF inoculated irCCaMK plants grown in competition with EV. The analysis of putative targets of selected miRNAs revealed an involvement in P starvation (miR399), phytohormone signaling (Nat-R-PN59, miR172, miR393) and defense (e.g. miR482, miR8667, Nat-R-PN-47).

CONCLUSIONS

Our study demonstrates (1) a large-scale reprograming of miRNAs induced by AMF colonization and (2) that the impaired AMF signaling due to CCaMK silencing and the resulting reduced competitive ability of irCCaMK plants play a role in modulating signal-dependent miRNA accumulation.

Differential expression of NBS-LRR-encoding genes in the root transcriptomes of two Solanum phureja genotypes with contrasting resistance to Globodera rostochiensis

BACKGROUND

The characterization of major resistance genes (R genes) in the potato remains an important task for molecular breeding. However, R genes are rapidly evolving and frequently occur in genomes as clusters with complex structures, and their precise mapping and identification are complicated and time consuming.

RESULTS

Comparative analysis of root transcriptomes of Solanum phureja genotypes with contrasting resistance to Globodera rostochiensis revealed a number of differentially expressed genes. However, compiling a list of candidate R genes for further segregation analysis was hampered by their scarce annotation. Nevertheless, combination of transcriptomic analysis with data on predicted potato NBS-LRR-encoding genes considerably improved the quality of the results and provided a reasonable number of candidate genes that provide S. phureja with strong resistance to the potato golden cyst nematode.

CONCLUSION

Combination of comparative analyses of tissue-specific transcriptomes in resistant and susceptible genotypes may be used as an approach for the rapid identification of candidate potato R genes for co-segregation analysis and may be used in parallel with more sophisticated studies based on genome resequencing.

Divergence and Conservative Evolution of XTNX Genes in Land Plants

The Toll-interleukin-1 receptor (TIR) and Nucleotide-binding site (NBS) domains are two major components of the TIR-NBS-leucine-rich repeat family plant disease resistance genes. Extensive functional and evolutionary studies have been performed on these genes; however, the characterization of a small group of genes that are composed of atypical TIR and NBS domains, namely XTNX genes, is limited. The present study investigated this specific gene family by conducting genome-wide analyses of 59 green plant genomes. A total of 143 XTNX genes were identified in 51 of the 52 land plant genomes, whereas no XTNX gene was detected in any green algae genomes, which indicated that XTNX genes originated upon emergence of land plants. Phylogenetic analysis revealed that the ancestral XTNX gene underwent two rounds of ancient duplications in land plants, which resulted in the formation of clades I/II and clades IIa/IIb successively. Although clades I and IIb have evolved conservatively in angiosperms, the motif composition difference and sequence divergence at the amino acid level suggest that functional divergence may have occurred since the separation of the two clades. In contrast, several features of the clade IIa genes, including the absence in the majority of dicots, the long branches in the tree, the frequent loss of ancestral motifs, and the loss of expression in all detected tissues of Zea mays, all suggest that the genes in this lineage might have undergone pseudogenization. This study highlights that XTNX genes are a gene family originated anciently in land plants and underwent specific conservative pattern in evolution.