R proteins as fundamentals of plant innate immunity

The chloroplast-localized casein kinase II α subunit, CPCK2, negatively regulates plant innate immunity through promoting S-nitrosylation of SABP3

The casein kinase II (CK2) complex consists of catalytic (α) and regulatory (β) subunits and is highly conserved throughout eukaryotes. Plant CK2 plays critical roles in multiple physiological processes; however, its function in plant immunity remains obscure. In this study, we demonstrated that the unique chloroplast-localized CK2 α subunit (CPCK2) is a negative regulator of Arabidopsis thaliana innate immunity. cpck2 mutants displayed enhanced resistance against the fungal pathogen powdery mildew, Golovinomyces cichoracearum and the virulent bacterial pathogen, Pseudomonas syringae pv. tomato (Pto) DC3000. Moreover, the cpck2-1 mutant accumulated higher salicylic acid (SA) levels and mutations that disabled SA biosynthesis or signaling inhibited cpck2-1-mediated disease resistance. CPCK2 interacted with the chloroplast-localized carbonic anhydrase (CA), SA-binding protein 3 (SABP3), which was required for cpck2-mediated immunity. Significantly, CPCK2 phosphorylated SABP3, which promoted S-nitrosylation of this enzyme. It has previously been established that S-nitrosylation of SABP3 reduces both its SA binding function and its CA activity, which compromises the immune-related function of SABP3. Taken together, our results establish CPCK2 as a negative regulator of SA accumulation and associated immunity. Importantly, our findings unveil a mechanism by which CPCK2 negatively regulates plant immunity by promoting S-nitrosylation of SABP3 through phosphorylation, which provides the first example in plants of S-nitrosylation being promoted by cognate phosphorylation.

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.

Alternative Splicing, An Overlooked Defense Frontier of Plants with Respect to Bacterial Infection

Disease represents a major problem in sustainable agricultural development. Plants interact closely with various microorganisms during their development and in response to the prevailing environment. In particular, pathogenic microorganisms can cause plant diseases, affecting the fertility, yield, and longevity of plants. During the long coevolution of plants and their pathogens, plants have evolved both molecular pattern-triggered immunity (PTI) and effector-triggered immunity (ETI) signaling networks in order to regulate host cells in response to pathogen infestation. Additionally, in the postgenomic era, alternative splicing (AS) has become uncovered as one of the major drivers of proteome diversity, and abnormal RNA splicing is closely associated with bacterial infections. Currently, the complexity of host-bacteria interactions is a much studied area of research that has shown steady progress over the past decade. Although the development of high-throughput sequencing technologies and their application in transcriptomes have revolutionized our understanding of AS, many mechanisms related to host-bacteria interactions remain still unclear. To this end, this review summarizes the changes observed in AS during host-bacteria interactions and outlines potential therapeutics for bacterial diseases based on existing studies. In doing so, we hope to provide guidelines for plant disease management in agriculture.

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

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

Microgravity and evasion of plant innate immunity by human bacterial pathogens

Spaceflight microgravity and modeled-microgravity analogs (MMA) broadly alter gene expression and physiology in both pathogens and plants. Research elucidating plant and bacterial responses to normal gravity or microgravity has shown the involvement of both physiological and molecular mechanisms. Under true and simulated microgravity, plants display differential expression of pathogen-defense genes while human bacterial pathogens exhibit increased virulence, antibiotic resistance, stress tolerance, and reduced LD50 in animal hosts. Human bacterial pathogens including Salmonella enterica and E. coli act as cross-kingdom foodborne pathogens by evading and suppressing the innate immunity of plants for colonization of intracellular spaces. It is unknown if evasion and colonization of plants by human pathogens occurs under microgravity and if there is increased infection capability as demonstrated using animal hosts. Understanding the relationship between microgravity, plant immunity, and human pathogens could prevent potentially deadly outbreaks of foodborne disease during spaceflight. This review will summarize (1) alterations to the virulency of human pathogens under microgravity and MMA, (2) alterations to plant physiology and gene expression under microgravity and MMA, (3) suppression and evasion of plant immunity by human pathogens under normal gravity, (4) studies of plant-microbe interactions under microgravity and MMA. A conclusion suggests future study of interactions between plants and human pathogens under microgravity is beneficial to human safety, and an investment in humanity's long and short-term space travel goals.

Molecular insights into the responses of barley to yellow mosaic disease through transcriptome analysis

BACKGROUND

Barley (Hordeum vulgare L.) represents the fourth most essential cereal crop in the world, vulnerable to barley yellow mosaic virus (BaYMV) and/or barley mild mosaic virus (BaMMV), leading to the significant yield reduction. To gain a better understanding of the mechanisms regarding barley crop tolerance to virus infection, we employed a transcriptome sequencing approach and investigated global gene expression among three barley varieties under both infected and control conditions.

RESULTS

High-throughput sequencing outputs revealed massive genetic responses, reflected by the barley transcriptome after BaYMV and/or BaMMV infection. Significant enrichments in peptidase complex and protein processing in endoplasmic reticulum were clustered through Gene ontology and KEGG analysis. Many genes were identified as transcription factors, antioxidants, disease resistance genes and plant hormones and differentially expressed between infected and uninfected barley varieties. Importantly, general response genes, variety-specific and infection-specific genes were also discovered. Our results provide useful information for future barley breeding to resist BaYMV and BaMMV.

CONCLUSIONS

Our study elucidates transcriptomic adaptations in barley response to BaYMV/BaMMV infection through high-throughput sequencing technique. The analysis outcome from GO and KEGG pathways suggests that BaYMV disease induced regulations in multiple molecular-biology processes and signalling pathways. Moreover, critical DEGs involved in defence and stress tolerance mechanisms were displayed. Further functional investigations focusing on these DEGs contributes to understanding the molecular mechanisms of plant response to BaYMV disease infection, thereby offering precious genetic resources for breeding barley varieties resistant to BaYMV disease.

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.

Deep polygenic neural network for predicting and identifying yield-associated genes in Indonesian rice accessions

As the fourth most populous country in the world, Indonesia must increase the annual rice production rate to achieve national food security by 2050. One possible solution comes from the nanoscopic level: a genetic variant called Single Nucleotide Polymorphism (SNP), which can express significant yield-associated genes. The prior benchmark of this study utilized a statistical genetics model where no SNP position information and attention mechanism were involved. Hence, we developed a novel deep polygenic neural network, named the NucleoNet model, to address these obstacles. The NucleoNets were constructed with the combination of prominent components that include positional SNP encoding, the context vector, wide models, Elastic Net, and Shannon’s entropy loss. This polygenic modeling obtained up to 2.779 of Mean Squared Error (MSE) with 47.156% of Symmetric Mean Absolute Percentage Error (SMAPE), while revealing 15 new important SNPs. Furthermore, the NucleoNets reduced the MSE score up to 32.28% compared to the Ordinary Least Squares (OLS) model. Through the ablation study, we learned that the combination of Xavier distribution for weights initialization and Normal distribution for biases initialization sparked more various important SNPs throughout 12 chromosomes. Our findings confirmed that the NucleoNet model was successfully outperformed the OLS model and identified important SNPs to Indonesian rice yields.

