RNA-Seq Analysis of Differential Gene Expression Responding to Different Rhizobium Strains in Soybean (Glycine max) Roots

Comparative genomic and transcriptomic analyses provide new insight into symbiotic host specificity

Host specificity plays important roles in expanding the host range of rhizobia, while the genetic information responsible for host specificity remains largely unexplored. In this report, the roots of four symbiotic systems with notable different symbiotic phenotypes and the control were studied at four different post-inoculation time points by RNA sequencning (RNA-seq). The differentially expressed genes (DEGs) were divided into "found only in soybean or Lotus," "only expressed in soybean or Lotus," and "expressed in both hosts" according to the comparative genomic analysis. The distributions of enriched function ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways vary significantly in different symbiotic systems. Host specific genes account for the majority of the DEGs involved in response to stimulus, associated with plant-pathogen interaction pathways, and encoding resistance (R) proteins, the symbiotic nitrogen fixation (SNF) proteins and the target proteins in the SNF-related modules. Our findings provided molecular candidates for better understanding the mechanisms of symbiotic host-specificity.

RNA-Seq analysis of mung bean (Vigna radiata L.) roots shows differential gene expression and predicts regulatory pathways responding to taxonomically different rhizobia

Symbiotic interaction among legume and rhizobia is a complex phenomenon which results in the formation of nitrogen-fixing nodules. Mung bean is promiscuous host however expression profile of this important legume plant in response to rhizobial infection was particularly lacking and urgently needed. We have demonstrated the pattern of gene expression of mung bean roots inoculated with two symbionts Bradyrhizobium yuanmingense Vr50 and Sinorhizobium (Ensifer) aridi Vr33 and non-inoculated control (CK). The RNA-Seq data analyzed at two growth stages i.e., 1-3 h and 10-16 days post inoculation revealed significantly higher number of differentially expressed genes (DEGs) at nodulation stage. The DEGs encoding receptor kinases identified at early stage might be involved in perception of Nod factors produced by different rhizobia. At nodulation stage important genes involved in plant hormone signal transduction, nitrogen and sulfur metabolism were identified. KEGG pathway enrichment analysis showed that metabolic pathways were most prominent in both groups (Group 1: Vr33 vs CK; Group 2: Vr50 vs CK), followed by biosynthesis of secondary metabolites, plant hormone signal transduction and biosynthesis of amino acids. Furthermore, DEGs involved in cell communication and plant hormone signal transduction were found to be different among two symbiotic systems while DEGs involved in carbon, nitrogen and sulfur metabolism were similar but their expression varied in response to two rhizobial strains. This study provides the first insight into the mechanisms underlying interactions of mung bean host with two taxonomically different symbionts (Bradyrhizobium and Sinorhizobium) and the candidate genes for better understanding the mechanisms of symbiotic host-specificity.

Soybean gene co-expression network analysis identifies two co-regulated gene modules associated with nodule formation and development

Gene co-expression network analysis is an efficient systems biology approach for the discovery of novel gene functions and trait-associated gene modules. To identify clusters of functionally related genes involved in soybean nodule formation and development, we performed a weighted gene co-expression network analysis. Two nodule-specific modules (NSM-1 and NSM-2, containing 304 and 203 genes, respectively) were identified. The NSM-1 gene promoters were significantly enriched in cis-binding elements for ERF, MYB, and C2H2-type zinc transcription factors, whereas NSM-2 gene promoters were enriched in cis-binding elements for TCP, bZIP, and bHLH transcription factors, suggesting a role of these regulatory factors in the transcriptional activation of nodule co-expressed genes. The co-expressed gene modules included genes with potential novel roles in nodulation, including those involved in xylem development, transmembrane transport, the ethylene signalling pathway, cytoskeleton organization, cytokinesis and regulation of the cell cycle, regulation of meristem initiation and growth, transcriptional regulation, DNA methylation, and histone modifications. Functional analysis of two co-expressed genes using TILLING mutants provided novel insight into the involvement of unsaturated fatty acid biosynthesis and folate metabolism in nodule formation and development. The identified gene co-expression modules provide valuable resources for further functional genomics studies to dissect the genetic basis of nodule formation and development in soybean.

Identifying the Soybean microRNAs Related to Phytophthora sojae Based on RNA Sequencing and Bioinformatics Analysis

