Nodulin 41, a novel late nodulin of common bean with peptidase activity

Chitin nanofibers promote rhizobial symbiotic nitrogen fixation in Lotus japonicus

Chitin, an N-acetyl-D-glucosamine polymer, has multiple functions in living organisms, including the induction of disease resistance and growth promotion in plants. In addition, chitin oligosaccharides (COs) are used as the backbone of the signaling molecule Nod factor secreted by soil bacteria rhizobia to establish a mutual symbiosis with leguminous plants. Nod factor perception triggers host plant responses for rhizobial symbiosis. In this study, the effects of chitins on rhizobial symbiosis were examined in the leguminous plants Lotus japonicus and soybean. Chitin nanofiber (CNF), retained with polymeric structures, and COs elicited calcium spiking in L. japonicus roots expressing a nuclear-localized cameleon reporter. Shoot growth and symbiotic nitrogen fixation were significantly increased by CNF but not COs in L.japonicus and soybean. However, treatments with chitin and cellulose nanofiber, structurally similar polymers to CNF, did not affect shoot growth and nitrogen fixation in L.japonicus. Transcriptome analysis also supported the specific effects of CNF on rhizobial symbiosis in L.japonicus. Although chitins comprise the same monosaccharides and nanofibers share similar physical properties, only CNF can promote rhizobial nitrogen fixation in leguminous plants. Taking the advantages on physical properties, CNF could be a promising material for improving legume yield by enhancing rhizobial symbiosis.

Genome-Wide Analyses of Aspartic Proteases on Potato Genome (Solanum tuberosum): Generating New Tools to Improve the Resistance of Plants to Abiotic Stress

Aspartic proteases are proteolytic enzymes widely distributed in living organisms and viruses. Although they have been extensively studied in many plant species, they are poorly described in potatoes. The present study aimed to identify and characterize S. tuberosum aspartic proteases. Gene structure, chromosome and protein domain organization, phylogeny, and subcellular predicted localization were analyzed and integrated with RNAseq data from different tissues, organs, and conditions focused on abiotic stress. Sixty-two aspartic protease genes were retrieved from the potato genome, distributed in 12 chromosomes. A high number of intronless genes and segmental and tandem duplications were detected. Phylogenetic analysis revealed eight StAP groups, named from StAPI to StAPVIII, that were differentiated into typical (StAPI), nucellin-like (StAPIIIa), and atypical aspartic proteases (StAPII, StAPIIIb to StAPVIII). RNAseq data analyses showed that gene expression was consistent with the presence of cis-acting regulatory elements on StAP promoter regions related to water deficit. The study presents the first identification and characterization of 62 aspartic protease genes and proteins on the potato genome and provides the baseline material for functional gene determinations and potato breeding programs, including gene editing mediated by CRISPR.

Dynamics of the cell fate specifications during female gametophyte development in Arabidopsis

The female gametophytes of angiosperms contain cells with distinct functions, such as those that enable reproduction via pollen tube attraction and fertilization. Although the female gametophyte undergoes unique developmental processes, such as several rounds of nuclear division without cell plate formation and final cellularization, it remains unknown when and how the cell fate is determined during development. Here, we visualized the living dynamics of female gametophyte development and performed transcriptome analysis of individual cell types to assess the cell fate specifications in Arabidopsis thaliana. We recorded time lapses of the nuclear dynamics and cell plate formation from the 1-nucleate stage to the 7-cell stage after cellularization using an in vitro ovule culture system. The movies showed that the nuclear division occurred along the micropylar–chalazal (distal–proximal) axis. During cellularization, the polar nuclei migrated while associating with the forming edge of the cell plate, and then, migrated toward each other to fuse linearly. We also tracked the gene expression dynamics and identified that the expression of MYB98pro::GFP–MYB98, a synergid-specific marker, was initiated just after cellularization in the synergid, egg, and central cells and was then restricted to the synergid cells. This indicated that cell fates are determined immediately after cellularization. Transcriptome analysis of the female gametophyte cells of the wild-type and myb98 mutant revealed that the myb98 synergid cells had egg cell–like gene expression profiles. Although in myb98, egg cell–specific gene expression was properly initiated in the egg cells only after cellularization, but subsequently expressed ectopically in one of the 2 synergid cells. These results, together with the various initiation timings of the egg cell–specific genes, suggest complex regulation of the individual gametophyte cells, such as cellularization-triggered fate initiation, MYB98-dependent fate maintenance, cell morphogenesis, and organelle positioning. Our system of live-cell imaging and cell type–specific gene expression analysis provides insights into the dynamics and mechanisms of cell fate specifications in the development of female gametophytes in plants.