Genome wide association study of MAGIC population reveals a novel QTL for salinity and sodicity tolerance in rice

The present study was conducted to identify the novel QTLs controlling salinity and sodicity tolerance using indica MAGIC rice population. Phenotyping was carried out in salinity (EC ~ 10 dS/m) and sodicity (pH ~ 9.8) at the seedling stage. Among 391 lines, 43 and 98 lines were found tolerant and moderately tolerant to salinity. For sodicity condition, 2 and 45 lines were showed tolerance and moderately tolerance at seedling stage. MAGIC population was genotyped with the help of genotyping by sequencing (GBS) and filtered 27041SNPs were used for genome wide marker trait association studies. With respect to salinity tolerance, 25 SNPs were distributed on chromosomes 1, 5, 11 and 12, whereas 18 SNPs were mapped on chromosomes 6, 4 and 11 with LOD value of > 3.25 to sodicity tolerance in rice. The candidate gene analysis detected twelve causal genes including SKC1 gene at Saltol region for salinity and six associated genes for sodic stress tolerance. The significant haplotypes responsible for core histone protein coding gene (LOC_Os12g25120) and three uncharacterized protein coding genes (LOC_Os01g20710, LOC_Os01g20870 and LOC_Os12g22020) were identified under saline stress. Likewise, five significant haplotypes coding for ribose 5-phosphate isomerise (LOC_Os04g24140), aspartyl protease (LOC_Os06g15760), aluminum-activated malate transporter (LOC_Os06g15779), OsFBX421-Fbox domain containing protein (LOC_Os11g32940) and one uncharacterized protein (LOC_Os11g32930) were detected for sodic stress tolerance. The identified novel SNPs could be the potential candidates for functional characterization. These candidate genes aid to further understanding of genetic mechanism on salinity and sodicity stress tolerance in rice. The tolerant line could be used in future breeding programme to enhance the salinity and sodicity tolerance in rice.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-022-01174-8.

RIN4 homologs from important crop species differentially regulate the Arabidopsis NB-LRR immune receptor, RPS2

KEY MESSAGE

RIN4 homologs from important crop species differ in their ability to prevent ectopic activity of the nucleotide binding-leucine rich repeat resistance protein, RPS2. Pathogens deploy virulence effectors to perturb host processes. Plants utilize intracellular resistance (R) proteins to recognize pathogen effectors either by direct interaction or indirectly via effector-mediated perturbations of host components. RPM1-INTERACTING PROTEIN4 (RIN4) is a plant immune regulator that mediates the indirect activation of multiple, independently evolved R-proteins by multiple, unrelated effector proteins. One of these, RPS2 (RESISTANT TO P. SYRINGAE2), is activated upon cleavage of Arabidopsis (At)RIN4 by the Pseudomonas syringae effector AvrRpt2. To gain insight into the AvrRpt2-RIN4-RPS2 defense-activation module, we compared the function of AtRIN4 with RIN4 homologs present in a diverse range of plant species. We selected seven homologs containing conserved features of AtRIN4, including two NOI (Nitrate induced) domains, each containing a predicted cleavage site for AvrRpt2, and a C-terminal palmitoylation site predicted to mediate membrane tethering of the proteins. Palmitoylation-mediated tethering of AtRIN4 to the plasma membrane and cleavage by AvrRpt2 are required for suppression and activation of RPS2, respectively. While all seven homologs are localized at the plasma membrane, only four suppress RPS2 when transiently expressed in Nicotiana benthamiana. All seven homologs are cleaved by AvrRpt2 and, for those homologs that are able to suppress RPS2, cleavage relieves suppression of RPS2. Further, we demonstrate that the membrane-tethered, C-terminal AvrRpt2-generated cleavage fragment is sufficient for the suppression of RPS2. Lastly, we show that the membrane localization of RPS2 is unaffected by its suppression or activation status.

The role of chaperone complex HSP90-SGT1-RAR1 as the associated machinery for hybrid inviability between Nicotiana gossei Domin and N. tabacum L

Two cultured cell lines (GTH4 and GTH4S) of a Nicotiana interspecific F1 hybrid (N. gossei × N. tabacum) were comparatively analyzed to find genetic factors related to hybrid inviability. Both cell lines proliferated at 37 °C, but after shifting to 26 °C, GTH4 started to die similar to the F1 hybrid seedlings, whereas GTH4S survived. As cell death requires de novo expression of genes and proteins, we compared expressed protein profiles between the two cell lines, and found that NgSGT1, a cochaperone of the chaperone complex (HSP90-SGT1-RAR1), was expressed in GTH4 but not in GTH4S. Agrobacterium-mediated transient expression of NgSGT1, but not NtSGT1, induced cell death in leaves of N. tabacum, suggesting its possible role in hybrid inviability. Cell death in N. tabacum was also induced by transient expression of NgRAR1, but not NtRAR1. In contrast, transient expression of any parental combinations of three components revealed that NgRAR1 promoted cell death, whereas NtRAR1 suppressed it in N. tabacum. A specific inhibitor of HSP90, geldanamycin, inhibited the progression of hypersensitive response-like cell death in GTH4 and leaf tissue after agroinfiltration. The present study suggested that components of the chaperone complex are involved in the inviability of Nicotiana interspecific hybrid.

QTL-seq analysis of powdery mildew resistance in a Korean cucumber inbred line

QTL mapping and RT-PCR analyses identified the CsGy5G015660 as a strong powdery mildew resistance candidate gene and natural variation of CsGy5G015660 allele was observed using 115 core germplasm. Powdery mildew (PM) is among the most serious fungal diseases encountered in the cultivation of cucurbits. The development of PM-resistant inbred lines is thus of considerable significance for cucumber breeding programs. In this study, we applied bulked segregant analysis combined with QTL-seq to identify PM resistance loci using F₂ population derived from a cross between two Korean cucumber inbred lines, PM-R (resistant) and PM-S (susceptible). Genome-wide SNP profiling using bulks of the two extreme phenotypes identified two QTLs on chromosomes 5 and 6, designated pm5.2 and pm6.1, respectively. The two PM resistance loci were validated using molecular marker-based classical QTL analysis: pm5.2 (30% R² at LOD 11) and pm6.1 (11% R² at LOD 3.2). Furthermore, reverse transcriptase-PCR analyses, using genes found to be polymorphic between PM-R and PM-S, were conducted to identify the candidate gene(s) responsible for PM resistance. We found that transcripts of the gene CsGy5G015660, encoding a putative leucine-rich repeat receptor-like serine/threonine-protein kinase (RPK2), showed specific accumulation in PM-R prior to the appearance of disease symptoms, and was accordingly considered a strong candidate gene for PM resistance. In addition, cleaved amplified polymorphic sequence markers from CsGy5G015660 were developed and used to screen 35 inbred lines. Natural variation in the CsGy5G015660 allele was also observed based on analysis of a core collection of 115 cucumber accessions. Our results provide new genetic insights for gaining a better understanding of the genetic basis of PM resistance in cucumber, and pave the way for further utilization in cucumber PM resistance breeding programs.

Glycine max NNL1 restricts symbiotic compatibility with widely distributed bradyrhizobia via root hair infection