Phytophthora root rot in soybeans is caused by a pathogen called Phytophthora sojae (P. sojae), which results in a significant decrease in soybean production within affected regions. MicroRNAs (miRNAs) are a class of small non-coding RNA molecules that play a key post-transcriptional regulatory role in eukaryotes. In this paper, the miRNAs that respond to P. sojae were analyzed from the gene level to complement the study of molecular resistance mechanisms in soybean. The study utilized high-throughput sequencing of soybean data to predict miRNAs that respond to P. sojae, analyze their specific functions, and verify regulatory relationships using qRT-PCR. The results showed that the miRNAs in soybean respond to P. sojae infection. MiRNAs can be transcribed independently, suggesting the presence of transcription factor binding sites in the promoter regions. Additionally, we performed an evolutionary analysis on conserved miRNAs that respond to P. sojae. Finally, we investigated the regulatory relationships among miRNAs, genes, and transcription factors, and identified five regulatory patterns. These findings lay the groundwork for future studies on the evolution of miRNAs responsive to P. sojae.

Current understanding of plant-microbe interaction through the lenses of multi-omics approaches and their benefits in sustainable agriculture

The success of sustainable agricultural practices has now become heavily dependent on the interactions between crop plants and their associated microbiome. Continuous advancement in high throughput sequencing platforms, omics-based approaches, and gene editing technologies has remarkably accelerated this area of research. It has enabled us to characterize the interactions of plants with associated microbial communities more comprehensively and accurately. Furthermore, the genomic and post-genomic era has significantly refined our perspective toward the complex mechanisms involved in those interactions, opening new avenues for efficiently deploying the knowledge in developing sustainable agricultural practices. This review focuses on our fundamental understanding of plant-microbe interactions and the contribution of existing multi-omics approaches, including those under active development and their tremendous success in unraveling different aspects of the complex network between plant hosts and microbes. In addition, we have also discussed the importance of sustainable and eco-friendly agriculture and the associated outstanding challenges ahead.

Phosphate-Solubilizing Pseudomonas sp. Strain WS32 Rhizosphere Colonization-Induced Expression Changes in Wheat Roots

Rhizosphere colonization is a pre-requisite for the favorable application of plant growth-promoting rhizobacteria (PGPR). Exchange and mutual recognition of signaling molecules occur frequently between plants and microbes. Here, the luciferase luxAB gene was electrotransformed into the phosphate-solubilizing strain Pseudomonas sp. WS32, a type of plant growth-promoting rhizobacterium with specific affinity for wheat. A labeled WS32 strain (WS32-L) was applied to determine the temporal and spatial traits of colonization within the wheat rhizosphere using rhizoboxes experimentation under natural condition. The effects of colonization on wheat root development and seedling growth were evaluated, and RNA sequencing (RNA-seq) was performed to explore the transcriptional changes that occur in wheat roots under WS32 colonization. The results showed that WS32-L could survive in the wheat rhizosphere for long periods and could expand into new zones following wheat root extension. Significant increases in seedling fresh and dry weight, root fresh and dry weight, root surface area, number of root tips, and phosphorus accumulation in the wheat leaves occurred in response to WS32 rhizosphere colonization. RNA-seq analysis showed that a total of 1485 genes in wheat roots were differentially expressed between the inoculated conditions and the uninoculated conditions. Most of the transcriptional changes occurred for genes annotated to the following functional categories: "phosphorus and other nutrient transport," "hormone metabolism and organic acid secretion," "flavonoid signal recognition," "membrane transport," and "transcription factor regulation." These results are therefore valuable to future studies focused on the molecular mechanisms underlying the growth-promoting activities of PGPR on their host plants.

A Nodule-Localized Small Heat Shock Protein GmHSP17.1 Confers Nodule Development and Nitrogen Fixation in Soybean

Small heat shock proteins (sHSPs) are ubiquitous proteins present in all organisms. The sHSPs are not only upregulated under heat shock as well as other stresses but also are expressed in unstressed cells, indicating quite diverse functions of sHSPs. However, there is little known about the role of sHSPs in nodulation and nitrogen fixation in soybean. In this study, we cloned a candidate protein of sHSP, GmHSP17.1, from proteome of nodule and analyzed its function in soybean nodulation. We found that GmHSP17.1 was a cytosolic protein and preferentially expressed during nodule development. An overexpression of GmHSP17.1 in composite transgenic plants showed increases in nodule number, fresh weight, nodule size, area of infection cells, and nitrogenase activity, and subsequently promoted the content of nitrogen and growth of soybean plants. While GmHSP17.1 RNA interference (RNAi) lines showed significantly impaired nodule development and nitrogen fixation efficiency. Through liquid chromatography-tandem mass spectrometry (LC-MS/MS), GmRIP1 was identified as the first potential target of GmHSP17.1, and was shown to be specifically expressed in soybean nodules. The interaction between GmHSP17.1 and GmRIP1 was further confirmed by yeast-two hybrid (Y2H), bimolecular fluorescence complementation (BiFC) in vivo and pull-down assay in vitro. Furthermore, peroxidase activity was markedly increased in GmHSP17.1 overexpressed nodules and decreased in RNAi lines. As a result, the reactive oxygen species (ROS) content greatly decreased in GmHSP17.1 overexpression lines and increased in suppression lines. Taken together, we conclude that GmHSP17.1 plays an important role in soybean nodulation through interacting with GmRIP1. Our results provide foundation for studying the mechanism of nitrogen fixation and for the genetics improvement of legume plants.

