RNA-Seq and Comparative Transcriptomic Analyses of Asian Soybean Rust Resistant and Susceptible Soybean Genotypes Provide Insights into Identifying Disease Resistance Genes.

Asian soybean rust (ASR), caused by Phakopsora pachyrhizi, is one of the most destructive foliar diseases that affect soybeans. Developing resistant cultivars is the most cost-effective, environmentally friendly, and easy strategy for controlling the disease. However, the current understanding of the mechanisms underlying soybean resistance to P. pachyrhizi remains limited, which poses a significant challenge in devising effective control strategies. In this study, comparative transcriptomic profiling using one resistant genotype and one susceptible genotype was performed under infected and control conditions to understand the regulatory network operating between soybean and P. pachyrhizi. RNA-Seq analysis identified a total of 6540 differentially expressed genes (DEGs), which were shared by all four genotypes. The DEGs are involved in defense responses, stress responses, stimulus responses, flavonoid metabolism, and biosynthesis after infection with P. pachyrhizi. A total of 25,377 genes were divided into 33 modules using weighted gene co-expression network analysis (WGCNA). Two modules were significantly associated with pathogen defense. The DEGs were mainly enriched in RNA processing, plant-type hypersensitive response, negative regulation of cell growth, and a programmed cell death process. In conclusion, these results will provide an important resource for mining resistant genes to P. pachyrhizi infection and valuable resources to potentially pyramid quantitative resistance loci for improving soybean germplasm.

Overexpression of GmMYB14 improves high-density yield and drought tolerance of soybean through regulating plant architecture mediated by the brassinosteroid pathway.

MYB transcription factors (TFs) have been reported to regulate the biosynthesis of secondary metabolites, as well as to mediate plant adaption to abiotic stresses, including drought. However, the roles of MYB TFs in regulating plant architecture and yield potential remain poorly understood. Here, we studied the roles of the dehydration-inducible GmMYB14 gene in regulating plant architecture, high-density yield and drought tolerance through the brassinosteroid (BR) pathway in soybean. GmMYB14 was shown to localize to nucleus and has a transactivation activity. Stable GmMYB14-overexpressing (GmMYB14-OX) transgenic soybean plants displayed a semi-dwarfism and compact plant architecture associated with decreased cell size, resulting in a decrease in plant height, internode length, leaf area, leaf petiole length and leaf petiole angle, and improved yield in high density under field conditions. Results of the transcriptome sequencing suggested the involvement of BRs in regulating GmMYB14-OX plant architecture. Indeed, GmMYB14-OX plants showed reduced endogenous BR contents, while exogenous application of brassinolide could partly rescue the phenotype of GmMYB14-OX plants. Furthermore, GmMYB14 was shown to directly bind to the promoter of GmBEN1 and up-regulate its expression, leading to reduced BR content in GmMYB14-OX plants. GmMYB14-OX plants also displayed improved drought tolerance under field conditions. GmBEN1 expression was also up-regulated in the leaves of GmMYB14-OX plants under polyethylene glycol treatment, indicating that the GmBEN1-mediated reduction in BR level under stress also contributed to drought/osmotic stress tolerance of the transgenic plants. Our findings provided a strategy for stably increasing high-density yield and drought tolerance in soybean using a single TF-encoding gene.

Genome-Wide Investigation and Expression Profiling Under Abiotic Stresses of a Soybean Unknown Function (DUF21) and Cystathionine-ß-Synthase (CBS) Domain-Containing Protein Family.

Cystathionine-ß-synthase (CBS) domain-containing proteins (CDCPs) constitute a large family in plants, and members of this family have been implicated in a variety of biological processes. However, the precise functions and the underlying mechanisms of most members of this family in plants remain to be elucidated. CBSDUF proteins belong to the CDCP superfamily, which contains one domain of unknown function (DUF21) and an N terminus that is adjacent to two intracellular CBS domains. In this study, a comprehensive genome database analysis of soybean was performed to investigate the role(s) of these CBSDUFs and to explore their nomenclature, classification, chromosomal distribution, exon-intron organization, protein structure, and phylogenetic relationships; the analysis identified a total of 18 putative CBSDUF genes. Using specific protein domains and phylogenetic analysis, the CBSDUF gene family was subdivided into eight groups. The soybean CBSDUF genes showed an uneven distribution on 12 chromosomes of Glycine max. RNA-seq transcriptome data from different tissues in public databases revealed tissue-specific and differential expression profiles of the GmCBSDUFs, and qPCR analysis revealed that certain groups of soybean CBSDUFs are likely involved in specific stress responses. In addition, GmCBSDUF3 transgenic Arabidopsis was subjected to phenotypic analysis under NaCl, PEG, and ABA stress treatments. The overexpression of GmCBSDUF3 could enhance tolerance to drought and salt stress in Arabidopsis. This study presents a first comprehensive look at soybean CBSDUF proteins and provides valuable resources for functionally elucidating this protein subgroup within the CBS domain-containing protein family.