Symbiosis between soybean (Glycine max) and rhizobia is essential for efficient nitrogen fixation. Rhizobial effectors secreted through the type-III secretion system are key for mediating the interactions between plants and rhizobia, but the molecular mechanism remains largely unknown. Here, our genome-wide association study for nodule number identified G. max Nodule Number Locus 1 (GmNNL1), which encodes a new R protein. GmNNL1 directly interacts with the nodulation outer protein P (NopP) effector from Bradyrhizobium USDA110 to trigger immunity and inhibit nodulation through root hair infection. The insertion of a 179 bp short interspersed nuclear element (SINE)-like transposon into GmNNL1 leads to the loss of function of GmNNL1, enabling bradyrhizobia to successfully nodulate soybeans through the root hair infection route and enhancing nitrogen fixation. Our findings provide important insights into the coevolution of soybean-bradyrhizobia compatibility and offer a way to design new legume-rhizobia interactions for efficient symbiotic nitrogen fixation.

NLRP3 Sensing of Diverse Inflammatory Stimuli Requires Distinct Structural Features

The NLRP3 inflammasome is central to host defense and implicated in various inflammatory diseases and conditions. While the favored paradigm of NLRP3 inflammasome activation stipulates a unifying signal intermediate that de-represses NLRP3, this view has not been tested. Further, structures within NLRP3 required for inflammasome activation are poorly defined. Here we demonstrate that while the NLRP3 LRRs are not auto-repressive and are not required for inflammasome activation by all agonists, distinct sequences within the NLRP3 LRRs positively and negatively modulate inflammasome activation by specific ligands. In addition, elements within the HD1/HD2 “hinge” of NLRP3 and the nucleotide-binding domain have contrasting functions depending upon the specific agonists. Further, while NLRP3 1–432 is minimally sufficient for inflammasome activation by all agonists tested, the pyrin, and linker domains (1–134) function cooperatively and are sufficient for inflammasome activation by certain agonists. Conserved cysteines 8 and 108 appear important for inflammasome activation by sterile, but not infectious insults. Our results define common and agonist-specific regions of NLRP3 that likely mediate ligand-specific responses, discount the hypothesis that NLRP3 inflammasome activation has a unified mechanism, and implicate NLRP3 as an integrator of agonist-specific, inflammasome activating signals.

In Arabidopsis hybrids and Hybrid Mimics, up-regulation of cell wall biogenesis is associated with the increased plant size

Hybrid breeding is of economic importance in agriculture for increasing yield, yet the basis of heterosis is not well understood. In Arabidopsis, crosses between different accessions produce hybrids with different levels of heterosis relative to parental phenotypes in biomass. In all hybrids, the advantage of the F1 hybrid in both phenotypic uniformity and yield gain is lost in the heterogeneous F2. F5/F6 Hybrid Mimics generated from a cross between C24 and Landsberg erecta (Ler) ecotypes demonstrated that the large plant phenotype of the F1 hybrids can be stabilized. Hybrid Mimic selection was applied to Wassilewskija (Ws)/Ler and Col/Ler hybrids. The two hybrids show different levels of heterosis. The Col/Ler hybrid generated F7 Hybrid Mimics with rosette diameter and fresh weight equivalent to the F1 hybrid at 30 DAS; F7 Ws/Ler Hybrid Mimics outperformed the F1 hybrid in both the rosette size and biomass. Transcriptome analysis revealed up-regulation of cell wall biosynthesis, and cell wall expansion genes could be a common pathway in increased size in the Arabidopsis hybrids and Hybrid Mimics. Intercross of two independent Hybrid Mimic lines can further increase the biomass gain. Our results encourage the use of Hybrid Mimics for breeding and for investigating the molecular basis of heterosis.

Comparative Genomic Analysis of Rice with Contrasting Photosynthesis and Grain Production under Salt Stress

Unfavourable environmental conditions, including soil salinity, lead to decreased rice (Oryza sativa L.) productivity, especially at the reproductive stage. In this study, we examined 30 rice varieties, which revealed significant differences in the photosynthetic performance responses under salt stress conditions during the reproductive stage, which ultimately affected yield components after recovery. In rice with a correlation between net photosynthetic rate (P_N) and intercellular CO₂ concentration (C_i) under salt stress, P_N was found to be negatively correlated with filled grain number after recovery. Applying stringent criteria, we identified 130,317 SNPs and 15,396 InDels between two “high-yield rice” varieties and two “low-yield rice” varieties with contrasting photosynthesis and grain yield characteristics. A total of 2089 genes containing high- and moderate-impact SNPs or InDels were evaluated by gene ontology (GO) enrichment analysis, resulting in over-represented terms in the apoptotic process and kinase activity. Among these genes, 262 were highly expressed in reproductive tissues, and most were annotated as receptor-like protein kinases. These findings highlight the importance of variations in signaling components in the genome and these loci can serve as potential genes in rice breeding to produce a variety with salt avoidance that leads to increased yield in saline soil.

Identification and expression profiling analysis of NBS-LRR genes involved in Fusarium oxysporum f.sp. conglutinans resistance in cabbage

As one of the most important resistance (R) gene families in plants, the NBS-LRR genes, encoding proteins with nucleotide-binding site (NBS) and leucine-rich repeat (LRR) domains, play significant roles in resisting pathogens. The published genomic data for cabbage (Brassica oleracea L.) provide valuable data to identify and characterize the genomic organization of cabbage NBS-LRR genes. Ultimately, we identified 105 TIR (N-terminal Toll/interleukin-1 receptor)-NBS-LRR (TNL) genes and 33 CC (coiled-coil)-NBS-LRR (CNL) genes. Further research indicated that 50.7% of the 138 NBS-LRR genes exist in 27 clusters and there are large differences among the gene structures and protein characteristics. Conserved motif and phylogenetic analysis showed that the structures of TNLs and CNLs were similar, with some differences. These NBS-LRRs are evolved under negative selection and mostly arose from whole-genome duplication events during evolution. Tissue-expression profiling of NBS-LRR genes revealed that 37.1% of the TNL genes are highly or specifically expressed in roots, especially the genes on chromosome 7 (76.5%). Digital gene expression and reverse transcription PCR analyses revealed the expression patterns of the NBS-LRR genes upon challenge by Fusarium oxysporum f.sp. conglutinans: nine genes were upregulated, and five were downregulated. The major resistance gene Foc1 probably works together with the other four genes in the same cluster to resist F. oxysporum infection.

PK-profiling method for identifying the expression of resistance-associated genes in partially resistant oats to crown rust

BACKGROUND

Protein kinases play a key role in plant cell homeostasis and the activation of defense mechanisms. Partial resistance to fungi in plants is interesting because of its durability. However, the variable number of minor loci associated with this type of resistance hampers the reliable identification of the full range of genes involved. The present work reports the technique of protein kinase (PK)-profiling for the identification of the PK genes induced in the partially resistant oats line MN841801-1 following exposure to the fungus Puccinia coronata. This is the first time this technique has been used with cDNA (complementary DNA) from a suppression subtractive hybridization library obtained after the hybridization of cDNAs from inoculated and mock-inoculated plants.

RESULTS