Transcriptomic analysis of Mesoamerican and Andean Phaseolus vulgaris accessions revealed mRNAs and lncRNAs associated with strain selectivity during symbiosis

Legume plants establish a nitrogen-fixing symbiosis with soil bacteria known as rhizobia. Compatibility between legumes and rhizobia is determined at species-specific level, but variations in the outcome of the symbiotic process are also influenced by the capacity of the plant to discriminate and select specific strains that are better partners. We compared the transcriptional response of two genetically diverse accessions of Phaseolus vulgaris from Mesoamerica and South Andes to Rhizobium etli strains that exhibit variable degrees of symbiotic affinities. Our results indicate that the plant genotype is the major determinant of the transcriptional reprogramming occurring in roots at early stages of the symbiotic interaction. Differentially expressed genes (DEGs) regulated in the Mesoamerican and the Andean accessions in response to specific strains are different, but they belong to the same functional categories. The common and strain-specific transcriptional responses to rhizobia involve distinct transcription factors  and cis-elements present in the promoters of DEGs in each accession, showing that diversification and domestication of common bean at different geographic regions influenced the evolution of symbiosis differently in each genetic pool. Quantitative PCR analysis validated our transcriptional datasets, which constitute a valuable source of coding and non-coding candidate genes to further unravel the molecular determinants governing the mechanisms by which plants select bacterial strains that produce a better symbiotic outcome.

Avoidance of detrimental defense responses in beneficial plant-microbe interactions

In the environment microbes interact with plants and provide them with benefits that include protection against biotic and abiotic stresses as well as improved nutrition. However, plants are also exposed to parasites and pathogens. To manage appropriate responses, evolution has resulted in improved tolerance of plants to beneficial microbes while keeping the ability to recognize detrimental ones and to develop defense responses. Here we review the mechanisms involved in these interactions. We also discuss how the interactions might be handled to improve crop resistance to pathogens without losing the ability to establish beneficial interactions.

Whole-Genome Sequencing of Bradyrhizobium diazoefficiens 113-2 and Comparative Genomic Analysis Provide Molecular Insights Into Species Specificity and Host Specificity

In the present study, we sequenced the complete genome of Bradyrhizobium diazoefficiens 113-2. The genomic characteristics of six selected rhizobial strains (two fast-growing rhizobia, two medium-slow-growing rhizobia and two slow-growing rhizobia) with four different legume hosts were analyzed by comparative genomic analysis. Genomes of B. diazoefficiens 113-2 and B. diazoefficiens USDA110 were found to share a large synteny blocks and a high ANI value, supporting 113-2 as a strain of B. diazoefficiens. 5,455 singletons and 11,656 clusters were identified among the six rhizobia genomes, and most of the pair-wise comparisons clusters were shared by the two genomes of strains in the same genus. Similar genus-specific gene numbers in the assigned COG functional terms were present in the two strains of the same genus, while the numbers were decreased with the increase of growth rate in most of the COG terms. KEGG pathway analysis of B. diazoefficiens 113-2 suggested that the rhizobial genes in ABC transporters and Two-Component system were mainly species-specific. Besides, the candidate genes related to secretion system and surface polysaccharides biosynthesis in the genomes of the six strains were explored and compared. 39 nodulation gene families, 12 nif gene families and 10 fix gene families in the genomes of these six strains were identified, and gene classes in most of gene families and the types and total gene numbers of gene families were substantially different among these six genomes. We also performed synteny analyses for above-mentioned nod, nif, and fix gene groupings, and selected NodW, NolK, NoeJ, NifB, FixK, and FixJ gene families to perform phylogeny analyses. Our results provided valuable molecular insights into species specificity and host specificity. The genetic information responsible for host specificity will play important roles in expanding the host range of rhizobia among legumes, which might provide new clues for the understanding of the genetic determinants of non-legume-rhizobium symbiosis.