The female gametophytes of angiosperms contain cells with distinct functions, such as those that enable reproduction via pollen tube attraction and fertilization. Live-cell imaging and transcriptome analysis of single female gametophyte cell reveal novel insights into the dynamics and mechanisms of cell fate specifications in the model plant Arabidopsis.

Dynamics of the cell fate specifications during female gametophyte development in Arabidopsis

The female gametophytes of angiosperms contain cells with distinct functions, such as those that enable reproduction via pollen tube attraction and fertilization. Although the female gametophyte undergoes unique developmental processes, such as several rounds of nuclear division without cell plate formation and final cellularization, it remains unknown when and how the cell fate is determined during development. Here, we visualized the living dynamics of female gametophyte development and performed transcriptome analysis of individual cell types to assess the cell fate specifications in Arabidopsis thaliana. We recorded time lapses of the nuclear dynamics and cell plate formation from the 1-nucleate stage to the 7-cell stage after cellularization using an in vitro ovule culture system. The movies showed that the nuclear division occurred along the micropylar-chalazal (distal-proximal) axis. During cellularization, the polar nuclei migrated while associating with the forming edge of the cell plate, and then, migrated toward each other to fuse linearly. We also tracked the gene expression dynamics and identified that the expression of MYB98pro::GFP-MYB98, a synergid-specific marker, was initiated just after cellularization in the synergid, egg, and central cells and was then restricted to the synergid cells. This indicated that cell fates are determined immediately after cellularization. Transcriptome analysis of the female gametophyte cells of the wild-type and myb98 mutant revealed that the myb98 synergid cells had egg cell-like gene expression profiles. Although in myb98, egg cell-specific gene expression was properly initiated in the egg cells only after cellularization, but subsequently expressed ectopically in one of the 2 synergid cells. These results, together with the various initiation timings of the egg cell-specific genes, suggest complex regulation of the individual gametophyte cells, such as cellularization-triggered fate initiation, MYB98-dependent fate maintenance, cell morphogenesis, and organelle positioning. Our system of live-cell imaging and cell type-specific gene expression analysis provides insights into the dynamics and mechanisms of cell fate specifications in the development of female gametophytes in plants.

The rhizobial autotransporter determines the symbiotic nitrogen fixation activity of Lotus japonicus in a host-specific manner

Leguminous plants establish endosymbiotic associations with rhizobia and form root nodules in which the rhizobia fix atmospheric nitrogen. The host plant and intracellular rhizobia strictly control this symbiotic nitrogen fixation. We recently reported a Lotus japonicus Fix^- mutant, apn1 (aspartic peptidase nodule-induced 1), that impairs symbiotic nitrogen fixation. APN1 encodes a nodule-specific aspartic peptidase involved in the Fix^- phenotype in a rhizobial strain-specific manner. This host-strain specificity implies that some molecular interactions between host plant APN1 and rhizobial factors are required, although the biological function of APN1 in nodules and the mechanisms governing the interactions are unknown. To clarify how rhizobial factors are involved in strain-specific nitrogen fixation, we explored transposon mutants of Mesorhizobium loti strain TONO, which normally form Fix^- nodules on apn1 roots, and identified TONO mutants that formed Fix⁺ nodules on apn1 The identified causal gene encodes an autotransporter, part of a protein secretion system of Gram-negative bacteria. Expression of the autotransporter gene in M. loti strain MAFF3030399, which normally forms Fix⁺ nodules on apn1 roots, resulted in Fix^- nodules. The autotransporter of TONO functions to secrete a part of its own protein (a passenger domain) into extracellular spaces, and the recombinant APN1 protein cleaved the passenger protein in vitro. The M. loti autotransporter showed the activity to induce the genes involved in nodule senescence in a dose-dependent manner. Therefore, we conclude that the nodule-specific aspartic peptidase, APN1, suppresses negative effects of the rhizobial autotransporter in order to maintain effective symbiotic nitrogen fixation in root nodules.

An atypical aspartic protease modulates lateral root development in Arabidopsis thaliana