Six degenerate primers based on the conserved domains of protein kinases were used in a PK-profiling assay including cDNA from mock-inoculated leaves and subtracted cDNA. Of the 75.7% of sequences cloned and sequenced that showed significant similarity to resistance genes, 76% were found to code for PKs. Translation and ClustalW2 alignment of each sequence cloned with the complete sequences of the most similar B. distachyon PKs allowed those of the partially resistant oat line to be deduced and characterized. Further, a phylogenetic study carried out after alignment of these B. distachyon PK sequences with the most similar protein sequences of related species also allowed to deduce different functions for the PK cloned. RT-qPCR (Reverse Transcription-quantitative PCR) was analyzed on nine representative sequences to validate the reliability of the employed PK-profiling method as a tool for identifying the expression of resistance-associated genes.

CONCLUSIONS

PK-profiling would appear to be a useful tool for the identification of the PKs expressed in oats after challenge by P. coronata, and perhaps other pathogens. Most of the PKs studied are related to receptor-like protein kinases expressed shortly after infection. This is in agreement with previous studies indicating a close relationship between partial resistance and the first layer of defense against pathogen used by plants.

Development of transcriptome based web genomic resources of yellow mosaic disease in Vigna mungo

Vigna mungo (Urdbean) is cultivated in the tropical and sub-tropical continental region of Asia. It is not only important source of dietary protein and nutritional elements, but also of immense value to human health due to medicinal properties. Yellow mosaic disease caused by Mungbean Yellow Mosaic India Virus is known to incur huge loss to crop, adversely affecting crop yield. Contrasting genotypes are ideal source for knowledge discovery of plant defence mechanism and associated candidate genes for varietal improvement. Whole genome sequence of this crop is yet to be completed. Moreover, genomic resources are also not freely accessible, thus available transcriptome data can be of immense use. V. mungo Transcriptome database, accessible at http://webtom.cabgrid.res.in/vmtdb/ has been developed using available data of two contrasting varieties viz., cv. VM84 (resistant) and cv. T9 (susceptible). De novo assembly was carried out using Trinity and CAP3. Out of total 240,945 unigenes, 165,894 (68.8%) showed similarity with known genes against NR database, and remaining 31.2% were found to be novel. We found 22,101 differentially expressed genes in all datasets, 44,335 putative genic SSR markers, 4105 SNPs and Indels, 64,964 transcriptional factor, 546 mature miRNA target prediction in 703 differentially expressed unigenes and 137 pathways. MAPK, salicylic acid-binding protein 2-like, pathogenesis-related protein and NBS-LRR domain were found which may play an important role in defence against pathogens. This is the first web genomic resource of V. mungo for future genome annotation as well as ready to use markers for future variety improvement program.

Characterization of CNL like protein fragment (CNL-LPF) from mature Lageneria siceraria seeds

Coiled coil domain-nucleotide binding site-leucine rich repeat (CC-NBS-LRR; CNL) proteins are highly conserved family of plant disease resistance proteins, remarkably comprise of coiled-coil domain, which plays significant role in plant innate immunity. The present study reports that moderately elicited oligomerization of plant CNL like protein fragment (CNL-LPF) in presence of ATP/Mg using various biophysical methods Circular dichroism (CD) results depicted a substantial increase in β-sheet structure content of CNL-LPF. ATP/Mg induced conformational change in protein was observed by increase in blue shift with extrinsic fluorescence measurement, which indicates the exposure of hydrophobic regions of CNL-LPF and leads to self-association i.e. oligomerization. Likewise, cluster of protein oligomer and alteration in protein surface morphology were observed in presence of ATP/Mg by Transmission electron microscopy (TEM) and Atomic force microscopy (AFM), respectively. Also, augmented antiproliferation of HT1376 cells (urinary bladder cancer cell lines) was observed by CNL-LPF in presence of ATP/Mg. In conclusion, the current study illustrates that extent of CNL-LPF oligomerization was enhanced in presence of ATP/Mg (as compared to its absence). Utilization of enhanced oligomerization property of CNL-LPF as an anti-proliferative agent needs more assessment.

Genome-wide identification and characterization of NB-ARC resistant genes in wheat (Triticum aestivum L.) and their expression during leaf rust infection

NB-ARC domain-containing resistance genes from the wheat genome were identified, characterized and localized on chromosome arms that displayed differential yet positive response during incompatible and compatible leaf rust interactions. Wheat (Triticum aestivum L.) is an important cereal crop; however, its production is affected severely by numerous diseases including rusts. An efficient, cost-effective and ecologically viable approach to control pathogens is through host resistance. In wheat, high numbers of resistance loci are present but only few have been identified and cloned. A comprehensive analysis of the NB-ARC-containing genes in complete wheat genome was accomplished in this study. Complete NB-ARC encoding genes were mined from the Ensembl Plants database to predict 604 NB-ARC containing sequences using the HMM approach. Genome-wide analysis of orthologous clusters in the NB-ARC-containing sequences of wheat and other members of the Poaceae family revealed maximum homology with Oryza sativa indica and Brachypodium distachyon. The identification of overlap between orthologous clusters enabled the elucidation of the function and evolution of resistance proteins. The distributions of the NB-ARC domain-containing sequences were found to be balanced among the three wheat sub-genomes. Wheat chromosome arms 4AL and 7BL had the most NB-ARC domain-containing contigs. The spatio-temporal expression profiling studies exemplified the positive role of these genes in resistant and susceptible wheat plants during incompatible and compatible interaction in response to the leaf rust pathogen Puccinia triticina. Two NB-ARC domain-containing sequences were modelled in silico, cloned and sequenced to analyze their fine structures. The data obtained in this study will augment isolation, characterization and application NB-ARC resistance genes in marker-assisted selection based breeding programs for improving rust resistance in wheat.

Combined de novo and genome guided assembly and annotation of the Pinus patula juvenile shoot transcriptome

Background

Pines are the most important tree species to the international forestry industry, covering 42 % of the global industrial forest plantation area. One of the most pressing threats to cultivation of some pine species is the pitch canker fungus, Fusarium circinatum, which can have devastating effects in both the field and nursery. Investigation of the Pinus-F. circinatum host-pathogen interaction is crucial for development of effective disease management strategies. As with many non-model organisms, investigation of host-pathogen interactions in pine species is hampered by limited genomic resources. This was partially alleviated through release of the 22 Gbp Pinus taeda v1.01 genome sequence (http://pinegenome.org/pinerefseq/) in 2014. Despite the fact that the fragmented state of the genome may hamper comprehensive transcriptome analysis, it is possible to leverage the inherent redundancy resulting from deep RNA sequencing with Illumina short reads to assemble transcripts in the absence of a completed reference sequence. These data can then be integrated with available genomic data to produce a comprehensive transcriptome resource. The aim of this study was to provide a foundation for gene expression analysis of disease response mechanisms in Pinus patula through transcriptome assembly.

Results

Eighteen de novo and two reference based assemblies were produced for P. patula shoot tissue. For this purpose three transcriptome assemblers, Trinity, Velvet/OASES and SOAPdenovo-Trans, were used to maximise diversity and completeness of assembled transcripts. Redundancy in the assembly was reduced using the EvidentialGene pipeline. The resulting 52 Mb P. patula v1.0 shoot transcriptome consists of 52 112 unigenes, 60 % of which could be functionally annotated.

Conclusions

The assembled transcriptome will serve as a major genomic resource for future investigation of P. patula and represents the largest gene catalogue produced to date for this species. Furthermore, this assembly can help detect gene-based genetic markers for P. patula and the comparative assembly workflow could be applied to generate similar resources for other non-model species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2277-7) contains supplementary material, which is available to authorized users.