Identification of soybean drought-tolerant genotypes and loci correlated with agronomic traits contributes new candidate genes for breeding.

KEY MESSAGE: Drought tolerance level of 136 soybean genotypes, the correlations among traits were evaluated, and several important drought-tolerant genotypes, traits, SNPs and genes were possibly useful for soybean genetic breeding. Drought is an adverse environmental factor affecting crops growth, development, and yield. Promising genotypes and genes with improved tolerance to drought are probably effective ways to alleviate the situation. In this study, our main task was to determine drought tolerance level of 136 soybean genotypes, the correlations among physiological and agronomic traits under drought, and drought-tolerant single nucleotide polymorphism (SNPs) and genes. In this study, twenty-six varieties were identified as excellent tolerant genotypes to stress among which S14, S93 and S135 with high drought-tolerant index (DTI > 1.3) and yield (Y > 300 kg). Fourteen varieties were identified as drought-sensitive genotypes, such as S25, S45 and S58, with low drought-tolerant index (DTI < 0.5). 422 SNPs and 302 genes correlated with seed number per plant (SNPP), maturity (M), number of seeds per pod (NSPP), node number of main stem (NNMS), Stem diameter (SD) and pull stem (PS) were detected under well-watered and drought conditions by genome-wide association study (GWAS). Among them, we found SNPs (Chr 3:1758920-1958934) between drought-tolerant and sensitive genotypes. 13 genes (Glyma.03G017800, Glyma.03G018000, Glyma.03G018200, Glyma.03G018400, Glyma.03G018500, Glyma.03G018600, Glyma.03G018700, Glyma.03G018800, Glyma.03G018900, Glyma.03G019000, Glyma.03G019100, Glyma.03G019200, Glyma.03G019300) correlated with NNMS were detected. By qRT-PCR, the expression level of Glyma.03G018000 and Glyma.03G018900 in drought-tolerant varieties was significantly increased, but low or no expression in sensitive varieties under drought stress. This study provides important drought-tolerant genotypes, traits, SNPs and potential genes, possibly useful for soybean genetic breeding.

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.

CRISPR/Cas9-mediated targeted mutagenesis of GmSPL9 genes alters plant architecture in soybean.

BACKGROUND: The plant architecture has significant effects on grain yield of various crops, including soybean (Glycine max), but the knowledge on optimization of plant architecture in order to increase yield potential is still limited. Recently, CRISPR/Cas9 system has revolutionized genome editing, and has been widely utilized to edit the genomes of a diverse range of crop plants. RESULTS: In the present study, we employed the CRISPR/Cas9 system to mutate four genes encoding SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) transcription factors of the SPL9 family in soybean. These four GmSPL9 genes are negatively regulated by GmmiR156b, a target for the improvement of soybean plant architecture and yields. The soybean Williams 82 was transformed with the binary CRISPR/Cas9 plasmid, assembled with four sgRNA expression cassettes driven by the Arabidopsis thaliana U3 or U6 promoter, targeting different sites of these four SPL9 genes via Agrobacterium tumefaciens-mediated transformation. A 1-bp deletion was detected in one target site of the GmSPL9a and one target site of the GmSPL9b, respectively, by DNA sequencing analysis of two T0-generation plants. T2-generation spl9a and spl9b homozygous single mutants exhibited no obvious phenotype changes; but the T2 double homozygous mutant spl9a/spl9b possessed shorter plastochron length. In T4 generation, higher-order mutant plants carrying various combinations of mutations showed increased node number on the main stem and branch number, consequently increased total node number per plants at different levels. In addition, the expression levels of the examined GmSPL9 genes were higher in the spl9b-1 single mutant than wild-type plants, which might suggest a feedback regulation on the expression of the investigated GmSPL9 genes in soybean. CONCLUSIONS: Our results showed that CRISPR/Cas9-mediated targeted mutagenesis of four GmSPL9 genes in different combinations altered plant architecture in soybean. The findings demonstrated that GmSPL9a, GmSPL9b, GmSPL9c and GmSPL9 function as redundant transcription factors in regulating plant architecture in soybean.

