Best practices for the execution, analysis, and data storage of plant single-cell/nucleus transcriptomics

Single-cell and single-nucleus RNA-sequencing technologies capture the expression of plant genes at an unprecedented resolution. Therefore, these technologies are gaining traction in plant molecular and developmental biology for elucidating the transcriptional changes across cell types in a specific tissue or organ, upon treatments, in response to biotic and abiotic stresses, or between genotypes. Despite the rapidly accelerating use of these technologies, collective and standardized experimental and analytical procedures to support the acquisition of high-quality data sets are still missing. In this commentary, we discuss common challenges associated with the use of single-cell transcriptomics in plants and propose general guidelines to improve reproducibility, quality, comparability, and interpretation and to make the data readily available to the community in this fast-developing field of research.

Identification and characterization of a temperature sensitive chlorotic soybean mutant

Screening a transposon-mutagenized soybean population led to the discovery of a recessively inherited chlorotic phenotype. This "vir1" phenotype results in smaller stature, weaker stems, and a smaller root system with smaller nodules. Genome sequencing identified 15 candidate genes with mutations likely to result in a loss of function. Amplicon sequencing of a segregating population was then used to narrow the list to a single candidate mutation, a single-base change in Glyma.07G102300 that disrupts splicing of the second intron. Single cell transcriptomic profiling indicates that this gene is expressed primarily in mesophyll cells and RNA sequencing data indicates it is upregulated in germinating seedlings by cold stress. Previous studies have shown that mutations to Os05g34040, the rice homolog of Glyma.07G102300, produced a chlorotic phenotype that was more pronounced in cool temperatures. Growing soybean vir1 mutants at lower temperatures also resulted in a more severe phenotype. In addition, transgenic expression of wild type Glyma.07G102300 in the knockout mutant of the Arabidopsis homolog At4930720 rescues the chlorotic phenotype, further supporting the hypothesis that the mutation in Glyma.07G102300 is causal of the vir1 phenotype.

The roles of a novel CDKB/KRP/FB3 cell cycle core complex in rice gametes and initiation of embryogenesis

The cell cycle controls division and proliferation of all eukaryotic cells and is tightly regulated at multiple checkpoints by complexes of core cell cycle proteins. Due to the difficulty in accessing female gametes and zygotes of flowering plants, little is known about the molecular mechanisms underlying embryogenesis initiation despite the crucial importance of this process for seed crops. In this study, we reveal three levels of factors involved in rice zygotic cell cycle control and characterize their functions and regulation. Protein-protein interaction studies, including within zygote cells, and in vitro biochemical analyses delineate a model of the zygotic cell cycle core complex for rice. In this model, CDKB1, a major regulator of plant mitosis, is a cyclin (CYCD5)-dependent kinase; its activity is coordinately inhibited by two cell cycle inhibitors, KRP4 and KRP5; and both KRPs are regulated via F-box protein 3 (FB3)-mediated proteolysis. Supporting their critical roles in controlling the rice zygotic cell cycle, mutations in KRP4, KRP5 and FB3 result in the compromised function of sperm cells and abnormal organization of female germ units, embryo and endosperm, thus significantly reducing seed-set rate. This work helps reveal regulatory mechanisms controlling the zygotic cell cycle toward seed formation in angiosperms.

Plant Single-Cell/Nucleus RNA-seq Workflow

Single-cell transcriptomics technologies allow researchers to investigate how individual cells, in complex multicellular organisms, differentially use their common genomic DNA. In plant biology, these technologies were recently applied to reveal the transcriptomes of various plant cells isolated from different organs and different species and in response to environmental stresses. These first studies support the potential of single-cell transcriptomics technology to decipher the biological function of plant cells, their developmental programs, cell-type-specific gene networks, programs controlling plant cell response to environmental stresses, etc. In this chapter, we provide information regarding the critical steps and important information to consider when developing an experimental design in plant single-cell biology. We also describe the current status of bioinformatics tools used to analyze single-cell RNA-seq datasets and how additional emerging technologies such as spatial transcriptomics and long-read sequencing technologies will provide additional information on the differential use of the genome by plant cells.

Cell-specific pathways recruited for symbiotic nodulation in the Medicago truncatula legume

Medicago truncatula is a model legume species that has been studied for decades to understand the symbiotic relationship between legumes and soil bacteria collectively named rhizobia. This symbiosis called nodulation is initiated in roots with the infection of root hair cells by the bacteria, as well as the initiation of nodule primordia from root cortical, endodermal, and pericycle cells, leading to the development of a new root organ, the nodule, where bacteria fix and assimilate the atmospheric dinitrogen for the benefit of the plant. Here, we report the isolation and use of the nuclei from mock and rhizobia-inoculated roots for the single nuclei RNA-seq (sNucRNA-seq) profiling to gain a deeper understanding of early responses to rhizobial infection in Medicago roots. A gene expression map of the Medicago root was generated, comprising 25 clusters, which were annotated as specific cell types using 119 Medicago marker genes and orthologs to Arabidopsis cell-type marker genes. A focus on root hair, cortex, endodermis, and pericycle cell types, showing the strongest differential regulation in response to a short-term (48 h) rhizobium inoculation, revealed not only known genes and functional pathways, validating the sNucRNA-seq approach, but also numerous novel genes and pathways, allowing a comprehensive analysis of early root symbiotic responses at a cell type-specific level.

Review: Challenges and perspectives in applying single nuclei RNA-seq technology in plant biology

Plant single-cell RNA-seq technology quantifies the abundance of plant transcripts at a single-cell resolution. Deciphering the transcriptomes of each plant cell, their regulation during plant cell development, and their response to environmental stresses will support the functional study of genes, the establishment of precise transcriptional programs, the prediction of more accurate gene regulatory networks, and, in the long term, the design of de novo gene pathways to enhance selected crop traits. In this review, we will discuss the opportunities, challenges, and problems, and share tentative solutions associated with the generation and analysis of plant single-cell transcriptomes. We will discuss the benefit and limitations of using plant protoplasts vs. nuclei to conduct single-cell RNA-seq experiments on various plant species and organs, the functional annotation of plant cell types based on their transcriptomic profile, the characterization of the dynamic regulation of the plant genes during cell development or in response to environmental stress, the need to characterize and integrate additional layers of -omics datasets to capture new molecular modalities at the single-cell level and reveal their causalities, the deposition and access to single-cell datasets, and the accessibility of this technology to plant scientists.

Cell-Type-Specific Profiling of the Arabidopsis thaliana Membrane Protein-Encoding Genes

