The Effects of Gene Duplication Modes on the Evolution of Regulatory Divergence in Wild and Cultivated Soybean

Genome-wide identification of glutamate receptor-like gene family in soybean

Glutamate receptor-like genes (GLRs) are essential in the growth and development of plants and many physiological and biochemical processes; however, related information in soybean is lacking. In this study, 105 GLRs, including 67 Glycine soja and 38 Glycine max GLRs, were identified and divided into two clades (Clades II and III) according to their phylogenetic relationships. GLR members in the same branch had a relatively conservative motif composition and genetic structure. Furthermore, the soybean GLR family mainly experienced purification selection during evolution. Cis-acting element analysis, gene ontology, and Kyoto Encyclopedia of Genes and Genomic annotations indicated the complexity of the gene regulation and functional diversity of the soybean GLR. Moreover, transcriptome data analysis showed that these GLRs had different expression profiles in different tissues, and Clade III members had higher and more common expression patterns. Additionally, the expression profiles under jasmonic acid treatment and salt stress indicate that the GLR participated in the jasmonic acid signaling pathway and plays a role in salt treatment. This study provides information for a comprehensive understanding of the soybean GLR family and a reference for further functional research and genetic improvement.

Flagellin and mannitol modulate callose biosynthesis and deposition in soybean seedlings

Callose is a polymer deposited on the cell wall and is necessary for plant growth and development. Callose is synthesized by genes from the glucan synthase-like family (GSL) and dynamically responds to various types of stress. Callose can inhibit pathogenic infection, in the case of biotic stresses, and maintain cell turgor and stiffen the plant cell wall in abiotic stresses. Here, we report the identification of 23 GSL genes (GmGSL) in the soybean genome. We performed phylogenetic analyses, gene structure prediction, duplication patterns, and expression profiles on several RNA-Seq libraries. Our analyses show that WGD/Segmental duplication contributed to expanding this gene family in soybean. Next, we analyzed the callose responses in soybean under abiotic and biotic stresses. The data show that callose is induced by both osmotic stress and flagellin 22 (flg22) and is related to the activity of β-1,3-glucanases. By using RT-qPCR, we evaluated the expression of GSL genes during the treatment of soybean roots with mannitol and flg22. The GmGSL23 gene was upregulated in seedlings treated with osmotic stress or flg22, showing the essential role of this gene in the soybean defense response to pathogenic organisms and osmotic stress. Our results provide an important understanding of the role of callose deposition and regulation of GSL genes in response to osmotic stress and flg22 infection in soybean seedlings.

Genome-wide characterization of the xyloglucan endotransglucosylase/hydrolase gene family in Solanum lycopersicum L. and gene expression analysis in response to arbuscular mycorrhizal symbiosis

Xyloglucan endotransglucosylase/hydrolases (XTHs) are a glycoside hydrolase protein family involved in the biosynthesis of xyloglucans, with essential roles in the regulation of plant cell wall extensibility. By taking advantage of the whole genome sequence in Solanum lycopersicum, 37 SlXTHs were identified in the present work. SlXTHs were classified into four subfamilies (ancestral, I/II, III-A, III-B) when aligned to XTHs of other plant species. Gene structure and conserved motifs showed similar compositions in each subfamily. Segmental duplication was the primary mechanism accounting for the expansion of SlXTH genes. In silico expression analysis showed that SlXTH genes exhibited differential expression in several tissues. GO analysis and 3D protein structure indicated that all 37 SlXTHs participate in cell wall biogenesis and xyloglucan metabolism. Promoter analysis revealed that some SlXTHs have MeJA- and stress-responsive elements. qRT-PCR expression analysis of nine SlXTHs in leaves and roots of mycorrhizal colonized vs. non-colonized plants showed that eight of these genes were differentially expressed in leaves and four in roots, suggesting that SlXTHs might play roles in plant defense induced by arbuscular mycorrhiza. Our results provide valuable insight into the function of XTHs in S. lycopersicum, in addition to the response of plants to mycorrhizal colonization.

Fitness and Rhizobacteria of F2, F3 Hybrids of Herbicide-Tolerant Transgenic Soybean and Wild Soybean