Few atypical aspartic proteases (APs) present in plants have been functionally studied to date despite having been implicated in developmental processes and stress responses. Here we characterize a novel atypical AP that we name Atypical Aspartic Protease in Roots 1 (ASPR1), denoting its expression in Arabidopsis roots. Recombinant ASPR1 produced by transient expression in Nicotiana benthamiana was active and displayed atypical properties, combining optimum acidic pH, partial sensitivity to pepstatin, pronounced sensitivity to redox agents, and unique specificity preferences resembling those of fungal APs. ASPR1 overexpression suppressed primary root growth and lateral root development, implying a previously unknown biological role for an AP. Quantitative comparison of wild-type and aspr1 root proteomes revealed deregulation of proteins associated with both reactive oxygen species and auxin homeostasis in the mutant. Together, our findings on ASPR1 reinforce the diverse pattern of enzymatic properties and biological roles of atypical APs and raise exciting questions on how these distinctive features impact functional specialization among these proteases.

Atypical and nucellin-like aspartic proteases: emerging players in plant developmental processes and stress responses

Members of the pepsin-like family (A1) of aspartic proteases (APs) are widely distributed in plants. A large number of genes encoding putative A1 APs are found in different plant genomes, the vast majority of which exhibit distinct features when compared with the so-called typical APs (and, therefore, grouped as atypical and nucellin-like APs). These features include the absence of the plant-specific insert; an unusually high number of cysteine residues; the nature of the amino acids preceding the first catalytic aspartate; and unexpected localizations. The over-representation of atypical and nucellin-like APs in plants is suggestive of greater diversification of protein functions and a more regulatory role for these APs, as compared with the housekeeping function generally attributed to typical APs. New functions have been uncovered for non-typical APs, with proposed roles in biotic and abiotic stress responses, chloroplast metabolism, and reproductive development, clearly suggesting functional specialization and tight regulation of activity. Furthermore, unusual enzymatic properties have also been documented for some of these proteases. Here, we give an overview of the current knowledge on the distinctive features and functions of both atypical and nucellin-like APs, and discuss this emerging pattern of functional complexity and specialization among plant pepsin-like proteases.

Uncovering Bax inhibitor-1 dual role in the legume-rhizobia symbiosis in common bean roots

Bax-inhibitor 1 (BI-1) is a cell death suppressor conserved in all eukaryotes that modulates cell death in response to abiotic stress and pathogen attack in plants. However, little is known about its role in the establishment of symbiotic interactions. Here, we demonstrate the functional relevance of an Arabidopsis thaliana BI-1 homolog (PvBI-1a) to symbiosis between the common bean (Phaseolus vulgaris) and Rhizobium tropici. We show that the changes in expression of PvBI-1a observed during early symbiosis resemble those of some defence response-related proteins. By using gain- and loss-of-function approaches, we demonstrate that the overexpression of PvBI-1a in the roots of common bean increases the number of rhizobial infection events (and therefore the final number of nodules per root), but induces the premature death of nodule cells, affecting their nitrogen fixation efficiency. Nodule morphological alterations are known to be associated with changes in the expression of genes tied to defence, autophagy, and vesicular trafficking. Results obtained in the present work suggest that BI-1 has a dual role in the regulation of programmed cell death during symbiosis, extending our understanding of its critical function in the modulation of host immunity while responding to beneficial microbes.

Loss-of-function of ASPARTIC PEPTIDASE NODULE-INDUCED 1 (APN1) in Lotus japonicus restricts efficient nitrogen-fixing symbiosis with specific Mesorhizobium loti strains

The nitrogen-fixing symbiosis of legumes and Rhizobium bacteria is established by complex interactions between the two symbiotic partners. Legume Fix^- mutants form apparently normal nodules with endosymbiotic rhizobia but fail to induce rhizobial nitrogen fixation. These mutants are useful for identifying the legume genes involved in the interactions essential for symbiotic nitrogen fixation. We describe here a Fix^- mutant of Lotus japonicus, apn1, which showed a very specific symbiotic phenotype. It formed ineffective nodules when inoculated with the Mesorhizobium loti strain TONO. In these nodules, infected cells disintegrated and successively became necrotic, indicating premature senescence typical of Fix^- mutants. However, it formed effective nodules when inoculated with the M. loti strain MAFF303099. Among nine different M. loti strains tested, four formed ineffective nodules and five formed effective nodules on apn1 roots. The identified causal gene, ASPARTIC PEPTIDASE NODULE-INDUCED 1 (LjAPN1), encodes a nepenthesin-type aspartic peptidase. The well characterized Arabidopsis aspartic peptidase CDR1 could complement the strain-specific Fix^- phenotype of apn1. LjAPN1 is a typical late nodulin; its gene expression was exclusively induced during nodule development. LjAPN1 was most abundantly expressed in the infected cells in the nodules. Our findings indicate that LjAPN1 is required for the development and persistence of functional (nitrogen-fixing) symbiosis in a rhizobial strain-dependent manner, and thus determines compatibility between M. loti and L. japonicus at the level of nitrogen fixation.