Membrane proteins work in large complexes to perceive and transduce external signals and to trigger a cellular response leading to the adaptation of the cells to their environment. Biochemical assays have been extensively used to reveal the interaction between membrane proteins. However, such analyses do not reveal the unique and complex composition of the membrane proteins of the different plant cell types. Here, we conducted a comprehensive analysis of the expression of Arabidopsis membrane proteins in the different cell types composing the root. Specifically, we analyzed the expression of genes encoding membrane proteins interacting in large complexes. We found that the transcriptional profiles of membrane protein-encoding genes differ between Arabidopsis root cell types. This result suggests that different cell types are characterized by specific sets of plasma membrane proteins, which are likely a reflection of their unique biological functions and interactions. To further explore the complexity of the Arabidopsis root cell membrane proteomes, we conducted a co-expression analysis of genes encoding interacting membrane proteins. This study confirmed previously reported interactions between membrane proteins, suggesting that the co-expression of genes at the single cell-type level can be used to support protein network predictions.

First Plant Cell Atlas symposium report

The Plant Cell Atlas (PCA) community hosted a virtual symposium on December 9 and 10, 2021 on single cell and spatial omics technologies. The conference gathered almost 500 academic, industry, and government leaders to identify the needs and directions of the PCA community and to explore how establishing a data synthesis center would address these needs and accelerate progress. This report details the presentations and discussions focused on the possibility of a data synthesis center for a PCA and the expected impacts of such a center on advancing science and technology globally. Community discussions focused on topics such as data analysis tools and annotation standards; computational expertise and cyber-infrastructure; modes of community organization and engagement; methods for ensuring a broad reach in the PCA community; recruitment, training, and nurturing of new talent; and the overall impact of the PCA initiative. These targeted discussions facilitated dialogue among the participants to gauge whether PCA might be a vehicle for formulating a data synthesis center. The conversations also explored how online tools can be leveraged to help broaden the reach of the PCA (i.e., online contests, virtual networking, and social media stakeholder engagement) and decrease costs of conducting research (e.g., virtual REU opportunities). Major recommendations for the future of the PCA included establishing standards, creating dashboards for easy and intuitive access to data, and engaging with a broad community of stakeholders. The discussions also identified the following as being essential to the PCA's success: identifying homologous cell-type markers and their biocuration, publishing datasets and computational pipelines, utilizing online tools for communication (such as Slack), and user-friendly data visualization and data sharing. In conclusion, the development of a data synthesis center will help the PCA community achieve these goals by providing a centralized repository for existing and new data, a platform for sharing tools, and new analytical approaches through collaborative, multidisciplinary efforts. A data synthesis center will help the PCA reach milestones, such as community-supported data evaluation metrics, accelerating plant research necessary for human and environmental health.

Isolating and cryopreserving pig skin cells for single-cell RNA sequencing study

The pig skin architecture and physiology are similar to those of humans. Thus, the pig model is very valuable for studying skin biology and testing therapeutics. The single-cell RNA sequencing (scRNA-seq) technology allows quantitatively analyzing cell types, compositions, states, signaling, and receptor-ligand interactome at single-cell resolution and at high throughput. scRNA-seq has been used to study mouse and human skins. However, studying pig skin with scRNA-seq is still rare. A critical step for successful scRNA-seq is to obtain high-quality single cells from the pig skin tissue. Here we report a robust method for isolating and cryopreserving pig skin single cells for scRNA-seq. We showed that pig skin could be efficiently dissociated into single cells with high cell viability using the Miltenyi Human Whole Skin Dissociation kit and the Miltenyi gentleMACS Dissociator. Furthermore, the obtained single cells could be cryopreserved using 90% FBS + 10% DMSO without causing additional cell death, cell aggregation, or changes in gene expression profiles. Using the developed protocol, we were able to identify all the major skin cell types. The protocol and results from this study are valuable for the skin research scientific community.

Vision, challenges and opportunities for a Plant Cell Atlas

With growing populations and pressing environmental problems, future economies will be increasingly plant-based. Now is the time to reimagine plant science as a critical component of fundamental science, agriculture, environmental stewardship, energy, technology and healthcare. This effort requires a conceptual and technological framework to identify and map all cell types, and to comprehensively annotate the localization and organization of molecules at cellular and tissue levels. This framework, called the Plant Cell Atlas (PCA), will be critical for understanding and engineering plant development, physiology and environmental responses. A workshop was convened to discuss the purpose and utility of such an initiative, resulting in a roadmap that acknowledges the current knowledge gaps and technical challenges, and underscores how the PCA initiative can help to overcome them.

Enhancing Our Understanding of Plant Cell-to-Cell Interactions Using Single-Cell Omics

Plants are composed of cells that physically interact and constantly adapt to their environment. To reveal the contribution of each plant cells to the biology of the entire organism, their molecular, morphological, and physiological attributes must be quantified and analyzed in the context of the morphology of the plant organs. The emergence of single-cell/nucleus omics technologies now allows plant biologists to access different modalities of individual cells including their epigenome and transcriptome to reveal the unique molecular properties of each cell composing the plant and their dynamic regulation during cell differentiation and in response to their environment. In this manuscript, we provide a perspective regarding the challenges and strategies to collect plant single-cell biological datasets and their analysis in the context of cellular interactions. As an example, we provide an analysis of the transcriptional regulation of the Arabidopsis genes controlling the differentiation of the root hair cells at the single-cell level. We also discuss the perspective of the use of spatial profiling to complement existing plant single-cell omics.

Single-nucleus RNA and ATAC sequencing reveals the impact of chromatin accessibility on gene expression in Arabidopsis roots at the single-cell level

Similar to other complex organisms, plants consist of diverse and specialized cell types. The gain of unique biological functions of these different cell types is the consequence of the establishment of cell-type-specific transcriptional programs. As a necessary step in gaining a deeper understanding of the regulatory mechanisms controlling plant gene expression, we report the use of single-nucleus RNA sequencing (sNucRNA-seq) and single-nucleus assay for transposase accessible chromatin sequencing (sNucATAC-seq) technologies on Arabidopsis roots. The comparison of our single-nucleus transcriptomes to the published protoplast transcriptomes validated the use of nuclei as biological entities to establish plant cell-type-specific transcriptomes. Furthermore, our sNucRNA-seq results uncovered the transcriptomes of additional cell subtypes not identified by single-cell RNA-seq. Similar to our transcriptomic approach, the sNucATAC-seq approach led to the distribution of the Arabidopsis nuclei into distinct clusters, suggesting the differential accessibility of chromatin between groups of cells according to their identity. To reveal the impact of chromatin accessibility on gene expression, we integrated sNucRNA-seq and sNucATAC-seq data and demonstrated that cell-type-specific marker genes display cell-type-specific patterns of chromatin accessibility. Our data suggest that the differential chromatin accessibility is a critical mechanism to regulate gene activity at the cell-type level.