The introduction of herbicide-tolerant (HT) transgenic soybeans (Glycine max (L.) Merr.) into farming systems raises great concern that transgenes may flow to endemic wild soybeans (Glycine soja Sieb. et Zucc.) via pollen, which may increase the ecological risks by increasing the fitness of hybrids under certain conditions and threaten the genetic diversity of wild soybean populations. In order to demonstrate the potential risk of gene flow from the HT soybean to the wild soybean, the fitness of F2 and F3 hybrids obtained from two wild soybean populations (HLJHRB-1, JSCZ) collected from China and the HT soybean was measured under farmland and wasteland soil conditions, as well as with or without weed competition. Compared with their wild progenitors, the F2 and F3 hybrids of HLJHRB-1 displayed a higher emergence rate, higher aboveground dry biomass, more pods and filled-seed plants, as well as better composite fitness under four planting conditions. The F2 and F3 hybrids of JSCZ also displayed a higher emergence rate, higher aboveground dry biomass, more pods, and more filled seeds per plant under mixed planting, whereas these characteristics were lower under pure planting conditions in wasteland and farmland soil. Therefore, the composite fitness of JSCZ hybrids was higher or lower depending on the planting conditions. Furthermore, the soil microbial communities of the F3 of HLJHRB-1, JSCZ, and the wild soybean were investigated with 16S rDNA sequencing, which showed that low alpha diversity of rhizobacteria was relative to high fitness, and Rhizobium played an important role in promoting F3 plant growth.

ABA-Dependent Regulation of Calcium-Dependent Protein Kinase Gene GmCDPK5 in Cultivated and Wild Soybeans

Calcium-dependent protein kinases (CDPKs) regulate plant development and stress responses. However, the interaction of these protein kinases with the abscisic acid (ABA) stress hormone signalling system has not been studied in detail. In Arabidopsis, AtCPK1 plays an important role in the acclimation of plants to environmental stresses. Phylogenetic and molecular analyses showed that, among 50 isoforms of Glycine max (L.) Merrill CDPKs, the GmCDPK27/GmCDPK48, GmCDPK5/GmCDPK24, and GmCDPK10/GmCDPK46 paralogous pairs were the isoforms most related to AtCDPK1. We investigated the expression of the corresponding six GmCDPKs genes during treatment with cold, heat, and salt stress. Wild soybean was the most resistant to stresses, and among the three cultivars studied (Sfera, Hodgson, and Hefeng25), Sfera was close to the wild type in terms of resistance. GmCDPK5 and GmCDPK10 were the most responsive to stress treatments, especially in wild soybean, compared with cultivars. Among the studied GmCDPK isoforms, only GmCDPK5 expression increased after treatment with abscisic acid (ABA) in a dose- and time-dependent manner. Targeted LC-MS/MS analysis of endogenous ABA levels showed that wild soybean and Sfera had nearly twice the ABA content of Hodgson and Hefeng25. An analysis of the expression of marker genes involved in ABA biosynthesis showed that GmNCED1-gene-encoding 9-cis-epoxycarotenoid dioxygenase 1 is induced to the greatest extent in wild soybean and Sfera under salt, cold, and heat exposure. Our data established a correlation between the induction of GmCDPK5 and ABA biosynthesis genes. GmCDPK5 is an interesting target for genetic and bioengineering purposes and can be used for genetic editing, overexpression, or as a marker gene in soybean varieties growing under unfavourable conditions.

Transcriptomic Variations and Network Hubs Controlling Seed Size and Weight During Maize Seed Development

To elucidate the mechanisms underlying seed development in maize, comprehensive RNA-seq analyses were conducted on Zhengdan1002 (ZD1002), Zhengdan958 (ZD958), and their parental lines during seven seed developmental stages. We found that gene expression levels were largely nonadditive in hybrids and that cis-only or trans × cis pattern played a large role in hybrid gene regulation during seed developmental stage. Weighted gene co-expression network (WGCNA) analysis showed that 36 modules were highly correlated (r = -0.90-0.92, p < 0.05) with kernel weight, length, and width during seed development. Forty-five transcription factors and 38 ribosomal protein genes were identified as major hub genes determining seed size/weight. We also described a network hub, Auxin Response Factor 12 of maize (ZmARF12), a member of a family of transcription factor that mediate gene expression in response to auxin, potentially links auxin signal pathways, cell division, and the size of the seeds. The ZmARF12 mutant exhibited larger seed size and higher grain weight. ZmARF12 transcription was negatively associated with cell division during seed development, which was confirmed by evaluating the yield of protoplasts that isolated from the kernels of the mutant and other inbred lines. Transient knock-down of ZmARF12 in maize plants facilitated cell expansion and division, whereas transient silencing of its potential interactor ZmIAA8 impaired cell division. ZmIAA8 expression was repressed in the ZmARF12 over-expressed protoplasts. The mutant phenotype and the genetics studies presented here illustrated evidence that ZmARF12 is a cell division repressor, and potentially determines the final seed size.

