Functional conservation and divergence of GmCHLI genes in polyploid soybean

Molecular, genetic, and genomic basis of seed size and yield characteristics in soybean

Soybean (Glycine max L. Merr.) is a crucial oilseed cash crop grown worldwide and consumed as oil, protein, and food by humans and feed by animals. Comparatively, soybean seed yield is lower than cereal crops, such as maize, rice, and wheat, and the demand for soybean production does not keep up with the increasing consumption level. Therefore, increasing soybean yield per unit area is the most crucial breeding objective and is challenging for the scientific community. Moreover, yield and associated traits are extensively researched in cereal crops, but little is known about soybeans' genetics, genomics, and molecular regulation of yield traits. Soybean seed yield is a complex quantitative trait governed by multiple genes. Understanding the genetic and molecular processes governing closely related attributes to seed yield is crucial to increasing soybean yield. Advances in sequencing technologies have made it possible to conduct functional genomic research to understand yield traits' genetic and molecular underpinnings. Here, we provide an overview of recent progress in the genetic regulation of seed size in soybean, molecular, genetics, and genomic bases of yield, and related key seed yield traits. In addition, phytohormones, such as auxin, gibberellins, cytokinins, and abscisic acid, regulate seed size and yield. Hence, we also highlight the implications of these factors, challenges in soybean yield, and seed trait improvement. The information reviewed in this study will help expand the knowledge base and may provide the way forward for developing high-yielding soybean cultivars for future food demands.

Characterization and fine mapping of a yellow leaf gene regulating chlorophyll biosynthesis and chloroplast development in cotton (Gossypium arboreum)

Chlorophyll biosynthesis and chloroplast development are essential for photosynthesis and plant growth. Gossypium arboreum, a valuable source of genetic variation for cotton improvement, remains poorly studied for the mechanisms regulating chlorophyll biosynthesis and chloroplast development. Here we created a G. arboreum etiolated leaf and stuntedness (els) mutant that displayed a distinct yellow color of leaves, bracts and stems throughout the whole growth, where chlorophyll accumulation in leaves was reduced and chloroplast development was delayed. The GaCHLH gene, which encodes the H subunit of magnesium chelatase (Mg-chelatase), was screened by MutMap and KASP analysis. Compared to GaCHLH, the gene Gachlh of the mutant had a single nucleotide transition (G to A) at 1549 bp, which causes the substitution of a glycine (G) by a serine (S) at the 517th amino acid, resulting in an abnormal secondary structure of the Gachlh protein. GaCHLH-silenced SXY1 and ZM24 plants exhibited a lower GaCHLH expression level, a lower chlorophyll content, and the yellow-leaf phenotype. Gachlh expression affected the expression of key genes in the tetrapyrrole pathway. GaCHLH and Gachlh were located in the chloroplasts and that alteration of the mutation site did not affect the final target position. The BiFC assay result indicated that Gachlh could not bind to GaCHLD properly, which prevented the assembly of Mg-chelatase and thus led to the failure of chlorophyll synthesis. In this study, the Gachlh gene of G. arboreum els was finely localized and identified for the first time, providing new insights into the chlorophyll biosynthesis pathway in cotton.

Systematic Characterization of the OSCA Family Members in Soybean and Validation of Their Functions in Osmotic Stress

Since we discovered OSCA1, a hyperosmolarity-gated calcium-permeable channel that acted as an osmosensor in Arabidopsis, the OSCA family has been identified genome-wide in several crops, but only a few OSCA members' functions have been experimentally demonstrated. Osmotic stress seriously restricts the yield and quality of soybean. Therefore, it is essential to decipher the molecular mechanism of how soybean responds to osmotic stress. Here, we first systematically studied and experimentally demonstrated the role of OSCA family members in the osmotic sensing of soybean. Phylogenetic relationships, gene structures, protein domains and structures analysis revealed that 20 GmOSCA members were divided into four clades, of which members in the same cluster may have more similar functions. In addition, GmOSCA members in clusters III and IV may be functionally redundant and diverged from those in clusters I and II. Based on the spatiotemporal expression patterns, GmOSCA1.6, GmOSCA2.1, GmOSCA2.6, and GmOSCA4.1 were extremely low expressed or possible pseudogenes. The remaining 16 GmOSCA genes were heterologously overexpressed in an Arabidopsis osca1 mutant, to explore their functions. Subcellular localization showed that most GmOSCA members could localize to the plasma membrane (PM). Among 16 GmOSCA genes, only overexpressing GmOSCA1.1, GmOSCA1.2, GmOSCA1.3, GmOSCA1.4, and GmOSCA1.5 in cluster I could fully complement the reduced hyperosmolality-induced [Ca²⁺]_i increase (OICI) in osca1. The expression profiles of GmOSCA genes against osmotic stress demonstrated that most GmOSCA genes, especially GmOSCA1.1, GmOSCA1.2, GmOSCA1.3, GmOSCA1.4, GmOSCA1.5, GmOSCA3.1, and GmOSCA3.2, strongly responded to osmotic stress. Moreover, overexpression of GmOSCA1.1, GmOSCA1.2, GmOSCA1.3, GmOSCA1.4, GmOSCA1.5, GmOSCA3.1, and GmOSCA3.2 rescued the drought-hypersensitive phenotype of osca1. Our findings provide important clues for further studies of GmOSCA-mediated calcium signaling in the osmotic sensing of soybean and contribute to improving soybean drought tolerance through genetic engineering and molecular breeding.