Isolation of Plant Root Nuclei for Single Cell RNA Sequencing

The characterization of the transcriptional similarities and differences existing between plant cells and cell types is important to better understand the biology of each cell composing the plant, to reveal new molecular mechanisms controlling gene activity, and to ultimately implement meaningful strategies to enhance plant cell biology. To gain a deeper understanding of the regulation of plant gene activity, the individual transcriptome of each plant cell needs to be established. Until recently, single cell approaches were mostly limited to bulk transcriptomic studies on selected cell types. Accessing specific cell types required the development of labor-intensive strategies. Recently, single cell sequencing strategies were successfully applied on isolated Arabidopsis thaliana root protoplasts. However, this strategy relies on the successful isolation of viable protoplasts upon the optimization of the enzymatic cocktails required to digest the cell wall and on the compatibility of fragile plant protoplasts with the use of microfluidic systems to generate single cell transcriptomic libraries. To overcome these difficulties, we present a simple and fast alternative strategy: the isolation and use of plant nuclei to access meaningful transcriptomic information from plant cells. This protocol was specifically developed to enable the use of the plant nuclei with 10× Genomics' Chromium technology partitions technology. Briefly, the plant nuclei are released from the root by chopping into a nuclei isolation buffer before purification by filtration then nuclei sorting. Upon sorting, the nuclei are resuspended in a low divalent ion buffer compatible with the Chromium technology in order to create single nuclei ribonucleic acid-sequencing libraries (sNucRNA-seq). © 2020 Wiley Periodicals LLC. Basic Protocol 1: Arabidopsis seed sterilization and planting Basic Protocol 2: Nuclei isolation from Arabidopsis roots Basic Protocol 3: Fluorescent-activated nuclei sorting (FANS) purification Support Protocol: Estimation of nuclei density using Countess II automated cell counter Alternate Protocol 1: Proper growth conditions for Medicago truncatula and Sorghum bicolor Alternate Protocol 2: Estimation of nuclei density using sNucRNA-seq technology.

Single cell RNA sequence analysis reveals novel clusters unique to post-infectious myocarditis in the mouse model of Coxsackievirus infection

Myocarditis is a common occurrence in children and adolescents that can lead to dilated cardiomyopathy. Since it is difficult to identify the triggers of this complex disease, mechanistic understanding of DCM pathogenesis may create avenues to develop new treatment strategies. We have been engaged in delineating the autoimmune pathways underlying the development of DCM using various mouse models of myocarditis that resemble human disease. One such model is Coxsackievirus B3 (CVB)-induced myocarditis in A/J mice that show two distinct phases of disease, but occur in continuum. While, the acute, infectious phase (up to ~10 days post-infection) is marked by viremia and infiltrations in the hearts, infectious virions become undetected during the chronic phase (~beyond 18 days post-infection), and animals develop DCM. The latter phenotype can occur over a period of ~60 to 90 days, and hitherto, it has been a challenge to determine the underlying mechanisms. To dissect this complex nature of disease, we recently performed single cell RNA sequence analysis on heart infiltrates obtained from CVB-infected mice. Using panels of 15 to 30 marker genes for various immune and non-immune cells, we identified 18 and 26 clusters respectively in healthy and infected hearts that also include few unknown clusters unique for each group. Likewise, transcriptome analyses led us to identify a total of 14,051 genes to be differentially expressed, and of these, 206 and 345 transcripts were unique to healthy and CVB-infected groups, respectively. The data suggest that uniquely expressed genes in the CVB-infected hearts may have a role in the development of post-infectious myocarditis leading to DCM.

Single cell RNA sequence analysis reveals novel clusters unique to post-infectious myocarditis in the mouse model of Coxsackievirus infection

Myocarditis is a common occurrence in children and adolescents that can lead to dilated cardiomyopathy. Since it is difficult to identify the triggers of this complex disease, mechanistic understanding of DCM pathogenesis may create avenues to develop new treatment strategies. We have been engaged in delineating the autoimmune pathways underlying the development of DCM using various mouse models of myocarditis that resemble human disease. One such model is Coxsackievirus B3 (CVB)-induced myocarditis in A/J mice that show two distinct phases of disease, but occur in continuum. While, the acute, infectious phase (up to ~10 days post-infection) is marked by viremia and infiltrations in the hearts, infectious virions become undetected during the chronic phase (~beyond 18 days post-infection), and animals develop DCM. The latter phenotype can occur over a period of ~60 to 90 days, and hitherto, it has been a challenge to determine the underlying mechanisms. To dissect this complex nature of disease, we recently performed single cell RNA sequence analysis on heart infiltrates obtained from CVB-infected mice. Using panels of 15 to 30 marker genes for various immune and non-immune cells, we identified 18 and 26 clusters respectively in healthy and infected hearts that also include few unknown clusters unique for each group. Likewise, transcriptome analyses led us to identify a total of 14,051 genes to be differentially expressed, and of these, 206 and 345 transcripts were unique to healthy and CVB-infected groups, respectively. The data suggest that uniquely expressed genes in the CVB-infected hearts may have a role in the development of post-infectious myocarditis leading to DCM.

Biological and Cellular Functions of the Microdomain-Associated FWL/CNR Protein Family in Plants

Membrane microdomains/nanodomains are sub-compartments of the plasma membrane enriched in sphingolipids and characterized by their unique protein composition. They play important roles in regulating plant development and plant-microbe interactions including mutualistic symbiotic interactions. Several protein families are associated with the microdomain fraction of biological membranes such as flotillins, prohibitins, and remorins. More recently, GmFWL1, a FWL/CNR protein exclusively expressed in the soybean nodule, was functionally characterized as a new microdomain-associated protein. Interestingly, GmFWL1 is homologous to the tomato FW2-2 protein, a major regulator of tomato fruit development. In this review, we summarize the knowledge gained about the biological, cellular, and physiological functions of members of the FWL/CNR family across various plant species. The role of the FWL/CNR proteins is also discussed within the scope of their evolution and transcriptional regulation.

Plant Hormones Differentially Control the Sub-Cellular Localization of Plasma Membrane Microdomains during the Early Stage of Soybean Nodulation