Genome-wide survey of soybean papain-like cysteine proteases and their expression analysis in root nodule symbiosis

BACKGROUND

Plant papain-like cysteine proteases (PLCPs) are a large class of proteolytic enzymes and play important roles in root nodule symbiosis (RNS), while the whole-genome studies of PLCP family genes in legume are quite limited, and the roles of Glycine max PLCPs (GmPLCPs) in nodulation, nodule development and senescence are not fully understood.

RESULTS

In the present study, we identified 97 GmPLCPs and performed a genome-wide survey to explore the expansion of soybean PLCP family genes and their relationships to RNS. Nineteen paralogous pairs of genomic segments, consisting of 77 GmPLCPs, formed by whole-genome duplication (WGD) events were identified, showing a high degree of complexity in duplication. Phylogenetic analysis among different species showed that the lineage differentiation of GmPLCPs occurred after family expansion, and large tandem repeat segment were specifically in soybean. The expression patterns of GmPLCPs in symbiosis-related tissues and nodules identified RNS-related GmPLCPs and provided insights into their putative symbiotic functions in soybean. The symbiotic function analyses showed that a RNS-related GmPLCP gene (Glyma.04G190700) really participate in nodulation and nodule development.

CONCLUSIONS

Our findings improved our understanding of the functional diversity of legume PLCP family genes, and provided insights into the putative roles of the legume PLCPs in nodulation, nodule development and senescence.

New tools for characterizing early brown stem rot disease resistance signaling in soybean

Brown stem rot (BSR) reduces soybean [Glycine max (L.) Merr.] yield by up to 38%. The BSR causal agent is Phialophora gregata f. sp. sojae, a slow-growing, necrotrophic fungus whose life cycle includes latent and pathogenic phases, each lasting several weeks. Brown stem rot foliar symptoms are often misdiagnosed as other soybean diseases or nutrient stress, making BSR resistance especially difficult to phenotype. To shed light on the genes and networks contributing to P. gregata resistance, we conducted RNA sequencing (RNA-seq) of a resistant genotype (PI 437970, Rbs3). Leaf, stem, and root tissues were collected 12, 24, and 36 h after stab inoculation with P. gregata, or mock infection, in the plant stem. By using multiple tissues and time points, we could see that leaves, stems, and roots use the same defense pathways. Our analyses suggest that P. gregata induces a biphasic defense response, with pathogen-associated molecular pattern (PAMP) triggered immunity observed in leaves at 12 and 24 h after infection (HAI) and effector triggered immunity detected at 36 h after infection in the stems. Gene networks associated with defense, photosynthesis, nutrient homeostasis, DNA replication, and growth are the hallmarks of resistance to P. gregata. While P. gregata is a slow-growing pathogen, our results demonstrate that pathogen recognition occurs hours after infection. By exploiting the genes and networks described here, we will be able to develop novel diagnostic tools to facilitate breeding and screening for BSR resistance.

The Significance of Flavonoids in the Process of Biological Nitrogen Fixation

Nitrogen is essential for the growth of plants. The ability of some plant species to obtain all or part of their requirement for nitrogen by interacting with microbial symbionts has conferred a major competitive advantage over those plants unable to do so. The function of certain flavonoids (a group of secondary metabolites produced by the plant phenylpropanoid pathway) within the process of biological nitrogen fixation carried out by Rhizobium spp. has been thoroughly researched. However, their significance to biological nitrogen fixation carried out during the actinorhizal and arbuscular mycorrhiza-Rhizobium-legume interaction remains unclear. This review catalogs and contextualizes the role of flavonoids in the three major types of root endosymbiosis responsible for biological nitrogen fixation. The importance of gaining an understanding of the molecular basis of endosymbiosis signaling, as well as the potential of and challenges facing modifying flavonoids either quantitatively and/or qualitatively are discussed, along with proposed strategies for both optimizing the process of nodulation and widening the plant species base, which can support nodulation.

Interactions of gene expression, alternative splicing, and DNA methylation in determining nodule identity