Identification and Comparative Analysis of CBS Domain-Containing Proteins in Soybean (Glycine max) and the Primary Function of GmCBS21 in Enhanced Tolerance to Low Nitrogen Stress.

Nitrogen is an important macronutrient required for plant growth, and is a limiting factor for crop productivity. Improving the nitrogen use efficiency (NUE) is therefore crucial. At present, the NUE mechanism is unclear and information on the genes associated with NUE in soybeans is lacking. cystathionine beta synthase (CBS) domain-containing proteins (CDCPs) may be implicated in abiotic stress tolerance in plants. We identified and classified a CBS domain-containing protein superfamily in soybean. A candidate gene for NUE, GmCBS21, was identified. GmCBS21 gene characteristics, the temporal expression pattern of the GmCBS21 gene, and the phenotype of GmCBS21 overexpression in transgenic Arabidopsis thaliana under low nitrogen stress were analyzed. The phenotypes suggested that the transgenic Arabidopsis thaliana seedlings performed better under the nitrogen-deficient condition. GmCBS21-overexpressing transgenic plants exhibit higher low nitrogen stress tolerance than WT plants, and this suggests its role in low nitrogen stress tolerance in plants. We conclude that GmCBS21 may serve as an excellent candidate for breeding crops with enhanced NUE and better yield.

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.

An pair of an atypical NLR encoding genes confer Asian soybean rust resistance in soybean.

Asian soybean rust (ASR), caused by Phakopsora pachyrhizi, is a devastating disease that is present in all major soybean-producing regions. The limited availability of resistant germplasm has resulted in a scarcity of commercial soybean cultivars that are resistant to the disease. To date, only the Chinese soybean landrace SX6907 has demonstrated an immune response to ASR. In this study, we present the isolation and characterization of Rpp6907-7 and Rpp6907-4, a gene pair that confer broad-spectrum resistance to ASR. Rpp6907-7 and Rpp6907-4 encode atypic nucleotide-binding leucine-rich repeat (NLR) proteins that are found to be required for NLR-mediated immunity. Genetic analysis shows that only Rpp6907-7 confers resistance, while Rpp6907-4 regulates Rpp6907-7 signaling activity by acting as a repressor in the absence of recognized effectors. Our work highlights the potential value of using Rpp6907 in developing resistant soybean cultivars.

Heterologous expression of the Glycine soja Kunitz-type protease inhibitor GsKTI improves resistance to drought stress and Helicoverpa armigera in transgenic Arabidopsis lines.

Kunitz-like protease inhibitors (KTIs) have been identified to play critical roles in insect defense, but evidence for their involvement in drought stress is sparse. The aim of this study was to identify and functionally characterize a Kunitz-like protease inhibitor, GsKTI, from the wild soybean (Glycine soja) variety ED059. Expression patterns suggest that drought stress and insect herbivory may induce GsKTI transcript levels. Transgenic Arabidopsis lines overexpressing GsKTI have been shown to exhibit enhanced drought tolerance by regulating the ABA signaling pathway and increasing xylem cell number. Transgenic Arabidopsis leaves overexpressing GsKTI interfered with insect digestion and thus had a negative effect on the growth of Helicoverpa armigera. It is concluded that GsKTI increases resistance to drought stress and insect attack in transgenic Arabidopsis lines.

Melatonin Attenuates the Urea-Induced Yields Improvement Through Remodeling Transcriptome and Rhizosphere Microbial Community Structure in Soybean.

Foliar application of nitrogen to enhance crop productivity has been widely used. Melatonin is an effective regulator in promoting plant growth. However, the effects of melatonin and the combination of melatonin and nitrogen on soybeans yields production remain largely unknown. In this study, a field experiment was conducted to evaluate the effects and mechanisms of spraying leaves with melatonin and urea on soybeans. Foliar application of urea significantly increased soybean yields and melatonin did not affect the yields, while combination of melatonin and urea significantly reduced the yields compared to the application of urea alone. A leaf transcriptional profile was then carried out to reveal the underlying mechanism and found that foliar spraying of urea specifically induced the expression of genes related to amino acid transport and nitrogen metabolism. However, foliar application of melatonin significantly changed the transcriptional pattern established by urea application and increased the expression of genes related to abiotic stress signaling pathways. The effects of melatonin and urea treatment on soil microbiome were also investigated. Neither melatonin nor urea application altered the soil microbial alpha diversity, but melatonin application changed rhizosphere microbial community structure, whereas the combination of melatonin and urea did not. Melatonin or urea application altered the abundance of certain taxa. The number of taxa changed by melatonin treatment was higher than urea treatment. Collectively, our results provide new and valuable insights into the effects of foliar application of melatonin to urea and further show that melatonin exerts strong antagonistic effects on urea-induced soybean yields, gene expression and certain soil microorganisms.