Phytohormones regulate the mutualistic symbiotic interaction between legumes and rhizobia, nitrogen-fixing soil bacteria, notably by controlling the formation of the infection thread in the root hair (RH). At the cellular level, the formation of the infection thread is promoted by the translocation of plasma membrane microdomains at the tip of the RH. We hypothesize that phytohormones regulate the translocation of plasma membrane microdomains to regulate infection thread formation. Accordingly, we treated with hormone and hormone inhibitors transgenic soybean roots expressing fusions between the Green Fluorescent Protein (GFP) and GmFWL1 or GmFLOT2/4, two microdomain-associated proteins translocated at the tip of the soybean RH in response to rhizobia. Auxin and cytokinin treatments are sufficient to trigger or inhibit the translocation of GmFWL1 and GmFLOT2/4 to the RH tip independently of the presence of rhizobia, respectively. Unexpectedly, the application of salicylic acid, a phytohormone regulating the plant defense system, also promotes the translocation of GmFWL1 and GmFLOT2/4 to the RH tip regardless of the presence of rhizobia. These results suggest that phytohormones are playing a central role in controlling the early stages of rhizobia infection by regulating the translocation of plasma membrane microdomains. They also support the concept of crosstalk of phytohormones to control nodulation.

Transcriptional Reprogramming of Legume Genomes: Perspective and Challenges Associated With Single-Cell and Single Cell-Type Approaches During Nodule Development

Transcriptomic approaches revealed thousands of genes differentially or specifically expressed during nodulation, a biological process resulting from the symbiosis between leguminous plant roots and rhizobia, atmospheric nitrogen-fixing symbiotic bacteria. Ultimately, nodulation will lead to the development of a new root organ, the nodule. Through functional genomic studies, plant transcriptomes have been used by scientists to reveal plant genes potentially controlling nodulation. However, it is important to acknowledge that the physiology, transcriptomic programs, and biochemical properties of the plant cells involved in nodulation are continuously regulated. They also differ between the different cell-types composing the nodules. To generate a more accurate picture of the transcriptome, epigenome, proteome, and metabolome of the cells infected by rhizobia and cells composing the nodule, there is a need to implement plant single-cell and single cell-types strategies and methods. Accessing such information would allow a better understanding of the infection of plant cells by rhizobia and will help understanding the complex interactions existing between rhizobia and the plant cells. In this mini-review, we are reporting the current knowledge on legume nodulation gained by plant scientists at the level of single cell-types, and provide perspectives on single cell/single cell-type approaches when applied to legume nodulation.

Enhancing Phenotyping and Molecular Analysis of Plant Root System Using Ultrasound Aeroponic Technology

Several plant growth systems are available to enhance the observation of the root system (e.g., hydroponic and aeroponic plant growth systems, use of transparent soils, etc.). This article describes the use of the ultrasound aeroponic system (USAS) to treat and to enhance access to the root systems of various model plant and crop species (e.g., Arabidopsis thaliana, Medicago truncatula, soybean, etc.). This system is also compatible with short-term (hr) and long-term (days/weeks) biotic and abiotic treatments of plants. Upon treatment, the ease of access to the plant root system facilitates phenotyping (e.g., analysis of root architecture, establishment of root light spectrum using remote sensing technology), microscopic, molecular, and biochemical experiments. In addition, to facilitate functional genomic studies, we combined the use of the USAS with the hairy root transformation system to grow and observe transgenic roots on composite legume plants. © 2018 by John Wiley & Sons, Inc.

Phosphate Deficiency Negatively Affects Early Steps of the Symbiosis between Common Bean and Rhizobia

Phosphate (Pi) deficiency reduces nodule formation and development in different legume species including common bean. Despite significant progress in the understanding of the genetic responses underlying the adaptation of nodules to Pi deficiency, it is still unclear whether this nutritional deficiency interferes with the molecular dialogue between legumes and rhizobia. If so, what part of the molecular dialogue is impaired? In this study, we provide evidence demonstrating that Pi deficiency negatively affects critical early molecular and physiological responses that are required for a successful symbiosis between common bean and rhizobia. We demonstrated that the infection thread formation and the expression of PvNSP2, PvNIN, and PvFLOT2, which are genes controlling the nodulation process were significantly reduced in Pi-deficient common bean seedlings. In addition, whole-genome transcriptional analysis revealed that the expression of hormones-related genes is compromised in Pi-deficient seedlings inoculated with rhizobia. Moreover, we showed that regardless of the presence or absence of rhizobia, the expression of PvRIC1 and PvRIC2, two genes participating in the autoregulation of nodule numbers, was higher in Pi-deficient seedlings compared to control seedlings. The data presented in this study provides a mechanistic model to better understand how Pi deficiency impacts the early steps of the symbiosis between common bean and rhizobia.

Plant Systems Biology at the Single-Cell Level

Our understanding of plant biology is increasingly being built upon studies using 'omics and system biology approaches performed at the level of the entire plant, organ, or tissue. Although these approaches open new avenues to better understand plant biology, they suffer from the cellular complexity of the analyzed sample. Recent methodological advances now allow plant scientists to overcome this limitation and enable biological analyses of single-cells or single-cell-types. Coupled with the development of bioinformatics and functional genomics resources, these studies provide opportunities for high-resolution systems analyses of plant phenomena. In this review, we describe the recent advances, current challenges, and future directions in exploring the biology of single-cells and single-cell-types to enhance our understanding of plant biology as a system.

Function of plasma membrane microdomain-associated proteins during legume nodulation

Plasma membrane microdomains are plasma membrane sub-compartments enriched in sphingolipids and sterols, and composed by a specific set of proteins. They are involved in recognizing signal molecules, transducing these signals, and controlling endocytosis and exocytosis processes. In a recent study, applying biochemical and microscopic methods, we characterized the soybean GmFWL1 protein, a major regulator of soybean nodulation, as a new membrane microdomain-associated protein. Interestingly, upon rhizobia inoculation of the soybean root system, GmFWL1 and one of its interacting partners, GmFLOT2/4, both translocate to the root hair cell tip, the primary site of interaction and infection between soybean and Rhizobium. The role of GmFWL1 as a plasma membrane microdomain-associated protein is also supported by immunoprecipitation assays performed on soybean nodules, which revealed 178 GmFWL1 protein partners including a large number of microdomain-associated proteins such as GmFLOT2/4. In this addendum, we provide additional information about the identity of the soybean proteins repetitively identified as GmFWL1 protein partners. Their function is discussed especially in regard to plant-microbe interactions and microbial symbiosis. This addendum will provide new insights in the role of plasma membrane microdomains in regulating legume nodulation.

Plant response to biotic stress: Is there a common epigenetic response during plant-pathogenic and symbiotic interactions?

Plants constantly interact with pathogenic and symbiotic microorganisms. Recent studies have revealed several regulatory mechanisms controlling these interactions. Among them, the plant defense system is activated not only in response to pathogenic, but also in response to symbiotic microbes. Interestingly, shortly after symbiotic microbial recognition, the plant defense system is suppressed to promote plant infection by symbionts. Research studies have demonstrated the influence of the plant epigenome in modulating both pathogenic and symbiotic plant-microbe interactions, thereby influencing plant survival, adaptation and evolution of the plant response to microbial infections. It is however unclear if plant pathogenic and symbiotic responses share similar epigenomic profiles or if epigenomic changes differentially regulate plant-microbe symbiosis and pathogenesis. In this mini-review, we provide an update of the current knowledge of epigenomic control on plant immune responses and symbiosis, with a special attention being paid to knowledge gap and potential strategies to fill-in the missing links.