Soybean nodulation is a highly controlled process that involves complex gene regulation at both transcriptional and post-transcriptional levels. In the present study, we profiled gene expression changes, alternative splicing events, and DNA methylation patterns during nodule formation, development, and senescence. The transcriptome data uncovered key transcription patterns of nodule development that included 9669 core genes and 7302 stage-specific genes. Alternative splicing analysis uncovered a total of 2323 genes that undergo alternative splicing events in at least one nodule developmental stage, with activation of exon skipping and repression of intron retention being the most common splicing events in nodules compared to roots. Approximately 40% of the differentially spliced genes were also differentially expressed at the same nodule developmental stage, implying a substantial association between gene expression and alternative splicing. Genome-wide-DNA methylation analysis revealed dynamic changes in nodule methylomes that were specific to each nodule stage, occurred in a sequence-specific manner, and impacted the expression of 1864 genes. An attractive hypothesis raised by our data is that increased DNA methylation may contribute to the efficiency of alternative splicing. Together, our results provide intriguing insights into the associations between gene expression, alternative splicing, and DNA methylation that may shape transcriptome complexity and proteome specificity in developing soybean nodules.

Plant transcriptome analysis reveals specific molecular interactions between alfalfa and its rhizobial symbionts below the species level

BACKGROUND

Leguminous plants alter patterns of gene expression in response to symbiotic colonization and infection by their cognate rhizobial bacteria, but the extent of the transcriptomic response has rarely been examined below the species level. Here we describe the identification of 12 rhizobial biotypes of Ensifer meliloti, which form nitrogen-fixing nodules in the roots of alfalfa (Medicago sativa L.), followed by a comparative RNA-seq analysis of four alfalfa cultivars each inoculated with two E. meliloti strains varying in symbiotic performance and phylogenetic relatedness.

RESULTS

Rhizobial biotypes were identified on the basis of their symbiotic performance, particularly shoot dry weight. Differentially expressed genes (DEGs) and metabolic pathways were determined by comparing the RNA-seq data with that of the uninoculated control plant. Significant differences were found between DEGs generated in each cultivar with the inoculation of two rhizobial strains in comparison (P < 0.01). A total of 8111 genes was differentially expressed, representing ~ 17.1% of the M. sativa genome. The proportion of DEGs ranges from 0.5 to 12.2% for each alfalfa cultivar. Interestingly, genes with predicted roles in flavonoid biosynthesis and plant-pathogen interaction (NBS-LRR) were identified as the most significant DEGs. Other DEGs include Medsa002106 and genes encoding nodulins and NCR peptides whose expression is specifically induced during the development of nitrogen-fixing nodules. More importantly, strong significant positive correlations were observed between plant transcriptomes (DEGs and KEGG pathways) and phylogenetic distances between the two rhizobial inoculants.

CONCLUSIONS

Alfalfa expresses significantly distinct sets of genes in response to infection by different rhizobial strains at the below-species levels (i.e. biotype or strain). Candidate genes underlying the specific interactions include Medsa002106 and those encoding nodulins and NCR peptides and proteins in the NBS-LRR family.

Dynamic gene expression changes in response to micronutrient, macronutrient, and multiple stress exposures in soybean

Preserving crop yield is critical for US soybean production and the global economy. Crop species have been selected for increased yield for thousands of years with individual lines selected for improved performance in unique environments, constraints not experienced by model species such as Arabidopsis. This selection likely resulted in novel stress adaptations, unique to crop species. Given that iron deficiency is a perennial problem in the soybean growing regions of the USA and phosphate deficiency looms as a limitation to global agricultural production, nutrient stress studies in crop species are critically important. In this study, we directly compared whole-genome expression responses of leaves and roots to iron (Fe) and phosphate (Pi) deficiency, representing a micronutrient and macronutrient, respectively. Conducting experiments side by side, we observed soybean responds to both nutrient deficiencies within 24 h. While soybean responds largely to -Fe deficiency, it responds strongly to Pi resupply. Though the timing of the responses was different, both nutrient stress signals used the same molecular pathways. Our study is the first to demonstrate the speed and diversity of the soybean stress response to multiple nutrient deficiencies. We also designed the study to examine gene expression changes in response to multiple stress events. We identified 865 and 3375 genes that either altered their direction of expression after a second stress exposure or were only differentially expressed after a second stress event. Understanding the molecular underpinnings of these responses in crop species could have major implications for improving stress tolerance and preserving yield.

Transcriptome analysis of soybean (Glycine max) root genes differentially expressed in rhizobial, arbuscular mycorrhizal, and dual symbiosis

