Micro RNA regulation of nodule zone-specific gene expression in soybean

Article history: Received 20210302 Received in revised form 20210312 Accepted 20210322 Available online 20210403 Nitrogen is a paramount important essential element for all living organisms. It has been found to bea crucial structural component of proteins, nucleic acids, enzymes and other cellular constituents which are inevitable for all forms of life. In the atmosphere, the percentage of nitrogen is very high (N2, 78%) compared to other inorganic gases. However, most organisms have practically no direct access to this nitrogen. While plants can not directly uptake nitrogen from atmosphere, they are capable of assimilating other forms of nitrogen, for example ammonium (NH4+) and nitrate (NO3-). For agricultural crop production, artificial fixation of nitrogen is heavily utilized and it is an expensive process that requires high temperatures (at least 400 °C) and pressures (around 200 atm). It has been conspicuously demonstrated that indiscriminate use of fertilizer hampers soil physical, chemical and micro biological properties and also a potential risk to environment e.g. water quality. Besides, chemically manufactured fertilizers are depleted from soils in various ways, for instance; denitrifying bacteria, volatilization, and leaching. Consequently, it results relatively poor availability of nitrogen to get into plants. On the flipside, only 1-2% of the nitrogen fixation in the world occurs through the natural process of lightening. Notably, microbial fixation is well characterized in diazotrophs for example; Rhizobia and Frankia, and blue-green algae. Against the backdrop, we are accentuated on an environmentally friendly and the most sustainable approach to increase productivity for legume and non-legume crops. Till today, the term biological nitrogen fixation (BNF) has received much attention as a sustainable alternative; this process facilitates atmospheric nitrogen to convert into ammonia by rhizobia in specialized plan organs termed “root nodules”. This review article seeks to better understand plant mechanisms involved in the development of root nodules in soybean.