Transcriptomic Analysis and the Expression of Disease-Resistant Genes in Oryza meyeriana under Native Condition

Oryza meyeriana (O. meyeriana), with a GG genome type (2n = 24), accumulated plentiful excellent characteristics with respect to resistance to many diseases such as rice shade and blast, even immunity to bacterial blight. It is very important to know if the diseases-resistant genes exist and express in this wild rice under native conditions. However, limited genomic or transcriptomic data of O. meyeriana are currently available. In this study, we present the first comprehensive characterization of the O. meyeriana transcriptome using RNA-seq and obtained 185,323 contigs with an average length of 1,692 bp and an N50 of 2,391 bp. Through differential expression analysis, it was found that there were most tissue-specifically expressed genes in roots, and next to stems and leaves. By similarity search against protein databases, 146,450 had at least a significant alignment to existed gene models. Comparison with the Oryza sativa (japonica-type Nipponbare and indica-type 93–11) genomes revealed that 13% of the O. meyeriana contigs had not been detected in O. sativa. Many diseases-resistant genes, such as bacterial blight resistant, blast resistant, rust resistant, fusarium resistant, cyst nematode resistant and downy mildew gene, were mined from the transcriptomic database. There are two kinds of rice bacterial blight-resistant genes (Xa1 and Xa26) differentially or specifically expressed in O. meyeriana. The 4 Xa1 contigs were all only expressed in root, while three of Xa26 contigs have the highest expression level in leaves, two of Xa26 contigs have the highest expression profile in stems and one of Xa26 contigs was expressed dominantly in roots. The transcriptomic database of O. meyeriana has been constructed and many diseases-resistant genes were found to express under native condition, which provides a foundation for future discovery of a number of novel genes and provides a basis for studying the molecular mechanisms associated with disease resistance in O. meyeriana.

Arabidopsis BNT1, an atypical TIR-NBS-LRR gene, acting as a regulator of the hormonal response to stress

During their life cycle, plants have to cope with fluctuating environmental conditions. The perception of the stressful environmental conditions induces a specific stress hormone signature specifying a proper response with an efficient fitness. By reverse genetics, we isolated and characterized a novel mutation in Arabidopsis, associated with environmental stress responses, that affects the At5g11250/BURNOUT1 (BNT1) gene which encode a Toll/Interleukin1 receptor-nucleotide binding site leucine-rich repeat (TIR-NBS-LRR) protein. The knock-out bnt1 mutants displayed, in the absence of stress conditions, a multitude of growth and development defects, suchas severe dwarfism, early senescence and flower sterility, similar to those observed in vitro in wild type plants upon different biotic and/or abiotic stresses. The disruption of BNT1 causes also a drastic increase of the jasmonic, salicylic and abscisic acids as well as ethylene levels. Which was consistent with the expression pattern observed in bnt1 showing an over representation of genes involved in the hormonal response to stress? Therefore, a defect in BNT1 forced the plant to engage in an exhausting general stress response, which produced frail, weakened and poorly adapted plants expressing "burnout" syndromes. Furthermore, by in vitro phenocopying experiments, physiological, chemical and molecular analyses, we propose that BNT1 could represent a molecular link between stress perception and specific hormonal signature.

Dissection of Two Complex Clusters of Resistance Genes in Lettuce (Lactuca sativa)

Of the over 50 phenotypic resistance genes mapped in lettuce, 25 colocalize to three major resistance clusters (MRC) on chromosomes 1, 2, and 4. Similarly, the majority of candidate resistance genes encoding nucleotide binding-leucine rich repeat (NLR) proteins genetically colocalize with phenotypic resistance loci. MRC1 and MRC4 span over 66 and 63 Mb containing 84 and 21 NLR-encoding genes, respectively, as well as 765 and 627 genes that are not related to NLR genes. Forward and reverse genetic approaches were applied to dissect MRC1 and MRC4. Transgenic lines exhibiting silencing were selected using silencing of β-glucuronidase as a reporter. Silencing of two of five NLR-encoding gene families resulted in abrogation of nine of 14 tested resistance phenotypes mapping to these two regions. At MRC1, members of the coiled coil-NLR-encoding RGC1 gene family were implicated in host and nonhost resistance through requirement for Dm5/8- and Dm45-mediated resistance to downy mildew caused by Bremia lactucae as well as the hypersensitive response to effectors AvrB, AvrRpm1, and AvrRpt2 of the nonpathogen Pseudomonas syringae. At MRC4, RGC12 family members, which encode toll interleukin receptor-NLR proteins, were implicated in Dm4-, Dm7-, Dm11-, and Dm44-mediated resistance to B. lactucae. Lesions were identified in the sequence of a candidate gene within dm7 loss-of-resistance mutant lines, confirming that RGC12G confers Dm7.

Chitosan-induced antiviral activity and innate immunity in plants

Immunity represents a trait common to all living organisms, and animals and plants share some similarities. Therefore, in susceptible host plants, complex defence machinery may be stimulated by elicitors. Among these, chitosan deserves particular attention because of its proved efficacy. This survey deals with the antiviral activity of chitosan, focusing on its perception by the plant cell and mechanism of action. Emphasis has been paid to benefits and limitations of this strategy in crop protection, as well as to the potential of chitosan as a promising agent in virus disease control.

Integrated signaling networks in plant responses to sedentary endoparasitic nematodes: a perspective

Sedentary plant endoparasitic nematodes can cause detrimental yield losses in crop plants making the study of detailed cellular, molecular, and whole plant responses to them a subject of importance. In response to invading nematodes and nematode-secreted effectors, plant susceptibility/resistance is mainly determined by the coordination of different signaling pathways including specific plant resistance genes or proteins, plant hormone synthesis and signaling pathways, as well as reactive oxygen signals that are generated in response to nematode attack. Crosstalk between various nematode resistance-related elements can be seen as an integrated signaling network regulated by transcription factors and small RNAs at the transcriptional, posttranscriptional, and/or translational levels. Ultimately, the outcome of this highly controlled signaling network determines the host plant susceptibility/resistance to nematodes.

Characterization of two distinct subfamilies of SUN-domain proteins in Arabidopsis and their interactions with the novel KASH-domain protein AtTIK