The GmFWL1 (FW2-2-like) nodulation gene encodes a plasma membrane microdomain-associated protein

The soybean gene GmFWL1 (FW2-2-like1) belongs to a plant-specific family that includes the tomato FW2-2 and the maize CNR1 genes, two regulators of plant development. In soybean, GmFWL1 is specifically expressed in root hair cells in response to rhizobia and in nodules. Silencing of GmFWL1 expression significantly reduced nodule numbers supporting its role during soybean nodulation. While the biological role of GmFWL1 has been described, its molecular function and, more generally, the molecular function of plant FW2-2-like proteins is unknown. In this study, we characterized the role of GmFWL1 as a membrane microdomain-associated protein. Specifically, using biochemical, molecular and cellular methods, our data show that GmFWL1 interacts with various proteins associated with membrane microdomains such as remorin, prohibitins and flotillins. Additionally, comparative genomics revealed that GmFWL1 interacts with GmFLOT2/4 (FLOTILLIN2/4), the soybean ortholog to Medicago truncatula FLOTILLIN4, a major regulator of the M. truncatula nodulation process. We also observed that, similarly to MtFLOT4 and GmFLOT2/4, GmFWL1 was localized at the tip of the soybean root hair cells in response to rhizobial inoculation supporting the early function of GmFWL1 in the rhizobium infection process.

A comparative genomic and transcriptomic analysis at the level of isolated root hair cells reveals new conserved root hair regulatory elements

KEY MESSAGE

A comparative transcriptomic and genomic analysis between Arabidopsis thaliana and Glycine max root hair genes reveals the evolution of the expression of plant genes after speciation and whole genome duplication. Our understanding of the conservation and divergence of the expression patterns of genes between plant species is limited by the quality of the genomic and transcriptomic resources available. Specifically, the transcriptomes generated from plant organs are the reflection of the contribution of the different cell types composing the samples weighted by their relative abundances in the sample. These contributions can vary between plant species leading to the generation of datasets which are difficult to compare. To gain a deeper understanding of the evolution of gene transcription in and between plant species, we performed a comparative transcriptomic and genomic analysis at the level of one single plant cell type, the root hair cell, and between two model plants: Arabidopsis (Arabidopsis thaliana) and soybean (Glycine max). These two species, which diverged 90 million years ago, were selected as models based on the large amount of genomic and root hair transcriptomic information currently available. Our analysis revealed in detail the transcriptional divergence and conservation between soybean paralogs (i.e., the soybean genome is the product of two successive whole genome duplications) and between Arabidopsis and soybean orthologs in this single plant cell type. Taking advantage of this evolutionary study, we combined bioinformatics, molecular, cellular and microscopic tools to characterize plant promoter sequences and the discovery of two root hair regulatory elements (RHE1 and RHE2) consistently and specifically active in plant root hair cells.

Comprehensive Comparative Genomic and Transcriptomic Analyses of the Legume Genes Controlling the Nodulation Process

Nitrogen is one of the most essential plant nutrients and one of the major factors limiting crop productivity. Having the goal to perform a more sustainable agriculture, there is a need to maximize biological nitrogen fixation, a feature of legumes. To enhance our understanding of the molecular mechanisms controlling the interaction between legumes and rhizobia, the symbiotic partner fixing and assimilating the atmospheric nitrogen for the plant, researchers took advantage of genetic and genomic resources developed across different legume models (e.g., Medicago truncatula, Lotus japonicus, Glycine max, and Phaseolus vulgaris) to identify key regulatory protein coding genes of the nodulation process. In this study, we are presenting the results of a comprehensive comparative genomic analysis to highlight orthologous and paralogous relationships between the legume genes controlling nodulation. Mining large transcriptomic datasets, we also identified several orthologous and paralogous genes characterized by the induction of their expression during nodulation across legume plant species. This comprehensive study prompts new insights into the evolution of the nodulation process in legume plant and will benefit the scientific community interested in the transfer of functional genomic information between species.

Identification and functional characterization of soybean root hair microRNAs expressed in response to Bradyrhizobium japonicum infection

Three soybean [Glycine max (L) Merr.] small RNA libraries were generated and sequenced using the Illumina platform to examine the role of miRNAs during soybean nodulation. The small RNA libraries were derived from root hairs inoculated with Bradyrhizobium japonicum (In_RH) or mock-inoculated with water (Un_RH), as well as from the comparable inoculated stripped root samples (i.e. inoculated roots with the root hairs removed). Sequencing of these libraries identified a total of 114 miRNAs, including 22 novel miRNAs. A comparison of miRNA abundance among the 114 miRNAs identified 66 miRNAs that were differentially expressed between root hairs and stripped roots, and 48 miRNAs that were differentially regulated in infected root hairs in response to B. japonicum when compared to uninfected root hairs (P ≤ 0.05). A parallel analysis of RNA ends (PARE) library was constructed and sequenced to reveal a total of 405 soybean miRNA targets, with most predicted to encode transcription factors or proteins involved in protein modification, protein degradation and hormone pathways. The roles of gma-miR4416 and gma-miR2606b during nodulation were further analysed. Ectopic expression of these two miRNAs in soybean roots resulted in significant changes in nodule numbers. miRNA target information suggested that gma-miR2606b regulates a Mannosyl-oligosaccharide 1, 2-alpha-mannosidase gene, while gma-miR4416 regulates the expression of a rhizobium-induced peroxidase 1 (RIP1)-like peroxidase gene, GmRIP1, during nodulation.

Xyloglucan, galactomannan, glucuronoxylan, and rhamnogalacturonan I do not have identical structures in soybean root and root hair cell walls

MAIN CONCLUSION

