Down-Regulating the Expression of 53 Soybean Transcription Factor Genes Uncovers a Role for SPEECHLESS in Initiating Stomatal Cell Lineages during Embryo Development

Stomatal density suppressor PagSDD1 is a "generalist" gene that promotes plant growth and improves water use efficiency

The serine protease SDD1 regulates stomatal density, but its potential impact on plant vegetative growth is unclear. Our study reveals a substantial upregulation of SDD1 in triploid poplar apical buds and leaves, suggesting its possible role in their growth regulation. We cloned PagSDD1 from poplar 84 K (Populus alba × P. glandulosa) and found that overexpression in poplar, soybean, and lettuce led to decreased leaf stomatal density. Furthermore, PagSDD1 represses PagEPF1, PagEPF2, PagEPFL9, PagSPCH, PagMUTE, and PagFAMA expression. In contrast, PagSDD1 promotes the expression of its receptors, PagTMM and PagERECTA. PagSDD1-OE poplars showed stronger drought tolerance than wild-type poplars. Simultaneously, PagSDD1-OE poplar, soybean, and lettuce had vegetative growth advantages. RNA sequencing revealed a significant upregulation of genes PagLHCB2.1 and PagGRF5, correlating positively with photosynthetic rate, and PagCYCA3;4 and PagEXPA8 linked to cell division and differentiation in PagSDD1-OE poplars. This increase promoted leaf photosynthesis, boosted auxin and cytokinin accumulation, and enhanced vegetative growth. SDD1 overexpression can increase the biomass of poplar, soybean, and lettuce by approximately 70, 176, and 155 %, respectively, and increase the water use efficiency of poplar leaves by over 52 %, which is of great value for the molecular design and breeding of plants with growth and water-saving target traits.

Identification and Expression Analysis of the Nucleotidyl Transferase Protein (NTP) Family in Soybean (Glycine max) under Various Abiotic Stresses

Nucleotidyl transferases (NTPs) are common transferases in eukaryotes and play a crucial role in nucleotide modifications at the 3' end of RNA. In plants, NTPs can regulate RNA stability by influencing 3' end modifications, which in turn affect plant growth, development, stress responses, and disease resistance. Although the functions of NTP family members have been extensively studied in Arabidopsis, rice, and maize, there is limited knowledge about NTP genes in soybeans. In this study, we identified 16 members of the NTP family in soybeans, including two subfamilies (G1 and G2) with distinct secondary structures, conserved motifs, and domain distributions at the protein level. Evolutionary analysis of genes in the NTP family across multiple species and gene collinearity analysis revealed a relatively conserved evolutionary pattern. Analysis of the tertiary structure of the proteins showed that NTPs have three conserved aspartic acids that bind together to form a possible active site. Tissue-specific expression analysis indicated that some NTP genes exhibit tissue-specific expression, likely due to their specific functions. Stress expression analysis showed significant differences in the expression levels of NTP genes under high salt, drought, and cold stress. Additionally, RNA-seq analysis of soybean plants subjected to salt and drought stress further confirmed the association of soybean NTP genes with abiotic stress responses. Subcellular localization experiments revealed that GmNTP2 and GmNTP14, which likely have similar functions to HESO1 and URT1, are located in the nucleus. These research findings provide a foundation for further investigations into the functions of NTP family genes in soybeans.

Genome-wide analysis of the SCAMPs gene family of soybean and functional identification of GmSCAMP5 in salt tolerance