SUN-domain proteins belong to a gene family including classical Cter-SUN and mid-SUN subfamilies differentiated by the position of the SUN domain within the protein. Although present in animal and plant species, mid-SUN proteins have so far remained poorly described. Here, we used a combination of genetics, yeast two-hybrid and in planta transient expression methods to better characterize the SUN family in Arabidopsis thaliana. First, we validated the mid-SUN protein subfamily as a monophyletic group conserved from yeast to plant. Arabidopsis Cter-SUN (AtSUN1 and AtSUN2) and mid-SUN (AtSUN3 and AtSUN4) proteins expressed as fluorescent protein fusions are membrane-associated and localize to the nuclear envelope (NE) and endoplasmic reticulum. However, only the Cter-SUN subfamily is enriched at the NE. We investigated interactions in and between members of the two subfamilies and identified the coiled-coil domain as necessary for mediating interactions. The functional significance of the mid-SUN subfamily was further confirmed in mutant plants as essential for early seed development and involved in nuclear morphology. Finally, we demonstrated that both subfamilies interact with the KASH domain of AtWIP1 and identified a new root-specific KASH-domain protein, AtTIK. AtTIK localizes to the NE and affects nuclear morphology. Our study indicates that Arabidopsis Cter-SUN and mid-SUN proteins are involved in a complex protein network at the nuclear membranes, reminiscent of the LInker of Nucleoskeleton and Cytoskeleton (LINC) complex found in other kingdoms.

Full-genome identification and characterization of NBS-encoding disease resistance genes in wheat

Host resistance is the most economical, effective and ecologically sustainable method of controlling diseases in crop plants. In bread wheat, despite the high number of resistance loci that have been cataloged to date, only few have been cloned, underlying the need for genomics-guided investigations capable of providing a prompt and acute knowledge on the identity of effective resistance genes that can be used in breeding programs. Proteins with a nucleotide-binding site (NBS) encoded by the major plant disease resistance (R) genes play an important role in the responses of plants to various pathogens. In this study, a comprehensive analysis of NBS-encoding genes within the whole wheat genome was performed, and the genome scale characterization of this gene family was established. From the recently published wheat genome sequence, we used a data mining and automatic prediction pipeline to identify 580 complete ORF candidate NBS-encoding genes and 1,099 partial-ORF ones. Among complete gene models, 464 were longer than 200 aa, among them 436 had less than 70 % of sequence identity to each other. This gene models set was deeply characterized. (1) First, we have analyzed domain architecture and identified, in addition to typical domain combinations, the presence of particular domains like signal peptides, zinc fingers, kinases, heavy-metal-associated and WRKY DNA-binding domains. (2) Functional and expression annotation via homology searches in protein and transcript databases, based on sufficient criteria, enabled identifying similar proteins for 60 % of the studied gene models and expression evidence for 13 % of them. (3) Shared orthologous groups were defined using NBS-domain proteins of rice and Brachypodium distachyon. (4) Finally, alignment of the 436 NBS-containing gene models to the full set of scaffolds from the IWGSC's wheat chromosome survey sequence enabled high-stringence anchoring to chromosome arms. The distribution of the R genes was found balanced on the three wheat sub-genomes. In contrast, at chromosome scale, 50 % of members of this gene family were localized on 6 of the 21 wheat chromosomes and ~22 % of them were localized on homeologous group 7. The results of this study provide a detailed analysis of the largest family of plant disease resistance genes in allohexaploid wheat. Some structural traits reported had not been previously identified and the genome-derived data were confronted with those stored in databases outlining the functional specialization of members of this family. The large reservoir of NBS-type genes presented and discussed will, firstly, form an important framework for marker-assisted improvement of resistance in wheat, and, secondly, open up new perspectives for a better understanding of the evolution dynamics of this gene family in grass species and in polyploid systems.

Characterization of novel wheat NBS domain-containing sequences and their utilization, in silico, for genome-scale R-gene mining

In crop improvement, the isolation, cloning and transfer of disease resistance genes (R-genes) is an ultimate goal usually starting from tentative R-gene analogs (RGAs) that are identified on the basis of their structure. For bread wheat, recent advances in genome sequencing are supporting the efforts of wheat geneticists worldwide. Among wheat R-genes, nucleotide-binding site (NBS)-encoding ones represent a major class. In this study, we have used a polymerase chain reaction-based approach to amplify and clone NBS-type RGAs from a bread wheat cultivar, 'Salambo 80.' Four novel complete ORF sequences showing similarities to previously reported R-genes/RGAs were used for in silico analyses. In a first step, where analyses were focused on the NBS domain, these sequences were phylogenetically assigned to two distinct groups: a first group close to leaf rust Lr21 resistance proteins; and a second one similar to cyst nematode resistance proteins. In a second step, sequences were used as initial seeds to walk up and downstream the NBS domain. This procedure enabled identifying 8 loci ranging in size between 2,115 and 7,653 bp. Ab initio gene prediction identified 8 gene models, among which two had complete ORFs. While GenBank survey confirmed the belonging of sequences to two groups, subsequent characterization using IWGSC genomic and proteomic data showed that the 8 gene models, reported in this study, were unique and their loci matched scaffolds on chromosome arms 1AS, 1BS, 4BS and 1DS. The gene model located on 1DS is a pseudo-Lr21 that was shown to have an NBS-LRR domain structure, while the potential association of the RGAs, here reported, is discussed. This study has produced novel R-gene-like loci and models in the wheat genome and provides the first steps toward further elucidation of their role in wheat disease resistance.

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.

Decoding the massive genome of loblolly pine using haploid DNA and novel assembly strategies

Background

The size and complexity of conifer genomes has, until now, prevented full genome sequencing and assembly. The large research community and economic importance of loblolly pine, Pinus taeda L., made it an early candidate for reference sequence determination.

Results

We develop a novel strategy to sequence the genome of loblolly pine that combines unique aspects of pine reproductive biology and genome assembly methodology. We use a whole genome shotgun approach relying primarily on next generation sequence generated from a single haploid seed megagametophyte from a loblolly pine tree, 20-1010, that has been used in industrial forest tree breeding. The resulting sequence and assembly was used to generate a draft genome spanning 23.2 Gbp and containing 20.1 Gbp with an N50 scaffold size of 66.9 kbp, making it a significant improvement over available conifer genomes. The long scaffold lengths allow the annotation of 50,172 gene models with intron lengths averaging over 2.7 kbp and sometimes exceeding 100 kbp in length. Analysis of orthologous gene sets identifies gene families that may be unique to conifers. We further characterize and expand the existing repeat library based on the de novo analysis of the repetitive content, estimated to encompass 82% of the genome.

Conclusions

In addition to its value as a resource for researchers and breeders, the loblolly pine genome sequence and assembly reported here demonstrates a novel approach to sequencing the large and complex genomes of this important group of plants that can now be widely applied.

Proteomics changes during the incompatible interaction between cowpea and Colletotrichum gloeosporioides (Penz.) Penz and Sacc

Anthracnose represents an important disease of cowpea [Vigna unguiculata L. (Walp.)] caused by the hemibiothrophic fungus Colletotrichum gloeosporioides that drastically reduces cowpea field production. In this study we investigated some biochemical aspects underlying the incompatible interaction between a resistant cowpea genotype and C. gloeosporioides using a proteomic approach. Analyses of two-dimensional gel electrophoresis patterns and protein identification indicate C. gloeosporioides infection-dependent cowpea leaf proteome changes associated with metabolism, photosynthesis, response to stress, oxidative burst and scavenging, defense signaling, and pathogenesis-related proteins. Moreover the C. gloeosporioides responsive proteins interaction network in cowpea revealed the interconnected modulation of key cellular processes involving particularly antioxidants proteins, photosynthetic apparatus forming proteins and proteins of the energetic metabolism that interact with each other suggesting that their expression changes are also important for resistance of cowpea to C. gloeosporioides.