Pathogenesis-related protein 1 (PR-1) genes in soybean: Genome-wide identification, structural analysis and expression profiling under multiple biotic and abiotic stresses

Plant pathogenesis-related (PR) proteins are a large group of proteins, classified in 17 families, that are induced by pathological conditions. Here, we characterized the soybean PR-1 (GmPR-1) gene repertoire at the sequence, structural and expression levels. We found 24 GmPR-1 genes, clustered in two phylogenetic groups. GmPR-1 genes are under strong purifying selection, particularly those that emerged by tandem duplications. GmPR-1 promoter regions are abundant in cis-regulatory elements associated with major stress-related transcription factor families, namely WRKY, ERF, HD-Zip, C2H2, NAC, and GATA. We observed that 23 GmPR-1 genes are induced by stress conditions or exclusively expressed upon stress. We explored 1972 transcriptome samples, including 26 stress conditions, revealing that most GmPR-1 genes are differentially expressed in a plethora of biotic and abiotic stresses. Our findings highlight stress-responsive GmPR-1 genes with potential biotechnological applications, such as the development of transgenic lines with increased resistance to biotic and abiotic stresses.

Distinct GmASMTs are involved in regulating transcription factors and signalling cross-talk across embryo development, biotic, and abiotic stress in soybean

N-Acetylserotonin O-methyltransferase (ASMT) is the final enzyme involved in melatonin biosynthesis. Identifying the expression of ASMT will reveal the regulatory role in the development and stress conditions in soybean. To identify and characterize ASMT in soybean (GmASMT), we employed genome-wide analysis, gene structure, cis-acting elements, gene expression, co-expression network analysis, and enzyme assay. We found seven pairs of segmental and tandem duplication pairs among the 44 identified GmASMTs by genome-wide analysis. Notably, co-expression network analysis reported that distinct GmASMTs are involved in various stress response. For example, GmASMT3, GmASMT44, GmASMT17, and GmASMT7 are involved in embryo development, heat, drought, aphid, and soybean cyst nematode infections, respectively. These distinct networks of GmASMTs were associated with transcription factors (NAC, MYB, WRKY, and ERF), stress signalling, isoflavone and secondary metabolites, calcium, and calmodulin proteins involved in stress regulation. Further, GmASMTs demonstrated auxin-like activities by regulating the genes involved in auxin transporter (WAT1 and NRT1/PTR) and auxin-responsive protein during developmental and biotic stress. The current study identified the key regulatory role of GmASMTs during development and stress. Hence GmASMT could be the primary target in genetic engineering for crop improvement under changing environmental conditions.

Genome-Wide Identification of Soybean ABC Transporters Relate to Aluminum Toxicity

ATP-binding cassette (ABC) transporter proteins are a gene super-family in plants and play vital roles in growth, development, and response to abiotic and biotic stresses. The ABC transporters have been identified in crop plants such as rice and buckwheat, but little is known about them in soybean. Soybean is an important oil crop and is one of the five major crops in the world. In this study, 255 ABC genes that putatively encode ABC transporters were identified from soybean through bioinformatics and then categorized into eight subfamilies, including 7 ABCAs, 52 ABCBs, 48 ABCCs, 5 ABCDs, 1 ABCEs, 10 ABCFs, 111 ABCGs, and 21 ABCIs. Their phylogenetic relationships, gene structure, and gene expression profiles were characterized. Segmental duplication was the main reason for the expansion of the GmABC genes. Ka/Ks analysis suggested that intense purifying selection was accompanied by the evolution of GmABC genes. The genome-wide collinearity of soybean with other species showed that GmABCs were relatively conserved and that collinear ABCs between species may have originated from the same ancestor. Gene expression analysis of GmABCs revealed the distinct expression pattern in different tissues and diverse developmental stages. The candidate genes GmABCB23, GmABCB25, GmABCB48, GmABCB52, GmABCI1, GmABCI5, and GmABCI13 were responsive to Al toxicity. This work on the GmABC gene family provides useful information for future studies on ABC transporters in soybean and potential targets for the cultivation of new germplasm resources of aluminum-tolerant soybean.