Chemical analyses and glycome profiling demonstrate differences in the structures of the xyloglucan, galactomannan, glucuronoxylan, and rhamnogalacturonan I isolated from soybean ( Glycine max ) roots and root hair cell walls. The root hair is a plant cell that extends only at its tip. All other root cells have the ability to grow in different directions (diffuse growth). Although both growth modes require controlled expansion of the cell wall, the types and structures of polysaccharides in the walls of diffuse and tip-growing cells from the same plant have not been determined. Soybean (Glycine max) is one of the few plants whose root hairs can be isolated in amounts sufficient for cell wall chemical characterization. Here, we describe the structural features of rhamnogalacturonan I, rhamnogalacturonan II, xyloglucan, glucomannan, and 4-O-methyl glucuronoxylan present in the cell walls of soybean root hairs and roots stripped of root hairs. Irrespective of cell type, rhamnogalacturonan II exists as a dimer that is cross-linked by a borate ester. Root hair rhamnogalacturonan I contains more neutral oligosaccharide side chains than its root counterpart. At least 90% of the glucuronic acid is 4-O-methylated in root glucuronoxylan. Only 50% of this glycose is 4-O-methylated in the root hair counterpart. Mono O-acetylated fucose-containing subunits account for at least 60% of the neutral xyloglucan from root and root hair walls. By contrast, a galacturonic acid-containing xyloglucan was detected only in root hair cell walls. Soybean homologs of the Arabidopsis xyloglucan-specific galacturonosyltransferase are highly expressed only in root hairs. A mannose-rich polysaccharide was also detected only in root hair cell walls. Our data demonstrate that the walls of tip-growing root hairs cells have structural features that distinguish them from the walls of other roots cells.

A Comparative Epigenomic Analysis of Polyploidy-Derived Genes in Soybean and Common Bean

Soybean (Glycine max) and common bean (Phaseolus vulgaris) share a paleopolyploidy (whole-genome duplication [WGD]) event, approximately 56.5 million years ago, followed by a genus Glycine-specific polyploidy, approximately 10 million years ago. Cytosine methylation is an epigenetic mark that plays an important role in the regulation of genes and transposable elements (TEs); however, the role of DNA methylation in the fate/evolution of genes following polyploidy and speciation has not been fully explored. Whole-genome bisulfite sequencing was used to produce nucleotide resolution methylomes for soybean and common bean. We found that, in soybean, CG body-methylated genes were abundant in WGD genes, which were, on average, more highly expressed than single-copy genes and had slower evolutionary rates than unmethylated genes, suggesting that WGD genes evolve more slowly than single-copy genes. CG body-methylated genes were also enriched in shared single-copy genes (single copy in both species) that may be responsible for the broad and high expression patterns of this class of genes. In addition, diverged methylation patterns in non-CG contexts between paralogs were due mostly to TEs in or near genes, suggesting a role for TEs and non-CG methylation in regulating gene expression post polyploidy. Reference methylomes for both soybean and common bean were constructed, providing resources for investigating epigenetic variation in legume crops. Also, the analysis of methylation patterns of duplicated and single-copy genes has provided insights into the functional consequences of polyploidy and epigenetic regulation in plant genomes.

Identification of microRNAs and their mRNA targets during soybean nodule development: functional analysis of the role of miR393j-3p in soybean nodulation

Plant microRNAs (miRNAs) play important regulatory roles in a number of developmental processes. The present work investigated the roles of miRNAs during nodule development in the crop legume soybean (Glycine max). Fifteen soybean small RNA libraries were sequenced from different stages of nodule development, including young nodules, mature nodules and senescent nodules. In order to identify the regulatory targets of the miRNAs, five parallel analysis of RNA ends (PARE) libraries were also sequenced from the same stages of nodule development. Sequencing identified 284 miRNAs, including 178 novel soybean miRNAs. Analysis of miRNA abundance identified 139 miRNAs whose expression was significantly regulated during nodule development, including 12 miRNAs whose expression changed > 10-fold. Analysis of the PARE libraries identified 533 miRNA targets, including three nodulation-related genes and eight nodule-specific genes. miR393j-3p was selected for detailed analysis as its expression was significantly regulated during nodule formation, and it targeted a nodulin gene, Early Nodulin 93 (ENOD93). Strong, ectopic expression of miR393j-3p, as well as RNAi silencing of ENOD93 expression, significantly reduced nodule formation. The data indicate that miR393j-3p regulation of ENOD93 mRNA abundance is a key control point for soybean nodule formation.

miR172 regulates soybean nodulation

Micro-RNAs (miRNAs) play a pivotal role in the control of gene expression and regulate plant developmental processes. miRNA 172 (miR172) is a conserved miRNA in plants reported to control the expression of genes involved in developmental phase transition, floral organ identity, and flowering time. However, the specific role of miR172 in legume nodulation is undefined. Ectopic expression of soybean miR172 resulted in an increase in nodule numbers in transgenic roots and an increase in the expression of both symbiotic leghemoglobin and nonsymbiotic hemoglobin. These nodules showed higher levels of nitrogenase activity. Further analysis revealed a complex regulatory circuit in which miR156 regulates miR172 expression and controls the level of an AP2 transcription factor. The latter, either directly or indirectly, controls the expression of nonsymbiotic hemoglobin, which is essential for regulating the levels of nodulation.

The fate of duplicated genes in a polyploid plant genome

Polyploidy is generally not tolerated in animals, but is widespread in plant genomes and may result in extensive genetic redundancy. The fate of duplicated genes is poorly understood, both functionally and evolutionarily. Soybean (Glycine max L.) has undergone two separate polyploidy events (13 and 59 million years ago) that have resulted in 75% of its genes being present in multiple copies. It therefore constitutes a good model to study the impact of whole-genome duplication on gene expression. Using RNA-seq, we tested the functional fate of a set of approximately 18 000 duplicated genes. Across seven tissues tested, approximately 50% of paralogs were differentially expressed and thus had undergone expression sub-functionalization. Based on gene ontology and expression data, our analysis also revealed that only a small proportion of the duplicated genes have been neo-functionalized or non-functionalized. In addition, duplicated genes were often found in collinear blocks, and several blocks of duplicated genes were co-regulated, suggesting some type of epigenetic or positional regulation. We also found that transcription factors and ribosomal protein genes were differentially expressed in many tissues, suggesting that the main consequence of polyploidy in soybean may be at the regulatory level.

Unleashing the potential of the root hair cell as a single plant cell type model in root systems biology

Plant root is an organ composed of multiple cell types with different functions. This multicellular complexity limits our understanding of root biology because -omics studies performed at the level of the entire root reflect the average responses of all cells composing the organ. To overcome this difficulty and allow a more comprehensive understanding of root cell biology, an approach is needed that would focus on one single cell type in the plant root. Because of its biological functions (i.e., uptake of water and various nutrients; primary site of infection by nitrogen-fixing bacteria in legumes), the root hair cell is an attractive single cell model to study root cell response to various stresses and treatments. To fully study their biology, we have recently optimized procedures in obtaining root hair cell samples. We culture the plants using an ultrasound aeroponic system maximizing root hair cell density on the entire root systems and allowing the homogeneous treatment of the root system. We then isolate the root hair cells in liquid nitrogen. Isolated root hair yields could be up to 800 to 1000~mg of plant cells from 60 root systems. Using soybean as a model, the purity of the root hair was assessed by comparing the expression level of genes previously identified as soybean root hair specific between preparations of isolated root hair cells and stripped roots, roots devoid in root hairs. Enlarging our tests to include other plant species, our results support the isolation of large quantities of highly purified root hair cells which is compatible with a systems biology approach.