The effect of salt damage on plants is mainly caused by the toxic effect of Na+. Studies showed that the secretory carrier membrane proteins were associated with the Na+ transport. However, the salt tolerance mechanism of secretory carrier protein (SCAMP) in soybean was yet to be defined. In this study, ten potential SCAMP genes distributed in seven soybean chromosomes were identified in the soybean genome. The phylogenetic tree of SCAMP domain sequences of several plants can divide SCAMPs into two groups. Most GmSCAMPs genes contained multiple Box4, MYB and MYC cis-elements indicated they may respond to abiotic stresses. We found that GmSCAMP1, GmSCAMP2 and GmSCAMP4 expressed in several tissues and GmSCAMP5 was significantly induced by salt stress. GmSCAMP5 showed the same expression patterns under NaCl treatment in salt-tolerant and salt-sensitive soybean varieties, but the induced time of GmSCAMP5 in salt-tolerant variety was earlier than that of salt-sensitive variety. To further study the effect of GmSCAMP5 on the salt tolerance of soybean plants, compared to GmSCAMP5-RNAi and EV-Control plants, GmSCAMP5-OE had less wilted leave and higher SPAD value. Compared to empty vector control, less trypan blue staining was observed in GmSCAMP5-OE leaves while more staining in GmSCAMP5-RNAi leaves. The Na+ of GmSCAMP5-RNAi plants leaves under NaCl stress were significantly higher than that in EV-Control plants, while significantly lower Na+ in GmSCAMP5-OE plants than in that EV-Control plants. The contents of leaves K+ of GmSCAMP5-RNAi, EV-Control, and GmSCAMP5-OE plants under NaCl stress were opposite to that of leaves Na+ content. Finally, salt stress-related genes NHX1, CLC1, TIP1, SOD1, and SOS1 in transformed hairy root changed significantly compared with the empty control. The research will provide novel information for study the molecular regulation mechanism of soybean salt tolerance.

Beyond skin-deep: targeting the plant surface for crop improvement

The above-ground plant surface is a well-adapted tissue layer that acts as an interface between the plant and its surrounding environment. As such, its primary role is to protect against desiccation and maintain the gaseous exchange required for photosynthesis. Further, this surface layer provides a barrier against pathogens and herbivory, while attracting pollinators and agents of seed dispersal. In the context of agriculture, the plant surface is strongly linked to post-harvest crop quality and yield. The epidermal layer contains several unique cell types adapted for these functions, while the non-lignified above-ground plant organs are covered by a hydrophobic cuticular membrane. This review aims to provide an overview of the latest understanding of the molecular mechanisms underlying crop cuticle and epidermal cell formation, with focus placed on genetic elements contributing towards quality, yield, drought tolerance, herbivory defence, pathogen resistance, pollinator attraction, and sterility, while highlighting the inter-relatedness of plant surface development and traits. Potential crop improvement strategies utilizing this knowledge are outlined in the context of the recent development of new breeding techniques.

Experimental validation of the mechanism of stomatal development diversification

Stomata are the structures responsible for gas exchange in plants. The established framework for stomatal development is based on the model plant Arabidopsis, but diverse patterns of stomatal development have been observed in other plant lineages and species. The molecular mechanisms behind these diversified patterns are still poorly understood. We recently proposed a model for the molecular mechanisms of the diversification of stomatal development based on the genus Callitriche (Plantaginaceae), according to which a temporal shift in the expression of key stomatal transcription factors SPEECHLESS and MUTE leads to changes in the behavior of meristemoids (stomatal precursor cells). In the present study, we genetically manipulated Arabidopsis to test this model. By altering the timing of MUTE expression, we successfully generated Arabidopsis plants with early differentiation or prolonged divisions of meristemoids, as predicted by the model. The epidermal morphology of the generated lines resembled that of species with prolonged or no meristemoid divisions. Thus, the evolutionary process can be reproduced by varying the SPEECHLESS to MUTE transition. We also observed unexpected phenotypes, which indicated the participation of additional factors in the evolution of the patterns observed in nature. This study provides novel experimental insights into the diversification of meristemoid behaviors.

Genome-Wide Tissue-Specific Genes Identification for Novel Tissue-Specific Promoters Discovery in Soybean

Promoters play a crucial role in controlling the spatial and temporal expression of genes at transcriptional levels in the process of higher plant growth and development. The spatial, efficient, and correct regulation of exogenous genes expression, as desired, is the key point in plant genetic engineering research. Constitutive promoters widely used in plant genetic transformation are limited because, sometimes, they may cause potential negative effects. This issue can be solved, to a certain extent, by using tissue-specific promoters. Compared with constitutive promoters, a few tissue-specific promoters have been isolated and applied. In this study, based on the transcriptome data, a total of 288 tissue-specific genes were collected, expressed in seven tissues, including the leaves, stems, flowers, pods, seeds, roots, and nodules of soybean (Glycine max). KEGG pathway enrichment analysis was carried out, and 52 metabolites were annotated. A total of 12 tissue-specific genes were selected via the transcription expression level and validated through real-time quantitative PCR, of which 10 genes showed tissue-specific expression. The 3-kb 5' upstream regions of ten genes were obtained as putative promoters. Further analysis showed that all the 10 promoters contained many tissue-specific cis-elements. These results demonstrate that high-throughput transcriptional data can be used as effective tools, providing a guide for high-throughput novel tissue-specific promoter discovery.