Introduction
Soybean (Glycine max) is one of the most important oil crops and a source of animal feed protein in the world. It has a salient feature to fix atmospheric nitrogen through symbioses with compatible rhizobia that yields to determinate type nodule (Oldroyd, Murray et al. 2011). Biological nitrogen fixation in soybean nodules reduces the use of chemical nitrogen fertilizers resulting in cost-savings to producers and minimizes environmental damage due to nitrogen run-off. A better understanding of how nodules form and function is important for selection or generation of soybean genotypes with better nitrogen fixation capacity. Soybean nodules originate from root cortex via de novo cell differentiation (Oldroyd 2013). Consequently, two major nodule development zones are formed for instance; the nodule primordium (Npr) in the middle and it is encircled by nodule parenchyma (Npa). At later time point, the Npr gives rise to Nfixation zone and the Npa holds vascular bundles. It is not clear what early signaling pathways driving the conspicuous development of the nodule zones. My research is aimed at filling this knowledge gap by illustrating the molecular signatures that paves the way to cellular differentiation in root nodule development in soybean. Based on initial evidence obtained by the Subramanian lab, we hypothesize that microRNAs (miRNAs) play important regulatory roles in spatio-temporal expression of their target genes during nodule developmental in soybean. For instance, the regulation of auxin sensitivity by miR160 has been found to be crucial for formation of nodule primordia and vasculature in the parenchyma (Marie Turner 2 2013). Against this backdrop, this review article focused on nuclear and cytoplasmic transcriptome as well as miRNA profiles of parenchyma and primordial tissues and determine the relative abundance and differentially expressed mRNAs and regulatory role of miRNAs in cell differentiation and nodule development.
Root nodule a sustainable alternative to fix atmospheric nitrogen Atmospheric nitrogen percentage is very high (N2, 78%) compared to other inorganic gases (Mary Elvira 1932). However, most of the organisms have practically no direct access to this nitrogen. Nevertheless, plants can not directly uptake nitrogen from atmosphere but they are capable of assimilating only very specific forms of nitrogen, for example ammonium (NH4+) and nitrate (NO3-) (Bytnerowicz and Fenn 1996, Peter M. Vitousek 1997) ( (Johnson 1996, Peter M. Vitousek 1997. Consequently, it results relatively poor availability of nitrogen to get into plants. On the flipside, over 90 % of the nitrogen fixation in the world occurs through the natural process of lightening and microorganisms. Furthermore, microbial fixation is well characterized in diazotrophs for example; Rhizobia and Frankia, and blue-green algae (Cheng 2008). It has been demonstrated that Bradyrhizobium strains substantially escalated soybean grain yield, and protein content up to 57% and 26%, respectively (Zimmer, Messmer et al. 2016). Against the backdrop, we are accentuated on an environmentally friendly and a sustainable approach to increase the productivity for legume and nonlegume crops. Literature mining depicted that biological nitrogen fixation in soybean nodules reduces the use of chemical nitrogen fertilizers resulting in cost-savings to producers and minimizes environmental damage due to nitrogen run-off.

Rhizobia infection leads to the root nodule development
In the natural environment, plants are continuously confronted with pathogenic and symbiotic microbes. Symbioses involves mutual exchange of diffusible signal molecules, first endophytic bacteria (rhizobia) are attracted by the plant root exudates flavonoids which are perceived and triggered the bacterial nodulation (nod) genes. Consequently, the bacteria synthesize specific lipochito-oligosaccharides, called nodulation (Nod) factors. This signal is perceived by the LysM receptor like kinase of host plant, it induces the root hair curling, and bacteria get access into the host epidermis through infection threads (ITs) and initiate cell division within the root cortex, leading to the progression of the root nodule meristem. In later stages of the interaction, bacteria are released from the infection threads into the plant cells, surrounded by membrane of plant origin. These bacteria multiply within the host cells and differentiate into the nitrogen fixing bacteroids (Udvardi and Day 1997) (Oldroyd 2013). Till now, integration of genetic and genomic approaches has revealed twenty-six genes to be involved in nodule development of Medicago truncatualaand Lotus japonicum (Kouchi, Imaizumi-Anraku et al. 2010). In addition, deep sequencing of the Medicago truncatularoot transcriptome has uncovered thousands of genes to be induced during Nod factor signaling and its resulting ethylene (ET) biosynthesis throughout the multiple development stages of indeterminate nodule (Larrainzar, Riely et al. 2015). Albeit the molecular mechanism of such regulation is not well understood. There has been a large-scale transcriptome analysis of B. japonicum-inoculated and mock-inoculated soybean root hairs. It has showed that a total of 1,973 soybean genes differentially expressed during root hair infection, particularly NFR5 and NIN genes (Libault, Farmer et al. 2010). Nevertheless, the signaling mechanisms directing the cellular differentiation of nodule are not known.

Soybean root nodule organogenesis
Soybean (Glycine max) has a genome size of 1.1 to1.5 Gb, it is partially diploidized tetraploid. It is one of the most important oil crops and a source of animal feed protein in the world (soybase.org/sb_about.php). It has a salient feature to fix atmospheric nitrogen through symbioses with compatible rhizobia that yields to determinate type nodule (Udvardi and Day 1997) (Oldroyd, Murray et al. 2011). Notwithstanding of the economic and environmental importance, there has been very few studies about quantitative trait loci (QTL) that controlling BNF traits, for instance nodule number, ration of nodule dry weight with nodule number, and shoot dry weight (SDW). It has been reported via composite interval mapping that approximately six QTLs bears very small effect on BNF traits (Santos, Geraldi et al. 2013). Besides, it has been demonstrated in earlier studies that nodules originate from root cortex via de novo cell differentiation into two different cell types, parenchymal and primordium (Celine Charon 1997) (Oldroyd&Downie 2008; Oldroyd 2013). In addition, early nodulin genes in legume for instance; Enod 40 gene reported to be expressed in root pericycle during the rhizobia infection and later it occupied in the dividing cortical cells (H. Kouchi and S. Hata 1993). Among the two major nodule development zones, the nodule primordium (Npr) in the middle which is encircled by nodule parenchyma (Npa). At later time point, the Npr gives rise to N-fixation zone and the Npa holds vascular bundles. Lately, a β-expansin gene, GmEXPB2 fused with GUS reporter gene which was observed to be preferentially expressed in nodule vascular trace and nodule vascular bundles. It indicated that GmEXPB2 might be crucial for nodule organogenesis. Over expression of GmEXPB2 contrast to suppressed transgenic lines found to be escalated nodule numb mass and nitrogenase activity. It further suggested that GmEXPB2 might have influenced over root architecture, nodule formation and development, and profoundly yielding to biological N2 fixation (Li, Zhao et al. 2015). Even though, it is not clear what early signaling pathways driving the conspicuous development of the nodule zones. Against the back drop, to understand the regulation of auxin sensitivity by miR160 which is believed to be crucial for the formation of nodule primordia (Marie Turner 2013).

Regulatory small RNAs biogenesis and its molecular functions
Regulatory small RNAs are ranged between 20 to 24 nucleotides which are ubiquitous elements of endogenous plant transcriptomics, a common response to exogenous viral infections and introduced double-stranded RNA (Axtell 2013). Three core enzymes families, for instance; RNAdependent RNA polymerase (RDR), Dicer like (DCL), and Argonaute (AGO) proteins paves the way of small RNA biogenesis and function in plants. Firstly, ribonuclease type III or DICERLIKE1 involves in the yield of a fold-back precursor RNA or primary miRNA (primiRNA) transcripts using an RNA templates in the nuclei. Later, the resulting miRNA-miRNA duplex which is originated in nucleus then translocated into cytoplasm. The guided miRNAmolecule is incorporated into ARGONAUTE (AGO) to Journal of Natural Products and Natural Products Synthesis 3 vascular trace and nodule vascular bundles. It indicated that GmEXPB2 might be crucial for nodule organogenesis. Over GmEXPB2 contrast to suppressed GmEXPB2 transgenic lines found to be escalated nodule number, nodule activity. It further suggested that over root architecture, nodule formation and development, and profoundly yielding to 2015). Even though, it is s driving the conspicuous ones. Against the back drop, to understand the regulation of auxin sensitivity by miR160 which formation of nodule primordia

Regulatory small RNAs biogenesis and its molecular
Regulatory small RNAs are ranged between 20 to 24 of endogenous plant transcriptomics, a common response to exogenous viral stranded RNA (Axtell 2013). Three core enzymes families, for instance; RNAdependent RNA polymerase (RDR), Dicer like (DCL), and Argonaute (AGO) ins paves the way of small RNA biogenesis and function in plants. Firstly, ribonuclease type III or DICERLIKE1 involves in back precursor RNA or primary miRNA (primiRNA) transcripts using an RNA templates in the nuclei. miRNA duplex which is originated in nucleus then translocated into cytoplasm. The guided miRNAmolecule is incorporated into ARGONAUTE (AGO) to form an active RISC complex to specific target RNAs that are complementary to the miRNA, and this proc follows up mRNA cleavage, represses the translation of the mRNAs or Chromatin modification. This phenomenon accentuated as an inhibition or silencing of the gene expression, which play a crucial role in the developmental process in plant and animal (Chapman and Carrington 2007) (Axtell 2013).  There are several miRNA families that are conserved across the vast evolutionary distances from flowering plants to mosses (Cuperus, Fahlgren et al. 2011). It has been observed in another study that miRNA, and its target pairing found to be stable for a prolonged periods of plant evolution. On the flip side, another group demonstrated that conserved plant miRNAs and their targets are to somehow flexible. For instance; miR159 is a highly conserved miRNA that targets not only a subset of MYB mRNAs but also observed to target a non MYB mRNA, SGN-U567133 (Buxdorf, Hendelman et al. 2010). A mutant tomato transgenic line (miR159-resistant line) showed higher level of the SGN-U567133 transcript and exhibited defects in leaf and flower development. This result suggests that miR159 involves in a post-transcriptional regulation. Additionally, it is found to be crucial for the normal tomato development.
Recently, the identification of miRNAs in the regulation of photoperiodic pathways in soybean have been reported through high throughput sequencing and qRT-PCR. Six libraries were constructed using Illumina Solexa, for instance; 0, 8, and 16 h under short day treatment, similar time points considered for the long the long day treatment. A total of 163 miRNAs families were reported which covered 318 plant miRNAs, and unclassified 81 novel predicted miRNAs. As expected, significant differences in abundance between short day and long day treatment was observed (Wenbin Li 2015). These findings provided evidence of miRNA in the regulation of flowering time that ultimately affects the seed yield and quality of soybean.
The complex regulatory network of miRNA-mRNA interactions during viral infection has been revealed via small RNA seq (sRNA), degradome seq, and genome-wide transcriptome analysis. There has been a total of 253 soybean miRNAs found to be two-folds abundance compared with mockinoculated control demonstrated through sRNA seq analysis. Among them 105 miRNAs were identified as potential targets of 125 transcripts that has been validated by degradome seq analyses. In addition, 2679 genes were detected via genome wide transcriptomic analysis. These genes have been differentially expressed during infection of soybean mosaic virus and among them 71 genes projected to induce in defense response (Hui Chen 2016). These findings suggested the regulatory role miRNA that governed the target gene expression during viral infection. Furthermore, the regulatory role of microRNAs (miRNAs) during Soybean-Bradyrhizobium japonicum mutualistic association was studied first by Subramanian et al. 2008. They sequenced approximately 350000 small RNAs of soybean root sample which were inoculated with B. japonicum. It helps to detect 20 conserved miRNAs loci based on the similarity to miRNAs in another plant species. In addition, 35 novel miRNAs were identified based on potential hairpin forming precursors in Soybean EST as well as shotgun genomic sequences (Subramanian, Fu et al. 2008). These findings advocated the potential role of miRNAs in the regulation of legumerhizobiumsymbiosis.
In another study, 120 hairpin-forming precursor genes have been identified in soybean by Turner et al. In addition, they reported three novel miRNAs for instance; miR160, miR164 and miR393 found to be involved in auxin signaling (Turner, Yu et al. 2012). Moreover, the plant hormone auxin is thought to have a pivotal role in nodule organogenesis in determinate and indeterminate type of nodule. It indicates a redundancy and diversity of miRNAs family members that governs the formation of root nodule.
It has been illustrated that auxin receptor gene family hushed by over expressed microRNA393. These plant roots found to be hypersensitive to auxin and yielded normal nodule. This observation advocated that only minimal/reduced auxin signaling is required for determinate nodule development. Likewise, overexpressed microRNA160 hushed a set of repressor auxin response transcription factor. These plant roots were hypersensitive to auxin and observed not to be reluctant in epidermal responses to rhizobia. Notably, it yielded to lower sized nodule primordium (Marie Turner 2013). This observation indicated that auxin hypersensitivity inhibits nodule organogenesis Organ specific expression of profile of miRNA and the potential targets were also studied. Two genes (Glyma10g10240 and Glyma17g05920) which were the target of miR169 but detected to be highly expressed in soybean nodule. Likewise, three potential targets of gma-new-miR13587 demonstrated to be highly expressed in the nodules than in the roots. As expected, gma-newmiR13587 found to be poorly expressed in the nodules than in the roots (Turner, Yu et al. 2012).
There was an inverse expression pattern observed in between roots and nodules. Li et al., studied the transgene expression of three novel miRNAs namely, miR482, miR1512, and miR1515 in Soybean. They noticed a significant increase of nodule numbers while root length and later root density were normal in all tested miRNA lines. As expected, there were differential expression of these miRNAs in supernodulating and nonnodulating soybean mutants. They reported that 6 novel miRNAs decoyed 22 predicted target genes. And it was estimated via real time polymerase chain reaction and qRT-PCR (Li, Deng et al. 2010). It advocates that miRNAs have the signatory roles in soybean nodule development.
Sequencing of small RNAs and Parallel analysis of RNA ends (PARE) libraries revealed to identify 284 nodule miRNAs, more than 500 target genes, and including 178 novel soybean miRNAs. It has been reported that ENOD93 only found to be expressed in nodule tissue not in other plant parts of Soybean. Ectopic expression of miR393j-3p and RNAi silencing approach to ENOD93 expression showed a significant reduction in nodule formation (Zhe Yan 2015). Therefore, this study showed a list of miRNAs and their potential target of nodulation genes. In the model legume (Medicago truncatula), 25 conserved miRNA families and 100 novel miRNA reads were detected by highthroughput sequencing. The expression of MtHAP2-1 (encodes a CCAAT binding transcription factor) to meristematic zones was restricted by miR169a which is found to be critical for the development of indeterminate type of nodule (Combier, Frugier et al. 2006). In another study, HDZIPIII transcripts were inhibited by overexpression of miR166, it dropped the number of symbiotic nodule and lateral root (Boualem, Laporte et al. 2008). To get insights into key genes of nodule zones, transcript profiles of specific cells/tissues were investigated at different time points from indeterminate nodules of M. truncatulausing laser capture micro dissection. It has been demonstrated from the comprehensive gene expression map that selected genes enriched in different cell/tissue types (Limpens, Moling et al. 2013). These findings indicated that organ specific gene expression could be controlled by the presence or absence of miRNAs.
Recently, Agrobacterium rhizogenesmediated hairy root transformation has been applied as tool for exploring cell type specific gene expression in tomato. Cell type or tissue specific promoter introduced into INTACT and TRAP constructs via gateway cloning technology to develop binary vectors. INTACT method used to capture biotin tagged nuceli from specific cell types and TRAP method used for profiling of mRNAs or foot printing of individual ribosomes (Ron 2014). TRAP methodology is not required tissue fixation or single cell suspension. It has been successfully used to date in organisms ranging from D. melanogaster to mice and human cultured cells. Multiple ribosomes or Polyribosomes (polysomes) are engaged in translation on a single mRNA. To evaluate the translation state of an mRNA, ribosomal subunits, ribosomes, and polysomes can be isolated from detergent-treated cell extracts (Heiman, Kulicke et al. 2014). In this study, we would perform polysome isolation deploying gene cassettes ENOD40p:HF-GFP-RPL18 for primordial tissues, and ENOD2p:HF-GFP-RPL18 for parenchymal tissues in Glycine max root nodules that express an epitope tagged version of ribosomal protein L18.
Over the last one decade, there has been several microarrays-based studies which characterized transcriptional variations deployed in nodule formation. It has been embedded with couple of shortcomings, for instance; relative late time points study, incomplete representation of plant genes,discrimination of close paralogs, and reduced sensitivity. Lately, next generation sequencingtechnology have widened the horizon of transcription analyses in different legume species to detectsymbiosis induced changes in late nodule developmental stages. Against this backdrop, we areaccentuated to reveal early transcriptional changes induced in determinate type of soybean noduleby Bradyrhizobium japonicum.
In determinate type of nodule, two major nodule development zones are formed for instance, the nodule primordium (Npr) in the middle and it is encircled by nodule parenchyma (Npa). At later time point, the Npr converted to Nfixation zone and the Npa contained vascular bundles. Of these facts, it is not clear what early signaling pathways driving the conspicuous development of thenodule zones. In this context, mechanisms regulate the distinct gene expression profiles in Npr andNpa cell types has not understood clearly. The proposed research study is aimed at filling this knowledge gap byillustrating the molecular signatures that paves the way to cellular differentiation in root noduledevelopment in soybean considering four different time points (5 dai, 7 dai, 10 dai& 14 dai).
The hypothesisis microRNAs(miRNAs) play important regulatory roles in spatio-temporal expression of their target genesduring nodule developmental in soybean. For example, a gradient of microRNA localizationbetween nodule primordium and parenchyma cells could result in distinct differentiation of thesecell types. To test this hypothesis, one has to obtain both cell type-specific miRNA andtranscriptome (miRNA target) profiles. Since, the majority of miRNA regulation occurs in thecytoplasm, we reasoned that comparison of nuclear and ribosomal transcriptome profiles wouldreveal genes whose expression is potentially regulated by post transcriptional mechanisms such asmiRNA cleavage. Combining this information with cell type-specific miRNA profiles, andto test the above hypothesis and identify key miRNA-target pairs important for nodule celldifferentiation. The use of translating ribosome affinity purification (TRAP) of nodule zonecells, namely from parenchyma and primordial tissues, to obtain cytoplasmic transcriptomes data.