Natural allelic variation of GmST05 controlling seed size and quality in soybean

Seed size is one of the most important agronomic traits determining the yield of crops. Cloning the key genes controlling seed size and pyramiding their elite alleles will facilitate yield improvement. To date, few genes controlling seed size have been identified in soybean, a major crop that provides half of the plant oil and one quarter of the plant protein globally. Here, through a genome-wide association study of over 1800 soybean accessions, we determined that natural allelic variation at GmST05 (Seed Thickness 05) predominantly controlled seed thickness and size in soybean germplasm. Further analyses suggested that the two major haplotypes of GmST05 differed significantly at the transcriptional level. Transgenic experiments demonstrated that GmST05 positively regulated seed size and influenced oil and protein contents, possibly by regulating the transcription of GmSWEET10a. Population genetic diversity analysis suggested that allelic variations of GmST05 were selected during geographical differentiation but have not been fixed. In summary, natural variation in GmST05 determines transcription levels and influences seed size and quality in soybean, making it an important gene resource for soybean molecular breeding.

Progress in soybean functional genomics over the past decade

Soybean is one of the most important oilseed and fodder crops. Benefiting from the efforts of soybean breeders and the development of breeding technology, large number of germplasm has been generated over the last 100 years. Nevertheless, soybean breeding needs to be accelerated to meet the needs of a growing world population, to promote sustainable agriculture and to address future environmental changes. The acceleration is highly reliant on the discoveries in gene functional studies. The release of the reference soybean genome in 2010 has significantly facilitated the advance in soybean functional genomics. Here, we review the research progress in soybean omics (genomics, transcriptomics, epigenomics and proteomics), germplasm development (germplasm resources and databases), gene discovery (genes that are responsible for important soybean traits including yield, flowering and maturity, seed quality, stress resistance, nodulation and domestication) and transformation technology during the past decade. At the end, we also briefly discuss current challenges and future directions.

Molecular Characterization of Mg-Chelatase CHLI Subunit in Pea (Pisum sativum L.)

As a rate-limiting enzyme for chlorophyll biosynthesis, Mg-chelatase is a promising target for improving photosynthetic efficiency. It consists of CHLH, CHLD, and CHLI subunits. In pea (Pisum sativum L.), two putative CHLI genes (PsCHLI1 and PsCHLI2) were revealed recently by the whole genome sequencing, but their molecular features are not fully characterized. In this study, PsCHLI1 and PsCHLI2 cDNAs were identified by PCR-based cloning and sequencing. Phylogenetic analysis showed that PsCHLIs were derived from an ancient duplication in legumes. Both PsCHLIs were more highly expressed in leaves than in other organs and downregulated by abscisic acid and heat treatments, while PsCHLI1 was more highly expressed than PsCHLI2. PsCHLI1 and PsCHLI2 encode 422- and 417-amino acid proteins, respectively, which shared 82% amino acid identity and were located in chloroplasts. Plants with a silenced PsCHLI1 closely resembled PsCHLI1 and PsCHLI2 double-silenced plants, as both exhibited yellow leaves with barely detectable Mg-chelatase activity and chlorophyll content. Furthermore, plants with a silenced PsCHLI2 showed no obvious phenotype. In addition, the N-terminal fragment of PsCHLI1 (PsCHLI1N, Val63-Cys191) and the middle fragment of PsCHLI1 (PsCHLI1M, Gly192-Ser336) mediated the formation of homodimers and the interaction with CHLD, respectively, while active PsCHLI1 was only achieved by combining PsCHLI1N, PsCHLI1M, and the C-terminal fragment of PsCHLI1 (Ser337-Ser422). Taken together, PsCHLI1 is the key CHLI subunit, and its peptide fragments are essential for maintaining Mg-chelatase activity, which can be used to improve photosynthetic efficiency by manipulating Mg-chelatase in pea.