Callose-mediated resistance to pathogenic intruders in plant defense-related papillae

Plants are exposed to a wide range of potential pathogens, which derive from diverse phyla. Therefore, plants have developed successful defense mechanisms during co-evolution with different pathogens. Besides many specialized defense mechanisms, the plant cell wall represents a first line of defense. It is actively reinforced through the deposition of cell wall appositions, so-called papillae, at sites of interaction with intruding microbial pathogens. The papilla is a complex structure that is formed between the plasma membrane and the inside of the plant cell wall. Even though the specific biochemical composition of papillae can vary between different plant species, some classes of compounds are commonly found which include phenolics, reactive oxygen species, cell wall proteins, and cell wall polymers. Among these polymers, the (1,3)-β-glucan callose is one of the most abundant and ubiquitous components. Whereas the function of most compounds could be directly linked with cell wall reinforcement or an anti-microbial effect, the role of callose has remained unclear. An evaluation of recent studies revealed that the timing of the different papilla-forming transport processes is a key factor for successful plant defense.

VIGS technology: an attractive tool for functional genomics studies in legumes

Legume species are among the most important crops worldwide. In recent years, six legume genomes have been completely sequenced, and there is now an urgent need for reverse-genetics tools to validate genes affecting yield and product quality. As most legumes are recalcitrant to stable genetic transformation, virus-induced gene silencing (VIGS) appears to be a powerful alternative technology for determining the function of unknown genes. VIGS technology is based on the property of plant viruses to trigger a defence mechanism related to post-transcriptional gene silencing (PTGS). Infection by a recombinant virus carrying a fragment of a plant target gene will induce homology-dependent silencing of the endogenous target gene. Several VIGS systems have been developed for legume species since 2004, including those based on Bean pod mottle virus, Pea early browning virus, and Apple latent spherical virus, and used in reverse-genetics studies of a wide variety of plant biological processes. In this work, we give an overview of the VIGS systems available for legumes, and present their successful applications in functional genomics studies. We also discuss the limitations of these VIGS systems and the future challenges to be faced in order to use VIGS to its full potential in legume species.

Identification of microRNAs potentially involved in male sterility of Brassica campestris ssp. chinensis using microRNA array and quantitative RT-PCR assays

microRNAs (miRNAs) are a class of newly identified, noncoding, small RNA molecules that negatively regulate gene expression. Many miRNAs are reportedly involved in plant growth, development and stress response processes. However, their roles in the sexual reproduction mechanisms in flowering plants remain unknown. Pollen development is an important process in the life cycle of a flowering plant, and it is closely related to the yield and quality of crop seeds. This study aimed to identify miRNAs involved in pollen development. A microarray assay was conducted using the known complementary sequences of plant miRNAs as probes on inflorescences of a sterile male line (Bcajh97-01A) and a fertile male line (Bcajh97-01B) of the Brassica campestris ssp. chinensis cv. 'Aijiaohuang' genic male sterility sister line system (Bcajh97-01A/B). The results showed that 44 miRNAs were differently expressed in the two lines. Of these, 15 had over 1.5-fold changes in their transcript levels, with 9 upregulated and 6 downregulated miRNAs in inflorescences of 'Bcajh97-01A' sterile line plants. We then focused on 3 of these 15 miRNAs (miR158, miR168 and miR172). Through computational methods, 13 family members were predicted for these 3 miRNAs and 22 genes were predicted to be their candidate target genes. By using 5' modified RACE, 2 target genes of miR168 and 5 target genes of miR172 were identified. Then, qRT-PCR was applied to verify the existence and expression patterns of the 3 miRNAs in the flower buds at five developmental stages. The results were generally consistent with those of the microarray. Thus, this study may give a valuable clue for further exploring the miRNA group that may function during pollen development.

Haplotype variability and identification of new functional alleles at the Rdg2a leaf stripe resistance gene locus

The barley Rdg2a locus confers resistance to the leaf stripe pathogen Pyrenophora graminea and, in the barley genotype Thibaut, it is composed of a gene family with three highly similar paralogs. Only one member of the gene family (called as Rdg2a) encoding for a CC-NB-LRR protein is able to confer resistance to the leaf stripe isolate Dg2. To study the genome evolution and diversity at the Rdg2a locus, sequences spanning the Rdg2a gene were compared in two barley cultivars, Thibaut and Morex, respectively, resistant and susceptible to leaf stripe. An overall high level of sequence conservation interrupted by several rearrangements that included three main deletions was observed in the Morex contig. The main deletion of 13,692 bp was most likely derived from unequal crossing over between Rdg2a paralogs leading to the generation of a chimeric Morex rdg2a gene which was not associated to detectable level of resistance toward leaf stripe. PCR-based analyses of genic and intergenic regions at the Rdg2a locus in 29 H. vulgare lines and one H. vulgare ssp. spontaneum accession indicated large haplotype variability in the cultivated barley gene pool suggesting rapid and recent divergence at this locus. Barley genotypes showing the same haplotype as Thibaut at the Rdg2a locus were selected for a Rdg2a allele mining through allele re-sequencing and two lines with polymorphic nucleotides leading to amino acid changes in the CC-NB and LRR encoding domains, respectively, were identified. Analysis of nucleotide diversity of the Rdg2a alleles revealed that the polymorphic sites were subjected to positive selection. Moreover, strong positively selected sites were located in the LRR encoding domain suggesting that both positive selection and divergence at homologous loci are possibly representing the molecular mechanism for the generation of high diversity at the Rdg2a locus in the barley gene pool.

Plant nucleotide binding site-leucine-rich repeat (NBS-LRR) genes: active guardians in host defense responses

The most represented group of resistance genes are those of the nucleotide binding site-leucine-rich repeat (NBS-LRR) class. These genes are very numerous in the plant genome, and they often occur in clusters at specific loci following gene duplication and amplification events. To date, hundreds of resistance genes and relatively few quantitative trait loci for plant resistance to pathogens have been mapped in different species, with some also cloned. When these NBS-LRR genes have been physically or genetically mapped, many cases have shown co-localization between resistance loci and NBS-LRR genes. This has allowed the identification of candidate genes for resistance, and the development of molecular markers linked to R genes. This review is focused on recent genomics studies that have described the abundance, distribution and evolution of NBS-LRR genes in plant genomes. Furthermore, in terms of their expression, NBS-LRR genes are under fine regulation by cis- and trans-acting elements. Recent findings have provided insights into the roles of alternative splicing, the ubiquitin/ proteasome system, and miRNAs and secondary siRNAs in the regulation of NBS-LRR gene expression at the post-transcriptional, post-translational and epigenetic levels. The possibility to use this knowledge for genetic improvement of plant resistance to pathogens is discussed.

The AvrB_AvrC domain of AvrXccC of Xanthomonas campestris pv. campestris is required to elicit plant defense responses and manipulate ABA homeostasis