Extensive embryonic patterning without cellular differentiation primes the plant epidermis for efficient post-embryonic stomatal activities

Plant leaves feature epidermal stomata that are organized in stereotyped patterns. How does the pattern originate? We provide transcriptomic, imaging, and genetic evidence that Arabidopsis embryos engage known stomatal fate and patterning factors to create regularly spaced stomatal precursor cells. Analysis of embryos from 36 plant species indicates that this trait is widespread among angiosperms. Embryonic stomatal patterning in Arabidopsis is established in three stages: first, broad SPEECHLESS (SPCH) expression; second, coalescence of SPCH and its targets into discrete domains; and third, one round of asymmetric division to create stomatal precursors. Lineage progression is then halted until after germination. We show that the embryonic stomatal pattern enables fast stomatal differentiation and photosynthetic activity upon germination, but it also guides the formation of additional stomata as the leaf expands. In addition, key stomatal regulators are prevented from driving the fate transitions they can induce after germination, identifying stage-specific layers of regulation that control lineage progression during embryogenesis.

Different Leaf Anatomical Responses to Water Deficit in Maize and Soybean

The stomata on leaf surfaces control gas exchange and water loss, closing during dry periods to conserve water. The distribution and size of stomatal complexes is determined by epidermal cell differentiation and expansion during leaf growth. Regulation of these processes in response to water deficit may result in stomatal anatomical plasticity as part of the plant acclimation to drought. We quantified the leaf anatomical plasticity under water-deficit conditions in maize and soybean over two experiments. Both species produced smaller leaves in response to the water deficit, partly due to the reductions in the stomata and pavement cell size, although this response was greater in soybean, which also produced thicker leaves under severe stress, whereas the maize leaf thickness did not change. The stomata and pavement cells were smaller with the reduced water availability in both species, resulting in higher stomatal densities. Stomatal development (measured as stomatal index, SI) was suppressed in both species at the lowest water availability, but to a greater extent in maize than in soybean. The result of these responses is that in maize leaves, the stomatal area fraction (f_gc) was consistently reduced in the plants grown under severe but not moderate water deficit, whereas the f_gc did not decrease in the water-stressed soybean leaves. The water deficit resulted in the reduced expression of one of two (maize) or three (soybean) SPEECHLESS orthologs, and the expression patterns were correlated with SI. The vein density (VD) increased in both species in response to the water deficit, although the effect was greater in soybean. This study establishes a mechanism of stomatal development plasticity that can be applied to other species and genotypes to develop or investigate stomatal development plasticity.

LEA13 and LEA30 Are Involved in Tolerance to Water Stress and Stomata Density in Arabidopsis thaliana

Late embryogenesis abundant (LEA) proteins are a large protein family that mainly function in protecting cells from abiotic stress, but these proteins are also involved in regulating plant growth and development. In this study, we performed a functional analysis of LEA13 and LEA30 from Arabidopsis thaliana. The results showed that the expression of both genes increased when plants were subjected to drought-stressed conditions. The insertional lines lea13 and lea30 were identified for each gene, and both had a T-DNA element in the regulatory region, which caused the genes to be downregulated. Moreover, lea13 and lea30 were more sensitive to drought stress due to their higher transpiration and stomatal spacing. Microarray analysis of the lea13 background showed that genes involved in hormone signaling, stomatal development, and abiotic stress responses were misregulated. Our results showed that LEA proteins are involved in drought tolerance and participate in stomatal density.

Transcriptome profiling reveals the spatial-temporal dynamics of gene expression essential for soybean seed development

BACKGROUND

Seeds are the economic basis of oilseed crops, especially soybeans, the most widely cultivated oilseed crop worldwide. Seed development is accompanied by a multitude of diverse cellular processes, and revealing the underlying regulatory activities is critical for seed improvement.

RESULTS