Soybean (Glycine max) roots establish associations with nodule-inducing rhizobia and arbuscular mycorrhizal (AM) fungi. Both rhizobia and AM fungi have been shown to affect the activity of and colonization by the other, and their interactions can be detected within host plants. Here, we report the transcription profiles of genes differentially expressed in soybean roots in the presence of rhizobial, AM, or rhizobial-AM dual symbiosis, compared with those in control (uninoculated) roots. Following inoculation, soybean plants were grown in a glasshouse for 6 weeks; thereafter their root transcriptomes were analyzed using an oligo DNA microarray. Among the four treatments, the root nodule number and host plant growth were highest in plants with dual symbiosis. We observed that the expression of 187, 441, and 548 host genes was up-regulated and 119, 1,439, and 1,298 host genes were down-regulated during rhizobial, AM, and dual symbiosis, respectively. The expression of 34 host genes was up-regulated in each of the three symbioses. These 34 genes encoded several membrane transporters, type 1 metallothionein, and transcription factors in the MYB and bHLH families. We identified 56 host genes that were specifically up-regulated during dual symbiosis. These genes encoded several nodulin proteins, phenylpropanoid metabolism-related proteins, and carbonic anhydrase. The nodulin genes up-regulated by the AM fungal colonization probably led to the observed increases in root nodule number and host plant growth. Some other nodulin genes were down-regulated specifically during AM symbiosis. Based on the results above, we suggest that the contribution of AM fungal colonization is crucial to biological N2-fixation and host growth in soybean with rhizobial-AM dual symbiosis.

Bradyrhizobium diazoefficiens USDA 110- Glycine max Interactome Provides Candidate Proteins Associated with Symbiosis

Although the legume-rhizobium symbiosis is a most-important biological process, there is a limited knowledge about the protein interaction network between host and symbiont. Using interolog- and domain-based approaches, we constructed an interspecies protein interactome containing 5115 protein-protein interactions between 2291 Glycine max and 290 Bradyrhizobium diazoefficiens USDA 110 proteins. The interactome was further validated by the expression pattern analysis in nodules, gene ontology term semantic similarity, co-expression analysis, and luciferase complementation image assay. In the G. max-B. diazoefficiens interactome, bacterial proteins are mainly ion channel and transporters of carbohydrates and cations, while G. max proteins are mainly involved in the processes of metabolism, signal transduction, and transport. We also identified the top 10 highly interacting proteins (hubs) for each species. Kyoto Encyclopedia of Genes and Genomes pathway analysis for each hub showed that a pair of 14-3-3 proteins (SGF14g and SGF14k) and 5 heat shock proteins in G. max are possibly involved in symbiosis, and 10 hubs in B. diazoefficiens may be important symbiotic effectors. Subnetwork analysis showed that 18 symbiosis-related soluble N-ethylmaleimide sensitive factor attachment protein receptor proteins may play roles in regulating bacterial ion channels, and SGF14g and SGF14k possibly regulate the rhizobium dicarboxylate transport protein DctA. The predicted interactome provide a valuable basis for understanding the molecular mechanism of nodulation in soybean.

Genome-wide analysis of genes encoding core components of the ubiquitin system in soybean (Glycine max) reveals a potential role for ubiquitination in host immunity against soybean cyst nematode

BACKGROUND

Ubiquitination is a major post-translational protein modification that regulates essentially all cellular and physiological pathways in eukaryotes. The ubiquitination process typically involves three distinct classes of enzymes, ubiquitin-activating enzyme (E1), ubiquitin-conjugating enzyme (E2) and ubiquitin ligase (E3). To date, a comprehensive identification and analysis of core components comprising of the whole soybean (Glycine max) ubiquitin system (UBS) has not been reported.

RESULTS

We performed a systematic, genome-wide analysis of genes that encode core members of the soybean UBS in this study. A total of 1431 genes were identified with high confidence to encode putative soybean UBS components, including 4 genes encoding E1s, 71 genes that encode the E2s, and 1356 genes encoding the E3-related components. Among the E3-encoding genes, 760 encode RING-type E3s, 124 encode U-box domain-containing E3s, and 472 encode F-box proteins. To find out whether the identified soybean UBS genes encode active enzymes, a set of genes were randomly selected and the enzymatic activities of their recombinant proteins were tested. Thioester assays indicated proteins encoded by the soybean E1 gene GmUBA1 and the majority of selected E2 genes are active E1 or E2 enzymes, respectively. Meanwhile, most of the purified RING and U-box domain-containing proteins displayed E3 activity in the in vitro ubiquitination assay. In addition, 1034 of the identified soybean UBS genes were found to express in at least one of 14 soybean tissues examined and the transcript level of 338 soybean USB genes were significantly changed after abiotic or biotic (Fusarium oxysporum and Rhizobium strains) stress treatment. Finally, the expression level of a large number of the identified soybean UBS-related genes was found significantly altered after soybean cyst nematode (SCN) treatment, suggesting the soybean UBS potentially plays an important role in soybean immunity against SCN.

CONCLUSIONS