Validation of reference genes for quantitative gene expression in the Lippia alba polyploid complex (Verbenaceae)

Lippia alba (Verbenaceae) is one of the most studied species of the genus Lippia, mainly due to its medicinal properties. The species was described as a polyploid complex with five cytotypes. The comparison of gene expression in species with several ploidal levels needs to be conducted carefully due to possible changes in gene regulation. Quantitative reverse transcription PCR (qRT-PCR) is a widely used method for transcript abundance analyses in plants. Besides being an extremely powerful technique, relative quantification by Real-Time quantitative PCR (RT-qPCR) needs the normalization with a stable reference gene. We evaluated the stability of nine candidate reference genes in Lippia alba with different ploidal levels using NormFinder, geNorm, and RefFinder software. The product of each primer showed a single peak in the melting curve. The R² value ranged from 0.998 to 1000 and primers efficiency ranged from 98.95% to 129%. The CIT gene came up as a stable housekeeping gene, being appropriate for studies in polyploid accessions of Lippia alba. Considering that polyploidy is widely documented in Angiosperms, the results can be used not only for further gene expression studies in L. alba but also as a possible reference gene for other polyploid complexes. Differential stability among different genes highlights the importance of the validation of reference genes used for RT-qPCR approach in polyploid studies.

Mutation of YL Results in a Yellow Leaf with Chloroplast RNA Editing Defect in Soybean

RNA editing plays a key role in organelle gene expression. Little is known about how RNA editing factors influence soybean plant development. Here, we report the isolation and characterization of a soybean yl (yellow leaf) mutant. The yl plants showed decreased chlorophyll accumulation, lower PS II activity, an impaired net photosynthesis rate, and an altered chloroplast ultrastructure. Fine mapping of YL uncovered a point mutation in Glyma.20G187000, which encodes a chloroplast-localized protein homologous to Arabidopsis thaliana (Arabidopsis) ORRM1. YL is mainly expressed in trifoliate leaves, and its deficiency affects the editing of multiple chloroplast RNA sites, leading to inferior photosynthesis in soybean. Taken together, these results demonstrate the importance of the soybean YL protein in chloroplast RNA editing and photosynthesis.

A class B heat shock factor selected for during soybean domestication contributes to salt tolerance by promoting flavonoid biosynthesis

Soybean (Glycine max) production is severely affected in unfavorable environments. Identification of the regulatory factors conferring stress tolerance would facilitate soybean breeding. In this study, through coexpression network analysis of salt-tolerant wild soybeans, together with molecular and genetic approaches, we revealed a previously unidentified function of a class B heat shock factor, HSFB2b, in soybean salt stress response. We showed that HSFB2b improves salt tolerance through the promotion of flavonoid accumulation by activating one subset of flavonoid biosynthesis-related genes and by inhibiting the repressor gene GmNAC2 to release another subset of genes in the flavonoid biosynthesis pathway. Moreover, four promoter haplotypes of HSFB2b were identified from wild and cultivated soybeans. Promoter haplotype II from salt-tolerant wild soybean Y20, with high promoter activity under salt stress, is probably selected for during domestication. Another promoter haplotype, III, from salt-tolerant wild soybean Y55, had the highest promoter activity under salt stress, had a low distribution frequency and may be subjected to the next wave of selection. Together, our results revealed the mechanism of HSFB2b in soybean salt stress tolerance. Its promoter variations were identified, and the haplotype with high activity may be adopted for breeding better soybean cultivars that are adapted to stress conditions.

Large-Scale Investigation of Soybean Gene Functions by Overexpressing a Full-Length Soybean cDNA Library in Arabidopsis