In this study, we profiled the transcriptomes of developing seeds at 20, 25, 30, and 40 days after flowering (DAF), as these stages represent critical time points of seed development from early to full development. We identified a set of highly abundant genes and highlighted the importance of these genes in supporting nutrient accumulation and transcriptional regulation for seed development. We identified 8925 differentially expressed genes (DEGs) that exhibited temporal expression patterns over the course and expression specificities in distinct tissues, including seeds and nonseed tissues (roots, stems, and leaves). Genes specific to nonseed tissues might have tissue-associated roles, with relatively low transcript abundance in developing seeds, suggesting their spatially supportive roles in seed development. Coexpression network analysis identified several underexplored genes in soybeans that bridge tissue-specific gene modules.

CONCLUSIONS

Our study provides a global view of gene activities and biological processes critical for seed formation in soybeans and prioritizes a set of genes for further study. The results of this study help to elucidate the mechanism controlling seed development and storage reserves.

Cerbantez-Bueno V, Verdugo-Perales M, R. Medina H, De Folter S, Acosta-García G. Arabidopsis cysteine-rich receptor-like protein kinase CRK33 affects stomatal density and drought tolerance

Cysteine-rich receptor-like protein kinases (CRKs) are transmembrane proteins containing two domains of unknown function 26 (DUF26) RLKs in their ectodomain. Despite that CRKs control important aspects of plant development, only few proteins have functionally been characterized. In this work, we analyzed the function of CRK33 by characterizing two insertional lines. The stomatal density and stomatal index were decreased in crk33-2 and crk33-3 plants in comparison to wild-type plants, correlating with a decreased transpiration in transgenic plants and a higher drought tolerance. Furthermore, photosynthesis and stomatal conductance changed. Finally, all four stomata cell fate genes were upregulated, especially the expression of TMM and SPCH in the mutant background, suggesting a role for CRK33 in stomatal spacing.

Genome-Wide Association Analysis for Phosphorus Use Efficiency Traits in Mungbean (Vigna radiata L. Wilczek) Using Genotyping by Sequencing Approach

Mungbean (Vigna radiata L. Wilczek) is an annual grain legume crop affected by low availability of phosphorus. Phosphorus deficiency mainly affects the growth and development of plants along with changes in root morphology and increase in root-to-shoot ratio. Deciphering the genetic basis of phosphorus use efficiency (PUE) traits can benefit our understanding of mungbean tolerance to low-phosphorus condition. To address this issue, 144 diverse mungbean genotypes were evaluated for 12 PUE traits under hydroponics with optimum- and low-phosphorus levels. The broad sense heritability of traits ranged from 0.63 to 0.92 and 0.58 to 0.92 under optimum- and low-phosphorus conditions, respectively. This study, reports for the first time such a large number of genome wide Single nucleotide polymorphisms (SNPs) (76,160) in mungbean. Further, genome wide association study was conducted using 55,634 SNPs obtained by genotyping-by-sequencing method. The results indicated that total 136 SNPs shared by both GLM and MLM models were associated with tested PUE traits under different phosphorus regimes. We have identified SNPs with highest p value (-log10(p)) for some traits like, TLA and RDW with p value (-log10(p)) of more than 6.0 at LP/OP and OP condition. We have identified nine SNPs (three for TLA and six for RDW trait) which was found to be present in chromosomes 8, 4, and 7. One SNP present in Vradi07g06230 gene contains zinc finger CCCH domain. In total, 71 protein coding genes were identified, of which 13 genes were found to be putative candidate genes controlling PUE by regulating nutrient uptake and root architectural development pathways in mungbean. Moreover, we identified three potential candidate genes VRADI11G08340, VRADI01G05520, and VRADI04G10750 with missense SNPs in coding sequence region, which results in significant variation in protein structure at tertiary level. The identified SNPs and candidate genes provide the essential information for genetic studies and marker-assisted breeding program for improving low-phosphorus tolerance in mungbean.

Does Molecular and Structural Evolution Shape the Speedy Grass Stomata?

