Background: GRAS proteins are crucial transcription factors, which are plant-specific and participate in various plant biological processes. Thanks to the rapid progress of the whole genome sequencing technologies, theGRASgene families in different plants have been broadly explored and studied. However, comprehensive research on the soybean (Glycine max)GRASgene family is relatively lagging.Results: In this study, 117Glycine maxGRASgenes (GmGRAS) were identified. Further phylogenetic analyses showed that theGmGRASgenes could be categorized into nine gene subfamilies: DELLA, HAM, LAS, LISCL, PAT1, SCL3, SCL4/7, SCR and SHR. Gene structure analyses turned out that theGmGRASgenes lacked introns and were relatively conserved. Conserved domains and motif patterns of theGmGRASmembers in the same subfamily or clade exhibited similarities. Notably, the expansion of theGmGRASgene family was driven both by gene tandem and segmental duplication events. Whereas, segmental duplications took the major role in generating newGmGRASgenes. Moreover, the synteny and evolutionary constraints analyses of the GRAS proteins among soybean and distinct species (two monocots and four dicots) provided more detailed evidence forGmGRASgene evolution.Cis-element analyses indicated that theGmGRASgenes may be responsive to diverse environmental stresses and regulate distinct biological processes. Besides, the expression patterns of theGmGRASgenes were varied in various tissues, during saline and dehydration stresses and during seed germination processes.Conclusions: We conducted a systematic investigation of theGRASgenes in soybean, which may be valuable in paving the way for futureGmGRASgene studies and soybean breeding.

The utilization of soybean heterosis is probably one of the potential approaches in future yield breakthrough as was the situation in rice breeding in China. Cytoplasmic male sterility (CMS) plays an important role in the production of hybrid seeds. However, the molecular mechanism of CMS in soybean remains unclear. The comparative transcriptome analysis between cytoplasmic male sterile line NJCMS1A and its near-isogenic maintainer NJCMS1B in soybean was conducted using Illumina sequencing technology. A total of 88,643 transcripts were produced in Illumina sequencing. Then 56,044 genes were obtained matching soybean reference genome. Three hundred and sixty five differentially expressed genes (DEGs) between NJCMS1A and NJCMS1B were screened by threshold, among which, 339 down-regulated and 26 up-regulated in NJCMS1A compared to in NJCMS1B. Gene Ontology (GO) annotation showed that 242 DEGs were annotated to 19 functional categories. Clusters of Orthologous Groups of proteins (COG) annotation showed that 265 DEGs were classified into 19 categories. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that 46 DEGs were assigned to 33 metabolic pathways. According to functional and metabolic pathway analysis combined with reported literatures, the relations between some key DEGs and the male sterility of NJCMS1A were discussed. qRT-PCR analysis validated that the gene expression pattern in RNA-Seq was reliable. Finally, enzyme activity assay showed that energy supply was decreased in NJCMS1A compared to in NJCMS1B. We concluded that the male sterility of NJCMS1A might be related to the disturbed functions and metabolism pathways of some key DEGs, such as DEGs involved in carbohydrate and energy metabolism, transcription factors, regulation of pollen development, elimination of reactive oxygen species (ROS), cellular signal transduction, and programmed cell death (PCD) etc. Future research will focus on cloning and transgenic function validation of possible candidate genes associated with soybean CMS.

The three-amino-acid-loop-extension (TALE) superfamily genes broadly existed in plants, which played important roles in plant growth, development and abiotic stress responses. In this study, we identified 68 Glycine max TALE (GmTALE) superfamily members. Phylogenetic analysis divided the GmTALE superfamily into the BEL1-like (BLH/BELL homeodomain) and the KNOX (KNOTTED-like homeodomain) subfamilies. Moreover, the KNOX subfamily could be further categorized into three clades (KNOX Class I, KNOX Class II and KNOX Class III). The GmTALE genes showed similarities in the gene structures in the same subfamily or clade, whose coding proteins exhibited analogous motif and conserved domain compositions. Besides, synteny analyses and evolutionary constraint evaluations of the TALE members among soybean and different species provided more clues for GmTALE superfamily evolution. The cis-element analyses in gene promoter regions and relevant gene expression profiling revealed different regulating roles of GmTALE genes during soybean plant development, saline and dehydration stresses. Genome-wide characterization, evolution, and expression profile analyses of GmTALE genes can pave the way for future gene functional research and facilitate their roles for applications in genetic improvement on soybean in saline and dehydration stresses.

Cytoplasmic male sterility (CMS) provides crucial breeding materials that facilitate hybrid seed production in various crops, and thus plays an important role in the study of hybrid vigor (heterosis), in plants. However, the CMS regulatory network in soybean remains unclear. MicroRNAs (miRNAs) play crucial roles in flower and pollen development by targeting genes that regulate their expression in plants. To identify the miRNAs and their targets that exist in the soybean CMS line NJCMS1A and its maintainer NJCMS1B, high-throughput sequencing and degradome analysis were conducted in this study. Two small RNA libraries were constructed from the flower buds of the soybean CMS line NJCMS1A and its maintainer NJCMS1B. A total of 105 new miRNAs present on the other arm of known pre-miRNAs, 23 new miRNA members, 158 novel miRNAs and 160 high-confidence soybean miRNAs were identified using high-throughput sequencing. Among the identified miRNAs, 101 differentially expressed miRNAs with greater than two-fold changes between NJCMS1A and NJCMS1B were discovered. The different expression levels of selected miRNAs were confirmed by stem-loop quantitative real-time PCR. A degradome analysis showed that 856 targets were predicted to be targeted by 296 miRNAs, including a squamosa promoter-binding protein-like transcription factor family protein, a pentatricopeptide repeat-containing protein, and an auxin response factor, which were previously shown to be involved in floral organ or anther development in plants. Additionally, some targets, including a MADS-box transcription factor, NADP-dependent isocitrate dehydrogenase and NADH-ubiquinone oxidoreductase 24 kDa subunit, were identified, and they may have some relationship with the programmed cell death, reactive oxygen species accumulation and energy deficiencies, which might lead to soybean male sterility. The present study is the first to use deep sequencing technology to identify miRNAs and their targets in the flower buds of the soybean CMS line NJCMS1A and its maintainer NJCMS1B. The results revealed that the miRNAs might participate in flower and pollen development, which could facilitate our understanding of the molecular mechanisms behind CMS in soybean.