Molecular breeding has become an important approach for crop improvement, and a prerequisite for molecular breeding is elucidation of the functions of genetic loci or genes. Soybean is one of the most important food and oil crops worldwide. However, due to the difficulty of genetic transformation in soybean, studies of its functional genomics lag far behind those of other crops such as rice, which severely impairs the progress of molecular improvement in soybean. Here, we describe an effective large-scale strategy to investigate the functions of soybean genes via overexpression of a full-length soybean cDNA library in Arabidopsis. The overexpression vector pJL12 was modified for use in the construction of a normalized full-length cDNA library. The constructed cDNA library showed good quality; repetitive clones represented approximately 4%, insertion fragments were approximately 2.2 kb, and the full-length rate was approximately 98%. This cDNA library was then overexpressed in Arabidopsis, and approximately 2000 transgenic lines were preliminarily obtained. Phenotypic analyses of the positive T₁ transgenic plants showed that more than 5% of the T₁ transgenic lines displayed abnormal developmental phenotypes, and approximately 1% of the transgenic lines exhibited potentially favorable traits. We randomly amplified 4 genes with obvious phenotypes (enlarged seeds, yellowish leaves, more branches, and dense siliques) and repeated the transgenic analyses in Arabidopsis. Subsequent phenotypic observation demonstrated that these phenotypes were indeed due to the overexpression of soybean genes. We believe our strategy represents an effective large-scale approach to investigate the functions of soybean genes and further reveal genes favorable for molecular improvement in soybean.

Fine mapping and candidate gene analysis of the virescent gene v 1 in Upland cotton (Gossypium hirsutum)

The young leaves of virescent mutants are yellowish and gradually turn green as the plants reach maturity. Understanding the genetic basis of virescent mutants can aid research of the regulatory mechanisms underlying chloroplast development and chlorophyll biosynthesis, as well as contribute to the application of virescent traits in crop breeding. In this study, fine mapping was employed, and a recessive gene (v ₁) from a virescent mutant of Upland cotton was narrowed to an 84.1-Kb region containing ten candidate genes. The GhChlI gene encodes the cotton Mg-chelatase I subunit (CHLI) and was identified as the candidate gene for the virescent mutation using gene annotation. BLAST analysis showed that the GhChlI gene has two copies, Gh_A10G0282 and Gh_D10G0283. Sequence analysis indicated that the coding region (CDS) of GhChlI is 1269 bp in length, with three predicted exons and one non-synonymous nucleotide mutation (G1082A) in the third exon of Gh_D10G0283, with an amino acid (AA) substitution of arginine (R) to lysine (K). GhChlI-silenced TM-1 plants exhibited a lower GhChlI expression level, a lower chlorophyll content, and the virescent phenotype. Analysis of upstream regulatory elements and expression levels of GhChlI showed that the expression quantity of GhChlI may be normal, and with the development of the true leaf, the increase in the Gh_A10G0282 dosage may partially make up for the deficiency of Gh_D10G0283 in the v ₁ mutant. Phylogenetic analysis and sequence alignment revealed that the protein sequence encoded by the third exon of GhChlI is highly conserved across diverse plant species, in which AA substitutions among the completely conserved residues frequently result in changes in leaf color in various species. These results suggest that the mutation (G1082A) within the GhChlI gene may cause a functional defect of the GhCHLI subunit and thus the virescent phenotype in the v₁ mutant. The GhChlI mutation not only provides a tool for understanding the associations of CHLI protein function and the chlorophyll biosynthesis pathway but also has implications for cotton breeding.

GmPHR25, a GmPHR member up-regulated by phosphate starvation, controls phosphate homeostasis in soybean

As an essential nutrient element, phosphorus (P) plays an important role in plant growth and development. Low P availability is a limiting factor for crop production, especially for legume crops (e.g. soybean), which require additional P to sustain nitrogen fixation through symbiotic associations with rhizobia. Although PHOSPHATE STARVATION RESPONSE 1 (PHR1) or PHR1-like is considered as a central regulator of phosphate (Pi) homeostasis in several plant species, it remains undefined in soybean. In this study, 35 GmPHR members were cloned from the soybean genome and expression patterns in soybean were assayed under nitrogen (N) and P deficiency conditions. GmPHR25, which is up-regulated in response to Pi starvation, was then overexpressed in soybean hairy roots in vitro and in vivo to investigate its functions. The results showed that overexpressing GmPHR25 increased Pi concentration in transgenic soybean hairy roots under normal conditions, accompanied with a significant decrease in hairy root growth. Furthermore, transcripts of 11 out of 14 high-affinity Pi transporter (GmPT) members as well as five other Pi starvation-responsive genes were significantly increased in soybean hairy roots with GmPHR25 overexpression. Taken together, this study suggests that GmPHR25 is a vital regulator in the P signaling network, and controls Pi homeostasis in soybean.