It has been increasingly important for breeding programs to be aimed at crops that are capable of coping with a changing climate, especially with regards to higher frequency and intensity of drought events. Grass stomatal complex has been proposed as an important factor that may enable grasses to adapt to water stress and variable climate conditions. There are many studies focusing on the stomatal morphology and development in the eudicot model plant Arabidopsis and monocot model plant Brachypodium. However, the comprehensive understanding of the distinction of stomatal structure and development between monocots and eudicots, especially between grasses and eudicots, are still less known at evolutionary and comparative genetic levels. Therefore, we employed the newly released version of the One Thousand Plant Transcriptome (OneKP) database and existing databases of green plant genome assemblies to explore the evolution of gene families that contributed to the formation of the unique structure and development of grass stomata. This review emphasizes the differential stomatal morphology, developmental mechanisms, and guard cell signaling in monocots and eudicots. We provide a summary of useful molecular evidences for the high water use efficiency of grass stomata that may offer new horizons for future success in breeding climate resilient crops.

The Tomato Genome Encodes SPCH, MUTE, and FAMA Candidates That Can Replace the Endogenous Functions of Their Arabidopsis Orthologs

Stomatal abundance determines the maximum potential for gas exchange between the plant and the atmosphere. In Arabidopsis, it is set during organ development through complex genetic networks linking epidermal differentiation programs with environmental response circuits. Three related bHLH transcription factors, SPCH, MUTE, and FAMA, act as positive drivers of stomata differentiation. Mutant alleles of some of these genes sustain different stomatal numbers in the mature organs and have potential to modify plant performance under different environmental conditions. However, knowledge about stomatal genes in dicotyledoneous crops is scarce. In this work, we identified the Solanum lycopersicum putative orthologs of these three master regulators and assessed their functional orthology by their ability to complement Arabidopsis loss-of-function mutants, the epidermal phenotypes elicited by their conditional overexpression, and the expression patterns of their promoter regions in Arabidopsis. Our results indicate that the tomato proteins are functionally equivalent to their Arabidopsis counterparts and that the tomato putative promoter regions display temporal and spatial expression domains similar to those reported for the Arabidopsis genes. In vivo tracking of tomato stomatal lineages in developing cotyledons revealed cell division and differentiation histories similar to those of Arabidopsis. Interestingly, the S. lycopersicum genome harbors a FAMA-like gene, expressed in leaves but functionally distinct from the true FAMA orthologue. Thus, the basic program for stomatal development in S. lycopersicum uses key conserved genetic determinants. This opens the possibility of modifying stomatal abundance in tomato through previously tested Arabidopsis alleles conferring altered stomata abundance phenotypes that correlate with physiological traits related to water status, leaf cooling, or photosynthesis.

Multiple transcriptional factors control stomata development in rice

Grass stomata can balance gas exchange and evaporation effectively in rapidly changing environments via their unique anatomical features. Although the key components of stomatal development in Arabidopsis have been largely elucidated over the past decade, the molecular mechanisms that govern stomatal development in grasses are poorly understood. Via the genome editing system and T-DNA insertion lines, the key transcriptional factors (TFs) regulating stomatal development in rice (Oryza sativa) were knocked out. A combination of genetic and biochemical assays subsequently revealed the functions of these TFs. OsSPCH/OsICE is essential for the initiation of stomatal lineage. OsMUTE/OsICE determines meristemoid to guard mother cell (GMC) transition. OsFAMA/OsICE influences subsidiary mother cell asymmetric division and mature stoma differentiation. OsFLP regulates the orientation of GMC symmetrical division. More importantly, we found that OsSCR/OsSHR controls the initiation of stomatal lineage cells and the formation of subsidiary cells. The transcription of OsSCR is activated by OsSPCH and OsMUTE. This study characterised the functions of master regulatory TFs that control each stomatal developmental stage in rice. Our findings are helpful for elucidating how various species reprogramme the molecular mechanisms to generate different stomatal types during evolution.

Leveraging RNA-Seq to Characterize Resistance to Brown Stem Rot and the Rbs3 Locus in Soybean

Brown stem rot, caused by the fungus Phialophora gregata, reduces soybean yield by up to 38%. Although three dominant resistance loci have been identified (Rbs1 to Rbs3), the gene networks responsible for pathogen recognition and defense remain unknown. Further, identification and characterization of resistant and susceptible germplasm remains difficult. We conducted RNA-Seq of infected and mock-infected leaf, stem, and root tissues of a resistant (PI 437970, Rbs3) and susceptible (Corsoy 79) genotype. Combining historical mapping data with genotype expression differences allowed us to identify a cluster of receptor-like proteins that are candidates for the Rbs3 resistance gene. Reads mapping to the Rbs3 locus were used to identify potential novel single-nucleotide polymorphisms within candidate genes that could improve phenotyping and breeding efficiency. Comparing responses to infection revealed little overlap in differential gene expression between genotypes or tissues. Gene networks associated with defense, DNA replication, and iron homeostasis are hallmarks of resistance to P. gregata. This novel research demonstrates the utility of combining contrasting genotypes, gene expression, and classical genetic studies to characterize complex disease resistance loci.