CRISPR/Cas9-Mediated Targeted Mutagenesis of GmUGT Enhanced Soybean Resistance Against Leaf-Chewing Insects Through Flavonoids Biosynthesis.

Leaf-chewing insects are important pests that cause yield loss and reduce seed quality in soybeans (Glycine max). Breeding soybean varieties that are resistant to leaf-chewing insects can minimize the need for insecticide use and reduce yield loss. The marker gene for QTL-M, Glyma.07g110300 (LOC100775351) that encodes a UDP-glycosyltransferase (UGT) is the major determinant of resistance against leaf-chewing insects in soybean; it exhibits a loss of function in insect-resistant soybean germplasms. In this study, Agrobacterium-mediated transformation introduced the CRISPR/Cas9 expression vector into the soybean cultivar Tianlong No. 1 to generate Glyma.07g110300-gene mutants. We obtained two novel types of mutations, a 33-bp deletion and a single-bp insertion in the GmUGT coding region, which resulted in an enhanced resistance to Helicoverpa armigera and Spodoptera litura. Additionally, overexpressing GmUGT produced soybean varieties that were more sensitive to H. armigera and S. litura. Both mutant and overexpressing lines exhibited no obvious phenotypic changes. The difference in metabolites and gene expression suggested that GmUGT is involved in imparting resistance to leaf-chewing insects by altering the flavonoid content and expression patterns of genes related to flavonoid biosynthesis and defense. Furthermore, ectopic expression of the GmUGT gene in the ugt72b1 mutant of Arabidopsis substantially rescued the phenotype of H. armigera resistance in the atugt72b1 mutant. Our study presents a strategy for increasing resistance against leaf-chewing insects in soybean through CRISPR/Cas9-mediated targeted mutagenesis of the UGT genes.

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.

A GATA Transcription Factor from Soybean (Glycine max) Regulates Chlorophyll Biosynthesis and Suppresses Growth in the Transgenic Arabidopsis thaliana.

Chlorophyll plays an essential role in photosynthetic light harvesting and energy transduction in green tissues of higher plants and is closely related to photosynthesis and crop yield. Identification of transcription factors (TFs) involved in regulating chlorophyll biosynthesis is still limited in soybean (Glycine max), and the previously identified GmGATA58 is suggested to potentially modulate chlorophyll and nitrogen metabolisms, but its complete function is still unknown. In this study, subcellular localization assay showed that GmGATA58 was localized in the nucleus. Histochemical GUS assay and qPCR assay indicated that GmGATA58 was mainly expressed in leaves and responded to nitrogen, light and phytohormone treatments. Overexpression of GmGATA58 in the Arabidopsis thaliana ortholog AtGATA21 (gnc) mutant complemented the greening defect, while overexpression in Arabidopsis wild-type led to increasing chlorophyll content in leaves through up-regulating the expression levels of the large of chlorophyll biosynthetic pathway genes, but suppressing plant growth and yield, although the net photosynthetic rate was slightly improved. Dual-luciferase reporter assay also supported that GmGATA58 activated the transcription activities of three promoters of key chlorophyll biosynthetic genes of soybean in transformed protoplast of Arabidopsis. It is concluded that GmGATA58 played an important role in regulating chlorophyll biosynthesis, but suppressed plant growth and yield in transgenic Arabidopsis.

Genome-Wide Analysis of the WOX Gene Family and Function Exploration of GmWOX18 in Soybean.