Our findings indicate the presence of a large and diverse number of core UBS proteins in the soybean genome, which suggests that target-specific modification by ubiquitin is a complex and important part of cellular and physiological regulation in soybean. We also revealed certain members of the soybean UBS may be involved in immunity against soybean cyst nematode (SCN). This study sets up an essential foundation for further functional characterization of the soybean UBS in various physiological processes, such as host immunity against SCN.

Metabolomics and Transcriptomics Identify Multiple Downstream Targets of Paraburkholderia phymatum σ54 During Symbiosis with Phaseolus vulgaris

RpoN (or σ⁵⁴) is the key sigma factor for the regulation of transcription of nitrogen fixation genes in diazotrophic bacteria, which include α- and β-rhizobia. Our previous studies showed that an rpoN mutant of the β-rhizobial strain Paraburkholderia phymatum STM815^T formed root nodules on Phaseolus vulgaris cv. Negro jamapa, which were unable to reduce atmospheric nitrogen into ammonia. In an effort to further characterize the RpoN regulon of P. phymatum, transcriptomics was combined with a powerful metabolomics approach. The metabolome of P. vulgaris root nodules infected by a P. phymatumrpoN Fix⁻ mutant revealed statistically significant metabolic changes compared to wild-type Fix⁺ nodules, including reduced amounts of chorismate and elevated levels of flavonoids. A transcriptome analysis on Fix⁻ and Fix⁺ nodules—combined with a search for RpoN binding sequences in promoter regions of regulated genes—confirmed the expected control of σ⁵⁴ on nitrogen fixation genes in nodules. The transcriptomic data also allowed us to identify additional target genes, whose differential expression was able to explain the observed metabolite changes in numerous cases. Moreover, the genes encoding the two-component regulatory system NtrBC were downregulated in root nodules induced by the rpoN mutant, and contained a putative RpoN binding motif in their promoter region, suggesting direct regulation. The construction and characterization of an ntrB mutant strain revealed impaired nitrogen assimilation in free-living conditions, as well as a noticeable symbiotic phenotype, as fewer but heavier nodules were formed on P. vulgaris roots.

Genome-Wide Identification and Classification of Soybean C2H2 Zinc Finger Proteins and Their Expression Analysis in Legume-Rhizobium Symbiosis

Root nodule symbiosis (RNS) is one of the most productive and economical systems for nitrogen fixation, and previous studies have shown that several nodule-specific C2H2-zinc finger proteins (ZFPs) play important roles in symbiosis establishment and nodule function. However, C2H2-ZFPs are the most widespread ZFPs in eukaryotes, and a great variation of structure and function exist among the family members. It remains largely unclear whether or not special types of C2H2-ZF genes participate in symbiosis, especially in soybean. In the present study, we performed a genome-wide survey of soybean C2H2-ZF genes, and 321 soybean C2H2-ZF genes were identified and classified into 11 clearly distinguishable subsets (Gm-t1-SF, Gm-t2-SF, Gm-1i-Q-SF, Gm-1i-M-SF, Gm-1i-Z-SF, Gm-1i-D-SF, Gm-2i-Q-SF, Gm-2i-M-SF, Gm-2i-Mix-SF, Gm-3i-SF, and Gm-4i-SF) based on the arrangements, numbers, and types of C2H2-ZF domains. Phylogenetic and gene ontology analyses were carried out to assess the conserved sequence and GO function among these subsets, and the results showed that the classification of soybean C2H2-ZFPs was reasonable. The expression profile of soybean C2H2-ZFPs in multiple tissues showed that nearly half of soybean C2H2-ZFPs within different subsets had expressions in nodules, including a clustering branch consisting of 11 Gm-1i-Q-SF genes specifically expressed in symbiotic-relative tissues. RNA-Seq was used to identify symbiosis-related soybean C2H2-ZFPs, and the expression pattern of the soybean C2H2-ZFPs in roots and nodules at different development stages showed that soybean C2H2-ZFPs mainly played roles in nodule development or nodule function rather than nodulation signal transduction, and nearly half of these genes had high expressions and/or different expression patterns during soybean nodule development, especially for the six clustering branches of genes consisting of different subsets of C2H2-ZFPs. Furthermore, the selected symbiosis-related soybean C2H2-ZFPs might function in legume-rhizobium symbiosis through regulating or interacting with other key proteins. Taken together, our findings provided useful information for the study on classification and conservative function of C2H2-ZFPs, and offered solid evidence for investigation of rhizobium symbiosis-related C2H2-ZFPs in soybean or other legumes.

RNA sequencing and analysis of three Lupinus nodulomes provide new insights into specific host-symbiont relationships with compatible and incompatible Bradyrhizobium strains