A bHLH Transcription Factor Regulates Bisbibenzyl Biosynthesis in the Liverwort Plagiochasma appendiculatum

Liverworts, a section of the bryophyte plants which pioneered the colonization of terrestrial habitats, produce cyclic bisbibenzyls as secondary metabolites. These compounds are generated via the phenylpropanoid pathway, similar to flavonoid biosynthesis, for which basic helix-loop-helix (bHLH) transcription factors have been identified as one of the important regulators in higher plants. Here, a bHLH gene homolog (PabHLH) was isolated from the liverwort species Plagiochasma appendiculatum and its contribution to bisbibenzyl biosynthesis was explored. Variation in the abundance of PabHLH transcript mirrored that of tissue bisbibenzyl content in three different liverwort tissues. A phylogenetic analysis based on the bHLH domain sequence suggested that the gene encodes a member of bHLH subgroup IIIf, which clusters proteins involved in flavonoid synthesis. The gene's transient expression in onion epidermal cells implied that its product localized to the nucleus, and a transactivation assays in yeast showed that it was able to activate transcription. In both callus and thallus, the overexpression of PabHLH boosted bisbibenzyl accumulation, while also up-regulating PaPAL, Pa4CL1, PaSTCS1 and two genes encoding P450 cytochromes, and its RNA interference (RNAi)-induced suppression down-regulated the same set of genes and reduced the accumulation of bisbibenzyls. The abundance of PaCHS and PaFNSI transcript was related to flavonoid accumulation in transgenic thallus. PabHLH represents a candidate for the metabolic engineering of bisbibenzyl content.

Genome-Wide Identification of Chalcone Reductase Gene Family in Soybean: Insight into Root-Specific GmCHRs and Phytophthora sojae Resistance

Soybean (Glycine max [L.] Merr) is one of the main grain legumes worldwide. Soybean farmers lose billions of dollars' worth of yield annually due to root and stem rot disease caused by the oomycete Phytophthora sojae. Many strategies have been developed to combat the disease, however, these methods have proven ineffective in the long term. A more cost effective and durable approach is to select a trait naturally found in soybean that can increase resistance. One such trait is the increased production of phytoalexin glyceollins in soybean. Glyceollins are isoflavonoids, synthesized via the legume-specific branch of general phenylpropanoid pathway. The first key enzyme exclusively involved in glyceollin synthesis is chalcone reductase (CHR) which coacts with chalcone synthase for the production of isoliquiritigenin, the precursor for glyceollin biosynthesis. Here we report the identification of 14 putative CHR genes in soybean where 11 of them are predicted to be functional. Our results show that GmCHRs display tissue-specific gene expression, and that only root-specific GmCHRs are induced upon P. sojae infection. Among 4 root-specific GmCHRs, GmCHR2A is located near a QTL that is linked to P. sojae resistance suggesting GmCHR2A as a novel locus for partial resistance that can be utilized for resistance breeding.

Inference of Transcription Regulatory Network in Low Phytic Acid Soybean Seeds

A dominant loss of function mutation in myo-inositol phosphate synthase (MIPS) gene and recessive loss of function mutations in two multidrug resistant protein type-ABC transporter genes not only reduce the seed phytic acid levels in soybean, but also affect the pathways associated with seed development, ultimately resulting in low emergence. To understand the regulatory mechanisms and identify key genes that intervene in the seed development process in low phytic acid crops, we performed computational inference of gene regulatory networks in low and normal phytic acid soybeans using a time course transcriptomic data and multiple network inference algorithms. We identified a set of putative candidate transcription factors and their regulatory interactions with genes that have functions in myo-inositol biosynthesis, auxin-ABA signaling, and seed dormancy. We evaluated the performance of our unsupervised network inference method by comparing the predicted regulatory network with published regulatory interactions in Arabidopsis. Some contrasting regulatory interactions were observed in low phytic acid mutants compared to non-mutant lines. These findings provide important hypotheses on expression regulation of myo-inositol metabolism and phytohormone signaling in developing low phytic acid soybeans. The computational pipeline used for unsupervised network learning in this study is provided as open source software and is freely available at https://lilabatvt.github.io/LPANetwork/.