Plant disease induced by Xanthomonas campestris pv. campestris depends on type III effectors but the molecular basis is poorly understood. Here, AvrXccC8004 was characterized, and it was found that the AvrB_AvrC domain was essential and sufficient to elicit defense responses in an Arabidopsis-resistant ecotype (Col-0). An upregulation of genes in responding to the AvrB_AvrC domain of AvrXccC8004 was shown in a profile of host gene expression. The molecular changes were correlated with morphological changes observed in phenotypic and ultrastructural characterizations. Interestingly, the abscisic acid (ABA)-signaling pathway was also a prominent target for the AvrB_AvrC domain of AvrXccC8004. The highly elicited NCED5, encoding a key enzyme of ABA biosynthesis, was increased in parallel with ABA levels in AvrXccC8004 transgenic plants. Consistently, the X. campestris pv. campestris 8004 ΔavrXccC mutant was severely impaired in the ability to manipulate the accumulation of ABA and induction of ABA-related genes in challenged leaves. Moreover, exogenous application of ABA also enhanced the susceptibility of Arabidopsis to the X. campestris pv. campestris strains. These results indicate that the AvrB_AvrC domain of AvrXccC8004 alone has the activity to manipulate ABA homeostasis, which plays an important role in regulating the interactions of X. campestris pv. campestris and Arabidopsis.

Comparative analysis of zinc finger proteins involved in plant disease resistance

A meta-analysis was performed to understand the role of zinc finger domains in proteins of resistance (R) genes cloned from different crops. We analyzed protein sequences of seventy R genes of various crops in which twenty six proteins were found to have zinc finger domains along with nucleotide binding sites - leucine rice repeats (NBS-LRR) domains. We identified thirty four zinc finger domains in the R proteins of nine crops and were grouped into 19 types of zinc fingers. The size of individual zinc finger domain within the R genes varied from 11 to 84 amino acids, whereas the size of proteins containing these domains varied from 263 to 1305 amino acids. The biophysical analysis revealed that molecular weight of Pi54 zinc finger was lowest whereas the highest one was found in rice Pib zinc finger named as Transposes Transcription Factor (TTF). The instability (R(2) =0.95) and the aliphatic (R(2) =0.94) indices profile of zinc finger domains follows the polynomial distribution pattern. The pairwise identity analysis showed that the Lin11, Isl-1 & Mec-3 (LIM) zinc finger domain of rice blast resistance protein pi21 have 12.3% similarity with the nuclear transcription factor, X-box binding-like 1 (NFX) type zinc finger domain of Pi54 protein. For the first time, we reported that Pi54 (Pi-k(h)-Tetep), a rice blast resistance (R) protein have a small zinc finger domain of NFX type located on the C-terminal in between NBS and LRR domains of the R-protein. Compositional analysis depicted by the helical wheel diagram revealed the presence of a hydrophobic region within this domain which might help in exposing the LRR region for a possible R-Avr interaction. This domain is unique among all other cloned plant disease resistance genes and might play an important role in broad-spectrum nature of rice blast resistance gene Pi54.

Mechanistic insight into pentatricopeptide repeat proteins as sequence-specific RNA-binding proteins for organellar RNAs in plants

The pentatricopeptide repeat (PPR) protein family is highly expanded in terrestrial plants. Arabidopsis contains 450 PPR genes, which represents 2% of the total protein-coding genes. PPR proteins are eukaryote-specific RNA-binding proteins implicated in multiple aspects of RNA metabolism of organellar genes. Most PPR proteins affect a single or small subset of gene(s), acting in a gene-specific manner. Studies over the last 10 years have revealed the significance of this protein family in coordinated gene expression in different compartments: the nucleus, chloroplast and mitochondrion. Here, we summarize recent studies addressing the mechanistic aspect of PPR proteins.

Genome Wide Analysis of Nucleotide-Binding Site Disease Resistance Genes in Brachypodium distachyon

Nucleotide-binding site (NBS) disease resistance genes play an important role in defending plants from a variety of pathogens and insect pests. Many R-genes have been identified in various plant species. However, little is known about the NBS-encoding genes in Brachypodium distachyon. In this study, using computational analysis of the B. distachyon genome, we identified 126 regular NBS-encoding genes and characterized them on the bases of structural diversity, conserved protein motifs, chromosomal locations, gene duplications, promoter region, and phylogenetic relationships. EST hits and full-length cDNA sequences (from Brachypodium database) of 126 R-like candidates supported their existence. Based on the occurrence of conserved protein motifs such as coiled-coil (CC), NBS, leucine-rich repeat (LRR), these regular NBS-LRR genes were classified into four subgroups: CC-NBS-LRR, NBS-LRR, CC-NBS, and X-NBS. Further expression analysis of the regular NBS-encoding genes in Brachypodium database revealed that these genes are expressed in a wide range of libraries, including those constructed from various developmental stages, tissue types, and drought challenged or nonchallenged tissue.

A novel cold-inducible gene from Pak-choi (Brassica campestris ssp. chinensis), BcWRKY46, enhances the cold, salt and dehydration stress tolerance in transgenic tobacco

WRKY TFs belong to one of the largest families of transcriptional regulators in plants and form integral parts of signaling webs that modulate many plant processes. BcWRKY46, a cDNA clone encoding a polypeptide of 284 amino acids and exhibited the structural features of group III of WRKY protein family, was isolated from the cold-treated leaves of Pak-choi (Brassica campestris ssp. chinensis Makino, syn. B. rapa ssp. chinensis) using the cDNA-AFLP technique. Expression of this gene was induced quickly and strongly in response to various environmental stresses, including low temperatures, ABA, salt and dehydration. Constitutive expression of BcWRKY46 in tobacco under the control of the CaMV35S promoter reduced the susceptibility of transgenic tobacco to freezing, ABA, salt and dehydration stresses. Our studies suggest that BcWRKY46 plays an important role in responding to ABA and abiotic stress.

Alternative splicing: enhancing ability to cope with stress via transcriptome plasticity

Alternative splicing is a mechanism for the regulation of gene expression that is widespread in higher eukaryotes. Genome-wide approaches, based on comparison of expressed and genomic sequences, on tiling arrays, and on next-generation sequencing, have provided growing knowledge of the extent, distribution and association of alternative splicing with stress-related genes in plants. The functional meaning of alternative splicing in response to stress has been defined for many genes, and in particular for those involved in the regulation of the stress responses, such as protein kinases, transcription factors, splicing regulators and pathogen-resistance genes. The production of proteins with diverse domain rearrangements from the same gene is the main alternative splicing mechanism for pathogen-resistance genes. The plant response to abiotic stress is also characterized by a second mechanism, which consists of the expression of alternative transcripts that are targeted to nonsense-mediated decay. These quantitatively regulate stress-related gene expression. Many alternative splicing events are well conserved among plant species, and also across kingdoms, especially those observed in response to stress, for genes encoding splicing regulators, and other classes of RNA-binding proteins. Nevertheless, non-conserved events indicate that alternative splicing represents an evolutionary strategy that rapidly increases genome plasticity and develops new gene functions, along with other mechanisms such as gene duplication. Finally, the study of the naturally occurring variability of alternative splicing and the identification of genomic regions involved in the regulation of alternative splicing in crops are proposed as strategies for selecting genotypes with superior performance under adverse environmental conditions.