14-3-3 proteins SGF14c and SGF14l play critical roles during soybean nodulation

The soybean (Glycine max) genome contains 18 members of the 14-3-3 protein family, but little is known about their association with specific phenotypes. Here, we report that the Glyma0529080 Soybean G-box Factor 14-3-3c (SGF14c) and Glyma08g12220 (SGF14l) genes, encoding 14-3-3 proteins, appear to play essential roles in soybean nodulation. Quantitative reverse transcription-polymerase chain reaction and western-immunoblot analyses showed that SGF14c mRNA and protein levels were specifically increased in abundance in nodulated soybean roots 10, 12, 16, and 20 d after inoculation with Bradyrhizobium japonicum. To investigate the role of SGF14c during soybean nodulation, RNA interference was employed to silence SGF14c expression in soybean roots using Agrobacterium rhizogenes-mediated root transformation. Due to the paleopolyploid nature of soybean, designing a specific RNA interference sequence that exclusively targeted SGF14c was not possible. Therefore, two highly similar paralogs (SGF14c and SGF14l) that have been shown to function as dimers were silenced. Transcriptomic and proteomic analyses showed that mRNA and protein levels were significantly reduced in the SGF14c/SGF14l-silenced roots, and these roots exhibited reduced numbers of mature nodules. In addition, SGF14c/SGF14l-silenced roots contained large numbers of arrested nodule primordia following B. japonicum inoculation. Transmission electron microscopy further revealed that the host cytoplasm and membranes, except the symbiosome membrane, were severely degraded in the failed nodules. Altogether, transcriptomic, proteomic, and cytological data suggest a critical role of one or both of these 14-3-3 proteins in early development stages of soybean nodules.

Quantitative phosphoproteomic analysis of soybean root hairs inoculated with Bradyrhizobium japonicum

Root hairs are single hair-forming cells on roots that function to increase root surface area, enhancing water and nutrient uptake. In leguminous plants, root hairs also play a critical role as the site of infection by symbiotic nitrogen fixing rhizobia, leading to the formation of a novel organ, the nodule. The initial steps in the rhizobia-root hair infection process are known to involve specific receptor kinases and subsequent kinase cascades. Here, we characterize the phosphoproteome of the root hairs and the corresponding stripped roots (i.e. roots from which root hairs were removed) during rhizobial colonization and infection to gain insight into the molecular mechanism of root hair cell biology. We chose soybean (Glycine max L.), one of the most important crop plants in the legume family, for this study because of its larger root size, which permits isolation of sufficient root hair material for phosphoproteomic analysis. Phosphopeptides derived from root hairs and stripped roots, mock inoculated or inoculated with the soybean-specific rhizobium Bradyrhizobium japonicum, were labeled with the isobaric tag eight-plex iTRAQ, enriched using Ni-NTA magnetic beads and subjected to nanoRPLC-MS/MS1 analysis using HCD and decision tree guided CID/ETD strategy. A total of 1625 unique phosphopeptides, spanning 1659 nonredundant phosphorylation sites, were detected from 1126 soybean phosphoproteins. Among them, 273 phosphopeptides corresponding to 240 phosphoproteins were found to be significantly regulated (>1.5-fold abundance change) in response to inoculation with B. japonicum. The data reveal unique features of the soybean root hair phosphoproteome, including root hair and stripped root-specific phosphorylation suggesting a complex network of kinase-substrate and phosphatase-substrate interactions in response to rhizobial inoculation.

Identification of soybean proteins from a single cell type: the root hair

Root hairs (RH) are a terminally differentiated single cell type, mainly involved in water and nutrient uptake from the soil. The soybean RH cell represents an excellent model for the study of single cell systems biology. In this study, we identified 5702 proteins, with at least two peptides, from soybean RH using an accurate mass and time tag approach, establishing a comprehensive proteome reference map of this single cell type. We also showed that trypsin is the most appropriate enzyme for soybean proteomic studies by performing an in silico digestion of the soybean proteome using different proteases. Although the majority of proteins identified in this study are involved in basal metabolism, the function of others are more related to RH formation/function and include proteins involved in nutrient uptake (transporters) or vesicular trafficking (cytoskeleton and ras-associated binding proteins). Interestingly, some of these proteins appear to be specifically detected in RH and constitute promising candidates for further studies to elucidate unique features of this single-cell model.

Soybean Knowledge Base (SoyKB): a web resource for soybean translational genomics

BACKGROUND

Soybean Knowledge Base (SoyKB) is a comprehensive all-inclusive web resource for soybean translational genomics. SoyKB is designed to handle the management and integration of soybean genomics, transcriptomics, proteomics and metabolomics data along with annotation of gene function and biological pathway. It contains information on four entities, namely genes, microRNAs, metabolites and single nucleotide polymorphisms (SNPs).

METHODS

SoyKB has many useful tools such as Affymetrix probe ID search, gene family search, multiple gene/metabolite search supporting co-expression analysis, and protein 3D structure viewer as well as download and upload capacity for experimental data and annotations. It has four tiers of registration, which control different levels of access to public and private data. It allows users of certain levels to share their expertise by adding comments to the data. It has a user-friendly web interface together with genome browser and pathway viewer, which display data in an intuitive manner to the soybean researchers, producers and consumers.

CONCLUSIONS

SoyKB addresses the increasing need of the soybean research community to have a one-stop-shop functional and translational omics web resource for information retrieval and analysis in a user-friendly way. SoyKB can be publicly accessed at http://soykb.org/.

Identification of quantitative trait loci controlling gene expression during the innate immunity response of soybean

Microbe-associated molecular pattern-triggered immunity (MTI) is an important component of the plant innate immunity response to invading pathogens. However, most of our knowledge of MTI comes from studies of model systems with relatively little work done with crop plants. In this work, we report on variation in both the microbe-associated molecular pattern-triggered oxidative burst and gene expression across four soybean (Glycine max) genotypes. Variation in MTI correlated with the level of pathogen resistance for each genotype. A quantitative trait locus analysis on these traits identified four loci that appeared to regulate gene expression during MTI in soybean. Likewise, we observed that both MTI variation and pathogen resistance were quantitatively inherited. The approach utilized in this study may have utility for identifying key resistance loci useful for developing improved soybean cultivars.

A dual-targeted soybean protein is involved in Bradyrhizobium japonicum infection of soybean root hair and cortical cells