Similarity between soybean and Arabidopsis seed methylomes and loss of non-CG methylation does not affect seed development

We profiled soybean and Arabidopsis methylomes from the globular stage through dormancy and germination to understand the role of methylation in seed formation. CHH methylation increases significantly during development throughout the entire seed, targets primarily transposable elements (TEs), is maintained during endoreduplication, and drops precipitously within the germinating seedling. By contrast, no significant global changes in CG- and CHG-context methylation occur during the same developmental period. An Arabidopsis ddcc mutant lacking CHH and CHG methylation does not affect seed development, germination, or major patterns of gene expression, implying that CHH and CHG methylation does not play a significant role in seed development or in regulating seed gene activity. By contrast, over 100 TEs are transcriptionally de-repressed in ddcc seeds, suggesting that the increase in CHH-context methylation may be a failsafe mechanism to reinforce transposon silencing. Many genes encoding important classes of seed proteins, such as storage proteins, oil biosynthesis enzymes, and transcription factors, reside in genomic regions devoid of methylation at any stage of seed development. Many other genes in these classes have similar methylation patterns, whether the genes are active or repressed. Our results suggest that methylation does not play a significant role in regulating large numbers of genes important for programming seed development in both soybean and Arabidopsis. We conclude that understanding the mechanisms controlling seed development will require determining how cis-regulatory elements and their cognate transcription factors are organized in genetic regulatory networks.

Dynamic Changes of Genome-Wide DNA Methylation during Soybean Seed Development

Seed development is programmed by expression of many genes in plants. Seed maturation is an important developmental process to soybean seed quality and yield. DNA methylation is a major epigenetic modification regulating gene expression. However, little is known about the dynamic nature of DNA methylation and its effects on gene expression during plant development. Through whole-genome bisulfite sequencing, we showed that DNA methylation went through dynamic changes during seed maturation. An average of 66% CG, 45% CHG and 9% CHH contexts was methylated in cotyledons. CHH methylation levels in cotyledons changed greatly from 6% at the early stage to 11% at the late stage. Transcribed genes were approximately two-fold more likely to be differentially methylated than non-transcribed genes. We identified 40, 66 and 2136 genes containing differentially methylated regions (DMRs) with negative correlation between their expression and methylation in the CG, CHG and CHH contexts, respectively. The majority of the DMR genes in the CHH context were transcriptionally down-regulated as seeds mature: 99% of them during early maturation were down-regulated, and preferentially associated with DNA replication and cell division. The results provide novel insights into the dynamic nature of DNA methylation and its relationship with gene regulation in seed development.

Identification of ZOUPI Orthologs in Soybean Potentially Involved in Endosperm Breakdown and Embryogenic Development

Soybean (Glycine max Merr.) is the world's most widely grown legume and provides an important source of protein and oil. Improvement of seed quality requires deep insights into the genetic regulation of seed development. The endosperm serves as a temporary source of nutrients that are transported from maternal to filial tissues, and it also generates signals for proper embryo formation. Endosperm cell death is associated with the processes of nutrient transfer and embryo expansion. The bHLH domain transcription factor AtZHOUPI (AtZOU) plays a key role in both the lysis of the transient endosperm and the formation of embryo cuticle in Arabidopsis thaliana. There are two copies of soybean GmZOU (GmZOU-1 and GmZOU-2), which fall into the same phylogenetic clade as AtZOU. These two copies share the same transcription orientation and are the result of tandem duplication. The expression of GmZOUs is limited to the endosperm, where it peaks during the heart embryo stage. When the exogenous GmZOU-1 and GmZOU-2 were expressed in the zou-4 mutant of Arabidopsis, only GmZOU-1 partially complemented the zou mutant phenotype, as indicated by endosperm breakdown and embryo cuticle formation in the transgenic lines. This research confirmed that the GmZOU-1 is a ZOU ortholog that may be responsible for endosperm breakdown and embryo cuticle formation in soybean.