Genome resequencing and transcriptome profiling reveal structural diversity and expression patterns of constitutive disease resistance genes in Huanglongbing-tolerant Poncirus trifoliata and its hybrids

Huanglongbing (HLB) is the most destructive bacterial disease of citrus worldwide. While most citrus varieties are susceptible to HLB, Poncirus trifoliata, a close relative of Citrus, and some of its hybrids with Citrus are tolerant to HLB. No specific HLB tolerance genes have been identified in P. trifoliata but recent studies have shown that constitutive disease resistance (CDR) genes were expressed at much higher levels in HLB-tolerant Poncirus hybrids and the expression of CDR genes was modulated by Candidatus Liberibacter asiaticus (CLas), the pathogen of HLB. The current study was undertaken to mine and characterize the CDR gene family in Citrus and Poncirus and to understand its association with HLB tolerance in Poncirus. We identified 17 CDR genes in two citrus genomes, deduced their structures, and investigated their phylogenetic relationships. We revealed that the expansion of the CDR family in Citrus seems to be due to segmental and tandem duplication events. Through genome resequencing and transcriptome sequencing, we identified eight CDR genes in the Poncirus genome (PtCDR1-PtCDR8). The number of SNPs was the highest in PtCDR2 and the lowest in PtCDR7. Most of the deletion and insertion events were observed in the UTR regions of Citrus and Poncirus CDR genes. PtCDR2 and PtCDR8 were in abundance in the leaf transcriptomes of two HLB-tolerant Poncirus genotypes and were also upregulated in HLB-tolerant, Poncirus hybrids as revealed by real-time PCR analysis. These two CDR genes seem to be good candidate genes for future studies of their role in citrus-CLas interactions.

Toxicogenomic Responses of the Model Legume Medicago truncatula to Aged Biosolids Containing a Mixture of Nanomaterials (TiO₂, Ag, and ZnO) from a Pilot Wastewater Treatment Plant

Toxicogenomic responses in Medicago truncatula A17 were monitored following exposure to biosolids-amended soils. Treatments included biosolids produced using a pilot wastewater treatment plant with either no metal introduced into the influent (control); bulk/ionic TiO2, ZnO, and AgNO3 added to influent (bulk/dissolved treatment); or Ag, ZnO, and TiO2 engineered nanomaterials added to influent (ENM treatment) and then added to soil, which was aged in the field for 6 months. In our companion study, we found inhibition of nodulation in the ENM but not in the bulk/dissolved treatment. Gene expression profiling revealed highly distinct profiles with more than 10-fold down-regulation in 239 genes in M. truncatula roots from the ENM treatment, while gene expression patterns were similar between bulk/dissolved and control treatments. In response to ENM exposure, many of the identified biological pathways, gene ontologies, and individual genes are associated with nitrogen metabolism, nodulation, metal homeostasis, and stress responses. Expression levels of nine genes were independently confirmed with qRT-PCR. Exposure to ENMs induced unique shifts in expression profiles and biological pathways compared with bulk/dissolved treatment, despite the lack of difference in bioavailable metal fractions, metal oxidation state, and coordination environment between ENM and bulk/dissolved biosolids. As populations of Sinorhizobium meliloti Rm2011 were similar in bulk/dissolved and ENM treatments, our results suggest that inhibition of nodulation in the ENM treatment was primarily due to phytotoxicity, likely caused by enhanced bioavailability of Zn ions.

Proteomic analysis of the soybean symbiosome identifies new symbiotic proteins