The symbiotic interaction between legumes and soil bacteria (e.g., soybean [Glycine max L.] and Bradyrhizobium japonicum]) leads to the development of a new root organ, the nodule, where bacteria differentiate into bacteroids that fix atmospheric nitrogen for assimilation by the plant host. In exchange, the host plant provides a steady carbon supply to the bacteroids. This carbon can be stored within the bacteroids in the form of poly-3-hydroxybutyrate granules. The formation of this symbiosis requires communication between both partners to regulate the balance between nitrogen fixation and carbon utilization. In the present study, we describe the soybean gene GmNMNa that is specifically expressed during the infection of soybean cells by B. japonicum. GmNMNa encodes a protein of unknown function. The GmNMNa protein was localized to the nucleolus and also to the mitochondria. Silencing of GmNMNa expression resulted in reduced nodulation, a reduction in the number of bacteroids per infected cell in the nodule, and a clear reduction in the accumulation of poly-3-hydroxybutyrate in the bacteroids. Our results highlight the role of the soybean GmNMNa gene in regulating symbiotic bacterial infection, potentially through the regulation of the accumulation of carbon reserves.

Enzymatic activity of the soybean ecto-apyrase GS52 is essential for stimulation of nodulation

Nitrogen is an essential nutrient for plant growth. In the Rhizobium-legume symbiosis, root nodules are the sites of bacterial nitrogen fixation, in which atmospheric nitrogen is converted into a form that plants can utilize. While recent studies suggested an important role for the soybean (Glycine max) ecto-apyrase GS52 in rhizobial root hair infection and root nodule formation, precisely how this protein impacts the nodulation process remains undetermined. In this study, the biochemical characteristics of the GS52 enzyme were investigated. Computer modeling of the GS52 apyrase structure identified key amino acid residues important for catalytic activity, which were subsequently mutagenized. Although the GS52 enzyme exhibited broad substrate specificity, its activity on pyrimidine nucleotides and diphosphate nucleotides was significantly higher than on ATP. This result was corroborated by structural modeling of GS52, which predicted a low specificity for the adenine base within the substrate-binding pocket of the enzyme. The wild-type enzyme and its inactive mutant forms were expressed in soybean roots in order to evaluate the importance of GS52 enzymatic activity for nodulation. The results indicated a clear correlation between GS52 enzymatic activity and nodule number. Altogether, our study indicates that the catalytic activity of the GS52 apyrase, likely acting on extracellular nucleotides, is critical for rhizobial infection and nodulation.

Root hair systems biology

Plant functional genomic studies have largely measured the response of whole plants, organs and tissues, resulting in the dilution of the signal from individual cells. Methods are needed where the full repertoire of functional genomic tools can be applied to a single plant cell. Root hair cells are an attractive model to study the biology of a single, differentiated cell type because of their ease of isolation, polar growth, and role in water and nutrient uptake, as well as being the site of infection by nitrogen-fixing bacteria. This review highlights the recent advances in our understanding of plant root hair biology and examines whether the root hair has potential as a model for plant cell systems biology.

Evolution of FW2.2-like (FWL) and PLAC8 genes in eukaryotes

The tomato FW2.2 quantitative trait locus, which regulates tomato fruit size, was genetically and physically mapped around 15 years ago. Subsequently, the FW2.2 gene was cloned and shown to contain a PLAC8 domain, originally identified in mammalian placental proteins. Data suggest that FW2.2 likely controls tomato cell size, perhaps by direct interaction with casein kinase II. Several FW2.2-like (FWL) genes have now been identified from a variety of plant species, but until recently only the tomato FW2.2 gene had been the subject of detailed investigation. Recently, soybean and maize FWL genes were identified and shown to have a role in plant organogenesis. It is now apparent that the FWL genes in plants are a large gene family, which is even larger given inclusion of genes for the various eukaryotic PLAC8-domain proteins. Although overall the protein sequence identity/similarity among the family members is relatively low, there is strong conservation of key domains, suggesting a conservation of the core biochemical function of these proteins. In this Addendum Article, we highlight the similarities and differences exiting between plant FWL genes and enlarge this comparison to the mammalian PLAC8 genes. These comparisons suggest the possible conservation of biological function for FWL proteins.

Soybean metabolites regulated in root hairs in response to the symbiotic bacterium Bradyrhizobium japonicum

Nodulation of soybean (Glycine max) root hairs by the nitrogen-fixing symbiotic bacterium Bradyrhizobium japonicum is a complex process coordinated by the mutual exchange of diffusible signal molecules. A metabolomic study was performed to identify small molecules produced in roots and root hairs during the rhizobial infection process. Metabolites extracted from roots and root hairs mock inoculated or inoculated with B. japonicum were analyzed by gas chromatography-mass spectrometry and ultraperformance liquid chromatography-quadrupole time of flight-mass spectrometry. These combined approaches identified 2,610 metabolites in root hairs. Of these, 166 were significantly regulated in response to B. japonicum inoculation, including various (iso)flavonoids, amino acids, fatty acids, carboxylic acids, and various carbohydrates. Trehalose was among the most strongly induced metabolites produced following inoculation. Subsequent metabolomic analyses of root hairs inoculated with a B. japonicum mutant defective in the trehalose synthase, trehalose 6-phosphate synthase, and maltooligosyltrehalose synthase genes showed that the trehalose detected in the inoculated root hairs was primarily of bacterial origin. Since trehalose is generally considered an osmoprotectant, these data suggest that B. japonicum likely experiences osmotic stress during the infection process, either on the root hair surface or within the infection thread.

An integrated transcriptome atlas of the crop model Glycine max, and its use in comparative analyses in plants

Soybean (Glycine max L.) is a major crop providing an important source of protein and oil, which can also be converted into biodiesel. A major milestone in soybean research was the recent sequencing of its genome. The sequence predicts 69,145 putative soybean genes, with 46,430 predicted with high confidence. In order to examine the expression of these genes, we utilized the Illumina Solexa platform to sequence cDNA derived from 14 conditions (tissues). The result is a searchable soybean gene expression atlas accessible through a browser (http://digbio.missouri.edu/soybean_atlas). The data provide experimental support for the transcription of 55,616 annotated genes and also demonstrate that 13,529 annotated soybean genes are putative pseudogenes, and 1736 currently unannotated sequences are transcribed. An analysis of this atlas reveals strong differences in gene expression patterns between different tissues, especially between root and aerial organs, but also reveals similarities between gene expression in other tissues, such as flower and leaf organs. In order to demonstrate the full utility of the atlas, we investigated the expression patterns of genes implicated in nodulation, and also transcription factors, using both the Solexa sequence data and large-scale qRT-PCR. The availability of the soybean gene expression atlas allowed a comparison with gene expression documented in the two model legume species, Medicago truncatula and Lotus japonicus, as well as data available for Arabidopsis thaliana, facilitating both basic and applied aspects of soybean research.