Soybean is one of the most important crops, providing large amounts of dietary proteins and edible oil, and is also an excellent model for studying evolution of duplicated genes. However, relative to the model plants Arabidopsis and rice, the present knowledge about soybean transcriptome is quite limited. In this study, we employed RNA-seq to investigate transcriptomes of 11 soybean tissues, for genome-wide discovery of truly expressed genes, and novel and alternative transcripts, as well as analyses of conservation and divergence of duplicated genes and their functional implications. We detected a total of 54,132 high-confidence expressed genes, and identified 6,718 novel transcriptional regions with a mean length of 372 bp. We also provided strong evidence for alternative splicing (AS) events for ~15.9% of the genes with two or more exons. Among them, 1,834 genes exhibited stage-dependent AS, and 202 genes had tissue-biased exon-skipping events. We further defined the conservation and divergence in expression patterns between duplicated gene pairs from recent whole genome duplications (WGDs); differentially expressed genes, tissue preferentially expressed genes, transcription factors and specific gene family members were identified for shoot apical meristem and flower development. Our results significantly improved soybean gene annotation, and also provide valuable resources for functional genomics and studies of the evolution of duplicated genes from WGDs in soybean.

The WRKY transcription factors (TFs) are one of the largest families of TFs in plants and play multiple roles in plant development and stress response. In the present study, GmWRKY16 encoding a WRKY transcription factor in soybean was functionally characterized in Arabidopsis. GmWRKY16 is a nuclear protein that contains a highly conserved WRKY domain and a C₂H₂ zinc-finger structure, and has the characteristics of transcriptional activation ability, presenting a constitutive expression pattern with relative expression levels of over fourfold in the old leaves, flowers, seeds and roots of soybean. The results of quantitative real time polymerase chain reaction (qRT-PCR) showed that GmWRKY16 could be induced by salt, alkali, ABA, drought and PEG-6000. As compared with the control, overexpression of GmWRKY16 in Arabidopsis increased the seed germination rate and root growth of seedlings in transgenic lines under higher concentrations of mannitol, NaCl and ABA. In the meantime, GmWRKY16 transgenic lines showed over 75% survival rate after rehydration and enhanced Arabidopsis tolerance to salt and drought with higher proline and lower MDA accumulation, less water loss of the detached leaves, and accumulated more endogenous ABA than the control under stress conditions. Further studies showed that AtWRKY8, KIN1, and RD29A were induced in GmWRKY16 transgenic plants under NaCl treatment. The expressions of the ABA biosynthesis gene (NCED3), signaling genes (ABI1, ABI2, ABI4, and ABI5), responsive genes (RD29A, COR15A, COR15B, and RD22) and stress-related marker genes (KIN1, LEA14, LEA76, and CER3) were regulated in transgenic lines under drought stress. In summary, these results suggest that GmWRKY16 as a WRKY TF may promote tolerance to drought and salt stresses through an ABA-mediated pathway.

Summary: The importance of rhizomicrobiome in plant development, nutrition acquisition and stress tolerance is unquestionable. Relevant plant genes corresponding to the above functions also regulate rhizomicrobiome construction. Deciphering the molecular regulatory network of plant‐microbe interactions could substantially contribute to improving crop yield and quality. Here, the plant gene‐related nutrient uptake, biotic and abiotic stress resistance, which may influence the composition and function of microbial communities, are discussed in this review. In turn, the influence of microbes on the expression of functional plant genes, and thereby plant growth and immunity, is also reviewed. Moreover, we have specifically paid attention to techniques and methods used to link plant functional genes and rhizomicrobiome. Finally, we propose to further explore the molecular mechanisms and signalling pathways of microbe‐host gene interactions, which could potentially be used for managing plant health in agricultural systems.

Background: MicroRNAs (miRNAs) play important regulatory roles in development and stress response in plants. Wild soybean (Glycine soja) has undergone long-term natural selection and may have evolved special mechanisms to survive stress conditions as a result. However, little information about miRNAs especially miRNAs responsive to aluminum (Al) stress is available in wild soybean. Results: Two small RNA libraries and two degradome libraries were constructed from the roots of Al-treated and Al-free G. soja seedlings. For miRNA identification, a total of 7,287,655 and 7,035,914 clean reads in Al-treated and Al-free small RNAs libraries, respectively, were generated, and 97 known miRNAs and 31 novel miRNAs were identified. In addition, 49 p3 or p5 strands of known miRNAs were found. Among all the identified miRNAs, the expressions of 30 miRNAs were responsive to Al stress. Through degradome sequencing, 86 genes were identified as targets of the known miRNAs and five genes were found to be the targets of the novel miRNAs obtained in this study. Gene ontology (GO) annotations of target transcripts indicated that 52 target genes cleaved by conserved miRNA families might play roles in the regulation of transcription. Additionally, some genes, such as those for the auxin response factor (ARF), domain-containing disease resistance protein (NB-ARC), leucine-rich repeat and toll/interleukin-1 receptor-like protein (LRR-TIR) domain protein, cation transporting ATPase, Myb transcription factors, and the no apical meristem (NAM) protein, that are known to be responsive to stress, were found to be cleaved under Al stress conditions. Conclusions: A number of miRNAs and their targets were detected in wild soybean. Some of them that were responsive to biotic and abiotic stresses were regulated by Al stress. These findings provide valuable information to understand the function of miRNAs in Al tolerance.

Sugarcane-soybean intercropping has been widely used to control disease and improve nutrition in the field. However, the response of the soil microbial community diversity and structure to intercropping is not well understood. Since microbial diversity corresponds to soil quality and plant health, a pot experiment was conducted with sugarcane intercropped with soybean. Rhizosphere soil was collected 40 days after sowing, and MiSeq sequencing was utilized to analyze the soil microbial community diversity and composition. Soil columns were used to assess the influence of intercropping on soil microbial activity (soil respiration and carbon-use efficiency: nitrogen-use efficiency ratio). PICRUSt and FUNGuild analysis were conducted to predict microbial functional profiling. Our results showed that intercropping decreased pH by approximately 8.9% and enhanced the soil organic carbon, dissolved organic carbon, and available nitrogen (N) by 5.5%, 13.4%, and 10.0%, respectively. These changes in physicochemical properties corresponded to increased microbial diversity and shifts in soil microbial communities. Microbial community correlated significantly (p < 0.05) with soil respiration rates and nutrient use efficiency. Furthermore, intercropping influenced microbial functions, such as carbon fixation pathways in prokaryotes, citrate cycle (TCA cycle) of bacteria and wood saprotrophs of fungi. These overrepresented functions might accelerate nutrient conversion and control phytopathogens in soil.