Legumes form a symbiosis with rhizobia in which the plant provides an energy source to the rhizobia bacteria that it uses to fix atmospheric nitrogen. This nitrogen is provided to the legume plant, allowing it to grow without the addition of nitrogen fertilizer. As part of the symbiosis, the bacteria in the infected cells of a new root organ, the nodule, are surrounded by a plant-derived membrane, the symbiosome membrane, which becomes the interface between the symbionts. Fractions containing the symbiosome membrane (SM) and material from the lumen of the symbiosome (peribacteroid space or PBS) were isolated from soybean root nodules and analyzed using nongel proteomic techniques. Bicarbonate stripping and chloroform-methanol extraction of isolated SM were used to reduce complexity of the samples and enrich for hydrophobic integral membrane proteins. One hundred and ninety-seven proteins were identified as components of the SM, with an additional fifteen proteins identified from peripheral membrane and PBS protein fractions. Proteins involved in a range of cellular processes such as metabolism, protein folding and degradation, membrane trafficking, and solute transport were identified. These included a number of proteins previously localized to the SM, such as aquaglyceroporin nodulin 26, sulfate transporters, remorin, and Rab7 homologs. Among the proteome were a number of putative transporters for compounds such as sulfate, calcium, hydrogen ions, peptide/dicarboxylate, and nitrate, as well as transporters for which the substrate is not easy to predict. Analysis of the promoter activity for six genes encoding putative SM proteins showed nodule specific expression, with five showing expression only in infected cells. Localization of two proteins was confirmed using GFP-fusion experiments. The data have been deposited to the ProteomeXchange with identifier PXD001132. This proteome will provide a rich resource for the study of the legume-rhizobium symbiosis.

An RNA-Seq based gene expression atlas of the common bean

BACKGROUND

Common bean (Phaseolus vulgaris) is grown throughout the world and comprises roughly 50% of the grain legumes consumed worldwide. Despite this, genetic resources for common beans have been lacking. Next generation sequencing, has facilitated our investigation of the gene expression profiles associated with biologically important traits in common bean. An increased understanding of gene expression in common bean will improve our understanding of gene expression patterns in other legume species.

RESULTS

Combining recently developed genomic resources for Phaseolus vulgaris, including predicted gene calls, with RNA-Seq technology, we measured the gene expression patterns from 24 samples collected from seven tissues at developmentally important stages and from three nitrogen treatments. Gene expression patterns throughout the plant were analyzed to better understand changes due to nodulation, seed development, and nitrogen utilization. We have identified 11,010 genes differentially expressed with a fold change ≥ 2 and a P-value < 0.05 between different tissues at the same time point, 15,752 genes differentially expressed within a tissue due to changes in development, and 2,315 genes expressed only in a single tissue. These analyses identified 2,970 genes with expression patterns that appear to be directly dependent on the source of available nitrogen. Finally, we have assembled this data in a publicly available database, The Phaseolus vulgaris Gene Expression Atlas (Pv GEA), http://plantgrn.noble.org/PvGEA/ . Using the website, researchers can query gene expression profiles of their gene of interest, search for genes expressed in different tissues, or download the dataset in a tabular form.

CONCLUSIONS

These data provide the basis for a gene expression atlas, which will facilitate functional genomic studies in common bean. Analysis of this dataset has identified genes important in regulating seed composition and has increased our understanding of nodulation and impact of the nitrogen source on assimilation and distribution throughout the plant.

Proteome reference maps of the Lotus japonicus nodule and root

Legume symbiosis with rhizobia results in the formation of a specialized organ, the root nodule, where atmospheric dinitrogen is reduced to ammonia. In Lotus japonicus (Lotus), several genes involved in nodule development or nodule function have been defined using biochemistry, genetic approaches, and high-throughput transcriptomics. We have employed proteomics to further understand nodule development. Two developmental stages representing nodules prior to nitrogen fixation (white) and mature nitrogen fixing nodules (red) were compared with roots. In addition, the proteome of a spontaneous nodule formation mutant (snf1) was determined. From nodules and roots, 780 and 790 protein spots from 2D gels were identified and approximately 45% of the corresponding unique gene accessions were common. Including a previous proteomics set from Lotus pod and seed, the common gene accessions were decreased to 7%. Interestingly, an indication of more pronounced PTMs in nodules than in roots was determined. Between the two nodule developmental stages, higher levels of pathogen-related 10 proteins, HSPs, and proteins involved in redox processes were found in white nodules, suggesting a higher stress level at this developmental stage. In contrast, protein spots corresponding to nodulins such as leghemoglobin, asparagine synthetase, sucrose synthase, and glutamine synthetase were prevalent in red nodules. The distinct biochemical state of nodules was further highlighted by the conspicuous presence of several nitrilases, ascorbate metabolic enzymes, and putative rhizobial effectors.