WUSCHEL-related homeobox (WOX) is a family of transcription factors that are unique to plants and is characterized by the presence of a homeodomain. The WOX transcription factor plays an important role in regulating plant growth and development and the response to abiotic stress. Soybean is one of the most important oil crops worldwide. In this study, based on the available genome data of soybean, the WOX gene family was identified by bioinformatics analysis. The chromosome distribution, gene and protein structures, phylogenetic relationship and gene expression patterns of this family were comprehensively compared. The results showed that a total of 33 putative WOX genes in the soybean genome were found and then designated as GmWOX1- GmWOX33, which were distributed across 19 chromosomes except chromosome 16. Multiple sequence analysis of the GmWOX gene family revealed a highly conserved homeodomain. Phylogenetic tree analysis showed that 33 WOX genes could be divided into three major clades (modern/WUS, intermediate and ancient) in soybean. Of these 33 WOX genes, some showed differential expression patterns in the tested tissues (leaves, pods, unopen and open flowers, nodules, seed, roots, root hairs, stems, shoot apical meristems and shoot tips). In addition, the expression profile and qRT-PCR analysis showed that most of the GmWOX genes responded to different abiotic stress treatments (cold and drought). According to the expression pattern of GmWOX genes in the high regeneration capacity soybean material P3, overexpression of GmWOX18 was selected for function analysis. The overexpression of GmWOX18 increased the regeneration ability of clustered buds. The results will provide valuable information for further studies on the roles of WOX genes in regulating soybean growth, development and responses to abiotic stress, as well as a basis for the functional identification and analysis of WOX genes in soybean.

GmSYP24, a putative syntaxin gene, confers osmotic/drought, salt stress tolerances and ABA signal pathway.

As major environment factors, drought or high salinity affect crop growth, development and yield. Transgenic approach is an effective way to improve abiotic stress tolerance of crops. In this study, we comparatively analyzed gene structures, genome location, and the evolution of syntaxin proteins containing late embryogenesis abundant (LEA2) domain. GmSYP24 was identified as a dehydration-responsive gene. Our study showed that the GmSYP24 protein was located on the cell membrane. The overexpression of GmSYP24 (GmSYP24ox) in soybean and heteroexpression of GmSYP24 (GmSYP24hx) in Arabidopsis exhibited insensitivity to osmotic/drought and high salinity. However, wild type soybean, Arabidopsis, and the mutant of GmSYP24 homologous gene of Arabidopsis were sensitive to the stresses. Under the abiotic stresses, transgenic soybean plants had greater water content and higher activities of POD, SOD compared with non-transgenic controls. And the leaf stomatal density and opening were reduced in transgenic Arabidopsis. The sensitivity to ABA was decreased during seed germination of GmSYP24ox and GmSYP24hx. GmSYP24hx induced up-regulation of ABA-responsive genes. GmSYP24ox alters the expression of some aquaporins under osmotic/drought, salt, or ABA treatment. These results demonstrated that GmSYP24 played an important role in osmotic/drought or salt tolerance in ABA signal pathway.

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.

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.

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

The root nodule symbiosis (RNS) between legume plants and rhizobia is the most efficient and productive source of nitrogen fixation, and has critical importance in agriculture and mesology. Soybean (Glycine max), one of the most important legume crops in the world, establishes a nitrogen-fixing symbiosis with different types of rhizobia, and the efficiency of symbiotic nitrogen fixation in soybean greatly depends on the symbiotic host-specificity. Although, it has been reported that rhizobia use surface polysaccharides, secretion proteins of the type-three secretion systems and nod factors to modulate host range, the host control of nodulation specificity remains poorly understood. In this report, the soybean roots of two symbiotic systems (Bradyrhizobium japonicum strain 113-2-soybean and Sinorhizobium fredii USDA205-soybean)with notable different nodulation phenotypes and the control were studied at five different post-inoculation time points (0.5, 7-24 h, 5, 16, and 21 day) by RNA-seq (Quantification). The results of qPCR analysis of 11 randomly-selected genes agreed with transcriptional profile data for 136 out of 165 (82.42%) data points and quality assessment showed that the sequencing library is of quality and reliable. Three comparisons (control vs. 113-2, control vs. USDA205 and USDA205 vs. 113-2) were made and the differentially expressed genes (DEGs) between them were analyzed. The number of DEGs at 16 days post-inoculation (dpi) was the highest in the three comparisons, and most of the DEGs in USDA205 vs. 113-2 were found at 16 dpi and 21 dpi. 44 go function terms in USDA205 vs. 113-2 were analyzed to evaluate the potential functions of the DEGs, and 10 important KEGG pathway enrichment terms were analyzed in the three comparisons. Some important genes induced in response to different strains (113-2 and USDA205) were identified and analyzed, and these genes primarily encoded soybean resistance proteins, NF-related proteins, nodulins and immunity defense proteins, as well as proteins involving flavonoids/flavone/flavonol biosynthesis and plant-pathogen interaction. Besides, 189 candidate genes are largely expressed in roots and\or nodules. The DEGs uncovered in this study provides molecular candidates for better understanding the mechanisms of symbiotic host-specificity and explaining the different symbiotic effects between soybean roots inoculated with different strains (113-2 and USDA205).