Nitrogen fixation in the legume root-nodule symbiosis has a critical importance in natural and agricultural ecosystems and depends on the proper choice of the symbiotic partners. However, the genetic determinism of symbiotic specificity remains unclear. To study this process, we inoculated three Lupinus species (L. albus, L. luteus, L. mariae-josephae), belonging to the under-investigated tribe of Genistoids, with two Bradyrhizobium strains (B. japonicum, B. valentinum) presenting contrasted degrees of symbiotic specificity depending on the host. We produced the first transcriptomes (RNA-Seq) from lupine nodules in a context of symbiotic specificity. For each lupine species, we compared gene expression between functional and non-functional interactions and determined differentially expressed (DE) genes. This revealed that L. luteus and L. mariae-josephae (nodulated by only one of the Bradyrhizobium strains) specific nodulomes were richest in DE genes than L. albus (nodulation with both microsymbionts, but non-functional with B. valentinum) and share a higher number of these genes between them than with L. albus. In addition, a functional analysis of DE genes highlighted the central role of the genetic pathways controlling infection and nodule organogenesis, hormones, secondary, carbon and nitrogen metabolisms, as well as the implication of plant defence in response to compatible or incompatible Bradyrhizobium strains.

Aberrant Meiotic Prophase I Leads to Genic Male Sterility in the Novel TE5A Mutant of Brassica napus

Genic male sterility (GMS) has already been extensively utilized for hybrid rapeseed production. TE5A is a novel thermo-sensitive dominant GMS line in Brassica napus, however, its mechanisms of GMS remain largely unclear. Histological and Transmission electron microscopy (TEM) analyses of anthers showed that the male gamete development of TE5A was arrested at meiosis prophase I. EdU uptake of S-phase meiocytes revealed that the TE5A mutant could accomplish DNA replication, however, chromosomal and fluorescence in situ hybridization (FISH) analyses of TE5A showed that homologous chromosomes could not pair, synapse, condense and form bivalents. We then analyzed the transcriptome differences between young floral buds of sterile plants and its near-isogenic fertile plants through RNA-Seq. A total of 3,841 differentially expressed genes (DEGs) were obtained, some of which were associated with homologous chromosome behavior and cell cycle control during meiosis. Dynamic expression changes of selected candidate DEGs were then analyzed at different anther developmental stages. The present study not only demonstrated that the TE5A mutant had defects in meiotic prophase I via detailed cytological analysis, but also provided a global insight into GMS-associated DEGs and elucidated the mechanisms of GMS in TE5A through RNA-Seq.

Search for Nodulation and Nodule Development-Related Cystatin Genes in the Genome of Soybean (Glycine max)

Nodulation, nodule development and senescence directly affects nitrogen fixation efficiency, and previous studies have shown that inhibition of some cysteine proteases delay nodule senescence, so their nature inhibitors, cystatin genes, are very important in nodulation, nodule development, and senescence. Although several cystatins are actively transcribed in soybean nodules, their exact roles and functional diversities in legume have not been well explored in genome-wide survey studies. In this report, we performed a genome-wide survey of cystatin family genes to explore their relationship to nodulation and nodule development in soybean and identified 20 cystatin genes that encode peptides with 97-245 amino acid residues, different isoelectric points (pI) and structure characteristics, and various putative plant regulatory elements in 3000 bp putative promoter fragments upstream of the 20 soybean cystatins in response to different abiotic/biotic stresses, hormone signals, and symbiosis signals. The expression profiles of these cystatin genes in soybean symbiosis with rhizobium strain Bradyrhizobium japonicum strain 113-2 revealed that 7 cystatin family genes play different roles in nodulation as well as nodule development and senescence. However, these genes were not root nodule symbiosis (RNS)-specific and did not encode special clade cystatin protein with structures related to nodulation and nodule development. Besides, only two of these soybean cystatins were not upregulated in symbiosis after ABA treatment. The functional analysis showed that a candidate gene Glyma.15G227500 (GmCYS16) was likely to play a positive role in soybean nodulation. Besides, evolutionary relationships analysis divided the cystatin genes from Arabidopsis thaliana, Nicotiana tabacum, rice, barley and four legume plants into three groups. Interestingly, Group A cystatins are special in legume plants, but only include one of the above-mentioned 7 cystatin genes related to nodulation and nodule development. Overall, our results provide useful information or clues for our understanding of the functional diversity of legume cystatin family proteins in soybean nodulation and nodule development and for finding nodule-specific cysteine proteases in soybean.